diff options
Diffstat (limited to 'ML_hydro/src/hydro_soundWave.cc')
| -rw-r--r-- | ML_hydro/src/hydro_soundWave.cc | 26 |
1 files changed, 14 insertions, 12 deletions
diff --git a/ML_hydro/src/hydro_soundWave.cc b/ML_hydro/src/hydro_soundWave.cc index 6b8c0be..64f5d15 100644 --- a/ML_hydro/src/hydro_soundWave.cc +++ b/ML_hydro/src/hydro_soundWave.cc @@ -28,8 +28,6 @@ static void hydro_soundWave_Body(cGH const * restrict const cctkGH, int const di DECLARE_CCTK_PARAMETERS; - /* Declare finite differencing variables */ - /* Include user-supplied include files */ /* Initialise finite differencing variables */ @@ -60,9 +58,9 @@ static void hydro_soundWave_Body(cGH const * restrict const cctkGH, int const di CCTK_REAL const p1o2dx = 0.5*INV(dx); CCTK_REAL const p1o2dy = 0.5*INV(dy); CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx)*INV(dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx)*INV(dz); - CCTK_REAL const p1o4dydz = 0.25*INV(dy)*INV(dz); + CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); + CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); + CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); CCTK_REAL const p1odx2 = INV(SQR(dx)); CCTK_REAL const p1ody2 = INV(SQR(dy)); CCTK_REAL const p1odz2 = INV(SQR(dz)); @@ -77,7 +75,7 @@ static void hydro_soundWave_Body(cGH const * restrict const cctkGH, int const di /* Loop over the grid points */ #pragma omp parallel - CCTK_LOOP3 (hydro_soundWave, + CCTK_LOOP3(hydro_soundWave, i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2], cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) { @@ -95,11 +93,11 @@ static void hydro_soundWave_Body(cGH const * restrict const cctkGH, int const di /* Calculate temporaries and grid functions */ CCTK_REAL rhoL = 1.; - CCTK_REAL vel1L = Sin(2*xL*Pi*INV(ToReal(L)))*ToReal(A); + CCTK_REAL vel1L = sin(2*xL*Pi*INV(ToReal(L)))*ToReal(A); - CCTK_REAL vel2L = Sin(2*xL*Pi*INV(ToReal(L)))*ToReal(A); + CCTK_REAL vel2L = sin(2*xL*Pi*INV(ToReal(L)))*ToReal(A); - CCTK_REAL vel3L = Sin(2*xL*Pi*INV(ToReal(L)))*ToReal(A); + CCTK_REAL vel3L = sin(2*xL*Pi*INV(ToReal(L)))*ToReal(A); CCTK_REAL epsL = 1.; @@ -110,7 +108,7 @@ static void hydro_soundWave_Body(cGH const * restrict const cctkGH, int const di vel2[index] = vel2L; vel3[index] = vel3L; } - CCTK_ENDLOOP3 (hydro_soundWave); + CCTK_ENDLOOP3(hydro_soundWave); } extern "C" void hydro_soundWave(CCTK_ARGUMENTS) @@ -129,11 +127,15 @@ extern "C" void hydro_soundWave(CCTK_ARGUMENTS) return; } - const char *groups[] = {"ML_hydro::eps_group","grid::coordinates","ML_hydro::rho_group","ML_hydro::vel_group"}; + const char *const groups[] = { + "ML_hydro::eps_group", + "grid::coordinates", + "ML_hydro::rho_group", + "ML_hydro::vel_group"}; GenericFD_AssertGroupStorage(cctkGH, "hydro_soundWave", 4, groups); - GenericFD_LoopOverEverything(cctkGH, &hydro_soundWave_Body); + GenericFD_LoopOverEverything(cctkGH, hydro_soundWave_Body); if (verbose > 1) { |