A code fragment for summation by parts. Experimental.

1 files changed, 25 insertions, 0 deletions

diff --git a/Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/sbp_calc_coeffs.h b/Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/sbp_calc_coeffs.h
new file mode 100644
index 0000000..092d36f
--- /dev/null
+++ b/Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/sbp_calc_coeffs.h
@@ -0,0 +1,25 @@
+int const deriv_options = -1;
+  
+/* (drho/dx)^i = D^i_j rho^j, with D^i_j = q(i,j) */
+/* imin and imax contain the first and last non-zero entry of the
+   two-dimensional array q */
+/* All indexing is according to Fortran conventions */
+  
+int const ni = cctk_lsh[0];
+int const nj = cctk_lsh[1];
+int const nk = cctk_lsh[2];
+  
+CCTK_INT  * restrict const sbp_imin = malloc (ni * sizeof *sbp_imin);
+CCTK_INT  * restrict const sbp_imax = malloc (ni * sizeof *sbp_imax);
+CCTK_REAL * restrict const qx   = malloc (ni * ni * sizeof *qx);
+Diff_coeff (cctkGH, 0, ni, sbp_imin, sbp_imax, qx, deriv_options);
+  
+CCTK_INT  * restrict const sbp_jmin = malloc (nj * sizeof *sbp_jmin);
+CCTK_INT  * restrict const sbp_jmax = malloc (nj * sizeof *sbp_jmax);
+CCTK_REAL * restrict const qy   = malloc (nj * nj * sizeof *qy);
+Diff_coeff (cctkGH, 1, nj, sbp_jmin, sbp_jmax, qy, deriv_options);
+  
+CCTK_INT  * restrict const sbp_kmin = malloc (nk * sizeof *sbp_kmin);
+CCTK_INT  * restrict const sbp_kmax = malloc (nk * sizeof *sbp_kmax);
+CCTK_REAL * restrict const qz   = malloc (nk * nk * sizeof *qz);
+Diff_coeff (cctkGH, 2, nk, sbp_kmin, sbp_kmax, qz, deriv_options);