aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authoreschnett <eschnett@57fe0bb3-ccba-405f-9b23-de0201f165b7>2013-06-13 18:56:08 +0000
committereschnett <eschnett@57fe0bb3-ccba-405f-9b23-de0201f165b7>2013-06-13 18:56:08 +0000
commit836dc2db3a4417b9730a19e0685df3bec6550178 (patch)
treef1e566c2d7a918db385dc602cb689e13ec01ee43
parentbc0686abb3264a6ae0d2adae5b9f4943a9b0611e (diff)
Initialise all grid points when grid functions are paddedHEADmaster
Use cctk_ash instead of cctk_lsh to determine total array size. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/HydroBase/trunk@67 57fe0bb3-ccba-405f-9b23-de0201f165b7
-rw-r--r--src/Initialisation.c12
1 files changed, 6 insertions, 6 deletions
diff --git a/src/Initialisation.c b/src/Initialisation.c
index 9c6e366..9aa31ec 100644
--- a/src/Initialisation.c
+++ b/src/Initialisation.c
@@ -9,7 +9,7 @@ void HydroBase_Zero (CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- int const np = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
+ int const np = cctk_ash[0] * cctk_ash[1] * cctk_ash[2];
#pragma omp parallel for
for (int i=0; i<np; ++i) {
@@ -116,7 +116,7 @@ void HydroBase_Y_e_one (CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- int const np = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
+ int const np = cctk_ash[0] * cctk_ash[1] * cctk_ash[2];
#pragma omp parallel for
for (int i=0; i<np; ++i) {
@@ -165,7 +165,7 @@ void HydroBase_Bvec_zero (CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- int const np = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
+ int const np = cctk_ash[0] * cctk_ash[1] * cctk_ash[2];
#pragma omp parallel for
for (int i=0; i<np; ++i) {
@@ -218,7 +218,7 @@ void HydroBase_Avec_zero (CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- int const np = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
+ int const np = cctk_ash[0] * cctk_ash[1] * cctk_ash[2];
#pragma omp parallel for
for (int i=0; i<np; ++i) {
@@ -271,7 +271,7 @@ void HydroBase_Aphi_zero (CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- int const np = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
+ int const np = cctk_ash[0] * cctk_ash[1] * cctk_ash[2];
#pragma omp parallel for
for (int i=0; i<np; ++i) {
@@ -318,7 +318,7 @@ void HydroBase_InitExcisionMask (CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- int const np = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
+ int const np = cctk_ash[0] * cctk_ash[1] * cctk_ash[2];
#pragma omp parallel for
for (int i=0; i<np; ++i) {