[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