In the first paragraph of the 2nd section of this article, it is stated that given a set of linearly independent $n$-qubit state vectors, Alice can perform the Gram-Schmidt procedure to obtain orthogonal quantum states.
But I know the Gram-Shmidt process is not a unitary operation, so is it possible to come up with an algorithm, using ancillary qubits or something?
I.e., if we have $M$ arbitrary, linearly independent $n$-qubit state vectors $\{ |x_1\rangle, |x_2\rangle, |x_3\rangle, \cdots , |x_M\rangle \}$ do we have an algorithm to perform Gram-Schmidt orthogonalization on this set to obtain a set of mutually orthogonal state vectors and can we construct the corresponding circuit?