diff options
author | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-05-02 20:59:32 +0000 |
---|---|---|
committer | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-05-02 20:59:32 +0000 |
commit | 74fb1e6ea34d6e03a35ff6c158f455c39904bf5a (patch) | |
tree | d8f9b95f30517e9bafd8c67301c7383bc8beb76e /src/GRHydro_Set_TP_GF.c | |
parent | 291e94d06b30046227fb075cbfa97b9656339d5a (diff) |
file/parameter string replacement from whisky to GRHydro
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@112 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Set_TP_GF.c')
-rw-r--r-- | src/GRHydro_Set_TP_GF.c | 89 |
1 files changed, 89 insertions, 0 deletions
diff --git a/src/GRHydro_Set_TP_GF.c b/src/GRHydro_Set_TP_GF.c new file mode 100644 index 0000000..6e9fd1a --- /dev/null +++ b/src/GRHydro_Set_TP_GF.c @@ -0,0 +1,89 @@ + /*@@ + @file Set_Trivial_Riemann_Problem_Grid_Function + @date Thu May 08 + @author Frank Loeffler + @desc + This routine sets the grid function for the trivial rieman problem + bits. This is only done for debugging purposes + + @enddesc + @@*/ + +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" + +#include "SpaceMask.h" + +void Set_Trivial_Riemann_Problem_Grid_Function(CCTK_ARGUMENTS) +{ + + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + int nx, ny, nz, i, j, k, point; + int type, state; + + nx = cctk_lsh[0]; + ny = cctk_lsh[1]; + nz = cctk_lsh[2]; + + if (flux_direction[0]==1) + { + type =SpaceMask_GetTypeBits("Hydro_RiemannProblemX"); + state=SpaceMask_GetStateBits("Hydro_RiemannProblemX", "trivial"); + for (k=0 ; k<nz ; k++) + for (j=0 ; j<ny ; j++) + for (i=0 ; i<nx ; i++) + { + point = CCTK_GFINDEX3D(cctkGH,i,j,k); + if (SpaceMask_CheckStateBits(space_mask, point, type, state)) + { + GRHydro_trivial_rp_gf_x[point]=1; + } + else + { + GRHydro_trivial_rp_gf_x[point]=0; + } + } + } + if (flux_direction[0]==2) + { + type =SpaceMask_GetTypeBits("Hydro_RiemannProblemY"); + state=SpaceMask_GetStateBits("Hydro_RiemannProblemY", "trivial"); + for (k=0 ; k<nz ; k++) + for (j=0 ; j<ny ; j++) + for (i=0 ; i<nx ; i++) + { + point = CCTK_GFINDEX3D(cctkGH,i,j,k); + if (SpaceMask_CheckStateBits(space_mask, point, type, state)) + { + GRHydro_trivial_rp_gf_y[point]=1; + } + else + { + GRHydro_trivial_rp_gf_y[point]=0; + } + } + } + if (flux_direction[0]==3) + { + type =SpaceMask_GetTypeBits("Hydro_RiemannProblemZ"); + state=SpaceMask_GetStateBits("Hydro_RiemannProblemZ", "trivial"); + for (k=0 ; k<nz ; k++) + for (j=0 ; j<ny ; j++) + for (i=0 ; i<nx ; i++) + { + point = CCTK_GFINDEX3D(cctkGH,i,j,k); + if (SpaceMask_CheckStateBits(space_mask, point, type, state)) + { + GRHydro_trivial_rp_gf_z[point]=1; + } + else + { + GRHydro_trivial_rp_gf_z[point]=0; + } + } + } + return; +} |