Step 1: Successful run#

Listing 307 Output when running Step 1.#
 1$ pylith step01a_gravity.cfg
 2
 3 >> software/pylith-debug/lib/python3.12/site-packages/pylith/apps/PyLithApp.py:76:main
 4 -- info (application-flow)
 5 -- Running on 1 process(es).
 6 >> src/cig/pylith/libsrc/pylith/utils/PetscOptions.cc:251:static void pylith::utils::_PetscOptions::write(pythia::journal::info_t&, const char*, const pylith::utils::PetscOptions&)
 7 -- info (application-flow)
 8 -- Setting PETSc options:
 9    ksp_atol = 1.0e-7
10    ksp_converged_reason = true
11    ksp_error_if_not_converged = true
12    ksp_gmres_restart = 100
13    ksp_guess_pod_size = 8
14    ksp_guess_type = pod
15    ksp_rtol = 1.0e-14
16    mg_fine_ksp_max_it = 5
17    pc_type = gamg
18    snes_atol = 5.0e-7
19    snes_converged_reason = true
20    snes_error_if_not_converged = true
21    snes_monitor = true
22    snes_rtol = 1.0e-14
23    ts_error_if_step_fails = true
24    ts_exact_final_time = matchstep
25    ts_monitor = true
26    ts_type = beuler
27    viewer_hdf5_collective = true
28
29 >> src/cig/pylith/libsrc/pylith/meshio/MeshIOPetsc.cc:204:virtual void pylith::meshio::MeshIOPetsc::_read()
30 -- info (application-flow)
31 -- Component 'meshiopetsc.reader': Reading finite-element mesh from 'mesh_tri.msh'.
32 >> src/cig/pylith/libsrc/pylith/meshio/MeshIO.cc:76:void pylith::meshio::MeshIO::read(pylith::topology::Mesh*, bool)
33 -- info (application-flow)
34 -- Component 'meshiopetsc.reader': Domain bounding box:
35    (-100000, 100000)
36    (-100000, 0)
37 >> src/cig/pylith/libsrc/pylith/problems/TimeDependent.cc:316:virtual void pylith::problems::TimeDependent::verifyConfiguration() const
38 -- info (application-flow)
39 -- Component 'timedependent.problem': Verifying problem configuration.
40 >> software/pylith-debug/lib/python3.12/site-packages/pylith/problems/Problem.py:238:_printInfo
41 -- info (application-flow)
42 -- Scales for nondimensionalization:
43    Length scale: 5000*m
44    Displacement scale: 1*m
45    Time scale: 3.15576e+09*s
46    Rigidity scale: 1e+10*m**-1*kg*s**-2
47    Temperature scale: 1*K
48 >> src/cig/pylith/libsrc/pylith/problems/TimeDependent.cc:342:virtual void pylith::problems::TimeDependent::initialize()
49 -- info (application-flow)
50 -- Component 'timedependent.problem': Initializing problem.
51 >> src/cig/pylith/libsrc/pylith/problems/TimeDependent.cc:473:void pylith::problems::TimeDependent::solve()
52 -- info (application-flow)
53 -- Component 'timedependent.problem': Solving equations.
540 TS dt 0.001 time 0.
55    0 SNES Function norm 3.007829881319e+03
56      Linear solve converged due to CONVERGED_ATOL iterations 11
57    1 SNES Function norm 3.524261285349e-07
58    Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 1
591 TS dt 0.001 time 0.001
60 >> software/pylith-debug/lib/python3.12/site-packages/pylith/problems/Problem.py:222:finalize
61 -- info (application-flow)
62 -- Finalizing problem.

Our simulation now runs without errors and the output looks correct.