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) {