diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/integrate.cc | 2 | ||||
| -rw-r--r-- | src/integrate.hh | 2 | ||||
| -rw-r--r-- | src/tests.cc | 1 |
3 files changed, 4 insertions, 1 deletions
diff --git a/src/integrate.cc b/src/integrate.cc index 90a27d9..1f64c29 100644 --- a/src/integrate.cc +++ b/src/integrate.cc @@ -48,7 +48,7 @@ CCTK_REAL Midpoint2DIntegral(CCTK_REAL const *f, int nx, int ny, CCTK_REAL hx, C return hx * hy * integrand_sum; } -static CCTK_REAL Trapezoidal2DIntegral(CCTK_REAL const *f, int nx, int ny, CCTK_REAL hx, CCTK_REAL hy) +CCTK_REAL Trapezoidal2DIntegral(CCTK_REAL const *f, int nx, int ny, CCTK_REAL hx, CCTK_REAL hy) { CCTK_REAL integrand_sum = 0.0; int ix = 0, iy = 0; diff --git a/src/integrate.hh b/src/integrate.hh index 168fd2f..85e2abb 100644 --- a/src/integrate.hh +++ b/src/integrate.hh @@ -6,6 +6,8 @@ CCTK_REAL Midpoint2DIntegral(CCTK_REAL const *f, int nx, int ny, CCTK_REAL hx, CCTK_REAL hy); +CCTK_REAL Trapezoidal2DIntegral(CCTK_REAL const *f, int nx, int ny, CCTK_REAL hx, CCTK_REAL hy); + CCTK_REAL Simpson2DIntegral(CCTK_REAL const *f, int nx, int ny, CCTK_REAL hx, CCTK_REAL hy); diff --git a/src/tests.cc b/src/tests.cc index ff5cc0a..09562c2 100644 --- a/src/tests.cc +++ b/src/tests.cc @@ -57,6 +57,7 @@ void Multipole_TestIntegrationConvergence(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; *test_simpson_convergence_order = integration_convergence_order(&Simpson2DIntegral); + *test_trapezoidal_convergence_order = integration_convergence_order(&Trapezoidal2DIntegral); } // void Multipole_TestIntegrate(CCTK_ARGUMENTS) |