Changed Update default solver settings for incompressible elasticity for Schur factorization from `full to lower. Set the default length scale to 1 km instead of 100 km for consistency with most quasi-static problems with faults. For simulations with faults, the length scale should be set to the discretization size on the fault; this is related to the preconditioner which links the length scale to the discretization size. Improved documentation for solver tolerances and adaptive time stepping. Removed manual marking of buried edges (use automatic marking) in examples and tests. Updated PETSc version to 3.25.4. Added examples/subduction-3d: Add expected Gmsh output for geometric operations in documentation. Added configure option for pip install without build isolation (`--without-build-isolation`). Fixed Optimized the check of cohesive cells, which greatly reduces the total time for inserting cohesive cells. examples/subduction-3d: Increased tolerance for slab surface creation. Fixed setting the mesh filename in parameter files using Cubit mesh. Corrected documentation for isotropic linear poroelasticity Jacobian terms (Jf3uu and Jf1eu; now matches code). Developer Changes Removed flag (--output-sync=target) in top-level Makefile.am that is not backwards compatible.