Use CarpetInterp2 if it is available, since this is more efficient

than CarpetInterp.


git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/PunctureTracker/trunk@11 a2659f00-0f4f-0410-9214-a4596bbb8a4f

1 files changed, 22 insertions, 8 deletions

diff --git a/src/puncture_tracker.c b/src/puncture_tracker.c
index 6f2773c..4e36c1f 100644
--- a/src/puncture_tracker.c
+++ b/src/puncture_tracker.c
@@ -105,6 +105,7 @@ PunctureTracker_Track (CCTK_ARGUMENTS)
   }
   
   // Interpolation parameter table
+  int const order = 4;
   int const param_table_handle = Util_TableCreateFromString ("order=4");
   if (param_table_handle < 0) {
     CCTK_WARN (CCTK_WARN_ALERT, "Can't create parameter table");
@@ -159,14 +160,27 @@ PunctureTracker_Track (CCTK_ARGUMENTS)
     output_arrays[2] = pt_betaz;
     
     // Interpolate
-    int const ierr = CCTK_InterpGridArrays
-      (cctkGH, dim,
-       operator_handle, param_table_handle, coordsys_handle,
-       num_points,
-       CCTK_VARIABLE_REAL,
-       interp_coords,
-       num_vars, input_array_indices,
-       num_vars, output_array_type_codes, output_arrays);
+    int ierr;
+    if (CCTK_IsFunctionAliased ("InterpGridArrays")) {
+      // TODO: use correct array types
+      // (CCTK_POINTER[] vs. CCTK_REAL[])
+      ierr = InterpGridArrays
+        (cctkGH, dim,
+         order,
+         num_points,
+         interp_coords,
+         num_vars, input_array_indices,
+         num_vars, output_arrays);
+    } else {
+      ierr = CCTK_InterpGridArrays
+        (cctkGH, dim,
+         operator_handle, param_table_handle, coordsys_handle,
+         num_points,
+         CCTK_VARIABLE_REAL,
+         interp_coords,
+         num_vars, input_array_indices,
+         num_vars, output_array_type_codes, output_arrays);
+    }
     if (ierr < 0) {
       CCTK_WARN (CCTK_WARN_ALERT, "Interpolation error");
       goto label_free_param_table;