| One probably needs multigrid elliptic solvers which are compatible with hyperbolic evolution using Adaptive Mesh Refinement |
| One might to solve elliptic equations at each time step for subsidary variables |
| However neither Grand Challenge nor (significant) part of numerical analysis community seemed to have experience in this and there appears no "off the shelf" solution. |