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/Extract.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/Extract.F')
| -rw-r--r-- | src/Extract.F | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/src/Extract.F b/src/Extract.F index 1ea5e86..ace99fb 100644 --- a/src/Extract.F +++ b/src/Extract.F @@ -378,7 +378,8 @@ c Do extraction at each radius extract_at_each_radius: DO WHILE (radius < r2) CALL D3_extract(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin, - & igrid,orig,myproc,interpolation_order,Nt,Np,all_modes,lmode, + & igrid,orig,myproc,interpolation_operator,interpolation_order, + & Nt,Np,all_modes,lmode, & mmode,x_1d,y_1d,z_1d,Dx,Dy,Dz,Psi_power,Psi,g00, & gxx,gxy,gxz,gyy,gyz,gzz,kxx,kxy,kxz,kyy,kyz,kzz, & radius,ADMmass,momentum,spin,mass,rsch, @@ -549,7 +550,8 @@ c Cannot use the conformal equation for ADM mass now END IF CALL D3_extract(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin, - & igrid,orig,myproc,interpolation_order,Nt,Np,all_modes, + & igrid,orig,myproc,interpolation_operator,interpolation_order, + & Nt,Np,all_modes, & lmode,mmode,x_1d,y_1d,z_1d, & Dx,Dy,Dz,Psi_power,Psi,g00,gxx,gxy,gxz,gyy,gyz,gzz,kxx,kxy,kxz, & kyy,kyz,kzz,radius,ADMmass,momentum,spin,mass,rsch, |