CarpetLib: Reorganise prolongation and restriction operators

Reorganise prolongation and restriction operators.  This is a major
implementation change.

Most operators are now written as C++ templates instead of as Fortran
77 code.  This simplifies the code, since C++ routines can be called
more easily, and they also have access to CarpetLib's high-level data
structures.

Previously, the operators combined temporal and spatial interpolation.
Now, time interpolation and space interpolation are handled
separately.  This may be less efficient, but simplifies the code
significantly, since there are now N+M instead of N*M routines, for N
time interpolation and M space interpolation methods.

Remove the minmod prolongation operator, which was previously
disabled.

Add support for cell centering, using a method described by Simon
Hern, and suggested for Carpet by Ian Hawke.

darcs-hash:20070112205812-dae7b-5329795aa698e7bbc3671b1504134885dd830238.gz

1 files changed, 1 insertions, 1 deletions

diff --git a/Carpet/CarpetLib/src/defs.hh b/Carpet/CarpetLib/src/defs.hh
index 21272b58e..3e2bc1734 100644
--- a/Carpet/CarpetLib/src/defs.hh
+++ b/Carpet/CarpetLib/src/defs.hh
@@ -53,7 +53,7 @@ typedef vect<vect<int,dim>,2>  i2vect;
 
 
 // A general type
-enum centering { vertex_centered, cell_centered };
+enum centering { error_centered, vertex_centered, cell_centered };