  The gap between a recently developed dynamical version of relaxed
magnetohydrodynamics (RxMHD) and ideal MHD (IMHD) is bridged by approximating
the zero-resistivity "Ideal" Ohm's Law (IOL) constraint using an augmented
Lagrangian method borrowed from optimization theory. The augmentation combines
a pointwise vector Lagrange multiplier method and global penalty function
method and can be used either for iterative enforcement of the IOL to arbitrary
accuracy, or for constructing a continuous sequence of magnetofluid dynamics
models running between RxMHD (no IOL) and weak IMHD (IOL almost everywhere).
This is illustrated by deriving dispersion relations for linear waves on an MHD
equilibrium.
