# Volcano, dike example on PyLith wiki **URL:** \ **Category:** PyLith\ **Tags:** examples\ **Created:** [December 19, 2018, 4:53pm UTC](https://community.geodynamics.org/t/volcano-dike-example-on-pylith-wiki/297 "2018-12-19T16:53:13Z")\ **Posts on this page:** 4\ **Page:** 1
**Author:** ![baagaard](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/baagaard/32/615_2.png) [@baagaard](https://community.geodynamics.org/u/baagaard)\ **Post date:** [December 19, 2018, 4:53pm UTC](https://community.geodynamics.org/t/volcano-dike-example-on-pylith-wiki/297/1 "2018-12-19T16:53:13Z")
I have updated the volcano / dike example on the [PyLith wiki](https://wiki.geodynamics.org/software:pylith:examples). In smoothing the mesh using CUBIT 15.3 (and presumably current versions of Trelis as well), the mesh topology became incorrect resulting in NaN for the norm of the residual when running PyLith. Simply multiple smoothing passes with a single pass creates a good quality mesh with correct topology. I also updated the field split solver settings in `pylithapp.cfg` for PyLith v2.2.1. Thanks to Josh Jones for reporting this issue. ---
**Author:** ![remy](https://avatars.discourse-cdn.com/v4/letter/r/71e660/32.png) [@remy](https://community.geodynamics.org/u/remy)\ **Post date:** [April 23, 2019, 1:23pm UTC](https://community.geodynamics.org/t/volcano-dike-example-on-pylith-wiki/297/2 "2019-04-23T13:23:28Z")
Dear Dr Aagard, I would like to use pylith for modeling ground displacement. I have tested the example two\_source\_dike\_slip.cfg. I have generated the mesh using trelis pro 16.5 but I have the following problem when I run pylith the following message appears : [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: [0]PETSC ERROR: KSPSolve has not converged due to Nan or Inf norm It seems that this error is similar to the one you mention. Thank you for helping me Best Regard Dominique ---
**Author:** ![baagaard](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/baagaard/32/615_2.png) [@baagaard](https://community.geodynamics.org/u/baagaard)\ **Post date:** [April 23, 2019, 5:34pm UTC](https://community.geodynamics.org/t/volcano-dike-example-on-pylith-wiki/297/3 "2019-04-23T17:34:50Z")
This could be similar. Please create a log file with everything from stdout and stderr (`pylith two_source_dike_slip.cfg >& run.log`) so that we can see the entire error message. Please attach the log file and the pylith\_parameters.json file to your post. ---
**Author:** ![remy](https://avatars.discourse-cdn.com/v4/letter/r/71e660/32.png) [@remy](https://community.geodynamics.org/u/remy)\ **Post date:** [April 24, 2019, 7:24am UTC](https://community.geodynamics.org/t/volcano-dike-example-on-pylith-wiki/297/4 "2019-04-24T07:24:25Z")
Dear Dr Aagard, I have resolved this problem by suppressing in the file mesh\_simple\_var.jou the lines corresponding to the smoothing of the mesh volume all … smooth volume all. Now it seems OK. Thank you very much for your reply. Dominique