diff options
Diffstat (limited to 'Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc')
-rw-r--r-- | Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc | 44 |
1 files changed, 22 insertions, 22 deletions
diff --git a/Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc b/Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc index 72c52f05e..0b1ff0751 100644 --- a/Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc +++ b/Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc @@ -6,7 +6,7 @@ #include <cctk.h> #include <cctk_Parameters.h> -#include "operator_prototypes.hh" +#include "operator_prototypes_3d.hh" #include "typeprops.hh" using namespace std; @@ -139,7 +139,7 @@ namespace CarpetLib { goto l81; // begin j loop - l80: + l80: j = 0; js = j0; jd = dstjoff; @@ -147,7 +147,7 @@ namespace CarpetLib { goto l801; // begin i loop - l800: + l800: i = 0; is = i0; id = dstioff; @@ -155,7 +155,7 @@ namespace CarpetLib { goto l8001; // kernel - l8000: + l8000: dst[DSTIND3(id,jd,kd)] = src[SRCIND3(is,js,ks)]; i = i+1; id = id+1; @@ -163,7 +163,7 @@ namespace CarpetLib { goto l900; // kernel - l8001: + l8001: dst[DSTIND3(id,jd,kd)] = + f1 * src[SRCIND3(is ,js,ks)] + f2 * src[SRCIND3(is+1,js,ks)]; @@ -174,14 +174,14 @@ namespace CarpetLib { goto l900; // end i loop - l900: + l900: j = j+1; jd = jd+1; if (j < regjext) goto l801; goto l90; // begin i loop - l801: + l801: i = 0; is = i0; id = dstioff; @@ -189,7 +189,7 @@ namespace CarpetLib { goto l8011; // kernel - l8010: + l8010: dst[DSTIND3(id,jd,kd)] = + f1 * src[SRCIND3(is,js ,ks)] + f2 * src[SRCIND3(is,js+1,ks)]; @@ -199,7 +199,7 @@ namespace CarpetLib { goto l901; // kernel - l8011: + l8011: dst[DSTIND3(id,jd,kd)] = + f1*f1 * src[SRCIND3(is ,js ,ks)] + f2*f1 * src[SRCIND3(is+1,js ,ks)] @@ -212,7 +212,7 @@ namespace CarpetLib { goto l901; // end i loop - l901: + l901: j = j+1; jd = jd+1; js = js+1; @@ -220,14 +220,14 @@ namespace CarpetLib { goto l90; // end j loop - l90: + l90: k = k+1; kd = kd+1; if (k < regkext) goto l81; goto l9; // begin j loop - l81: + l81: j = 0; js = j0; jd = dstjoff; @@ -235,7 +235,7 @@ namespace CarpetLib { goto l811; // begin i loop - l810: + l810: i = 0; is = i0; id = dstioff; @@ -243,7 +243,7 @@ namespace CarpetLib { goto l8101; // kernel - l8100: + l8100: dst[DSTIND3(id,jd,kd)] = + f1 * src[SRCIND3(is,js,ks )] + f2 * src[SRCIND3(is,js,ks+1)]; @@ -253,7 +253,7 @@ namespace CarpetLib { goto l910; // kernel - l8101: + l8101: dst[DSTIND3(id,jd,kd)] = + f1*f1 * src[SRCIND3(is ,js,ks )] + f2*f1 * src[SRCIND3(is+1,js,ks )] @@ -266,14 +266,14 @@ namespace CarpetLib { goto l910; // end i loop - l910: + l910: j = j+1; jd = jd+1; if (j < regjext) goto l811; goto l91; // begin i loop - l811: + l811: i = 0; is = i0; id = dstioff; @@ -281,7 +281,7 @@ namespace CarpetLib { goto l8111; // kernel - l8110: + l8110: dst[DSTIND3(id,jd,kd)] = + f1*f1 * src[SRCIND3(is,js ,ks )] + f2*f1 * src[SRCIND3(is,js+1,ks )] @@ -293,7 +293,7 @@ namespace CarpetLib { goto l911; // kernel - l8111: + l8111: { T const res1 = + f1*f1*f1 * src[SRCIND3(is ,js ,ks )] @@ -314,7 +314,7 @@ namespace CarpetLib { goto l911; // end i loop - l911: + l911: j = j+1; jd = jd+1; js = js+1; @@ -322,7 +322,7 @@ namespace CarpetLib { goto l91; // end j loop - l91: + l91: k = k+1; kd = kd+1; ks = ks+1; @@ -330,7 +330,7 @@ namespace CarpetLib { goto l9; // end k loop - l9:; + l9:; } |