diff options
author | herrmann <herrmann@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 2003-05-12 09:23:58 +0000 |
---|---|---|
committer | herrmann <herrmann@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 2003-05-12 09:23:58 +0000 |
commit | 7e155396eadfca19daeff52c0459d09e5574c980 (patch) | |
tree | 4891682caf72e98d0c43e6f6d91f8414da64e07b /src/D3_extract_int.F | |
parent | 98a1c9882513adf3daf0d65ba855cf590a5a1bf3 (diff) |
Added private STRING parameter Extract::interpolation_operator defaulting to
"Lagrange polynomial interpolation".
See LocalInterp for other interpolation_operator values.
At the moment LocalInterp will use
"Lagrange polynomial interpolation (tensor product)".
Changed the testsuites to use
"Lagrange polynomial interpolation (maximum degree)"
to get the same behaviour as before and the testsuites to pass again.
They didn't pass before, because the interpolator gives slightly different
answers for the "tensor product" version. The Moncrief Q's should be
zero for Schwarzschild. The result is for example 10^-6 rather than
10^-5 which was reported as a fail.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Extract/trunk@77 5301f0c2-dbc4-4cee-b2f5-8d7afba4d129
Diffstat (limited to 'src/D3_extract_int.F')
-rw-r--r-- | src/D3_extract_int.F | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/src/D3_extract_int.F b/src/D3_extract_int.F index 07d7dab..ce14b37 100644 --- a/src/D3_extract_int.F +++ b/src/D3_extract_int.F @@ -8,7 +8,7 @@ c ------------------------------------------------------------------ INTERFACE SUBROUTINE D3_extract(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin,igrid, - & origin,myproc,interpolation_order, + & origin,myproc,interpolation_operator,interpolation_order, & Nt,Np,all_modes,l,m,x,y,z,Dx,Dy,Dz,Psi_power,Psi, & g00,gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz, & eta,ADMmass,momentum,spin,mass,rsch,Qodd,Qeven, @@ -36,6 +36,9 @@ c ------------------------------------------------------------------ CCTK_REAL,INTENT(OUT) :: & ADMmass(2),mass,rsch,Qodd(:,:,:),Qeven(:,:,:),dtaudt, & momentum(3),spin(3) + CCTK_STRING,INTENT(IN) :: + & interpolation_operator + END SUBROUTINE |