Implement test for all MHD operators based on PSYDAC
Test all MHD operators based on PSYDAC using a for-loop. This new version will replace the old one in a future merge request.
Commits:
- Ignore vertical lines when computing y-axis limits.
- Implement test for psydac MHD operators, with major bug fixes.
- Fix older code in psydac MHD operators test.
- In polar spline vis., first write data, then read it back for plotting.
- Convert DomainType (analytical or numerical mapping) into an enum.
- Update line plot style.
- Special handling of X1, X2 MHD operators which give 3x scalar outputs.
- Provide individual functions to assemble MHD operators.
- Restore
MHD_ops.__init__()
after rebase. - Restore flattening returned DOF.
- Revert my own
MHD_ops.assemble_XX()
implementations. - Do not return
rhs
duringdot()
, due to special handling of X1 X2. - Resolve naming conflict between projectors and operators
P1
,P2
. - Rearrange statements in different versions of
psydac_lin_ops
tests. - Seed random number generation (fix mpirun assertion error).
- Move asserts into a function.
- Renamed psydac lin ops test file.