86 void eigen33(const Matrix& mat, Matrix& evecs, Vector& evals)  in eigen33()  argument
137     evecs.setIdentity();  in eigen33()
144     evecs.col (2) = tmp.row (0).cross (tmp.row (1)).normalized ();  in eigen33()
148     evecs.col(1) = tmp.row (0).cross(tmp.row (1));  in eigen33()
149     Scalar n1 = evecs.col(1).norm();  in eigen33()
151       evecs.col(1) = evecs.col(2).unitOrthogonal();  in eigen33()
153       evecs.col(1) /= n1;  in eigen33()
156     evecs.col(1) = evecs.col(2).cross(evecs.col(1).cross(evecs.col(2))).normalized();  in eigen33()
157     evecs.col(0) = evecs.col(2).cross(evecs.col(1));  in eigen33()
184   Mat evecs;  in main()  local
186   BENCH(t, tries, rep, eigen33(A,evecs,evals));  in main()