[petsc-users] arkimex rejecting all dt in petsc 3.6.0
Italo Tasso
italo at tasso.com.br
Thu Jun 18 05:48:29 CDT 2015
I just upgraded to 3.6.0 and my code stopped working. All dt are rejected.
I used the same configure line, same code, same everything.
With 3.5.4 I get:
0 TS dt 1e-06 time 0
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 6.107056905987e-03
2 SNES Function norm 1.483881932064e-10
3 SNES Function norm 9.122873794272e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.790429171165e+04
1 SNES Function norm 7.289068747803e-04
2 SNES Function norm 8.227639633330e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.154356516184e+05
1 SNES Function norm 2.309925413255e-03
2 SNES Function norm 6.382141981406e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 accepted t=0 + 1.000e-06
wlte=0.000654 family='arkimex' scheme=0:'3' dt=1.000e-05
With 3.6.0 I get:
0 TS dt 1e-06 time 0
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 6.107056925316e-03
2 SNES Function norm 1.519319591792e-10
3 SNES Function norm 9.070104116945e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.790429171165e+04
1 SNES Function norm 6.942541792651e-04
2 SNES Function norm 8.458781909516e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.154356516184e+05
1 SNES Function norm 2.287202942961e-03
2 SNES Function norm 6.585201377573e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 1.000e-06 wlte=
324 family='arkimex' scheme=0:'3' dt=1.311e-07
Any ideas? I attached the full output.
Options I use:
-ts_view -ts_type arkimex -ts_arkimex_fully_implicit -ts_adapt_monitor
-ts_monitor -snes_monitor -snes_converged_reason -ksp_type preonly -pc_type
lu -pc_factor_mat_solver_package mumps -snes_rtol 0 -snes_atol 1e-10
-snes_stol 0
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150618/4af9d52e/attachment.html>
-------------- next part --------------
0 TS dt 1e-06 time 0
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 6.107056905987e-03
2 SNES Function norm 1.483881932064e-10
3 SNES Function norm 9.122873794272e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.790429171165e+04
1 SNES Function norm 7.289068747803e-04
2 SNES Function norm 8.227639633330e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.154356516184e+05
1 SNES Function norm 2.309925413255e-03
2 SNES Function norm 6.382141981406e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 accepted t=0 + 1.000e-06 wlte=0.000654 family='arkimex' scheme=0:'3' dt=1.000e-05
1 TS dt 1e-05 time 1e-06
TS Object: 1 MPI processes
type: arkimex
maximum steps=1000000000
maximum time=1e-06
total number of nonlinear solver iterations=7
total number of nonlinear solve failures=0
total number of linear solver iterations=7
total number of rejected steps=0
ARK IMEX 3
Stiff abscissa ct = 0.000000 0.871733 0.600000 1.000000
Stiffly accurate: yes
Explicit first stage: yes
FSAL property: yes
Nonstiff abscissa c = 0.000000 0.871733 0.600000 1.000000
TSAdapt Object: 1 MPI processes
type: basic
number of candidates 1
Basic: clip fastest decrease 0.1, fastest increase 10
Basic: safety factor 0.9, extra factor after step rejection 0.5
SNES Object: 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=0, absolute=1e-10, solution=0
total number of linear solver iterations=2
total number of function evaluations=3
SNESLineSearch Object: 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
matrix ordering: nd
factor fill ratio given 0, needed 0
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=30603, cols=30603
package used to perform factorization: mumps
total: nonzeros=6.03922e+06, allocated nonzeros=6.03922e+06
total number of mallocs used during MatSetValues calls =0
MUMPS run parameters:
SYM (matrix type): 0
PAR (host participation): 1
ICNTL(1) (output for error): 6
ICNTL(2) (output of diagnostic msg): 0
ICNTL(3) (output for global info): 0
ICNTL(4) (level of printing): 0
ICNTL(5) (input mat struct): 0
ICNTL(6) (matrix prescaling): 7
ICNTL(7) (sequentia matrix ordering):7
ICNTL(8) (scalling strategy): 77
ICNTL(10) (max num of refinements): 0
ICNTL(11) (error analysis): 0
ICNTL(12) (efficiency control): 1
ICNTL(13) (efficiency control): 0
ICNTL(14) (percentage of estimated workspace increase): 20
ICNTL(18) (input mat struct): 0
ICNTL(19) (Shur complement info): 0
ICNTL(20) (rhs sparse pattern): 0
ICNTL(21) (somumpstion struct): 0
ICNTL(22) (in-core/out-of-core facility): 0
ICNTL(23) (max size of memory can be allocated locally):0
ICNTL(24) (detection of null pivot rows): 0
ICNTL(25) (computation of a null space basis): 0
ICNTL(26) (Schur options for rhs or solution): 0
ICNTL(27) (experimental parameter): -8
ICNTL(28) (use parallel or sequential ordering): 1
ICNTL(29) (parallel ordering): 0
ICNTL(30) (user-specified set of entries in inv(A)): 0
ICNTL(31) (factors is discarded in the solve phase): 0
ICNTL(33) (compute determinant): 0
CNTL(1) (relative pivoting threshold): 0.01
CNTL(2) (stopping criterion of refinement): 1.49012e-08
CNTL(3) (absomumpste pivoting threshold): 0
CNTL(4) (vamumpse of static pivoting): -1
CNTL(5) (fixation for null pivots): 0
RINFO(1) (local estimated flops for the elimination after analysis):
[0] 1.03709e+09
RINFO(2) (local estimated flops for the assembly after factorization):
[0] 1.07576e+07
RINFO(3) (local estimated flops for the elimination after factorization):
[0] 1.03736e+09
INFO(15) (estimated size of (in MB) MUMPS internal data for running numerical factorization):
[0] 73
INFO(16) (size of (in MB) MUMPS internal data used during numerical factorization):
[0] 73
INFO(23) (num of pivots eliminated on this processor after factorization):
[0] 30603
RINFOG(1) (global estimated flops for the elimination after analysis): 1.03709e+09
RINFOG(2) (global estimated flops for the assembly after factorization): 1.07576e+07
RINFOG(3) (global estimated flops for the elimination after factorization): 1.03736e+09
(RINFOG(12) RINFOG(13))*2^INFOG(34) (determinant): (0,0)*(2^0)
INFOG(3) (estimated real workspace for factors on all processors after analysis): 6039225
INFOG(4) (estimated integer workspace for factors on all processors after analysis): 352584
INFOG(5) (estimated maximum front size in the complete tree): 456
INFOG(6) (number of nodes in the complete tree): 2491
INFOG(7) (ordering option effectively use after analysis): 5
INFOG(8) (structural symmetry in percent of the permuted matrix after analysis): 100
INFOG(9) (total real/complex workspace to store the matrix factors after factorization): 6040215
INFOG(10) (total integer space store the matrix factors after factorization): 352602
INFOG(11) (order of largest frontal matrix after factorization): 456
INFOG(12) (number of off-diagonal pivots): 1207
INFOG(13) (number of delayed pivots after factorization): 9
INFOG(14) (number of memory compress after factorization): 0
INFOG(15) (number of steps of iterative refinement after solution): 0
INFOG(16) (estimated size (in MB) of all MUMPS internal data for factorization after analysis: value on the most memory consuming processor): 73
INFOG(17) (estimated size of all MUMPS internal data for factorization after analysis: sum over all processors): 73
INFOG(18) (size of all MUMPS internal data allocated during factorization: value on the most memory consuming processor): 73
INFOG(19) (size of all MUMPS internal data allocated during factorization: sum over all processors): 73
INFOG(20) (estimated number of entries in the factors): 6039225
INFOG(21) (size in MB of memory effectively used during factorization - value on the most memory consuming processor): 63
INFOG(22) (size in MB of memory effectively used during factorization - sum over all processors): 63
INFOG(23) (after analysis: value of ICNTL(6) effectively used): 0
INFOG(24) (after analysis: value of ICNTL(12) effectively used): 1
INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
INFOG(28) (after factorization: number of null pivots encountered): 0
INFOG(29) (after factorization: effective number of entries in the factors (sum over all processors)): 6040215
INFOG(30, 31) (after solution: size in Mbytes of memory used during solution phase): 67, 67
INFOG(32) (after analysis: type of analysis done): 1
INFOG(33) (value used for ICNTL(8)): 7
INFOG(34) (exponent of the determinant if determinant is requested): 0
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=30603, cols=30603, bs=3
total: nonzeros=815409, allocated nonzeros=815409
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 10201 nodes, limit used is 5
CONVERGED_TIME at time 1e-06 after 1 steps
-------------- next part --------------
0 TS dt 1e-06 time 0
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 6.107056925316e-03
2 SNES Function norm 1.519319591792e-10
3 SNES Function norm 9.070104116945e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.790429171165e+04
1 SNES Function norm 6.942541792651e-04
2 SNES Function norm 8.458781909516e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.154356516184e+05
1 SNES Function norm 2.287202942961e-03
2 SNES Function norm 6.585201377573e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 1.000e-06 wlte= 324 family='arkimex' scheme=0:'3' dt=1.311e-07
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 1.030893399688e-02
2 SNES Function norm 1.227244086383e-10
3 SNES Function norm 8.680479900786e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.850186237744e+04
1 SNES Function norm 4.435591601555e-04
2 SNES Function norm 8.317684980141e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.165153426209e+05
1 SNES Function norm 1.813148503621e-03
2 SNES Function norm 5.613694043108e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 1.311e-07 wlte= 320 family='arkimex' scheme=0:'3' dt=1.724e-08
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 4.263053096085e-02
2 SNES Function norm 6.252919451896e-10
3 SNES Function norm 9.058384623510e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.858128091677e+04
1 SNES Function norm 3.963115544432e-04
2 SNES Function norm 8.278281256019e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166564177803e+05
1 SNES Function norm 9.631229966553e-04
2 SNES Function norm 5.708329543623e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 1.724e-08 wlte= 320 family='arkimex' scheme=0:'3' dt=2.269e-09
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 1.079124564076e-01
2 SNES Function norm 1.130538110070e-09
3 SNES Function norm 8.991479075319e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.859174767183e+04
1 SNES Function norm 4.321261456929e-04
2 SNES Function norm 8.218414643337e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166749679782e+05
1 SNES Function norm 6.070076055571e-04
2 SNES Function norm 5.569380923621e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 2.269e-09 wlte= 320 family='arkimex' scheme=0:'3' dt=2.987e-10
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 2.312307641993e-02
2 SNES Function norm 2.026987092062e-10
3 SNES Function norm 8.926723703962e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.859312557607e+04
1 SNES Function norm 4.828786679744e-04
2 SNES Function norm 8.341615291368e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166774092974e+05
1 SNES Function norm 1.835403825645e-03
2 SNES Function norm 5.577036922069e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 2.987e-10 wlte= 320 family='arkimex' scheme=0:'3' dt=3.931e-11
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 2.351722737234e-02
2 SNES Function norm 2.529720083595e-10
3 SNES Function norm 8.892352147112e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.859330694478e+04
1 SNES Function norm 4.885408570662e-04
2 SNES Function norm 8.327396926987e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166777306269e+05
1 SNES Function norm 1.804733305674e-03
2 SNES Function norm 5.388931877150e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 3.931e-11 wlte= 320 family='arkimex' scheme=0:'3' dt=5.174e-12
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 9.516392283907e-02
2 SNES Function norm 1.084610044563e-09
3 SNES Function norm 9.117232285273e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.859333081724e+04
1 SNES Function norm 1.951960746471e-04
2 SNES Function norm 8.022315535248e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166777729214e+05
1 SNES Function norm 1.707739879535e-03
2 SNES Function norm 6.069805583058e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 5.174e-12 wlte= 320 family='arkimex' scheme=0:'3' dt=6.810e-13
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 9.238561718262e-02
2 SNES Function norm 1.201780471186e-09
3 SNES Function norm 8.749890804482e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.859333395939e+04
1 SNES Function norm 4.010933191742e-04
2 SNES Function norm 8.258990222529e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166777784888e+05
1 SNES Function norm 1.833489380639e-03
2 SNES Function norm 5.339548154287e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 6.810e-13 wlte= 320 family='arkimex' scheme=0:'3' dt=8.964e-14
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 7.725363551151e-02
2 SNES Function norm 3.035001046151e-10
3 SNES Function norm 9.104321363901e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 3
0 SNES Function norm 7.859333437279e+04
1 SNES Function norm 4.404284856037e-04
2 SNES Function norm 8.005734930179e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166777792253e+05
1 SNES Function norm 9.764901758589e-04
2 SNES Function norm 5.365314630014e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 8.964e-14 wlte= 320 family='arkimex' scheme=0:'3' dt=1.180e-14
0 SNES Function norm 2.549981005316e+05
1 SNES Function norm 1.267027453118e-05
2 SNES Function norm 8.908505647078e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 7.859333442583e+04
1 SNES Function norm 4.668956513888e-04
2 SNES Function norm 8.348862687250e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
0 SNES Function norm 1.166777793501e+05
1 SNES Function norm 9.432068723061e-04
2 SNES Function norm 5.390404010769e-11
Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 2
TSAdapt 'basic': step 0 rejected t=0 + 1.180e-14 wlte= 320 family='arkimex' scheme=0:'3' dt=1.553e-15
DIVERGED_STEP_REJECTED at time 0 after 0 steps
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR:
[0]PETSC ERROR: TSStep has failed due to DIVERGED_STEP_REJECTED
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.6.0, Jun, 09, 2015
[0]PETSC ERROR: ./planb on a arch-linux2-c-opt named localhost.localdomain by root Thu Jun 18 07:37:39 2015
[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-fblaslapack --download-mpich --download-mumps --download-scalapack --download-parmetis --download-metis --with-debugging=0 COPTFLAGS="-Ofast -march=native -mtune=native" CXXOPTFLAGS="-Ofast -march=native -mtune=native" FOPTFLAGS="-Ofast -march=native -mtune=native"
[0]PETSC ERROR: #1 TSStep() line 3107 in /root/petsc/petsc-3.6.0/src/ts/interface/ts.c
[0]PETSC ERROR: #2 TSSolve() line 3282 in /root/petsc/petsc-3.6.0/src/ts/interface/ts.c
More information about the petsc-users
mailing list