Dear Sir,
I would like to ask a question regarding my use of poroelastic simulations in Pylith 4.1.3. When I use the no_faultzone_permeability3.spatialdb file for all materials, the simulation runs successfully. However, when I use the no_faultzone_permeability.spatialdb file for all materials, an error occurs, displaying the following message:
[0]: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Residual norm computed by GMRES recursion formula 16442.8 is far from the computed residual norm 41861.1 at restart, residual norm at start of cycle 97030.1
Could you please advise on how to resolve this issue? I found a reference in the manual mentioning an “edge” issue, but my simulation only involves blocks.
Additionally, I would like to know how to output stress, pressure, and other values on the fault surface in Pylith 4.1.3. The syntax seems to differ significantly from that in Pylith 2, and I have not yet found any relevant examples. I apologize for the inconvenience.
file.zip (1.6 MB)