Step 6: Error 9#

Error Message#

Listing 292 Error message 9 when running Step 6.#
 1$ pylith step06_twofaults.cfg
 2
 3[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
 4[0]PETSC ERROR: Error in external library
 5[0]PETSC ERROR: Could not find values for initiation_time at (  -21860.3  -27621) in spatial database 'Fault rupture for main fault'.
 6[0]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!
 7[0]PETSC ERROR:   Option left: name:-ksp_atol value: 1.0e-12 source: code
 8[0]PETSC ERROR:   Option left: name:-ksp_converged_reason (no value) source: code
 9[0]PETSC ERROR:   Option left: name:-ksp_error_if_not_converged (no value) source: code
10[0]PETSC ERROR:   Option left: name:-ksp_guess_pod_size value: 8 source: code
11[0]PETSC ERROR:   Option left: name:-ksp_guess_type value: pod source: code
12[0]PETSC ERROR:   Option left: name:-ksp_rtol value: 1.0e-12 source: code
13[0]PETSC ERROR:   Option left: name:-mg_fine_pc_type value: vpbjacobi source: code
14[0]PETSC ERROR:   Option left: name:-pc_type value: gamg source: code
15[0]PETSC ERROR:   Option left: name:-snes_atol value: 1.0e-9 source: code
16[0]PETSC ERROR:   Option left: name:-snes_converged_reason (no value) source: code
17[0]PETSC ERROR:   Option left: name:-snes_error_if_not_converged (no value) source: code
18[0]PETSC ERROR:   Option left: name:-snes_max_it value: 1 source: command line
19[0]PETSC ERROR:   Option left: name:-snes_monitor (no value) source: code
20[0]PETSC ERROR:   Option left: name:-snes_rtol value: 1.0e-12 source: code
21[0]PETSC ERROR:   Option left: name:-ts_error_if_step_fails (no value) source: code
22[0]PETSC ERROR:   Option left: name:-ts_monitor (no value) source: code
23[0]PETSC ERROR:   Option left: name:-ts_type value: beuler source: code
24[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.
25[0]PETSC ERROR: Petsc Development GIT revision: v3.21.2-167-g4fed2113cae  GIT Date: 2024-05-31 10:11:14 -0400
26[0]PETSC ERROR: /software/unix/py3.12-venv/pylith-debug/bin/mpinemesis on a arch-clang-15.0_debug named IGSKCI164LM006 by baagaard Wed Jun  5 15:14:26 2024
27[0]PETSC ERROR: Configure options --PETSC_ARCH=arch-clang-15.0_debug --with-debugging=1 --with-clanguage=c --with-mpi-compilers=1 --with-shared-libraries=1 --with-64-bit-points=1 --with-large-file-io=1 --with-lgrind=0 --download-chaco=1 --download-parmetis=1 --download-metis=1 --download-triangle --download-ml=1 --download-superlu=1 --with-fc=0 --download-f2cblaslapack --with-hdf5=1 --with-hdf5-include=/software/unix/hdf5-1.14/clang-15.0/include --with-hdf5-lib=/software/unix/hdf5-1.14/clang-15.0/lib/libhdf5.dylib --with-zlib=1
28[0]PETSC ERROR: #1 static PetscErrorCode pylith::topology::FieldQuery::queryDBPointFn(PylithInt, PylithReal, const PylithReal *, PylithInt, PylithScalar *, void *)() at /src/cig/pylith/libsrc/pylith/topology/FieldQuery.cc:316
29[0]PETSC ERROR: #2 DMProjectPoint_Func_Private() at /software/unix/petsc-dev/src/dm/impls/plex/plexproject.c:128
30[0]PETSC ERROR: #3 DMProjectPoint_Private() at /software/unix/petsc-dev/src/dm/impls/plex/plexproject.c:483
31[0]PETSC ERROR: #4 DMProjectLocal_Generic_Plex() at /software/unix/petsc-dev/src/dm/impls/plex/plexproject.c:1003
32[0]PETSC ERROR: #5 DMProjectFunctionLocal_Plex() at /software/unix/petsc-dev/src/dm/impls/plex/plexproject.c:1034
33[0]PETSC ERROR: #6 DMProjectFunctionLocal() at /software/unix/petsc-dev/src/dm/interface/dm.c:8183
34[0]PETSC ERROR: #7 void pylith::topology::FieldQuery::queryDB()() at /src/cig/pylith/libsrc/pylith/topology/FieldQuery.cc:223
35Fatal error. Calling MPI_Abort() to abort PyLith application.
36Traceback (most recent call last):
37  File "/software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/apps/PetscApplication.py", line 55, in onComputeNodes
38    self.main(*args, **kwds)
39  File "/software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/apps/PyLithApp.py", line 103, in main
40    self.problem.initialize()
41  File "/software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/problems/Problem.py", line 187, in initialize
42    ModuleProblem.initialize(self)
43  File "/software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/problems/problems.py", line 165, in initialize
44    return _problems.Problem_initialize(self)
45           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
46RuntimeError: Error detected while in PETSc function.
47Abort(-1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0
48/software/unix/py3.12-venv/pylith-debug/bin/nemesis: mpiexec: exit 255
49/software/unix/py3.12-venv/pylith-debug/bin/pylith: /software/unix/py3.12-venv/pylith-debug/bin/nemesis: exit 1

Troubleshooting Strategy#

With the linear interpolation we get an error about not being able to find an initiation time for a point. This suggests there are errors in our spatial database file related to interpolation. We examine the header and data points for errors. We notice that our points lie along a line (data dimension is 1), but our header has data-dim=2.

Resolution#

Listing 293 Correct error in fault_slip.spatialdb.#
# Error
data-dim = 2
    
# Correct
data-dim = 1