The prolongation method can now be controlled from the tags table defined in the interface.ccl. The prolongation method defaults to "Lagrange" and the other options are "None" and "TVD" (which I'm now sure doesn't work properly).

The prolongation method can now be controlled from the tags table defined in the interface.ccl. The prolongation method defaults to "Lagrange" and the other options are "None" and "TVD" (which I'm now sure doesn't work properly).
As an example, from CartGrid3D:

REAL coordinates type=GF tags='Prolongation="None"'
{
  x, y, z, r
}

switches off prolongation for the coordinates.

Also, a one off opportunity to marvel at C++ code that combines the worst features of Fortran with the worst features of C...

darcs-hash:20031015051401-58737-870c23e54e58164a45ffe3438e393587b1cc670b.gz

2 files changed, 263 insertions, 121 deletions

diff --git a/Carpet/CarpetLib/src/data.cc b/Carpet/CarpetLib/src/data.cc
index da9b54b2b..f676c5806 100644
--- a/Carpet/CarpetLib/src/data.cc
+++ b/Carpet/CarpetLib/src/data.cc
@@ -1,4 +1,4 @@
-// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/data.cc,v 1.30 2003/10/14 16:39:16 schnetter Exp $
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/data.cc,v 1.31 2003/10/15 07:14:01 hawke Exp $
 
 #include <assert.h>
 #include <limits.h>
@@ -17,6 +17,9 @@
 
 #include "data.hh"
 
+#include "util_ErrorCodes.h"
+#include "util_Table.h"
+
 using namespace std;
 
 
@@ -318,6 +321,14 @@ extern "C" {
      const int srcbbox[3][3],
      const int dstbbox[3][3],
      const int regbbox[3][3]);
