aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAnton Khirnov <anton@khirnov.net>2014-12-01 13:18:53 +0100
committerAnton Khirnov <anton@khirnov.net>2014-12-01 13:18:53 +0100
commit67ca1a0ee0419c305dff9151c3f93205c553bd00 (patch)
tree43fdf6a1f0f7bf7b3b7a5788a98b4031deb1bf3d
parentdb77d3ae71ed866ea58ab0bafb9a1d8dba59981c (diff)
Drop commented out cruft.
-rw-r--r--brill_data.c93
1 files changed, 0 insertions, 93 deletions
diff --git a/brill_data.c b/brill_data.c
index 72a44a3..1f6378f 100644
--- a/brill_data.c
+++ b/brill_data.c
@@ -803,96 +803,3 @@ fail:
return ret;
}
-
-#if 0
-void brill_data(CCTK_ARGUMENTS)
-{
- static BrillDataContext *prev_bd;
-
- DECLARE_CCTK_ARGUMENTS;
- DECLARE_CCTK_PARAMETERS;
-
- BrillDataContext *bd;
-
- long double *basis_val_r, *basis_val_z;
-
- int64_t grid_size = CCTK_GFINDEX3D(cctkGH,
- cctk_lsh[0] - 1,
- cctk_lsh[1] - 1,
- cctk_lsh[2] - 1) + 1;
-
- /* on the first run, solve the equation for ψ */
- if (!prev_bd) {
- bd = init_bd(cctkGH, amplitude, eppley_n, basis_order_r, basis_order_z,
- colloc_grid_offset_r, colloc_grid_offset_z,
- scale_factor, overdet);
-
- brill_solve(bd);
-
- if (export_coeffs)
- memcpy(brill_coeffs, bd->psi_coeffs->data, sizeof(*brill_coeffs) * bd->nb_coeffs);
-
- prev_bd = bd;
- } else
- bd = prev_bd;
-
- memset(kxx, 0, sizeof(*kxx) * grid_size);
- memset(kyy, 0, sizeof(*kyy) * grid_size);
- memset(kzz, 0, sizeof(*kzz) * grid_size);
- memset(kxy, 0, sizeof(*kxy) * grid_size);
- memset(kxz, 0, sizeof(*kxz) * grid_size);
- memset(kyz, 0, sizeof(*kyz) * grid_size);
- memset(gxz, 0, sizeof(*kxz) * grid_size);
- memset(gyz, 0, sizeof(*kyz) * grid_size);
-
- /* precompute the basis values on the grid points */
- basis_val_r = malloc(sizeof(*basis_val_r) * bd->nb_coeffs_rho * cctk_lsh[1] * cctk_lsh[0]);
- basis_val_z = malloc(sizeof(*basis_val_z) * bd->nb_coeffs_z * cctk_lsh[2]);
-
- for (int i = 0; i < cctk_lsh[2]; i++) {
- CCTK_REAL zz = z[CCTK_GFINDEX3D(cctkGH, 0, 0, i)];
- for (int j = 0; j < bd->nb_coeffs_z; j++)
- basis_val_z[i * bd->nb_coeffs_z + j] = bd->basis->eval(zz, j, bd->basis_scale_factor_rho);
- }
-
- for (int j = 0; j < cctk_lsh[1]; j++)
- for (int i = 0; i < cctk_lsh[0]; i++) {
- CCTK_REAL xx = x[CCTK_GFINDEX3D(cctkGH, i, j, 0)];
- CCTK_REAL yy = y[CCTK_GFINDEX3D(cctkGH, i, j, 0)];
- CCTK_REAL r = sqrt(SQR(xx) + SQR(yy));
-
- for (int k = 0; k < bd->nb_coeffs_rho; k++)
- basis_val_r[(j * cctk_lsh[0] + i) * bd->nb_coeffs_rho + k] = bd->basis->eval(r, k, bd->basis_scale_factor_rho);
- }
-
-#pragma omp parallel for
- for (int k = 0; k < cctk_lsh[2]; k++)
- for (int j = 0; j < cctk_lsh[1]; j++)
- for (int i = 0; i < cctk_lsh[0]; i++) {
- int index = CCTK_GFINDEX3D(cctkGH, i, j, k);
- long double xx = x[index], yy = y[index], zz = z[index];
-
- long double r2 = SQR(xx) + SQR(yy);
- long double r = sqrtl(r2);
-
- long double q = bd->q_func(bd, r, zz);
- long double e2q = expl(2 * q);
-
- long double psi = 1.0, psi4;
-
- for (int n = 0; n < bd->nb_coeffs_z; n++)
- for (int m = 0; m < bd->nb_coeffs_rho; m++)
- psi += bd->psi_coeffs->data[n * bd->nb_coeffs_z + m] * basis_val_r[(j * cctk_lsh[0] + i) * bd->nb_coeffs_rho + m] * basis_val_z[k * bd->nb_coeffs_z + n];
-
- psi4 = SQR(SQR(psi));
-
- if (r2 < EPS)
- r2 = EPS;
-
- gxx[index] = psi4 * (e2q + (1 - e2q) * SQR(yy) / r2);
- gyy[index] = psi4 * (e2q + (1 - e2q) * SQR(xx) / r2);
- gzz[index] = psi4 * e2q;
- gxy[index] = psi4 * (-(1.0 - e2q) * xx * yy / r2);
- }
-}
-#endif