aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc
diff options
context:
space:
mode:
Diffstat (limited to 'Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc')
-rw-r--r--Carpet/CarpetLib/src/prolongate_3d_o1_rf2.cc44
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:;
}