+  void CCTK_FCALL CCTK_FNAME(prolongate_3d_real8_minmod)
+    (const CCTK_REAL8* src,
+     const int& srciext, const int& srcjext, const int& srckext,
+     CCTK_REAL8* dst,
+     const int& dstiext, const int& dstjext, const int& dstkext,
+     const int srcbbox[3][3],
+     const int dstbbox[3][3],
+     const int regbbox[3][3]);
   void CCTK_FCALL CCTK_FNAME(prolongate_3d_real8_o5)
     (const CCTK_REAL8* src,
      const int& srciext, const int& srcjext, const int& srckext,
@@ -345,6 +356,15 @@ extern "C" {
      const int srcbbox[3][3],
      const int dstbbox[3][3],
      const int regbbox[3][3]);
+  void CCTK_FCALL CCTK_FNAME(prolongate_3d_real8_2tl_minmod)
+    (const CCTK_REAL8* src1, const CCTK_REAL8& t1,
+     const CCTK_REAL8* src2, const CCTK_REAL8& t2,
+     const int& srciext, const int& srcjext, const int& srckext,
+     CCTK_REAL8* dst, const CCTK_REAL8& t,
+     const int& dstiext, const int& dstjext, const int& dstkext,
+     const int srcbbox[3][3],
+     const int dstbbox[3][3],
+     const int regbbox[3][3]);
   void CCTK_FCALL CCTK_FNAME(prolongate_3d_real8_2tl_o5)
     (const CCTK_REAL8* src1, const CCTK_REAL8& t1,
      const CCTK_REAL8* src2, const CCTK_REAL8& t2,
@@ -375,6 +395,16 @@ extern "C" {
      const int srcbbox[3][3],
      const int dstbbox[3][3],
      const int regbbox[3][3]);
+  void CCTK_FCALL CCTK_FNAME(prolongate_3d_real8_3tl_minmod)
+    (const CCTK_REAL8* src1, const CCTK_REAL8& t1,
+     const CCTK_REAL8* src2, const CCTK_REAL8& t2,
+     const CCTK_REAL8* src3, const CCTK_REAL8& t3,
+     const int& srciext, const int& srcjext, const int& srckext,
+     CCTK_REAL8* dst, const CCTK_REAL8& t,
+     const int& dstiext, const int& dstjext, const int& dstkext,
+     const int srcbbox[3][3],
+     const int dstbbox[3][3],
+     const int regbbox[3][3]);
   void CCTK_FCALL CCTK_FNAME(prolongate_3d_real8_3tl_o5)
     (const CCTK_REAL8* src1, const CCTK_REAL8& t1,
      const CCTK_REAL8* src2, const CCTK_REAL8& t2,
@@ -438,6 +468,70 @@ void data<CCTK_REAL8,3>
     regbbox[1][d] = box.upper()[d];
     regbbox[2][d] = box.stride()[d];
   }
+
+  const int varindex = (gsrcs[0])->var_index();
+  const int groupindex = CCTK_GroupIndexFromVarI(varindex);
+  const int group_tags_table = CCTK_GroupTagsTableI(groupindex);
+  assert(group_tags_table >= 0);
+  
+  int typecode = -1;
+  int keysize = -1;
+  if (! Util_TableQueryValueInfo(group_tags_table,
+                                 &typecode,
+                                 &keysize,
+                                 "Prolongation"))
+  {
+    CCTK_VWarn(4, __LINE__, __FILE__, CCTK_THORNSTRING,
+               "Tags table for group %s does not contain 'Prolongation'"
+               " tag", CCTK_GroupName(groupindex));
+  }
+  else
+  {  
+    assert(typecode == CCTK_VARIABLE_CHAR);
+  }
+
+#define PROLONG_NONE 0
+#define PROLONG_LAGRANGE 1
+#define PROLONG_TVD 2
+
+  int prolong_method;
+  if (keysize < 0) { // No key in table - default to Lagrange.
+    prolong_method = 1;
+  }
+  else {
+    char prolong_string[keysize+10];
+    const int error = Util_TableGetString(group_tags_table,
+                                          keysize+10,
+                                          prolong_string,
+                                          "Prolongation");
+    if (error < 0) {
+      CCTK_VWarn(0, __LINE__, __FILE__, CCTK_THORNSTRING,
+                 "Error code %d getting prolongation method"
+                 " from tags table for group %s.",
+                 error, CCTK_GroupName(groupindex));
+    }
+    if (CCTK_Equals(prolong_string, "None")){
+      prolong_method = PROLONG_NONE;
+    }
+    else if (CCTK_Equals(prolong_string, "Lagrange")){
+      prolong_method = PROLONG_LAGRANGE;
+    }
+    else if (CCTK_Equals(prolong_string, "TVD")){
+      prolong_method = PROLONG_TVD;
+    }
+    else {
+      CCTK_VWarn(0, __LINE__, __FILE__, CCTK_THORNSTRING,
+                 "Prolongation method %s for group %s not recognized.\n"
+                 "Is this is a typo? Known values are:\n"
+                 "None\nLagrange\nTVD",
+                 prolong_method, CCTK_GroupName(groupindex));
+    }
+      
+  }
+
+//   CCTK_VInfo(CCTK_THORNSTRING,
+//              "Variable %d is %s with method %d", 
+//              varindex, CCTK_VarName(varindex), prolong_method);
   
   assert (all(dext.stride() == box.stride()));
   if (all(sext.stride() < dext.stride())) {
@@ -451,132 +545,176 @@ void data<CCTK_REAL8,3>
        srcbbox, dstbbox, regbbox);
     
   } else if (all(sext.stride() > dext.stride())) {
-    
-    switch (order_time) {
-      
-    case 0:
-      assert (srcs.size()>=1);
-      switch (order_space) {
-      case 0:
-      case 1:
-	CCTK_FNAME(prolongate_3d_real8)
-	  ((const CCTK_REAL8*)srcs[0]->storage(),
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(),
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      case 2:
-      case 3:
-	CCTK_FNAME(prolongate_3d_real8_o3)
-	  ((const CCTK_REAL8*)srcs[0]->storage(),
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(),
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      case 4:
-      case 5:
-	CCTK_FNAME(prolongate_3d_real8_o5)
-	  ((const CCTK_REAL8*)srcs[0]->storage(),
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(),
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      default:
-	assert (0);
-      }
-      break;
-      
-    case 1:
-      assert (srcs.size()>=2);
-      switch (order_space) {
-      case 0:
-      case 1:
-	CCTK_FNAME(prolongate_3d_real8_2tl)
-	  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
-	   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(), time,
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      case 2:
-      case 3:
-	CCTK_FNAME(prolongate_3d_real8_2tl_o3)
-	  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
-	   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(), time,
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      case 4:
-      case 5:
-	CCTK_FNAME(prolongate_3d_real8_2tl_o5)
-	  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
-	   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(), time,
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      default:
-	assert (0);
-      }
-      break;
-      
-    case 2:
-      assert (srcs.size()>=3);
-      switch (order_space) {
-      case 0:
-      case 1:
-	CCTK_FNAME(prolongate_3d_real8_3tl)
-	  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
-	   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
-	   (const CCTK_REAL8*)srcs[2]->storage(), times[2],
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(), time,
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      case 2:
-      case 3:
-	CCTK_FNAME(prolongate_3d_real8_3tl_o3)
-	  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
-	   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
-	   (const CCTK_REAL8*)srcs[2]->storage(), times[2],
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(), time,
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
-      case 4:
-      case 5:
-	CCTK_FNAME(prolongate_3d_real8_3tl_o5)
-	  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
-	   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
-	   (const CCTK_REAL8*)srcs[2]->storage(), times[2],
-	   srcshp[0], srcshp[1], srcshp[2],
-	   (CCTK_REAL8*)storage(), time,
-	   dstshp[0], dstshp[1], dstshp[2],
-	   srcbbox, dstbbox, regbbox);
-	break;
+
+    switch (prolong_method) {
+
+      case PROLONG_NONE:  // PROLONG_NONE: NOTHING
+        break;
+      case PROLONG_LAGRANGE:  // PROLONG_LAGRANGE: DEFAULT
+        switch (order_time) {
+          
+          case 0:
+            assert (srcs.size()>=1);
+            switch (order_space) {
+              case 0:
+              case 1:
+                CCTK_FNAME(prolongate_3d_real8)
+                  ((const CCTK_REAL8*)srcs[0]->storage(),
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(),
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              case 2:
+              case 3:
+                CCTK_FNAME(prolongate_3d_real8_o3)
+                  ((const CCTK_REAL8*)srcs[0]->storage(),
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(),
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              case 4:
+              case 5:
+                CCTK_FNAME(prolongate_3d_real8_o5)
+                  ((const CCTK_REAL8*)srcs[0]->storage(),
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(),
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              default:
+                assert (0);
+            }
+            break;
+            
+          case 1:
+            assert (srcs.size()>=2);
+            switch (order_space) {
+              case 0:
+              case 1:
+                CCTK_FNAME(prolongate_3d_real8_2tl)
+                  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+                   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(), time,
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              case 2:
+              case 3:
+                CCTK_FNAME(prolongate_3d_real8_2tl_o3)
+                  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+                   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(), time,
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              case 4:
+              case 5:
+                CCTK_FNAME(prolongate_3d_real8_2tl_o5)
+                  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+                   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(), time,
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              default:
+                assert (0);
+            }
+            break;
+            
+          case 2:
+            assert (srcs.size()>=3);
+            switch (order_space) {
+              case 0:
+              case 1:
+                CCTK_FNAME(prolongate_3d_real8_3tl)
+                  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+                   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+                   (const CCTK_REAL8*)srcs[2]->storage(), times[2],
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(), time,
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              case 2:
+              case 3:
+                CCTK_FNAME(prolongate_3d_real8_3tl_o3)
+                  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+                   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+                   (const CCTK_REAL8*)srcs[2]->storage(), times[2],
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(), time,
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              case 4:
+              case 5:
+                CCTK_FNAME(prolongate_3d_real8_3tl_o5)
+                  ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+                   (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+                   (const CCTK_REAL8*)srcs[2]->storage(), times[2],
+                   srcshp[0], srcshp[1], srcshp[2],
+                   (CCTK_REAL8*)storage(), time,
+                   dstshp[0], dstshp[1], dstshp[2],
+                   srcbbox, dstbbox, regbbox);
+                break;
+              default:
+                assert (0);
+            }
+            break;
+            
+          default:
+            assert (0);
+        } // switch (order_time)
+        break;
+      case PROLONG_TVD: // PROLONG_TVD
+        switch (order_time) {
+          case 0: 
+            CCTK_FNAME(prolongate_3d_real8_minmod)
+              ((const CCTK_REAL8*)srcs[0]->storage(),
+               srcshp[0], srcshp[1], srcshp[2],
+               (CCTK_REAL8*)storage(),
+               dstshp[0], dstshp[1], dstshp[2],
+               srcbbox, dstbbox, regbbox);
+            break;
+         case 1:
+            CCTK_FNAME(prolongate_3d_real8_2tl_minmod)
+              ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+               (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+               srcshp[0], srcshp[1], srcshp[2],
+               (CCTK_REAL8*)storage(), time,
+               dstshp[0], dstshp[1], dstshp[2],
+               srcbbox, dstbbox, regbbox);
+          case 2: 
+            CCTK_FNAME(prolongate_3d_real8_3tl_minmod)
+              ((const CCTK_REAL8*)srcs[0]->storage(), times[0],
+               (const CCTK_REAL8*)srcs[1]->storage(), times[1],
+               (const CCTK_REAL8*)srcs[2]->storage(), times[2],
+               srcshp[0], srcshp[1], srcshp[2],
+               (CCTK_REAL8*)storage(), time,
+               dstshp[0], dstshp[1], dstshp[2],
+               srcbbox, dstbbox, regbbox);
+            break;
+          default:
+            assert (0);
+        }
+        break;
       default:
-	assert (0);
-      }
-      break;
-      
-    default:
-      assert (0);
-    }
+        assert(0);
+    } // switch (prolong_method)
+    
     
   } else {
     assert (0);
   }
 }
 
+  
+
 
 
 // Output
diff --git a/Carpet/CarpetLib/src/gdata.hh b/Carpet/CarpetLib/src/gdata.hh
index 75b12319e..edc08c8e7 100644
--- a/Carpet/CarpetLib/src/gdata.hh
+++ b/Carpet/CarpetLib/src/gdata.hh
@@ -1,4 +1,4 @@
-// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/gdata.hh,v 1.18 2003/10/14 16:39:16 schnetter Exp $
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/gdata.hh,v 1.19 2003/10/15 07:14:01 hawke Exp $
 
 #ifndef GDATA_HH
 #define GDATA_HH
@@ -67,6 +67,10 @@ public:
   virtual void free () = 0;
   
   // Accessors
+
+  int var_index () const {
+    return varindex;
+  }
   
   bool has_storage () const {
     return _has_storage;