[Commits] [svn:einsteintoolkit] HydroBase/trunk/src/ (Rev. 67)

[schnetter at cct.lsu.edu]

User: eschnett
Date: 2013/06/13 01:56 PM

Modified:
 /trunk/src/
  Initialisation.c

Log:
 Initialise all grid points when grid functions are padded
 
 Use cctk_ash instead of cctk_lsh to determine total array size.

File Changes:

Directory: /trunk/src/
======================

File [modified]: Initialisation.c
Delta lines: +6 -6
===================================================================
--- trunk/src/Initialisation.c	2013-05-27 22:32:12 UTC (rev 66)
+++ trunk/src/Initialisation.c	2013-06-13 18:56:08 UTC (rev 67)
@@ -9,7 +9,7 @@
   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 @@
   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 @@
   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 @@
   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 @@
   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 @@
   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) {