Non-symmetric magnetohydrostatic equilibria: a multigrid approach