Jacobi diagonalization with cyclic sweeps
----------------------
n=5
A random symmetric matrix A:
     0.134      0.847      0.764      0.255      0.495 
     0.847      0.449      0.652      0.789      0.094 
     0.764      0.652      0.028      0.836      0.433 
     0.255      0.789      0.836      0.762      0.002 
     0.495      0.094      0.433      0.002      0.445 
sweeps: 5
matrix A after diagonalization (upper triangle should be zeroed):
     0.134      0.000      0.000      0.000      0.000 
     0.847      0.449     -0.000     -0.000     -0.000 
     0.764      0.652      0.028      0.000      0.000 
     0.255      0.789      0.836      0.762      0.000 
     0.495      0.094      0.433      0.002      0.445 
V^T*A*V (should be diagonal):
    -0.898      0.000      0.000      0.000      0.000 
     0.000     -0.556      0.000     -0.000     -0.000 
     0.000      0.000      0.037      0.000      0.000 
     0.000     -0.000     -0.000      0.684      0.000 
     0.000     -0.000      0.000      0.000      2.554 
eigenvalues (should be eqal the diagonal elements of V^T*A*V):
   -0.8984    -0.5561    0.03691     0.6836      2.554