Switching to Piola transform for simulation_type 2, for the RHS. This need not be built numerically. Instead, I'm using Sage for CAS solutions.