aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authortradke <tradke@89daf98e-ef62-4674-b946-b8ff9de2216c>2003-05-06 08:13:10 +0000
committertradke <tradke@89daf98e-ef62-4674-b946-b8ff9de2216c>2003-05-06 08:13:10 +0000
commitb96fccc02e13fc603395c61d6cd44783268e3caf (patch)
treed4c303a796bf2e4dd43bd9d6cf3d073a15702401
parentac8c40a006b522c390906b67616af535a42c0824 (diff)
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
-rw-r--r--src/AHFinder_dis.F5
-rw-r--r--src/AHFinder_gau.F7
-rw-r--r--src/AHFinder_int.F6
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