geodynamics/pylith v4.0.0

Changes in user parameters

• Changed name of fault Lagrange multiplier field for solution component in Python from `lagrange_fault` to `lagrange_multiplier_fault` to match name of solution field in C++.
• Removed support for importing meshes from LaGriT.

Other changes

• Change in fault tractions are now included in the fault `data_fields` for prescribed slip.
• Fault and boundary orientation directions are now included in the `info_fields` for simulation output.
• State variables are now included in the default `data_fields` for simulation output.
• The default solver settings use the PETSc proper orthogonal decomposition (POD) methodology for initial guess of solutions to improve convergence.
• Add demonstration of `pylith_powerlaw_gendb` in Step 8 of `examples/reverse-2d`.
• Add demonstration of using poroelasticity with porosity as a state variable to `examples/magma-2d`.
• Switched from CppUnit to Catch2 for the C++ testing framework.
• Improve integration with VSCode for testing and debugging (see Developer Guide)
• Bug fixes
  • Fix errors in KinSrcTimeHistory.py
  • Fix creation of PETSc label for edges when importing Gmsh files. This fixes creation of faults with buried edges for 3D meshes imported from Gmsh.
  • Add containers for solution fields for poroelasticity with faults.
• Update PETSc to version 3.20.2.
• Update Python requirement to version 3.8 or later.
• Update Pyre requirement to version 1.1.0 or later.
• Update SpatialData requirement to version 3.1.0 or later.

Refer to the PyLith webpage on geodynamics.org for more information.