Removed old interpolator API code.

git-svn-id: http://svn.cactuscode.org/arrangements/CactusPUGH/PUGHInterp/trunk@79 1c20744c-e24a-42ec-9533-f5004cb800e5

1 files changed, 2 insertions, 196 deletions

diff --git a/src/Startup.c b/src/Startup.c
index 69be436..f1849e0 100644
--- a/src/Startup.c
+++ b/src/Startup.c
@@ -32,48 +32,6 @@ void PUGHInterp_Startup (void);
 #ifdef CCTK_MPI
 static void *SetupGH (tFleshConfig *config, int convergence_level, cGH *GH);
 #endif
-static int InterpGV_1stOrder (cGH *GH,
-                              const char *coord_system,
-                              int num_points,
-                              int num_in_array_indices,
-                              int num_out_arrays,
-                              const void *const interp_coord_arrays[],
-                              const int interp_coord_array_types[],
-                              const int in_array_indices[],
-                              void *const out_arrays[],
-                              const int out_array_types[]);
-static int InterpGV_2ndOrder (cGH *GH,
-                              const char *coord_system,
-                              int num_points,
-                              int num_in_array_indices,
-                              int num_out_arrays,
-                              const void *const interp_coord_arrays[],
-                              const int interp_coord_array_types[],
-                              const int in_array_indices[],
-                              void *const out_arrays[],
-                              const int out_array_types[]);
-static int InterpGV_3rdOrder (cGH *GH,
-                              const char *coord_system,
-                              int num_points,
-                              int num_in_array_indices,
-                              int num_out_arrays,
-                              const void *const interp_coord_arrays[],
-                              const int interp_coord_array_types[],
-                              const int in_array_indices[],
-                              void *const out_arrays[],
-                              const int out_array_types[]);
-static int PUGHInterp_InterpGV (cGH *GH,
-                                int order,
-                                const char *local_interpolator,
-                                const char *coord_system_name,
-                                int N_points,
-                                int N_input_arrays,
-                                int N_output_arrays,
-                                const void *const interp_coord_arrays[],
-                                const int interp_coord_array_types[],
-                                const int input_array_indices[],
-                                void *const output_arrays[],
-                                const int output_array_types[]);
 
 
  /*@@
@@ -82,12 +40,11 @@ static int PUGHInterp_InterpGV (cGH *GH,
    @author    Thomas Radke
    @desc
               The startup registration routine for PUGHInterp.
-              Registers the PUGHInterp GH extensions setup routine
-              and the interpolation operators with the flesh.
+              Registers the PUGHInterp GH extensions setup routine with the
+              flesh and overloads the CCTK_InterpGridArrays() routine.
    @enddesc
    @calls     CCTK_OverloadInterpGridArrays
               CCTK_RegisterGHExtensionSetupGH
-              CCTK_InterpRegisterOperatorGV
 @@*/
 void PUGHInterp_Startup (void)
 {
@@ -97,13 +54,6 @@ void PUGHInterp_Startup (void)
   CCTK_RegisterGHExtensionSetupGH (CCTK_RegisterGHExtension ("PUGHInterp"),
                                    SetupGH);
 #endif
-
-  CCTK_InterpRegisterOperatorGV (InterpGV_1stOrder,
-                                 "first-order uniform cartesian");
-  CCTK_InterpRegisterOperatorGV (InterpGV_2ndOrder,
-                                 "second-order uniform cartesian");
-  CCTK_InterpRegisterOperatorGV (InterpGV_3rdOrder,
-                                 "third-order uniform cartesian");
 }
 
 
@@ -152,147 +102,3 @@ static void *SetupGH (tFleshConfig *config, int convergence_level, cGH *GH)
   return (myGH);
 }
 #endif /* CCTK_MPI */
