29
30
31
32#include "implicit_f.inc"
33
34
35
36#include "mvsiz_p.inc"
37
38
39
40 INTEGER JFT, JLT
41
43 . gg(*),cc(mvsiz,3,3),
cg(mvsiz,3,3),g33(mvsiz,3,3)
44
45
46
47 INTEGER I,J,K
48
50 . g_1(mvsiz)
51
52
53 DO i=jft,jlt
54 g_1(i) =one/gg(i)
55 END DO
56
57 DO j=1,3
58 DO k=j,3
59 DO i=jft,jlt
60 cc(i,j,k) = cc(i,j,k)*g_1(i)
61 g33(i,j,k) = g33(i,j,k)*g_1(i)
62 cc(i,k,j) = cc(i,j,k)
63 g33(i,k,j) = g33(i,j,k)
64 ENDDO
65 END DO
66 END DO
67
68 DO j=1,3
69 DO k=1,3
70 DO i=jft,jlt
71 cg(i,j,k) =
cg(i,j,k)*g_1(i)
72 END DO
73 END DO
74 END DO
75
76 RETURN
subroutine cg(dim, mat, rhs, sol, max_iter, tol)