diff options
Diffstat (limited to 'src/WaveToy.cc')
-rw-r--r-- | src/WaveToy.cc | 27 |
1 files changed, 24 insertions, 3 deletions
diff --git a/src/WaveToy.cc b/src/WaveToy.cc index 66e3be3..bd2237d 100644 --- a/src/WaveToy.cc +++ b/src/WaveToy.cc @@ -115,6 +115,10 @@ extern "C" void WaveToyCXX_Boundaries(CCTK_ARGUMENTS) int ierr=-1; int sw[3]; + CCTK_REAL finf = 1.0; + + int npow = 1; + /* Set the stencil width */ sw[0]=1; sw[1]=1; @@ -122,14 +126,31 @@ extern "C" void WaveToyCXX_Boundaries(CCTK_ARGUMENTS) ierr = CartSymGN(cctkGH,"wavetoy::scalarevolve"); - if (CCTK_Equals(bound,"flat")) + if (CCTK_EQUALS(bound,"flat")) { ierr = BndFlatVN(cctkGH,sw,"wavetoy::phi"); } + else if (CCTK_Equals(bound,"static")) + { + ierr = BndStaticVN(cctkGH,sw,"wavetoy::phi"); + } else if (CCTK_Equals(bound,"radiation")) { - ierr = BndRadiativeVN(cctkGH,sw,0,1,"wavetoy::phi","wavetoy::phi"); - } + ierr = BndRadiativeVN(cctkGH,sw,0.0,1.0,"wavetoy::phi", + "wavetoy::phi"); + } + else if (CCTK_Equals(bound,"robin")) + { + ierr = BndRobinVN(cctkGH,sw,finf,npow,"wavetoy::phi"); + } + else if (CCTK_Equals(bound,"zero")) + { + ierr = BndScalarVN(cctkGH,sw,0.0,"wavetoy::phi"); + } + else if (! CCTK_Equals(bound,"none")) + { + CCTK_WARN(0,"Boundary condition not recognized"); + } if (ierr < 0) { |