50 arma::mat SymmetricOrth(const arma::mat & Svec, const arma::vec & Sval) {  in SymmetricOrth()  argument
60   return Svec*arma::diagmat(Sinvh)*trans(Svec);  in SymmetricOrth()
68   arma::mat Svec;  in SymmetricOrth()  local
69   eig_sym_ordered(Sval,Svec,S);  in SymmetricOrth()
72   return SymmetricOrth(Svec,Sval);  in SymmetricOrth()
75 arma::mat CanonicalOrth(const arma::mat & Svec, const arma::vec & Sval, double cutoff) {  in CanonicalOrth()  argument
77   const size_t Nbf=Svec.n_rows;  in CanonicalOrth()
89     Sinvh.col(i)=Svec.col(Ndep+i)/sqrt(Sval(Ndep+i));  in CanonicalOrth()
106   arma::mat Svec;  in CanonicalOrth()  local
109   eig_sym_ordered(Sval,Svec,S);  in CanonicalOrth()
110   return CanonicalOrth(Svec,Sval,cutoff);  in CanonicalOrth()
148   arma::mat Svec;  in BasOrth()  local
150   eig_sym_ordered(Sval,Svec,S);  in BasOrth()
167     return SymmetricOrth(Svec,Sval);  in BasOrth()
172     return CanonicalOrth(Svec,Sval,tol);  in BasOrth()
221   arma::mat Svec;  in S_half_invhalf()  local
224   eig_sym_ordered(Sval,Svec,S);  in S_half_invhalf()
239     Sinvh+=Svec.col(icol)*arma::trans(Svec.col(icol))/sqrt(Sval(icol));  in S_half_invhalf()
240     Shalf+=Svec.col(icol)*arma::trans(Svec.col(icol))*sqrt(Sval(icol));  in S_half_invhalf()