From 836dc2db3a4417b9730a19e0685df3bec6550178 Mon Sep 17 00:00:00 2001
From: eschnett <eschnett@57fe0bb3-ccba-405f-9b23-de0201f165b7>
Date: Thu, 13 Jun 2013 18:56:08 +0000
Subject: Initialise all grid points when grid functions are padded

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
---
 src/Initialisation.c | 12 ++++++------
 1 file changed, 6 insertions(+), 6 deletions(-)

(limited to 'src')

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) {
-- 
cgit v1.2.3