Update thorn ADMAnalysis to make it work better with mesh refinement.

We applied similar changes to ADMConstraints some time ago.

Set prolongation types of all grid functions.

Allow the Ricci tensor to have 3 time levels, so that its boundaries
can be interpolated.

Add parameters ricci_persist, ricci_timelevels, and
ricci_prolongation_type to specify how the boundaries are treated.
The default is to do nothing on mesh refinement boundaries, keeping
the old behaviour.

Select Ricci tensor boundary conditions in level mode in a new
routine.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMAnalysis/trunk@52 80bd93c7-81bc-453a-9e3f-619c0b4f6fe4

1 files changed, 11 insertions, 4 deletions

diff --git a/src/Ricci.c b/src/Ricci.c
index adb4693..04f8b18 100644
--- a/src/Ricci.c
+++ b/src/Ricci.c
@@ -68,7 +68,7 @@ void ADMAnalysis_Ricci(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
 
   CCTK_INT lsh[3];
-  CCTK_INT i,j,k, d, ijk, di, dj, dk, err;
+  CCTK_INT i,j,k, d, ijk, di, dj, dk;
 
 #include "CactusEinstein/ADMMacros/src/macro/RICCI_declare.h"
 
@@ -104,6 +104,16 @@ void ADMAnalysis_Ricci(CCTK_ARGUMENTS)
                     Ricci11, Ricci12, Ricci13, Ricci22, Ricci23, Ricci33,
                     Ricci, detg);
 
+  return;
+  
+}
+
+void ADMAnalysis_Ricci_Boundaries(CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS;
+
+  CCTK_INT err;
+
   /* Apply Flat Boundary Condition */
   err = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, 
                                   "ADMAnalysis::ricci_scalar", "Flat");
@@ -118,9 +128,6 @@ void ADMAnalysis_Ricci(CCTK_ARGUMENTS)
    *  value of the Ricci tensor is going to be used for something other
    *  than output, then the boundaries should be handled more properly.
    */
-
-  return;
-  
 }
 
 /********************************************************************