diff options
Diffstat (limited to 'src/qms.c')
-rw-r--r-- | src/qms.c | 36 |
1 files changed, 25 insertions, 11 deletions
@@ -232,7 +232,7 @@ static int context_init(cGH *cctkGH) qms->max_radius = 60.0; ret = qms_solver_init(&qms->solver, cctkGH, basis_order_r, basis_order_z, - outer_bound, filter_power, 0.0); + outer_bound, filter_power, 0.0, ccz4); if (ret < 0) return ret; @@ -460,16 +460,30 @@ void qms_init(CCTK_ARGUMENTS) if (!qms_context) context_init(cctkGH); - double *Kdot11 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot11"); - double *Kdot22 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot22"); - double *Kdot33 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot33"); - double *Kdot12 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot12"); - double *Kdot13 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot13"); - double *Kdot23 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot23"); - - double *Xtdot1 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Xtdot1"); - double *Xtdot2 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Xtdot2"); - double *Xtdot3 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Xtdot3"); + double *Kdot11, *Kdot22, *Kdot33, *Kdot12, *Kdot13, *Kdot23, *Xtdot1, *Xtdot2, *Xtdot3; + if (ccz4) { + Kdot11 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Kdot11"); + Kdot22 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Kdot22"); + Kdot33 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Kdot33"); + Kdot12 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Kdot12"); + Kdot13 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Kdot13"); + Kdot23 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Kdot23"); + + Xtdot1 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Xtdot1"); + Xtdot2 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Xtdot2"); + Xtdot3 = CCTK_VarDataPtr(cctkGH, 0, "ML_CCZ4::Xtdot3"); + } else { + Kdot11 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot11"); + Kdot22 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot22"); + Kdot33 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot33"); + Kdot12 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot12"); + Kdot13 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot13"); + Kdot23 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Kdot23"); + + Xtdot1 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Xtdot1"); + Xtdot2 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Xtdot2"); + Xtdot3 = CCTK_VarDataPtr(cctkGH, 0, "ML_BSSN::Xtdot3"); + } for (int k = 0; k < cctk_lsh[2]; k++) for (int j = 0; j < cctk_lsh[1]; j++) |