class JACOBI

Flattened version is here


jacobi(aa :MAT, out d :VEC, out v:MAT) :INT
**** Compute all eigenvalues and eigenvectors of a real symmetric matrix a[N][N]. On output elements of a above the disgonal are destroyed. d[N] returns the eigenvalues of a. v[N][N] is a matrix whose columns contain, on output, the normalized eigenvectors of a. nrot returns the number of Jacobi rotations that were required.


rotate(a :MAT, i, j, k, l :INT, inout h, inout g, s, tau :FLT)
**** Do rotations required by Jacobi Transformation

The Sather Home Page