CCTK_INT FUNCTION DriverInterpolate(
  CCTK_POINTER_TO_CONST IN cctkGH,
  CCTK_INT IN N_dims,
  CCTK_INT IN local_interp_handle,
  CCTK_INT IN param_table_handle,
  CCTK_INT IN coord_system_handle,
  CCTK_INT IN N_interp_points,
  CCTK_INT IN interp_coords_type,
  CCTK_POINTER_TO_CONST ARRAY IN interp_coords,
  CCTK_INT IN N_input_arrays,
  CCTK_INT ARRAY IN input_array_indices,
  CCTK_INT IN N_output_arrays,
  CCTK_INT ARRAY IN output_array_types,
  CCTK_POINTER ARRAY IN output_arrays)
REQUIRES FUNCTION DriverInterpolate

#include "util_Table.h"

  vector<vector<CCTK_REAL> > results2(nvars);

  vector<CCTK_REAL *> resultptrs2(nvars);

    results2.at(n).resize(npoints);
    resultptrs2.at(n) = results2.at(n).data();

  CCTK_INT const N_dims = 3;
  CCTK_INT const N_input = all_varinds.size();
  const void *interp_coords[N_dims];
  interp_coords[0] = coords[0].data();
  interp_coords[1] = coords[1].data();
  interp_coords[2] = coords[2].data();
  void * const *output_array = (void * const *)&resultptrs2[0];
  /* DriverInterpolate arguments that aren't currently used */
  CCTK_INT const coord_system_handle = 0;
  CCTK_INT const interp_coords_type_code = 0;
  CCTK_INT const output_array_type_codes[1] = {0};
  int interp_handle = 0;
  /* Table generation */
  
  int param_table_handle;
  
  param_table_handle = Util_TableCreate(UTIL_TABLE_FLAGS_DEFAULT);
  if (param_table_handle < 0)
          CCTK_ERROR("Can't create parameter table!");
  if (Util_TableSetInt(param_table_handle, 1, "order") < 0)
          CCTK_ERROR("Can't set order in parameter table!");
  if (Util_TableSetIntArray(param_table_handle,
                            nvars, varinds.data(),
                            "operand_indices") < 0)
          CCTK_ERROR("Can't set operand_indices array in parameter table!");
  if (Util_TableSetIntArray(param_table_handle,
                            nvars, operations.data(),
                            "operation_codes") < 0)
          CCTK_ERROR("Can't set operation_codes array in parameter table!");
  int ierr = DriverInterpolate(cctkGH, N_dims, interp_handle, param_table_handle,
                               coord_system_handle, npoints, interp_coords_type_code,
                               interp_coords, N_input, all_varinds.data(), 
                               nvars, output_array_type_codes, output_array);

        assert(chop(results2.at(v).at(n) - coords[d].at(n)) == 0);
        assert(chop(results.at(v).at(n) - results2.at(v).at(n)) == 0);

        assert(chop(results2.at(v).at(n) - (op - 1 == d ? 1 : 0)) == 0);
        assert(chop(results2.at(v).at(n) - results.at(v).at(n)) == 0);

        assert(chop(results2.at(v).at(n)) == 0);