diff options
Diffstat (limited to 'Auxiliary/Cactus/KrancNumericalTools/Perturb/src/id_perturb.c')
-rw-r--r-- | Auxiliary/Cactus/KrancNumericalTools/Perturb/src/id_perturb.c | 49 |
1 files changed, 0 insertions, 49 deletions
diff --git a/Auxiliary/Cactus/KrancNumericalTools/Perturb/src/id_perturb.c b/Auxiliary/Cactus/KrancNumericalTools/Perturb/src/id_perturb.c deleted file mode 100644 index 07bd106..0000000 --- a/Auxiliary/Cactus/KrancNumericalTools/Perturb/src/id_perturb.c +++ /dev/null @@ -1,49 +0,0 @@ -/* - $Header$ -*/ - -#include <stdlib.h> - -#include "perturb.h" - -#include "cctk_Arguments.h" -#include "cctk_Parameters.h" - -void -add_perturb_to_var (int idx, const char* optstring, void* cctkGH) -{ - DECLARE_CCTK_PARAMETERS; - int i, j, k, ijk; - CCTK_REAL* data; - cGH* GH = cctkGH; - int type; - - data = (CCTK_REAL*) CCTK_VarDataPtrI(GH, 0, idx); - - for (k=1; k< GH->cctk_lsh[2]-1; ++k) - { - for (j=1; j< GH->cctk_lsh[1]-1; ++j) - { - for (i=1; i< GH->cctk_lsh[0]-1; ++i) - { - ijk = CCTK_GFINDEX3D(GH, i, j, k); - - data[ijk] += amplitude * pow(-1, i/period) * pow(-1, j/period) * pow(-1, k/period); - } - } - } -} - - -void -id_perturb(CCTK_ARGUMENTS) -{ - DECLARE_CCTK_ARGUMENTS - DECLARE_CCTK_PARAMETERS - - if (CCTK_TraverseString(id_vars, add_perturb_to_var, cctkGH, - CCTK_GROUP_OR_VAR) < 0) - { - CCTK_WARN (1, "Failed to parse 'IDRandom::noisy_id_vars' parameter"); - } -} |