diff options
Diffstat (limited to 'src/sor_flat.F')
-rw-r--r-- | src/sor_flat.F | 18 |
1 files changed, 11 insertions, 7 deletions
diff --git a/src/sor_flat.F b/src/sor_flat.F index e36514a..ac78d52 100644 --- a/src/sor_flat.F +++ b/src/sor_flat.F @@ -11,6 +11,7 @@ #include "cctk_arguments.h" #include "cctk_parameters.h" +#include "CactusElliptic/EllBase/src/EllBase.h" /*@@ @routine sor @@ -23,7 +24,8 @@ @enddesc @@*/ - subroutine sor_flat_core3d(_CCTK_FARGUMENTS, + subroutine sor_flat_core3d( + & ierr,_CCTK_FARGUMENTS, $ Mlinear_lsh,Mlinear, $ Nsource_lsh,Nsource, $ var,var_idx, @@ -33,8 +35,9 @@ _DECLARE_CCTK_FARGUMENTS DECLARE_CCTK_PARAMETERS - INTEGER CCTK_Equals - + DECLARE_CCTK_FUNCTIONS + + integer ierr INTEGER Mlinear_lsh(3) CCTK_REAL Mlinear(Mlinear_lsh(1),Mlinear_lsh(2),Mlinear_lsh(3)) @@ -87,11 +90,13 @@ c Coeeficients for the solver: 19 point stencil... logical cheb, const, none, verb integer Mlinear_storage,Nsource_storage - INTEGER sum_handle,ierr + INTEGER sum_handle c stencil size INTEGER sw(3) + ierr = ELL_SUCCESS + tol = AbsTol(1) c Get the reduction handel for the sum operation @@ -248,9 +253,8 @@ c Apply octant Symmetries enddo - write (*,*) "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!" - write (*,*) "!! WARNING: SOR SOLVER DID NOT CONVERGE !!" - write (*,*) "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!" + call CCTK_INFO("SOR solver did not converge") + ierr = ELL_NOCONVERGENCE 123 continue |