Switch to using new version of Lagrange polynomial interpolation

which has separate tensor-product and maximum-degree interpolation
operators, with the former now the default.

This should fix the problem Yosef Zlowcher (sp?) found with
the interpolant not being continuous in multiple dimensions.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@145 df1f8a13-aa1d-4dd4-9681-27ded5b42416

1 files changed, 36 insertions, 18 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/InterpLocalUniform.h b/src/GeneralizedPolynomial-Uniform/InterpLocalUniform.h
index 8d9c551..d4821c4 100644
--- a/src/GeneralizedPolynomial-Uniform/InterpLocalUniform.h
+++ b/src/GeneralizedPolynomial-Uniform/InterpLocalUniform.h
@@ -192,24 +192,42 @@ struct	Jacobian_info
  */
 
 /* functions in InterpLocalUniform.c */
-int LocalInterp_ILU_Lagrange(int N_dims,
-			     int param_table_handle,
-			     /***** coordinate system *****/
-			     const CCTK_REAL coord_origin[],
-			     const CCTK_REAL coord_delta[],
-			     /***** interpolation points *****/
-			     int N_interp_points,
-			     int interp_coords_type_code,
-			     const void *const interp_coords[],
-			     /***** input arrays *****/
-			     int N_input_arrays,
-			     const CCTK_INT input_array_dims[],
-			     const CCTK_INT input_array_type_codes[],
-			     const void *const input_arrays[],
-			     /***** output arrays *****/
-			     int N_output_arrays,
-			     const CCTK_INT output_array_type_codes[],
-			     void *const output_arrays[]);
+int LocalInterp_ILU_Lagrange_TP(int N_dims,
+				int param_table_handle,
+				/***** coordinate system *****/
+				const CCTK_REAL coord_origin[],
+				const CCTK_REAL coord_delta[],
+				/***** interpolation points *****/
+				int N_interp_points,
+				int interp_coords_type_code,
+				const void *const interp_coords[],
+				/***** input arrays *****/
+				int N_input_arrays,
+				const CCTK_INT input_array_dims[],
+				const CCTK_INT input_array_type_codes[],
+				const void *const input_arrays[],
+				/***** output arrays *****/
+				int N_output_arrays,
+				const CCTK_INT output_array_type_codes[],
+				void *const output_arrays[]);
+int LocalInterp_ILU_Lagrange_MD(int N_dims,
+				int param_table_handle,
+				/***** coordinate system *****/
+				const CCTK_REAL coord_origin[],
+				const CCTK_REAL coord_delta[],
+				/***** interpolation points *****/
+				int N_interp_points,
+				int interp_coords_type_code,
+				const void *const interp_coords[],
+				/***** input arrays *****/
+				int N_input_arrays,
+				const CCTK_INT input_array_dims[],
+				const CCTK_INT input_array_type_codes[],
+				const void *const input_arrays[],
+				/***** output arrays *****/
+				int N_output_arrays,
+				const CCTK_INT output_array_type_codes[],
+				void *const output_arrays[]);
 int LocalInterp_ILU_Hermite(int N_dims,
 			    int param_table_handle,
 			    /***** coordinate system *****/