Step 1: Successful run#
1$ pylith step01a_gravity.cfg
2
3# Output
4>> software/pylith-debug/lib/python3.12/site-packages/pylith/apps/PyLithApp.py:76:main
5-- info (application-flow)
6-- Running on 1 process(es).
7>> 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&)
8-- info (application-flow)
9-- Setting PETSc options:
10 ksp_atol = 1.0e-7
11 ksp_converged_reason = true
12 ksp_error_if_not_converged = true
13 ksp_gmres_restart = 100
14 ksp_guess_pod_size = 8
15 ksp_guess_type = pod
16 ksp_rtol = 1.0e-14
17 mg_fine_ksp_max_it = 5
18 pc_type = gamg
19 snes_atol = 5.0e-7
20 snes_converged_reason = true
21 snes_error_if_not_converged = true
22 snes_monitor = true
23 snes_rtol = 1.0e-14
24 ts_error_if_step_fails = true
25 ts_exact_final_time = matchstep
26 ts_monitor = true
27 ts_type = beuler
28 viewer_hdf5_collective = true
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.