| ... | @@ -947,9 +947,9 @@ divergence-free setup: | 
... | @@ -947,9 +947,9 @@ divergence-free setup: | 
| 
 | 
components by exact integration of given analytical expressions. For a
 | 
 | 
components by exact integration of given analytical expressions. For a
 | 
| 
 | 
grid cell (ix,iy,iz) on superblock `g` the discretized components would then
 | 
 | 
grid cell (ix,iy,iz) on superblock `g` the discretized components would then
 | 
| 
 | 
read 
 | 
 | 
read 
 | 
| 
 | 
$$\mathtt{g->bx[iz][iy][ix]}=\frac{1}{\delta \mathcal{A}_x}\int
 | 
 | 
`g->bx[iz][iy][ix]`$=\frac{1}{\delta \mathcal{A}_x}\int
 | 
| 
 | 
\limits_\mathtt{g->y[iy]}^\mathtt{g->y[iy+1]}
 | 
 | 
\limits_\mathtt{g->y[iy]}^\mathtt{g->y[iy+1]}
 | 
| 
 | 
\int\limits_\mathtt{g->z[iz]}^\mathtt{g->z[iz+1]} B_x(\mathtt{g->x[ix]},y,z)h_yh_zdydz$$
 | 
 | 
\int\limits_\mathtt{g->z[iz]}^\mathtt{g->z[iz+1]} B_x(\mathtt{g->x[ix]},y,z)h_yh_zdydz$
 | 
| 
 | 
$$\mathtt{g->by[iz][iy][ix]}=\frac{1}{\delta \mathcal{A}_y}\int
 | 
 | 
$$\mathtt{g->by[iz][iy][ix]}=\frac{1}{\delta \mathcal{A}_y}\int
 | 
| 
 | 
\limits_\mathtt{g->x[ix]}^\mathtt{g->x[ix+1]}
 | 
 | 
\limits_\mathtt{g->x[ix]}^\mathtt{g->x[ix+1]}
 | 
| 
 | 
\int\limits_\mathtt{g->z[iz]}^\mathtt{g->z[iz+1]} B_y(x,\mathtt{g->y[iy]},z)h_zdxdz$$
 | 
 | 
\int\limits_\mathtt{g->z[iz]}^\mathtt{g->z[iz+1]} B_y(x,\mathtt{g->y[iy]},z)h_zdxdz$$
 | 
| ... | @@ -971,7 +971,7 @@ integration over cell faces turns out to be too complicated. | 
... | @@ -971,7 +971,7 @@ integration over cell faces turns out to be too complicated. | 
| 
 | 
The face-averaged magnetic field components are then obtained from the
 | 
 | 
The face-averaged magnetic field components are then obtained from the
 | 
| 
 | 
integral form of $\mathbf{B}=\nabla\times \mathbf{A}$:
 | 
 | 
integral form of $\mathbf{B}=\nabla\times \mathbf{A}$:
 | 
| 
 | 
 | 
 | 
 | 
| 
 | 
$\mathtt{g->bx[iz][iy][ix]}=\left[h_y\Delta_{iy} (h_{zy}\hat{A}_z)
 | 
 | 
`g->bx[iz][iy][ix]`$=\left[h_y\Delta_{iy} (h_{zy}\hat{A}_z)
 | 
| 
 | 
-h_y\Delta_{iz} \hat{A}_y\right]/\delta \mathcal{A}_x$
 | 
 | 
-h_y\Delta_{iz} \hat{A}_y\right]/\delta \mathcal{A}_x$
 | 
| 
 | 
 | 
 | 
 | 
| 
 | 
`g->by[iz][iy][ix]`$=\left[\Delta_{iz} \hat{A}_x
 | 
 | 
`g->by[iz][iy][ix]`$=\left[\Delta_{iz} \hat{A}_x
 | 
| ... |  | ... |  |