Fixes cubic solver. It was previously wrong.
Rewrites MHD jump solver replacing quartic solve with cubic by recognizing known monominal in quartic.