WaveToy Initial data in C++ (only gaussian initial data here for easyness)

git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWaveCXX/trunk@2 7ec00dc9-1e7a-42ad-bf73-a299ee168e72

2 files changed, 82 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;
+}
+
diff --git a/src/make.code.defn b/src/make.code.defn
new file mode 100644
index 0000000..4335306
--- /dev/null
+++ b/src/make.code.defn
@@ -0,0 +1,9 @@
+# Main make.code.defn file for thorn IDScalarWave
+# $Header$
+
+# Source files in this directory
+SRCS = InitialData.cc 
+
+# Subdirectories containing source files
+SUBDIRS = 
+