∅:D[
6.40969] → [
7.34480:34504]
B:BD[
7.34480] → [
7.34480:34504]
B:BD[
7.34504] → [
2.13194:13273]
∅:D[
2.13273] → [
7.34576:34600]
∅:D[
6.41051] → [
7.34576:34600]
B:BD[
7.34576] → [
7.34576:34600]
B:BD[
7.34600] → [
2.13274:13353]
∅:D[
2.13353] → [
7.34672:34696]
∅:D[
6.41133] → [
7.34672:34696]
B:BD[
7.34672] → [
7.34672:34696]
B:BD[
7.34696] → [
2.13354:13433]
∅:D[
2.13433] → [
6.41215:41286]
B:BD[
6.41215] → [
6.41215:41286]
∅:D[
6.41286] → [
7.34837:34869]
B:BD[
7.34837] → [
7.34837:34869]
B:BD[
7.34869] → [
6.41287:41326]
CCTK_REAL ddx_phi =
(phi[p.idx - p.di] - 2 * phi[p.idx] + phi[p.idx + p.di]) / pow(dx, 2);
CCTK_REAL ddy_phi =
(phi[p.idx - p.dj] - 2 * phi[p.idx] + phi[p.idx + p.dj]) / pow(dy, 2);
CCTK_REAL ddz_phi =
(phi[p.idx - p.dk] - 2 * phi[p.idx] + phi[p.idx + p.dk]) / pow(dz, 2);
CCTK_REAL psi_n = psi[p.idx] + dt * (ddx_phi + ddy_phi + ddz_phi);
CCTK_REAL dt_phi_p = psi_n;
CCTK_REAL dt_phi_m = psi_p[p.idx];
// CCTK_REAL ddx_phi =
// (phi[p.idx - p.di] - 2 * phi[p.idx] + phi[p.idx + p.di]) / pow(dx,
// 2);
// CCTK_REAL ddy_phi =
// (phi[p.idx - p.dj] - 2 * phi[p.idx] + phi[p.idx + p.dj]) / pow(dy,
// 2);
// CCTK_REAL ddz_phi =
// (phi[p.idx - p.dk] - 2 * phi[p.idx] + phi[p.idx + p.dk]) / pow(dz,
// 2);
// CCTK_REAL psi_n = psi[p.idx] + dt * (ddx_phi + ddy_phi + ddz_phi);
// CCTK_REAL dt_phi_p = psi_n;
// CCTK_REAL dt_phi_m = psi_p[p.idx];
CCTK_REAL dt_phi_p = psi[p.idx]; // we cheat
CCTK_REAL dt_phi_m = psi[p.idx];