219 |
I cg2d_b, |
I cg2d_b, |
220 |
U cg2d_x, |
U cg2d_x, |
221 |
O firstResidual, |
O firstResidual, |
222 |
|
O minResidualSq, |
223 |
O lastResidual, |
O lastResidual, |
224 |
U numIters, |
U numIters, |
225 |
|
O nIterMin, |
226 |
I myThid ) |
I myThid ) |
227 |
IMPLICIT NONE |
IMPLICIT NONE |
228 |
#include "SIZE.h" |
#include "SIZE.h" |
232 |
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) |
233 |
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) |
234 |
Real*8 firstResidual |
Real*8 firstResidual |
235 |
|
Real*8 minResidualSq |
236 |
Real*8 lastResidual |
Real*8 lastResidual |
237 |
INTEGER numIters |
INTEGER numIters |
238 |
|
INTEGER nIterMin |
239 |
INTEGER myThid |
INTEGER myThid |
240 |
C dummy self dependence (nontrivial so mfef90 doesn't kill it) |
C dummy self dependence (nontrivial so mfef90 doesn't kill it) |
241 |
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) |