72 |
I cg2d_b, |
I cg2d_b, |
73 |
U cg2d_x, |
U cg2d_x, |
74 |
O firstResidual, |
O firstResidual, |
75 |
|
O minResidualSq, |
76 |
O lastResidual, |
O lastResidual, |
77 |
U numIters, |
U numIters, |
78 |
|
O nIterMin, |
79 |
I myThid ) |
I myThid ) |
80 |
IMPLICIT NONE |
IMPLICIT NONE |
81 |
#include "SIZE.h" |
#include "SIZE.h" |
85 |
Real*8 cg2d_b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
Real*8 cg2d_b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
86 |
Real*8 cg2d_x(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
Real*8 cg2d_x(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
87 |
Real*8 firstResidual |
Real*8 firstResidual |
88 |
|
Real*8 minResidualSq |
89 |
Real*8 lastResidual |
Real*8 lastResidual |
90 |
INTEGER numIters |
INTEGER numIters |
91 |
|
INTEGER nIterMin |
92 |
INTEGER myThid |
INTEGER myThid |
93 |
C dummy self dependence (nontrivial so mfef90 doesn't kill it) |
C dummy self dependence (nontrivial so mfef90 doesn't kill it) |
94 |
cg2d_x(1,1,1,1)=2*cg2d_b(1,1,1,1) |
cg2d_x(1,1,1,1)=2*cg2d_b(1,1,1,1) |