[CIG-SHORT] some confusion about modeling

Matthew Knepley knepley at rice.edu
Thu Feb 9 08:20:22 PST 2017


On Thu, Feb 9, 2017 at 3:41 AM, fuqiang shi <shifuqiang121 at gmail.com> wrote:

> Dear,Sir,
>
>   I am a new learner of cubit+pylith. I am sorry to disturb you for some
> unexpected running errors, and I have tried all my best to solve them,
> but no answer.
> The error is totally like this:
> Maybe it have some relation with the solver, but I haven't do any modify
> on the solution parameter of the provided example.  I only modify the model
> slightly and my model is also attached (s. attachment).
> Your consideration and suggestions upon my request will be highly
> appreciated. I am looking forward to your favorable reply at your earliest
> convenience.
>

I responded to this question already. I said

This means the linear solver did not converge. For your initial tests I
recommend using a direct solver. This
will help you uncover inconsistencies in the specification, such as bad
boundary conditions. For problems with
no fault

  --petsc.pc_type lu

or with a fault use the settings in share/settings/solver_fault_exact.cfg.

When reporting a solver problem, always run with

  --petsc.ksp_view --petsc.ksp_pre_view --petsc.ksp_monitor

and send all the output.

  Thanks,

     Matt


>
>
> Sincerely,
>
>                                                                  Fuqiang
> Shi
>
> Linear solve did not converge due to DIVERGED_ITS iterations 100
> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [0]PETSC ERROR:
> [0]PETSC ERROR: KSPSolve has not converged
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [0]PETSC ERROR: Petsc Development GIT revision: v3.7.2-876-g7e45237  GIT
> Date: 2016-06-16 04:41:50 -0500
> [0]PETSC ERROR: /home/shifuqiang/pylith/bin/mpinemesis on a arch-pylith
> named tecto-30 by shifuqiang Thu Feb  9 10:10:23 2017
> [0]PETSC ERROR: Configure options --prefix=/home/brad/pylith-binary/dist
> --with-c2html=0 --with-x=0 --with-clanguage=C --with-mpicompilers=1
> --with-debugging=0 --with-shared-libraries=1 --with-64-bit-points=1
> --with-large-file-io=1 --download-chaco=1 --download-ml=1
> --download-f2cblaslapack=1 --with-hwloc=0 --with-ssl=0 --with-x=0
> --with-c2html=0 --with-lgrind=0 --with-hdf5=1 --with-hdf5-dir=/home/brad/pylith-binary/dist
> --LIBS=-lz --with-fc=0 CPPFLAGS="-I/home/brad/pylith-binary/dist/include
> " LDFLAGS="-L/home/brad/pylith-binary/dist/lib " CFLAGS="-g -O2"
> CXXFLAGS="-g -O2 -DMPICH_IGNORE_CXX_SEEK" FCFLAGS=
> PETSC_DIR=/home/brad/pylith-binary/build/petsc-pylith
> PETSC_ARCH=arch-pylith
> [0]PETSC ERROR: #1 KSPSolve() line 847 in /home/brad/pylith-binary/build
> /petsc-pylith/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: #2 void pylith::problems::SolverLinear
> ::solve(pylith::topology::Field*, pylith::topology::Jacobian*, const
> pylith::topology::Field&)() line 132 in ../../../pylith-2.1.4/libsrc/p
> ylith/problems/SolverLinear.cc
> Fatal error. Calling MPI_Abort() to abort PyLith application.
> Traceback (most recent call last):
>   File "/home/shifuqiang/pylith/lib/python2.7/site-packages/pylith/apps/PetscApplication.py",
> line 74, in onComputeNodes
>     self.main(*args, **kwds)
>   File "/home/shifuqiang/pylith/lib/python2.7/site-packages/pylith/apps/PyLithApp.py",
> line 135, in main
>     self.problem.run(self)
>   File "/home/shifuqiang/pylith/lib/python2.7/site-packages/pylith/problems/TimeDependent.py",
> line 155, in run
>     self.formulation.step(t, dt)
>   File "/home/shifuqiang/pylith/lib/python2.7/site-packages/pylith/problems/Implicit.py",
> line 212, in step
>     self.solver.solve(dispIncr, self.jacobian, residual)
>   File "/home/shifuqiang/pylith/lib/python2.7/site-packages/pylith/problems/problems.py",
> line 166, in solve
>     def solve(self, *args): return _problems.SolverLinear_solve(self,
> *args)
> RuntimeError: Error detected while in PETSc function.
> application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0
> /home/shifuqiang/pylith/bin/nemesis: mpirun: exit 255
> /home/shifuqiang/pylith/bin/pylith: /home/shifuqiang/pylith/bin/nemesis:
> exit 1
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170209/0a38eed3/attachment.html>


More information about the CIG-SHORT mailing list