Hello,
ASPECT 3.1.0-pre is not solving timestep 0 of the Dannberg et al. [2019], mid-ocean ridge melting model. It produces solutions, but when it gets down to the last refinement level the BlockSchurPreconditioner fails. I’ve tried relaxing the Linear solver tolerance to 1e-6, and reducing the global refinement from 6 to 5, but it ends up crashing in the same solver block (in timestep 1 or 2) . I’ve attached my input file (adapted from “mid_ocean_ridge_adaptive_porosity.prm”), std out, and error output. There was no solver_history.txt file produced. A few of the original input parameters have changed or gone away (noted by “G.Ito”) so I wonder if those contribute to the issues.
Thanks in advance,
Garrett
input_Lst1e-6.prm (11.3 KB)
log_Lst1e-6.out.txt (9.4 KB)
log_Lst1e-6.err.txt (20.2 KB)