diff options
Diffstat (limited to 'src/InitialData.cc')
-rw-r--r-- | src/InitialData.cc | 73 |
1 files changed, 73 insertions, 0 deletions
diff --git a/src/InitialData.cc b/src/InitialData.cc new file mode 100644 index 0000000..c8d790d --- /dev/null +++ b/src/InitialData.cc @@ -0,0 +1,73 @@ + /*@@ + @file InitialData.cc + @date + @author Werner Benger + @desc + Initial data for the 3D Wave Equation + Derived from Tom Goodale + @enddesc + @@*/ + +#include "cctk.h" +#include "cctk_Flesh.h" +#include "cctk_parameters.h" +#include "cctk_Groups.h" +#include "cctk_arguments.h" +#include "cctk_Comm.h" + +#include "CactusBase/CartGrid3D/src/Symmetry.h" + +#include <math.h> + +inline double sqr(double val) +{ + return val*val; +} + + + /*@@ + @routine IDScalarWave_InitialData + @date + @author Tom Goodale + @desc + Set up initial data for the wave equation + @enddesc + @calls + @calledby + @history + + @endhistory + +@@*/ + +extern "C" void IDScalarWaveCXX_InitialData(CCTK_CARGUMENTS) +{ + DECLARE_CCTK_CARGUMENTS + DECLARE_CCTK_PARAMETERS + + for(int k=0; k<cctk_lsh[2]; k++) + for(int j=0; j<cctk_lsh[1]; j++) + for(int i=0; i<cctk_lsh[0]; i++) + { + int index = CCTK_GFINDEX3D(cctkGH,i,j,k); + + double X = x[index], Y = y[index], Z = z[index]; + double R = sqrt(X*X + Y*Y + Z*Z); + + phi_old[index] = phi[index] = + amplitude*exp( - sqr( (R - radius) / sigma ) ); + } + + // + // Apply symmetry boundary conditions + // + ApplySymmetry(cctkGH,"wavetoy::scalarevolve"); + + // + // Synchronise + // + CCTK_SyncGroup(cctkGH,"wavetoy::scalarevolve"); + + return; +} + |