-
-
- /*@@
-   @routine   InterpGV_NthOrder
-   @date      Wed 14 Feb 2001
-   @author    Thomas Radke
-   @desc
-              Wrappers for the different interpolation operators
-              registered for first/second/third order interpolation.
-              These wrappers just call the common interpolation routine
-              passing all arguments plus the interpolation order.
-   @enddesc
-
-   @returntype int
-   @returndesc
-               the return code of the common interpolation routine
-   @endreturndesc
-@@*/
-static int InterpGV_1stOrder (cGH *GH,
-                              const char *coord_system,
-                              int num_points,
-                              int num_in_array_indices,
-                              int num_out_arrays,
-                              const void *const interp_coord_arrays[],
-                              const int interp_coord_array_types[],
-                              const int in_array_indices[],
-                              void *const out_arrays[],
-                              const int out_array_types[])
-{
-  return (PUGHInterp_InterpGV (GH, 1, "first-order uniform cartesian",
-                               coord_system, num_points,
-                               num_in_array_indices, num_out_arrays,
-                               interp_coord_arrays, interp_coord_array_types,
-                               in_array_indices, out_arrays, out_array_types));
-}
-
-
-static int InterpGV_2ndOrder (cGH *GH,
-                              const char *coord_system,
-                              int num_points,
-                              int num_in_array_indices,
-                              int num_out_arrays,
-                              const void *const interp_coord_arrays[],
-                              const int interp_coord_array_types[],
-                              const int in_array_indices[],
-                              void *const out_arrays[],
-                              const int out_array_types[])
-{
-  return (PUGHInterp_InterpGV (GH, 2, "second-order uniform cartesian",
-                               coord_system, num_points,
-                               num_in_array_indices, num_out_arrays,
-                               interp_coord_arrays, interp_coord_array_types,
-                               in_array_indices, out_arrays, out_array_types));
-}
-
-
-static int InterpGV_3rdOrder (cGH *GH,
-                              const char *coord_system,
-                              int num_points,
-                              int num_in_array_indices,
-                              int num_out_arrays,
-                              const void *const interp_coord_arrays[],
-                              const int interp_coord_array_types[],
-                              const int in_array_indices[],
-                              void *const out_arrays[],
-                              const int out_array_types[])
-{
-  return (PUGHInterp_InterpGV (GH, 3, "third-order uniform cartesian",
-                               coord_system, num_points,
-                               num_in_array_indices, num_out_arrays,
-                               interp_coord_arrays, interp_coord_array_types,
-                               in_array_indices, out_arrays, out_array_types));
-}
-
-
-static int PUGHInterp_InterpGV (cGH *GH,
-                                int order,
-                                const char *local_interpolator,
-                                const char *coord_system_name,
-                                int N_points,
-                                int N_input_arrays,
-                                int N_output_arrays,
-                                const void *const interp_coord_arrays[],
-                                const int interp_coord_array_types[],
-                                const int input_array_indices[],
-                                void *const output_arrays[],
-                                const int output_array_types[])
-{
-  int i, N_dims, retval;
-  int local_interp_handle, param_table_handle, coord_system_handle;
-  char table_string[64];
-  CCTK_INT *_input_array_indices, *_output_array_types;
-  union
-  {
-    const cGH *const_GH;
-          cGH *non_const_GH;
-  } _GH;
-
-
-  _GH.non_const_GH = GH;
-
-  /* get dimensionality of the coordinate system */
-  N_dims = CCTK_CoordSystemDim (coord_system_name);
-  if (N_dims <= 0)
-  {
-    CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING,
-                "Cannot get dimensions of coordinate system '%s'",
-                coord_system_name);
-    return (-1);
-  }
-  coord_system_handle = CCTK_CoordSystemHandle (coord_system_name);
-  local_interp_handle = CCTK_InterpHandle (local_interpolator);
-
-  /* turn native datatype arrays into CCTK_ datatypes */
-  _input_array_indices = malloc (2 * N_dims * sizeof (CCTK_INT));
-  _output_array_types = _input_array_indices + N_dims;
-  for (i = 0; i < N_dims; i++)
-  {
-    _input_array_indices[i] = input_array_indices[i];
-    _output_array_types[i] = output_array_types[i];
-  }
-
-  /* set up the parameter table (only supported option is "order") */
-  sprintf (table_string, "order = %d", order);
-  param_table_handle = Util_TableCreateFromString (table_string);
-
-  retval = PUGHInterp_InterpGridArrays (_GH.const_GH, N_dims,
-                                        local_interp_handle, param_table_handle,
-                                        coord_system_handle,
-                                        N_points,
-                                        interp_coord_array_types[0],
-                                        interp_coord_arrays,
-                                        N_input_arrays,
-                                        _input_array_indices,
-                                        N_output_arrays,
-                                        _output_array_types,
-                                        output_arrays);
-
-  /* release resources */
-  Util_TableDestroy (param_table_handle);
-  free (_input_array_indices);
-
-  return (retval);
-}