aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetLib/src/prolongate_3d_cc_o2_rf2.cc
diff options
context:
space:
mode:
Diffstat (limited to 'Carpet/CarpetLib/src/prolongate_3d_cc_o2_rf2.cc')
-rw-r--r--Carpet/CarpetLib/src/prolongate_3d_cc_o2_rf2.cc545
1 files changed, 0 insertions, 545 deletions
diff --git a/Carpet/CarpetLib/src/prolongate_3d_cc_o2_rf2.cc b/Carpet/CarpetLib/src/prolongate_3d_cc_o2_rf2.cc
deleted file mode 100644
index e6372f861..000000000
--- a/Carpet/CarpetLib/src/prolongate_3d_cc_o2_rf2.cc
+++ /dev/null
@@ -1,545 +0,0 @@
-#include <algorithm>
-#include <cassert>
-#include <cmath>
-#include <cstdlib>
-
-#include <cctk.h>
-#include <cctk_Parameters.h>
-
-#include "operator_prototypes_3d.hh"
-#include "typeprops.hh"
-
-using namespace std;
-
-
-
-namespace CarpetLib {
-
-
-
-#define SRCIND3(i,j,k) \
- index3 (i, j, k, \
- srciext, srcjext, srckext)
-#define DSTIND3(i,j,k) \
- index3 (i, j, k, \
- dstiext, dstjext, dstkext)
-
-
-
- template <typename T>
- void
- prolongate_3d_cc_o2_rf2 (T const * restrict const src,
- ivect3 const & restrict srcext,
- T * restrict const dst,
- ivect3 const & restrict dstext,
- ibbox3 const & restrict srcbbox,
- ibbox3 const & restrict dstbbox,
- ibbox3 const & restrict regbbox)
- {
- typedef typename typeprops<T>::real RT;
-
-
-
- if (any (srcbbox.stride() <= regbbox.stride() or
- dstbbox.stride() != regbbox.stride()))
- {
- CCTK_WARN (0, "Internal error: strides disagree");
- }
-
- if (any (srcbbox.stride() != reffact2 * dstbbox.stride())) {
- CCTK_WARN (0, "Internal error: source strides are not twice the destination strides");
- }
-
- if (any (dstbbox.stride() % 2 != 0)) {
- CCTK_WARN (0, "Internal error: destination strides are not even");
- }
-
- // This could be handled, but is likely to point to an error
- // elsewhere
- if (regbbox.empty()) {
- CCTK_WARN (0, "Internal error: region extent is empty");
- }
-
-
-
- ivect3 const regext = regbbox.shape() / regbbox.stride();
- assert (all ((regbbox.lower() - srcbbox.lower() + regbbox.stride() / 2) % regbbox.stride() == 0));
- ivect3 const srcoff = (regbbox.lower() - srcbbox.lower() + regbbox.stride() / 2) / regbbox.stride();
- assert (all ((regbbox.lower() - dstbbox.lower()) % regbbox.stride() == 0));
- ivect3 const dstoff = (regbbox.lower() - dstbbox.lower()) / regbbox.stride();
-
-
-
- bvect3 const needoffsetlo = srcoff % reffact2 != 0;
- bvect3 const needoffsethi = (srcoff + regext - 1) % reffact2 != 0;
- ivect3 const offsetlo = either (needoffsetlo, 2, 1);
- ivect3 const offsethi = either (needoffsethi, 2, 1);
-
-
-
- if (not regbbox.expand(offsetlo, offsethi).is_contained_in(srcbbox) or
- not regbbox .is_contained_in(dstbbox))
- {
- CCTK_WARN (0, "Internal error: region extent is not contained in array extent");
- }
-
- if (any (srcext != srcbbox.shape() / srcbbox.stride() or
- dstext != dstbbox.shape() / dstbbox.stride()))
- {
- CCTK_WARN (0, "Internal error: array sizes don't agree with bounding boxes");
- }
-
-
-
- size_t const srciext = srcext[0];
- size_t const srcjext = srcext[1];
- size_t const srckext = srcext[2];
-
- size_t const dstiext = dstext[0];
- size_t const dstjext = dstext[1];
- size_t const dstkext = dstext[2];
-
- size_t const regiext = regext[0];
- size_t const regjext = regext[1];
- size_t const regkext = regext[2];
-
- size_t const srcioff = srcoff[0];
- size_t const srcjoff = srcoff[1];
- size_t const srckoff = srcoff[2];
-
- size_t const dstioff = dstoff[0];
- size_t const dstjoff = dstoff[1];
- size_t const dstkoff = dstoff[2];
-
-
-
- size_t const fi = srcioff % 2;
- size_t const fj = srcjoff % 2;
- size_t const fk = srckoff % 2;
-
- size_t const i0 = srcioff / 2;
- size_t const j0 = srcjoff / 2;
- size_t const k0 = srckoff / 2;
-
- RT const one = 1;
-
- RT const f1 = 5*one/32;
- RT const f2 = 30*one/32;
- RT const f3 = -3*one/32;
-
-
-
- // Loop over fine region
- // Label scheme: l 8 fk fj fi
-
- size_t is, js, ks;
- size_t id, jd, kd;
- size_t i, j, k;
-
- // begin k loop
- k = 0;
- ks = k0;
- kd = dstkoff;
- if (fk == 0) goto l80;
- goto l81;
-
- // begin j loop
- l80:
- j = 0;
- js = j0;
- jd = dstjoff;
- if (fj == 0) goto l800;
- goto l801;
-
- // begin i loop
- l800:
- i = 0;
- is = i0;
- id = dstioff;
- if (fi == 0) goto l8000;
- goto l8001;
-
- // kernel
- l8000:
- dst[DSTIND3(id,jd,kd)] =
- + f1*f1*f1 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f1*f1 * src[SRCIND3(is ,js-1,ks-1)]
- + f3*f1*f1 * src[SRCIND3(is+1,js-1,ks-1)]
- + f1*f2*f1 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks-1)]
- + f3*f2*f1 * src[SRCIND3(is+1,js ,ks-1)]
- + f1*f3*f1 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f3*f1 * src[SRCIND3(is ,js+1,ks-1)]
- + f3*f3*f1 * src[SRCIND3(is+1,js+1,ks-1)]
- + f1*f1*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js-1,ks )]
- + f3*f1*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f1*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f3*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f1*f3*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js+1,ks )]
- + f3*f3*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f1*f1*f3 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f1*f3 * src[SRCIND3(is ,js-1,ks+1)]
- + f3*f1*f3 * src[SRCIND3(is+1,js-1,ks+1)]
- + f1*f2*f3 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks+1)]
- + f3*f2*f3 * src[SRCIND3(is+1,js ,ks+1)]
- + f1*f3*f3 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f3*f3 * src[SRCIND3(is ,js+1,ks+1)]
- + f3*f3*f3 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- if (i < regiext) goto l8001;
- goto l900;
-
- // kernel
- l8001:
- dst[DSTIND3(id,jd,kd)] =
- + f3*f1*f1 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f1*f1 * src[SRCIND3(is ,js-1,ks-1)]
- + f1*f1*f1 * src[SRCIND3(is+1,js-1,ks-1)]
- + f3*f2*f1 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks-1)]
- + f1*f2*f1 * src[SRCIND3(is+1,js ,ks-1)]
- + f3*f3*f1 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f3*f1 * src[SRCIND3(is ,js+1,ks-1)]
- + f1*f3*f1 * src[SRCIND3(is+1,js+1,ks-1)]
- + f3*f1*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js-1,ks )]
- + f1*f1*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f3*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f1*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f3*f3*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js+1,ks )]
- + f1*f3*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f3*f1*f3 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f1*f3 * src[SRCIND3(is ,js-1,ks+1)]
- + f1*f1*f3 * src[SRCIND3(is+1,js-1,ks+1)]
- + f3*f2*f3 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks+1)]
- + f1*f2*f3 * src[SRCIND3(is+1,js ,ks+1)]
- + f3*f3*f3 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f3*f3 * src[SRCIND3(is ,js+1,ks+1)]
- + f1*f3*f3 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- is = is+1;
- if (i < regiext) goto l8000;
- goto l900;
-
- // end i loop
- l900:
- j = j+1;
- jd = jd+1;
- if (j < regjext) goto l801;
- goto l90;
-
- // begin i loop
- l801:
- i = 0;
- is = i0;
- id = dstioff;
- if (fi == 0) goto l8010;
- goto l8011;
-
- // kernel
- l8010:
- dst[DSTIND3(id,jd,kd)] =
- + f1*f3*f1 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f3*f1 * src[SRCIND3(is ,js-1,ks-1)]
- + f3*f3*f1 * src[SRCIND3(is+1,js-1,ks-1)]
- + f1*f2*f1 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks-1)]
- + f3*f2*f1 * src[SRCIND3(is+1,js ,ks-1)]
- + f1*f1*f1 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f1*f1 * src[SRCIND3(is ,js+1,ks-1)]
- + f3*f1*f1 * src[SRCIND3(is+1,js+1,ks-1)]
- + f1*f3*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js-1,ks )]
- + f3*f3*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f1*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f3*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f1*f1*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js+1,ks )]
- + f3*f1*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f1*f3*f3 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f3*f3 * src[SRCIND3(is ,js-1,ks+1)]
- + f3*f3*f3 * src[SRCIND3(is+1,js-1,ks+1)]
- + f1*f2*f3 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks+1)]
- + f3*f2*f3 * src[SRCIND3(is+1,js ,ks+1)]
- + f1*f1*f3 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f1*f3 * src[SRCIND3(is ,js+1,ks+1)]
- + f3*f1*f3 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- if (i < regiext) goto l8011;
- goto l901;
-
- // kernel
- l8011:
- dst[DSTIND3(id,jd,kd)] =
- + f3*f3*f1 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f3*f1 * src[SRCIND3(is ,js-1,ks-1)]
- + f1*f3*f1 * src[SRCIND3(is+1,js-1,ks-1)]
- + f3*f2*f1 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks-1)]
- + f1*f2*f1 * src[SRCIND3(is+1,js ,ks-1)]
- + f3*f1*f1 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f1*f1 * src[SRCIND3(is ,js+1,ks-1)]
- + f1*f1*f1 * src[SRCIND3(is+1,js+1,ks-1)]
- + f3*f3*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js-1,ks )]
- + f1*f3*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f3*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f1*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f3*f1*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js+1,ks )]
- + f1*f1*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f3*f3*f3 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f3*f3 * src[SRCIND3(is ,js-1,ks+1)]
- + f1*f3*f3 * src[SRCIND3(is+1,js-1,ks+1)]
- + f3*f2*f3 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks+1)]
- + f1*f2*f3 * src[SRCIND3(is+1,js ,ks+1)]
- + f3*f1*f3 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f1*f3 * src[SRCIND3(is ,js+1,ks+1)]
- + f1*f1*f3 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- is = is+1;
- if (i < regiext) goto l8010;
- goto l901;
-
- // end i loop
- l901:
- j = j+1;
- jd = jd+1;
- js = js+1;
- if (j < regjext) goto l800;
- goto l90;
-
- // end j loop
- l90:
- k = k+1;
- kd = kd+1;
- if (k < regkext) goto l81;
- goto l9;
-
- // begin j loop
- l81:
- j = 0;
- js = j0;
- jd = dstjoff;
- if (fj == 0) goto l810;
- goto l811;
-
- // begin i loop
- l810:
- i = 0;
- is = i0;
- id = dstioff;
- if (fi == 0) goto l8100;
- goto l8101;
-
- // kernel
- l8100:
- dst[DSTIND3(id,jd,kd)] =
- + f1*f1*f3 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f1*f3 * src[SRCIND3(is ,js-1,ks-1)]
- + f3*f1*f3 * src[SRCIND3(is+1,js-1,ks-1)]
- + f1*f2*f3 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks-1)]
- + f3*f2*f3 * src[SRCIND3(is+1,js ,ks-1)]
- + f1*f3*f3 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f3*f3 * src[SRCIND3(is ,js+1,ks-1)]
- + f3*f3*f3 * src[SRCIND3(is+1,js+1,ks-1)]
- + f1*f1*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js-1,ks )]
- + f3*f1*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f1*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f3*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f1*f3*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js+1,ks )]
- + f3*f3*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f1*f1*f1 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f1*f1 * src[SRCIND3(is ,js-1,ks+1)]
- + f3*f1*f1 * src[SRCIND3(is+1,js-1,ks+1)]
- + f1*f2*f1 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks+1)]
- + f3*f2*f1 * src[SRCIND3(is+1,js ,ks+1)]
- + f1*f3*f1 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f3*f1 * src[SRCIND3(is ,js+1,ks+1)]
- + f3*f3*f1 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- if (i < regiext) goto l8101;
- goto l910;
-
- // kernel
- l8101:
- dst[DSTIND3(id,jd,kd)] =
- + f3*f1*f3 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f1*f3 * src[SRCIND3(is ,js-1,ks-1)]
- + f1*f1*f3 * src[SRCIND3(is+1,js-1,ks-1)]
- + f3*f2*f3 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks-1)]
- + f1*f2*f3 * src[SRCIND3(is+1,js ,ks-1)]
- + f3*f3*f3 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f3*f3 * src[SRCIND3(is ,js+1,ks-1)]
- + f1*f3*f3 * src[SRCIND3(is+1,js+1,ks-1)]
- + f3*f1*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js-1,ks )]
- + f1*f1*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f3*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f1*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f3*f3*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js+1,ks )]
- + f1*f3*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f3*f1*f1 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f1*f1 * src[SRCIND3(is ,js-1,ks+1)]
- + f1*f1*f1 * src[SRCIND3(is+1,js-1,ks+1)]
- + f3*f2*f1 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks+1)]
- + f1*f2*f1 * src[SRCIND3(is+1,js ,ks+1)]
- + f3*f3*f1 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f3*f1 * src[SRCIND3(is ,js+1,ks+1)]
- + f1*f3*f1 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- is = is+1;
- if (i < regiext) goto l8100;
- goto l910;
-
- // end i loop
- l910:
- j = j+1;
- jd = jd+1;
- if (j < regjext) goto l811;
- goto l91;
-
- // begin i loop
- l811:
- i = 0;
- is = i0;
- id = dstioff;
- if (fi == 0) goto l8110;
- goto l8111;
-
- // kernel
- l8110:
- dst[DSTIND3(id,jd,kd)] =
- + f1*f3*f3 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f3*f3 * src[SRCIND3(is ,js-1,ks-1)]
- + f3*f3*f3 * src[SRCIND3(is+1,js-1,ks-1)]
- + f1*f2*f3 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks-1)]
- + f3*f2*f3 * src[SRCIND3(is+1,js ,ks-1)]
- + f1*f1*f3 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f1*f3 * src[SRCIND3(is ,js+1,ks-1)]
- + f3*f1*f3 * src[SRCIND3(is+1,js+1,ks-1)]
- + f1*f3*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js-1,ks )]
- + f3*f3*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f1*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f3*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f1*f1*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js+1,ks )]
- + f3*f1*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f1*f3*f1 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f3*f1 * src[SRCIND3(is ,js-1,ks+1)]
- + f3*f3*f1 * src[SRCIND3(is+1,js-1,ks+1)]
- + f1*f2*f1 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks+1)]
- + f3*f2*f1 * src[SRCIND3(is+1,js ,ks+1)]
- + f1*f1*f1 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f1*f1 * src[SRCIND3(is ,js+1,ks+1)]
- + f3*f1*f1 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- if (i < regiext) goto l8111;
- goto l911;
-
- // kernel
- l8111:
- dst[DSTIND3(id,jd,kd)] =
- + f3*f3*f3 * src[SRCIND3(is-1,js-1,ks-1)]
- + f2*f3*f3 * src[SRCIND3(is ,js-1,ks-1)]
- + f1*f3*f3 * src[SRCIND3(is+1,js-1,ks-1)]
- + f3*f2*f3 * src[SRCIND3(is-1,js ,ks-1)]
- + f2*f2*f3 * src[SRCIND3(is ,js ,ks-1)]
- + f1*f2*f3 * src[SRCIND3(is+1,js ,ks-1)]
- + f3*f1*f3 * src[SRCIND3(is-1,js+1,ks-1)]
- + f2*f1*f3 * src[SRCIND3(is ,js+1,ks-1)]
- + f1*f1*f3 * src[SRCIND3(is+1,js+1,ks-1)]
- + f3*f3*f2 * src[SRCIND3(is-1,js-1,ks )]
- + f2*f3*f2 * src[SRCIND3(is ,js-1,ks )]
- + f1*f3*f2 * src[SRCIND3(is+1,js-1,ks )]
- + f3*f2*f2 * src[SRCIND3(is-1,js ,ks )]
- + f2*f2*f2 * src[SRCIND3(is ,js ,ks )]
- + f1*f2*f2 * src[SRCIND3(is+1,js ,ks )]
- + f3*f1*f2 * src[SRCIND3(is-1,js+1,ks )]
- + f2*f1*f2 * src[SRCIND3(is ,js+1,ks )]
- + f1*f1*f2 * src[SRCIND3(is+1,js+1,ks )]
- + f3*f3*f1 * src[SRCIND3(is-1,js-1,ks+1)]
- + f2*f3*f1 * src[SRCIND3(is ,js-1,ks+1)]
- + f1*f3*f1 * src[SRCIND3(is+1,js-1,ks+1)]
- + f3*f2*f1 * src[SRCIND3(is-1,js ,ks+1)]
- + f2*f2*f1 * src[SRCIND3(is ,js ,ks+1)]
- + f1*f2*f1 * src[SRCIND3(is+1,js ,ks+1)]
- + f3*f1*f1 * src[SRCIND3(is-1,js+1,ks+1)]
- + f2*f1*f1 * src[SRCIND3(is ,js+1,ks+1)]
- + f1*f1*f1 * src[SRCIND3(is+1,js+1,ks+1)];
- i = i+1;
- id = id+1;
- is = is+1;
- if (i < regiext) goto l8110;
- goto l911;
-
- // end i loop
- l911:
- j = j+1;
- jd = jd+1;
- js = js+1;
- if (j < regjext) goto l810;
- goto l91;
-
- // end j loop
- l91:
- k = k+1;
- kd = kd+1;
- ks = ks+1;
- if (k < regkext) goto l80;
- goto l9;
-
- // end k loop
- l9:;
-
- }
-
-
-
-#define TYPECASE(N,T) \
- template \
- void \
- prolongate_3d_cc_o2_rf2 (T const * restrict const src, \
- ivect3 const & restrict srcext, \
- T * restrict const dst, \
- ivect3 const & restrict dstext, \
- ibbox3 const & restrict srcbbox, \
- ibbox3 const & restrict dstbbox, \
- ibbox3 const & restrict regbbox);
-#include "typecase.hh"
-#undef TYPECASE
-
-
-
-} // CarpetLib
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-29 19:16:56 +0000