From b96fccc02e13fc603395c61d6cd44783268e3caf Mon Sep 17 00:00:00 2001 From: tradke Date: Tue, 6 May 2003 08:13:10 +0000 Subject: Get the name for the interpolation operator to use from the new parameter AHFinder::interpolation_operator. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@334 89daf98e-ef62-4674-b946-b8ff9de2216c --- src/AHFinder_dis.F | 5 ++++- src/AHFinder_gau.F | 7 +++++-- src/AHFinder_int.F | 6 ++++-- 3 files changed, 13 insertions(+), 5 deletions(-) diff --git a/src/AHFinder_dis.F b/src/AHFinder_dis.F index e44be98..44faf91 100644 --- a/src/AHFinder_dis.F +++ b/src/AHFinder_dis.F @@ -25,6 +25,8 @@ integer np,npoints,ierror integer param_table_handle,interp_handle,coord_system_handle character(30) options_string + character(128) operator + CCTK_INT nchars CCTK_POINTER, dimension(3) :: interp_coords CCTK_INT, dimension(6) :: in_array_indices @@ -495,7 +497,8 @@ call CCTK_WARN(0,"Cannot create parameter table for interpolator") endif - call CCTK_InterpHandle (interp_handle,"Lagrange polynomial interpolation") + call CCTK_FortranString (nchars, interpolation_operator, operator) + call CCTK_InterpHandle (interp_handle, operator) if (interp_handle.lt.0) then call CCTK_WARN(0,"Cannot get handle for interpolation ! Forgot to activate an implementation providing interpolation operators ??") endif diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F index ad79044..42a8acb 100644 --- a/src/AHFinder_gau.F +++ b/src/AHFinder_gau.F @@ -34,6 +34,8 @@ integer npoints integer param_table_handle,interp_handle,coord_system_handle,sum_handle integer ierror + character(len=128) :: operator + CCTK_INT nchars CCTK_INT rerror,error1,error2 @@ -67,7 +69,7 @@ CCTK_REAL :: m1p, m2p, x1p, x2p, y1p, y2p, z1p, z2p CCTK_REAL :: r12, r22, sigma CCTK_REAL, dimension(nx,ny,nz) :: fac - CCTK_INT :: use_att, apower, nchars + CCTK_INT :: use_att, apower logical :: use_rot_att logical :: str_comp logical :: gaussf_exists @@ -685,7 +687,8 @@ call CCTK_WARN(0,"Cannot create parameter table for interpolator") endif - call CCTK_InterpHandle (interp_handle,"Lagrange polynomial interpolation") + call CCTK_FortranString (nchars, interpolation_operator, operator) + call CCTK_InterpHandle (interp_handle, operator) if (interp_handle .lt. 0) then call CCTK_WARN(0,"Cannot get handle for interpolation ! Forgot to activate an implementation providing interpolation operators ??") endif diff --git a/src/AHFinder_int.F b/src/AHFinder_int.F index 060834a..e3c9833 100644 --- a/src/AHFinder_int.F +++ b/src/AHFinder_int.F @@ -37,7 +37,8 @@ integer interp_handle,coord_system_handle integer reduction_handle,param_table_handle character(30) options_string - CCTK_INT red_tmp + character(128) operator + CCTK_INT red_tmp, nchars CCTK_POINTER, dimension(3) :: interp_coords CCTK_INT, dimension(9) :: in_array_indices @@ -624,7 +625,8 @@ call CCTK_WARN(0,"Cannot create parameter table for interpolator") endif - call CCTK_InterpHandle (interp_handle,"Lagrange polynomial interpolation") + call CCTK_FortranString (nchars, interpolation_operator, operator) + call CCTK_InterpHandle (interp_handle, operator) if (interp_handle .lt. 0) then call CCTK_WARN(0,"Cannot get handle for interpolation ! Forgot to activate an implementation providing interpolation operators ??") endif -- cgit v1.2.3