[CIG-SHORT] Nonlinear solve did not converge
Marcelo Contreras
marcontk at gmail.com
Wed Oct 22 10:12:54 PDT 2014
Dear Brad,
I´m following some examples that include laws friction for a 3D subduction model (interseismic). Using petsc parameter from example step13.cfg (Slip weakening) there is no convergence (linear) . Now with the folllowing parameters I have solution (see figura1.png attached)
# ejemplo
snes_view = true
ksp_monitor_true_residual = true
fs_pc_type = fieldsplit
fs_pc_fieldsplit_schur_precondition = user
fs_pc_fieldsplit_real_diagonal = true
fs_pc_fieldsplit_type = schur
fs_pc_fieldsplit_schur_factorization_type = full
fs_fieldsplit_0_ksp_type = preonly
fs_fieldsplit_0_pc_type = ml
fs_fieldsplit_1_ksp_type = gmres
fs_fieldsplit_1_ksp_rtol = 1.0e-05
fs_fieldsplit_1_pc_type = jacobi
# Convergence parameters.
ksp_rtol = 1.0e-20
ksp_atol = 1.0e-09
ksp_max_it = 300
ksp_gmres_restart = 50
The problem stars when I activate "solver = pylith.problems.SolverNonlinear” the error is " Nonlinear solve did not converge”. Petsc parameters are: (Using parameters in step13.cfg didn’t work)
#--- NonLinear ----------
snes_monitor_short = true
snes_view = true
snes_type = ngmres
npc_snes_max_it = 1
npc_levels_snes_linesearch_type = basic
npc_snes_type_coarse_snes_converged_reason = true
npc_snes_type_levels_snes_type = newtonls
npc_snes_type_levels_snes_max_it = 10
npc_snes_type_levels_snes_linesearch_type = basic
npc_snes_type_levels_snes_max_linear_solve_fail = 30
npc_snes_type_levels_ksp_max_it = 20
npc_snes_type_levels_snes_converged_reason = true
npc_snes_type_coarse_snes_linesearch_type = basic
# Nonlinear solver monitoring options.
snes_rtol = 1.0e-20
snes_atol = 1.0e-09
snes_max_it = 10
#snes_monitor = true
snes_linesearch_monitor = true
snes_view = true
snes_converged_reason = true
Do you can help me with some idea?
best regards
Marcelo.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/utils/PetscManager.py:64:initialize
-- petsc(info)
-- Initialized PETSc.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/meshio/MeshIOObj.py:55:read
-- meshiocubit(info)
-- Reading finite-element mesh
>> meshio/MeshIOCubit.cc:142:<unknown>
-- meshiocubit(info)
-- Reading 87059 vertices.
>> meshio/MeshIOCubit.cc:197:<unknown>
-- meshiocubit(info)
-- Reading 473405 cells in 5 blocks.
>> meshio/MeshIOCubit.cc:257:<unknown>
-- meshiocubit(info)
-- Found 13 node sets.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'ALLNODES' with id 1 containing 87059 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'BOTTOM_MANTLE_BORDER' with id 2 containing 837 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'BOTTOM_SLAB_BORDER' with id 3 containing 480 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'EAST_BORDER' with id 4 containing 1539 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'FAULT_SEISMO' with id 5 containing 2135 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'FAULT_SLABBOT' with id 6 containing 10045 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'FAULT_SLABTOP' with id 7 containing 5742 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'GROUNDSURF' with id 8 containing 15822 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'LOFZ' with id 9 containing 293 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'NORTH_BORDER' with id 10 containing 1154 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'SOUTH_BORDER' with id 11 containing 839 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'WEST_MANTLE_BORDER' with id 12 containing 148 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'WEST_SLAB_BORDER' with id 13 containing 402 nodes.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:92:preinitialize
-- timedependent(info)
-- Pre-initializing problem.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:390:_setupMaterials
-- implicit(info)
-- Pre-initializing materials.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x105646990> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'C_CRUST_SUP'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065e8e40> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'C_CRUST_INF'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065f00c0> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'C_MANTLE'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065f0300> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'SLAB'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065f0540> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'O_MANTLE'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:420:_setupBC
-- implicit(info)
-- Pre-initializing boundary conditions.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'EAST_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'WEST_MANTLE_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'WEST_SLAB_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'BOTTOM_MANTLE_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:459:_setupInterfaces
-- implicit(info)
-- Pre-initializing interior interfaces.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x1065f0bd0> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'FAULT_SEISMO' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:108:preinitialize
-- faultcohesivekin(info)
-- Pre-initializing fault 'FAULT_SLABTOP'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x106601420> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'FAULT_SLABTOP' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:108:preinitialize
-- faultcohesivekin(info)
-- Pre-initializing fault 'FAULT_SLABBOT'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x1066016f0> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'FAULT_SLABBOT' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x1066019c0> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'LOFZ' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:188:preinitialize
-- implicit(info)
-- Pre-initializing output.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Problem.py:150:verifyConfiguration
-- timedependent(info)
-- Verifying compatibility of problem configuration.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:117:initialize
-- timedependent(info)
-- Initializing problem.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:500:_initialize
-- implicit(info)
-- Initializing integrators.
WARNING: No initial state given for friction model 'Slip weakening'. Using default value of zero.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:160:initialize
-- faultcohesivekin(info)
-- Initializing fault 'FAULT_SLABTOP'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:160:initialize
-- faultcohesivekin(info)
-- Initializing fault 'FAULT_SLABBOT'.
WARNING: No initial state given for friction model 'Slip weakening'. Using default value of zero.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:510:_initialize
-- implicit(info)
-- Initializing constraints.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:516:_initialize
-- implicit(info)
-- Setting up solution output.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:525:_initialize
-- implicit(info)
-- Creating solution field.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:136:initialize
-- implicit(info)
-- Creating other fields.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:159:initialize
-- implicit(info)
-- Creating Jacobian matrix.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:169:initialize
-- implicit(info)
-- Initializing solver.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:131:run
-- timedependent(info)
-- Solving problem.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:137:run
-- timedependent(info)
-- Preparing for prestep with elastic behavior.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:266:prestepElastic
-- implicit(info)
-- Setting constraints.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:582:_reformJacobian
-- implicit(info)
-- Integrating Jacobian operator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:152:run
-- timedependent(info)
-- Computing prestep with elastic behavior.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:607:_reformResidual
-- implicit(info)
-- Integrating residual term in operator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:215:step
-- implicit(info)
-- Solving equations.
0 SNES Function norm < 1.e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 0
SNES Object: 1 MPI processes
type: ngmres
Number of stored past updates: 30
Residual selection: gammaA=2e+00, gammaC=2e+00
Difference restart: epsilonB=1e-01, deltaB=9e-01
maximum iterations=10, maximum function evaluations=30000
tolerances: relative=1e-20, absolute=1e-09, solution=1e-08
total number of linear solver iterations=0
total number of function evaluations=1
SNESLineSearch Object: 1 MPI processes
type: shell
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=1
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:158:run
-- timedependent(info)
-- Finishing prestep with elastic behavior.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:250:poststep
-- implicit(info)
-- Writing solution fields.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:176:run
-- timedependent(info)
-- Main time loop, current time is t=0*s
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:187:run
-- timedependent(info)
-- Preparing to advance solution from time t=0*s to t=3.15576e+07*s.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:186:prestep
-- implicit(info)
-- Setting constraints.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:193:run
-- timedependent(info)
-- Advancing solution from t=0*s to t=3.15576e+07*s.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:607:_reformResidual
-- implicit(info)
-- Integrating residual term in operator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:215:step
-- implicit(info)
-- Solving equations.
0 SNES Function norm 0.948594
Line search: gnorm after quadratic fit 1.773100006763e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 2.002776491982e+00 lambda=1.0000000000000002e-02
1 SNES Function norm 2.00278
Line search: gnorm after quadratic fit 4.312556301433e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 5.540166796547e+00 lambda=1.0000000000000002e-02
2 SNES Function norm 0.569893
Line search: gnorm after quadratic fit 2.243738258592e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 1.991117051814e+00 lambda=1.0000000000000002e-02
3 SNES Function norm 0.425075
Line search: gnorm after quadratic fit 2.147718836488e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 2.049829198737e+00 lambda=1.0000000000000002e-02
4 SNES Function norm 0.391639
Line search: gnorm after quadratic fit 1.608131254064e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 1.636970057178e+00 lambda=1.0000000000000002e-02
5 SNES Function norm 1.63697
Line search: gnorm after quadratic fit 5.396512956049e+02
Line search: Cubic step no good, shrinking lambda, current gnorm 5.278748490904e+01 lambda=1.0000000000000002e-02
6 SNES Function norm 52.7875
Line search: gnorm after quadratic fit 2.322268614406e+04
Line search: Cubic step no good, shrinking lambda, current gnorm 2.275760046175e+03 lambda=1.0000000000000002e-02
7 SNES Function norm 2275.76
Line search: gnorm after quadratic fit 1.071505794820e+06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.051316283060e+05 lambda=1.0000000000000002e-02
8 SNES Function norm 105132
Line search: gnorm after quadratic fit 5.228582697605e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.135201186511e+06 lambda=1.0000000000000002e-02
9 SNES Function norm 5.1352e+06
Line search: gnorm after quadratic fit 2.686941666753e+09
Line search: Cubic step no good, shrinking lambda, current gnorm 2.641281837342e+08 lambda=1.0000000000000002e-02
10 SNES Function norm 2.64128e+08
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 10
SNES Object: 1 MPI processes
type: ngmres
Number of stored past updates: 30
Residual selection: gammaA=2e+00, gammaC=2e+00
Difference restart: epsilonB=1e-01, deltaB=9e-01
maximum iterations=10, maximum function evaluations=30000
tolerances: relative=1e-20, absolute=1e-09, solution=1e-08
total number of linear solver iterations=0
total number of function evaluations=51
SNESLineSearch Object: 1 MPI processes
type: shell
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=1
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:200:run
-- timedependent(info)
-- Finishing advancing solution from t=0*s to t=3.15576e+07*s.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:250:poststep
-- implicit(info)
-- Writing solution fields.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:300:finalize
-- implicit(info)
-- Formulation finalize.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/utils/PetscManager.py:75:finalize
-- petsc(info)
-- Finalizing PETSc.
marcont at OSXmate~/pylith190/modelos/test_viscoelastic_05: pylith inter_3D.cfg
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/utils/PetscManager.py:64:initialize
-- petsc(info)
-- Initialized PETSc.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/meshio/MeshIOObj.py:55:read
-- meshiocubit(info)
-- Reading finite-element mesh
>> meshio/MeshIOCubit.cc:142:<unknown>
-- meshiocubit(info)
-- Reading 87059 vertices.
>> meshio/MeshIOCubit.cc:197:<unknown>
-- meshiocubit(info)
-- Reading 473405 cells in 5 blocks.
>> meshio/MeshIOCubit.cc:257:<unknown>
-- meshiocubit(info)
-- Found 13 node sets.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'ALLNODES' with id 1 containing 87059 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'BOTTOM_MANTLE_BORDER' with id 2 containing 837 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'BOTTOM_SLAB_BORDER' with id 3 containing 480 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'EAST_BORDER' with id 4 containing 1539 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'FAULT_SEISMO' with id 5 containing 2135 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'FAULT_SLABBOT' with id 6 containing 10045 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'FAULT_SLABTOP' with id 7 containing 5742 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'GROUNDSURF' with id 8 containing 15822 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'LOFZ' with id 9 containing 293 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'NORTH_BORDER' with id 10 containing 1154 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'SOUTH_BORDER' with id 11 containing 839 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'WEST_MANTLE_BORDER' with id 12 containing 148 nodes.
>> meshio/MeshIOCubit.cc:285:<unknown>
-- meshiocubit(info)
-- Reading node set 'WEST_SLAB_BORDER' with id 13 containing 402 nodes.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:92:preinitialize
-- timedependent(info)
-- Pre-initializing problem.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:390:_setupMaterials
-- implicit(info)
-- Pre-initializing materials.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x105646990> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'C_CRUST_SUP'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065e8e40> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'C_CRUST_INF'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065f00c0> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'C_MANTLE'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065f0300> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'SLAB'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTet3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTet3D *' at 0x1065f0540> >
-- Vertices:
-- [[-1. -1. -1.]
[ 1. -1. -1.]
[-1. 1. -1.]
[-1. -1. 1.]]
-- Quad pts:
-- [[-0.5 -0.5 -0.5]]
-- Quad wts:
-- [ 1.33333333]
-- Basis fns @ quad pts ):
-- [[ 0.25 0.25 0.25 0.25]]
-- Basis fn derivatives @ quad pts:
-- [[[ -5.00000000e-01 -5.00000000e-01 -5.00000000e-01]
[ 5.00000000e-01 0.00000000e+00 1.38777878e-17]
[ 0.00000000e+00 5.00000000e-01 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 5.00000000e-01]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:404:_setupMaterials
-- implicit(info)
-- Added elasticity integrator for material 'O_MANTLE'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:420:_setupBC
-- implicit(info)
-- Pre-initializing boundary conditions.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'EAST_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'WEST_MANTLE_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'WEST_SLAB_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:439:_setupBC
-- implicit(info)
-- Added boundary condition 'BOTTOM_MANTLE_BORDER' as a constraint.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:459:_setupInterfaces
-- implicit(info)
-- Pre-initializing interior interfaces.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x1065f0bd0> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'FAULT_SEISMO' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:108:preinitialize
-- faultcohesivekin(info)
-- Pre-initializing fault 'FAULT_SLABTOP'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x106601420> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'FAULT_SLABTOP' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:108:preinitialize
-- faultcohesivekin(info)
-- Pre-initializing fault 'FAULT_SLABBOT'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x1066016f0> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'FAULT_SLABBOT' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/feassemble/FIATSimplex.py:152:initialize
-- fiatsimplex(info)
-- Cell geometry:
-- <pylith.feassemble.CellGeometry.GeometryTri3D; proxy of <Swig Object of type 'pylith::feassemble::GeometryTri3D *' at 0x1066019c0> >
-- Vertices:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad pts:
-- [[-1. -1.]
[ 1. -1.]
[-1. 1.]]
-- Quad wts:
-- [ 0.66666667 0.66666667 0.66666667]
-- Basis fns @ quad pts ):
-- [[ 1.00000000e+00 2.77555756e-17 0.00000000e+00]
[ 2.77555756e-17 1.00000000e+00 0.00000000e+00]
[ 5.55111512e-17 5.55111512e-17 1.00000000e+00]]
-- Basis fn derivatives @ quad pts:
-- [[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]
[[-0.5 -0.5]
[ 0.5 0. ]
[ 0. 0.5]]]
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:468:_setupInterfaces
-- implicit(info)
-- Added interface condition 'LOFZ' as an integrator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:188:preinitialize
-- implicit(info)
-- Pre-initializing output.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Problem.py:150:verifyConfiguration
-- timedependent(info)
-- Verifying compatibility of problem configuration.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:117:initialize
-- timedependent(info)
-- Initializing problem.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:500:_initialize
-- implicit(info)
-- Initializing integrators.
WARNING: No initial state given for friction model 'Slip weakening'. Using default value of zero.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:160:initialize
-- faultcohesivekin(info)
-- Initializing fault 'FAULT_SLABTOP'.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/FaultCohesiveKin.py:160:initialize
-- faultcohesivekin(info)
-- Initializing fault 'FAULT_SLABBOT'.
WARNING: No initial state given for friction model 'Slip weakening'. Using default value of zero.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:510:_initialize
-- implicit(info)
-- Initializing constraints.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:516:_initialize
-- implicit(info)
-- Setting up solution output.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:525:_initialize
-- implicit(info)
-- Creating solution field.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:136:initialize
-- implicit(info)
-- Creating other fields.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:159:initialize
-- implicit(info)
-- Creating Jacobian matrix.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:169:initialize
-- implicit(info)
-- Initializing solver.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:131:run
-- timedependent(info)
-- Solving problem.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:137:run
-- timedependent(info)
-- Preparing for prestep with elastic behavior.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:266:prestepElastic
-- implicit(info)
-- Setting constraints.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:582:_reformJacobian
-- implicit(info)
-- Integrating Jacobian operator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:152:run
-- timedependent(info)
-- Computing prestep with elastic behavior.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:607:_reformResidual
-- implicit(info)
-- Integrating residual term in operator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:215:step
-- implicit(info)
-- Solving equations.
0 SNES Function norm < 1.e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 0
SNES Object: 1 MPI processes
type: ngmres
Number of stored past updates: 30
Residual selection: gammaA=2e+00, gammaC=2e+00
Difference restart: epsilonB=1e-01, deltaB=9e-01
maximum iterations=10, maximum function evaluations=30000
tolerances: relative=1e-20, absolute=1e-09, solution=1e-08
total number of linear solver iterations=0
total number of function evaluations=1
SNESLineSearch Object: 1 MPI processes
type: shell
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=1
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:158:run
-- timedependent(info)
-- Finishing prestep with elastic behavior.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:250:poststep
-- implicit(info)
-- Writing solution fields.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:176:run
-- timedependent(info)
-- Main time loop, current time is t=0*s
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:187:run
-- timedependent(info)
-- Preparing to advance solution from time t=0*s to t=3.15576e+07*s.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:186:prestep
-- implicit(info)
-- Setting constraints.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:193:run
-- timedependent(info)
-- Advancing solution from t=0*s to t=3.15576e+07*s.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:607:_reformResidual
-- implicit(info)
-- Integrating residual term in operator.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:215:step
-- implicit(info)
-- Solving equations.
0 SNES Function norm 0.948594
Line search: gnorm after quadratic fit 1.773100006763e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 2.002776491982e+00 lambda=1.0000000000000002e-02
1 SNES Function norm 2.00278
Line search: gnorm after quadratic fit 4.312556301433e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 5.540166796547e+00 lambda=1.0000000000000002e-02
2 SNES Function norm 0.569893
Line search: gnorm after quadratic fit 2.243738258592e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 1.991117051814e+00 lambda=1.0000000000000002e-02
3 SNES Function norm 0.425075
Line search: gnorm after quadratic fit 2.147718836488e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 2.049829198737e+00 lambda=1.0000000000000002e-02
4 SNES Function norm 0.391639
Line search: gnorm after quadratic fit 1.608131254064e+01
Line search: Cubic step no good, shrinking lambda, current gnorm 1.636970057178e+00 lambda=1.0000000000000002e-02
5 SNES Function norm 1.63697
Line search: gnorm after quadratic fit 5.396512956049e+02
Line search: Cubic step no good, shrinking lambda, current gnorm 5.278748490904e+01 lambda=1.0000000000000002e-02
6 SNES Function norm 52.7875
Line search: gnorm after quadratic fit 2.322268614406e+04
Line search: Cubic step no good, shrinking lambda, current gnorm 2.275760046175e+03 lambda=1.0000000000000002e-02
7 SNES Function norm 2275.76
Line search: gnorm after quadratic fit 1.071505794820e+06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.051316283060e+05 lambda=1.0000000000000002e-02
8 SNES Function norm 105132
Line search: gnorm after quadratic fit 5.228582697605e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.135201186511e+06 lambda=1.0000000000000002e-02
9 SNES Function norm 5.1352e+06
Line search: gnorm after quadratic fit 2.686941666753e+09
Line search: Cubic step no good, shrinking lambda, current gnorm 2.641281837342e+08 lambda=1.0000000000000002e-02
10 SNES Function norm 2.64128e+08
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 10
SNES Object: 1 MPI processes
type: ngmres
Number of stored past updates: 30
Residual selection: gammaA=2e+00, gammaC=2e+00
Difference restart: epsilonB=1e-01, deltaB=9e-01
maximum iterations=10, maximum function evaluations=30000
tolerances: relative=1e-20, absolute=1e-09, solution=1e-08
total number of linear solver iterations=0
total number of function evaluations=51
SNESLineSearch Object: 1 MPI processes
type: shell
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=1
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:200:run
-- timedependent(info)
-- Finishing advancing solution from t=0*s to t=3.15576e+07*s.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Implicit.py:250:poststep
-- implicit(info)
-- Writing solution fields.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/faults/Fault.py:189:writeData
-- faultcohesivekin(info)
-- Writing fault data.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/problems/Formulation.py:300:finalize
-- implicit(info)
-- Formulation finalize.
>> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/utils/PetscManager.py:75:finalize
-- petsc(info)
-- Finalizing PETSc.
marcont at OSXmate~/pylith190/modelos/test_viscoelastic_05:
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20141022/ba2336dd/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: figura1.png
Type: image/png
Size: 237037 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20141022/ba2336dd/attachment-0001.png>
More information about the CIG-SHORT
mailing list