implement all regime well balanced 3 wave solver for MHD