Open iljah opened 8 months ago
Initializing magnetic field using one face-centered value might create too large divergence when AMR is involved. Investigate using e.g. https://www.boost.org/doc/libs/1_84_0/libs/math/doc/html/math_toolkit/gauss.html to initialize B0 and B1 more accurately in order to decrease div(B).
Quick test of using average of B values at 4 points on face decreased max div(B) by 4x in 2d wave test.
Initializing magnetic field using one face-centered value might create too large divergence when AMR is involved. Investigate using e.g. https://www.boost.org/doc/libs/1_84_0/libs/math/doc/html/math_toolkit/gauss.html to initialize B0 and B1 more accurately in order to decrease div(B).