From 67ca1a0ee0419c305dff9151c3f93205c553bd00 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Mon, 1 Dec 2014 13:18:53 +0100 Subject: Drop commented out cruft. --- brill_data.c | 93 ------------------------------------------------------------ 1 file changed, 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 -- cgit v1.2.3