CarpetLib: Add higher order interpolation operators

1 files changed, 30 insertions, 0 deletions

diff --git a/Carpet/CarpetLib/src/data.cc b/Carpet/CarpetLib/src/data.cc
index da1f8ab56..0629a60ef 100644
--- a/Carpet/CarpetLib/src/data.cc
+++ b/Carpet/CarpetLib/src/data.cc
@@ -599,6 +599,36 @@ transfer_prolongate (data const * const src,
                         this->extent(),
                         box);
       break;
+    case 7:
+      call_operator<T> (& prolongate_3d_o7_rf2,
+                        static_cast <T const *> (src->storage()),
+                        src->shape(),
+                        static_cast <T *> (this->storage()),
+                        this->shape(),
+                        src->extent(),
+                        this->extent(),
+                        box);
+      break;
+    case 9:
+      call_operator<T> (& prolongate_3d_o9_rf2,
+                        static_cast <T const *> (src->storage()),
+                        src->shape(),
+                        static_cast <T *> (this->storage()),
+                        this->shape(),
+                        src->extent(),
+                        this->extent(),
+                        box);
+      break;
+    case 11:
+      call_operator<T> (& prolongate_3d_o11_rf2,
+                        static_cast <T const *> (src->storage()),
+                        src->shape(),
+                        static_cast <T *> (this->storage()),
+                        this->shape(),
+                        src->extent(),
+                        this->extent(),
+                        box);
+      break;
     default:
       CCTK_WARN (CCTK_WARN_ABORT,
                  "There is no stencil for op=\"LAGRANGE\" with order_space not in {1, 3, 5}");