96 void eigen33(const Matrix& mat, Matrix& evecs, Vector& evals)  in eigen33()  argument
144     evecs.setIdentity();  in eigen33()
151     evecs.col (2) = tmp.row (0).cross (tmp.row (1)).normalized ();  in eigen33()
155     evecs.col(1) = tmp.row (0).cross(tmp.row (1));  in eigen33()
156     Scalar n1 = evecs.col(1).norm();  in eigen33()
158       evecs.col(1) = evecs.col(2).unitOrthogonal();  in eigen33()
160       evecs.col(1) /= n1;  in eigen33()
163     evecs.col(1) = evecs.col(2).cross(evecs.col(1).cross(evecs.col(2))).normalized();  in eigen33()
164     evecs.col(0) = evecs.col(2).cross(evecs.col(1));  in eigen33()
185   Mat evecs;  in main()  local
187   BENCH(t, tries, rep, eigen33(A,evecs,evals));  in main()
193     if(evecs.col(k).dot(eig.eigenvectors().col(k))<0)  in main()
194       evecs.col(k) = -evecs.col(k);  in main()
195   std::cerr << evecs - eig.eigenvectors() << "\n\n";  in main()