3 |
! 12:53:09 lmdzadmin Exp $ |
! 12:53:09 lmdzadmin Exp $ |
4 |
|
|
5 |
SUBROUTINE jacobi(a, n, np, d, v, nrot) |
SUBROUTINE jacobi(a, n, np, d, v, nrot) |
6 |
PARAMETER (nmax=400) |
implicit none |
7 |
DIMENSION a(np, np), d(np), v(np, np), b(nmax), z(nmax) |
integer, PARAMETER:: nmax=400 |
8 |
|
integer n, np |
9 |
|
real a(np, np), d(np), v(np, np), b(nmax), z(nmax) |
10 |
|
integer ip, iq, nrot, i, j |
11 |
|
real sm, g, h, t, theta, c, s, tau, tresh |
12 |
|
|
13 |
IF (n>nmax) THEN |
IF (n>nmax) THEN |
14 |
PRINT *, 'n, nmax=', n, nmax |
PRINT *, 'n, nmax=', n, nmax |