Even with spatial adaptivity, the number of mesh data points can be unnecessarily large if care is not
taken as to where refinement is done. This paper proposes an adaptive scheme that generates the mesh
by recursive bisection, allowing mesh blocks to be arbitrarily anisotropic to allow for fine structures in
some directions without over-refining in those directions that suffice with less refinement.