diff --git a/src/fft/ppm_fft_exec_3d_vec_bc2r_xy.f b/src/fft/ppm_fft_exec_3d_vec_bc2r_xy.f index de26a4962fea9cf7ccd936aaadda67665b365006..1f7823fd9c155192da11ba75b8ebe15d374ef4f1 100644 --- a/src/fft/ppm_fft_exec_3d_vec_bc2r_xy.f +++ b/src/fft/ppm_fft_exec_3d_vec_bc2r_xy.f @@ -117,3 +117,5 @@ RETURN END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_exec_3d_vec_c2c_z.f b/src/fft/ppm_fft_exec_3d_vec_c2c_z.f index ca01f76938fae0c0d8e6873477c0a14114fc1deb..fa87eb52fd2ac4b173447516e2bae34797f4deb6 100644 --- a/src/fft/ppm_fft_exec_3d_vec_c2c_z.f +++ b/src/fft/ppm_fft_exec_3d_vec_c2c_z.f @@ -122,3 +122,5 @@ RETURN END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_exec_3d_vec_fr2c_xy.f b/src/fft/ppm_fft_exec_3d_vec_fr2c_xy.f index 82ea0e3f128a82ba575754156c83c3fb34c4cd5b..b34ce856755749f8f2b57ea1797f3da2d4ef9409 100644 --- a/src/fft/ppm_fft_exec_3d_vec_fr2c_xy.f +++ b/src/fft/ppm_fft_exec_3d_vec_fr2c_xy.f @@ -101,3 +101,5 @@ RETURN END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_normalize_c.f b/src/fft/ppm_fft_normalize_c.f index d8f228fe0963ddfa9cd649a88000b61d8bd604f6..b95a358aafc2e8f09454ee78e57835703481f72c 100644 --- a/src/fft/ppm_fft_normalize_c.f +++ b/src/fft/ppm_fft_normalize_c.f @@ -125,3 +125,5 @@ END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_normalize_r.f b/src/fft/ppm_fft_normalize_r.f index 425ec12b10cb076ee3d918c95c96a5f7563e92ec..3109d8d9dd47442d018a27fa751ecc37f0d0a030 100644 --- a/src/fft/ppm_fft_normalize_r.f +++ b/src/fft/ppm_fft_normalize_r.f @@ -126,3 +126,5 @@ END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_plan_3d_vec_bc2c_z.f b/src/fft/ppm_fft_plan_3d_vec_bc2c_z.f index 9dd2ba22e1b6d53f99ca6d39ec3d4c43da1a3597..fdc8df86ed9f6e1b3d8cf34fca0321fcab8265a4 100644 --- a/src/fft/ppm_fft_plan_3d_vec_bc2c_z.f +++ b/src/fft/ppm_fft_plan_3d_vec_bc2c_z.f @@ -147,3 +147,5 @@ RETURN END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_plan_3d_vec_bc2r_xy.f b/src/fft/ppm_fft_plan_3d_vec_bc2r_xy.f index f83a32e45b28073063920967fe31edefd8715b03..d26221a13c0970d95c594cdaa2f36c00f57bf2c5 100644 --- a/src/fft/ppm_fft_plan_3d_vec_bc2r_xy.f +++ b/src/fft/ppm_fft_plan_3d_vec_bc2r_xy.f @@ -152,3 +152,5 @@ RETURN END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_plan_3d_vec_fc2c_z.f b/src/fft/ppm_fft_plan_3d_vec_fc2c_z.f index 262517640343fc57a631f074a5900f57d601d50c..f9f92e2d4435d763e6f4a0f9f988e38186404878 100644 --- a/src/fft/ppm_fft_plan_3d_vec_fc2c_z.f +++ b/src/fft/ppm_fft_plan_3d_vec_fc2c_z.f @@ -147,3 +147,5 @@ RETURN END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/fft/ppm_fft_plan_3d_vec_fr2c_xy.f b/src/fft/ppm_fft_plan_3d_vec_fr2c_xy.f index d1eab5b743a67e205d5818a915c43d7e1abc401f..95f61fcec4a1b68e515d2d282aa76cd05fed3bef 100644 --- a/src/fft/ppm_fft_plan_3d_vec_fr2c_xy.f +++ b/src/fft/ppm_fft_plan_3d_vec_fr2c_xy.f @@ -151,3 +151,5 @@ RETURN END SUBROUTINE __ROUTINE +#undef __ROUTINE +#undef __PREC diff --git a/src/poisson/ppm_poisson_init_predef.f b/src/poisson/ppm_poisson_init_predef.f index 104afef750e3e4f802c45cc75961e2e56677cebb..77e2d98a5146c14467a3543918fcbf1928552959 100644 --- a/src/poisson/ppm_poisson_init_predef.f +++ b/src/poisson/ppm_poisson_init_predef.f @@ -7,14 +7,14 @@ !------------------------------------------------------------------------- SUBROUTINE __ROUTINE(topoid,meshid,ppmpoisson,fieldin,fieldout,green,info& &,bc,derive) - !!! Routine to initialise Green's function solution of the Poisson - !!! equation. green is the flag defining which Green's function to use: + !!! Routine to initialise Greens function solution of the Poisson + !!! equation. green is the flag defining which Greens function to use: !!! * ppm_poisson_grn_pois_per - Poisson equation, periodic boundaries !!! * ppm_poisson_grn_pois_fre - Poisson equation, freespace boundaries (not implemented) !!! * ppm_poisson_grn_reprojec - Do vorticity reprojection to kill divergence - !!! Eventually the routine should be overloaded to accept custom Green's + !!! Eventually the routine should be overloaded to accept custom Greens !!! functions such that more general convolutions can be performed. - !!! green should be expanded to include more buildin Green's functions. + !!! green should be expanded to include more buildin Greens functions. !!! !!! The routine should accept an optional flag to toggle deallocation of !!! work arrays between calls to ppm_poisson_solve @@ -53,7 +53,7 @@ !!!ppm_poisson_grn_pois_fre - Poisson equation, freespace boundaries (not implemented) !!!ppm_poisson_grn_reprojec - Do vorticity reprojection to kill divergence !!! - !!!Eventually this should also accept custom Green's function + !!!Eventually this should also accept custom Greens function INTEGER, INTENT(OUT) :: info INTEGER,INTENT(IN),OPTIONAL :: bc !!!boundary condition for the convolution. Can be on of the following: @@ -93,7 +93,7 @@ INTEGER ,DIMENSION(__DIM) :: indl,indu INTEGER,PARAMETER :: MK = __PREC REAL(__PREC),PARAMETER :: PI=ACOS(-1.0_MK) !@ use ppm pi - !factor for the Green's function, including FFT normalization + !factor for the Greens function, including FFT normalization REAL(__PREC) :: normfac INTEGER :: i,j,k INTEGER :: kx,ky,kz @@ -512,7 +512,7 @@ ! then FFTed in those directions, mapped to z-pencils, FFTed in z and ! finally copied to ppmpoisson%fldgrnc. The real xy slabs have already ! been setup for FFTs etc so they offer a convenient container for the - ! FFTing the Green's function instead of setting up the whole apparatus + ! FFTing the Greens function instead of setting up the whole apparatus ! for this one-time affair. ! These loops must run over the padded(extended) domain thus %ndataxy ! \nabla \Psi = -\omega @@ -521,7 +521,7 @@ !------------------------------------------------------------------------- ELSE IF (green .EQ. ppm_poisson_grn_pois_fre) THEN !----------------------------------------------------------------------- - ! First initialise the real Green's function + ! First initialise the real Greens function !@alternatively this could come from as input !----------------------------------------------------------------------- !@write(*,*) 'what the fuck?' @@ -532,7 +532,7 @@ & REAL((ppmpoisson%nmxy(1))*(ppmpoisson%nmxy(2))*(ppmpoisson%nmxy(3)),MK))*dx*dy*dz !& REAL((ppmpoisson%nmxy(1))*(ppmpoisson%nmxy(2))*(ppmpoisson%nmxy(3)),MK)& !this is the correct normalization to bring one field back and forth. !remembering FFT normalization of ALL points: !vertex - !!& REAL((ppmpoisson%nmxy(1))*(ppmpoisson%nmxy(2))*(ppmpoisson%nmxy(3)),MK)& !this should be correct normalization. When back and forth transforming the green's function is correct + !!& REAL((ppmpoisson%nmxy(1))*(ppmpoisson%nmxy(2))*(ppmpoisson%nmxy(3)),MK)& !this should be correct normalization. When back and forth transforming the greens function is correct !!& *REAL((ppmpoisson%nmxy(1))*(ppmpoisson%nmxy(2))*(ppmpoisson%nmxy(3))/8,MK)) !this line is probably not necessary !!!& *REAL((ppmpoisson%nmxy(1))*(ppmpoisson%nmxy(2))*(ppmpoisson%nmxy(3)),MK)) !this line is probably not necessary !@write(*,*) ppmpoisson%nmxy, 'johannes' diff --git a/src/poisson/ppm_poisson_solve.f b/src/poisson/ppm_poisson_solve.f index 5cbd1dc70ef375ff557327d064ee9a65e6046b55..e9e03823d54b4b8e137d383731dd5f636fd74ea0 100644 --- a/src/poisson/ppm_poisson_solve.f +++ b/src/poisson/ppm_poisson_solve.f @@ -9,10 +9,10 @@ !------------------------------------------------------------------------- SUBROUTINE __ROUTINE(topoid,meshid,ppmpoisson,fieldin,fieldout,gstw,info,& & tmpcase) - !!! Routine to perform the Green's function solution of the Poisson + !!! Routine to perform the Greens function solution of the Poisson !!! equation. All settings are defined in ppm_poisson_initdef and stored !!! in the ppmpoisson plan. The tmpcase argument allows the use of a - !!! different Green's function or operation than initialised. This is + !!! different Greens function or operation than initialised. This is !!! particularly useful for vorticity reprojection !!! (ppm_poisson_grn_reprojec). !!! @@ -21,9 +21,6 @@ USE ppm_module_map_field USE ppm_module_map_field_global USE ppm_module_map - USE ppm_module_typedef !@ - USE ppm_module_data !@ - USE ppm_module_finalize !@ IMPLICIT NONE include 'mpif.h' @@ -62,10 +59,6 @@ REAL(__PREC) :: kx,ky,kz REAL(__PREC) :: phix,phiy,phiz REAL(__PREC) :: normfac - TYPE(ppm_t_equi_mesh), POINTER :: mesh => NULL() - TYPE(ppm_t_equi_mesh), POINTER :: target_mesh => NULL() - TYPE(ppm_t_topo), POINTER :: topo => NULL() - TYPE(ppm_t_topo), POINTER :: target_topo => NULL() #ifndef __NOPE INTEGER :: trank !@ @@ -353,7 +346,7 @@ trank =0 #endif !----------------------------------------------------------------------- - ! Apply the periodic Green's function + ! Apply the periodic Greens function !----------------------------------------------------------------------- IF (presentcase .EQ. ppm_poisson_grn_pois_per) THEN DO isub=1,ppmpoisson%nsublistz @@ -372,7 +365,7 @@ trank =0 ENDDO ENDDO !----------------------------------------------------------------------- - ! Apply the free-space Green's function + ! Apply the free-space Greens function !----------------------------------------------------------------------- ELSE IF (presentcase .EQ. ppm_poisson_grn_pois_fre) THEN DO isub=1,ppmpoisson%nsublistz @@ -646,20 +639,6 @@ trank =0 CALL ppm_write(ppm_rank, 'ppm_poisson_solve','Failed to push vector field.',info2) GOTO 9999 ENDIF - topo => ppm_topo(topoid)%t!@ - mesh => topo%mesh(meshid) !@ - target_topo => ppm_topo(ppmpoisson%topoidxy)%t !@ - target_mesh => target_topo%mesh(ppmpoisson%meshidxy) !@ - DO isub=1,topo%nsublist - isubl=topo%isublist(isub) - !@write(*,*) 'johannestest rank', ppm_rank,'istart', mesh%istart(:,isubl), & - !@'nnodes',mesh%nnodes(:,isubl), mesh%nm - ENDDO - DO isub=1,ppmpoisson%nsublistxy - isubl=ppmpoisson%isublistxy(isub) - !@write(*,*) 'johannestestxy rank', ppm_rank,'istart', target_mesh%istart(:,isubl), & - !@'nnodes', target_mesh%nnodes(:,isubl), target_mesh%nm - ENDDO !Send CALL ppm_map_field_send(info)