[CIG-SHORT] Pylith's error including `dtStable > 0.0' failed

hkj at mail.ustc.edu.cn hkj at mail.ustc.edu.cn
Fri Aug 17 02:45:54 PDT 2018


Dear Sir,


When I running pylith using my own FEM model, the error including" `dtStable > 0.0' failed"  appears very frequently, and I found that this can be avoided by using fewer nodes like


from
   pylith postseis.cfg mat_viscoelastic.cfg solver_fieldsplit.cfg --nodes=20
to
   pylith postseis.cfg mat_viscoelastic.cfg solver_fieldsplit.cfg --nodes=10


But after I rebuilt my FEM model by using a larger area, this error still occurs even if I use only 1 node,
I also chose dt from dt = 1*year to dt = 0.001*year , but it didn't work.




The .cfg file is attached, and the time formulation part is 


```
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.problem.formulation.time_step]
# Define the total time for the simulation and the time step size.
total_time = 2*year
dt = 0.001*year
```


And the output error is


``` [OUTPUT]
.......
.......
.......
-- Creating Jacobian matrix.
 >> /home/khuang/Pylith/lib/python2.7/site-packages/pylith/problems/Implicit.py:163:initialize
 -- implicit(info)
 -- Initializing solver.
 >> /home/khuang/Pylith/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:132:run
 -- timedependent(info)
 -- Solving problem.
 >> /home/khuang/Pylith/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:138:run
 -- timedependent(info)
 -- Preparing for prestep with elastic behavior.
mpinemesis: ../../../pylith-2.2.1/libsrc/pylith/materials/ElasticMaterial.cc:405: virtual PylithScalar pylith::materials::ElasticMaterial::stableTimeStepImplicit(const pylith::topology::Mesh&, pylith::topology::Field*): Assertion `dtStable > 0.0' failed.
[0]0:Return code = 0, signaled with Aborted
```


Can you tell me how to fix this error?


Sincerely
Kejing Huang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180817/3a5ca78c/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: postseis.cfg
Type: application/octet-stream
Size: 6039 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180817/3a5ca78c/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pylithapp.cfg
Type: application/octet-stream
Size: 7690 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180817/3a5ca78c/attachment-0001.obj>


More information about the CIG-SHORT mailing list