diff options
Diffstat (limited to 'src/call_derivs.c')
-rw-r--r-- | src/call_derivs.c | 32 |
1 files changed, 20 insertions, 12 deletions
diff --git a/src/call_derivs.c b/src/call_derivs.c index 992e028..aeab3d7 100644 --- a/src/call_derivs.c +++ b/src/call_derivs.c @@ -8,7 +8,6 @@ #include "stencil.h" - void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_REAL *var, CCTK_REAL *dvar, const CCTK_INT table_handle ) @@ -33,6 +32,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_4_2)(const CCTK_REAL *var, @@ -42,6 +42,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_4_3)(const CCTK_REAL *var, @@ -51,6 +52,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_4_3_opt)(const CCTK_REAL *var, @@ -60,6 +62,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_6_3)(const CCTK_REAL *var, @@ -69,6 +72,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_6_3_opt)(const CCTK_REAL *var, @@ -78,6 +82,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_6_5)(const CCTK_REAL *var, @@ -87,6 +92,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_6_5_opt)(const CCTK_REAL *var, @@ -96,6 +102,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_8_4)(const CCTK_REAL *var, @@ -105,6 +112,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_INT *dir, const CCTK_INT *bb, const CCTK_INT *gsize, + const CCTK_INT *offset, const CCTK_REAL *delta, CCTK_REAL *dvar); void CCTK_FCALL CCTK_FNAME(deriv_gf_8_4_opt)(const CCTK_REAL *var, @@ -225,19 +233,19 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, if ( CCTK_Equals(operator_type,"Minimal Bandwidth") ) { switch(loc_order) { case 2: { - CCTK_FNAME(deriv_gf_2_1)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_2_1)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 4: { - CCTK_FNAME(deriv_gf_4_2)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_4_2)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 6: { - CCTK_FNAME(deriv_gf_6_3)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_6_3)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 8: { - CCTK_FNAME(deriv_gf_8_4)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_8_4)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } default: @@ -246,15 +254,15 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, } else { switch(loc_order) { case 2: { - CCTK_FNAME(deriv_gf_2_1)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_2_1)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 4: { - CCTK_FNAME(deriv_gf_4_2)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_4_2)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 6: { - CCTK_FNAME(deriv_gf_6_3_opt)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_6_3_opt)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 8: { @@ -269,11 +277,11 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, if ( CCTK_Equals(operator_type,"Minimal Bandwidth") ) { switch(loc_order) { case 4: { - CCTK_FNAME(deriv_gf_4_3)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_4_3)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 6: { - CCTK_FNAME(deriv_gf_6_5)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_6_5)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } default: @@ -282,11 +290,11 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, } else { switch(loc_order) { case 4: { - CCTK_FNAME(deriv_gf_4_3_opt)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_4_3_opt)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } case 6: { - CCTK_FNAME(deriv_gf_6_5_opt)(var,&ni,&nj,&nk,&dir,bb,&gsize,&delta,dvar); + CCTK_FNAME(deriv_gf_6_5_opt)(var,&ni,&nj,&nk,&dir,bb,&gsize,loc_offset,&delta,dvar); break; } default: |