search

idaholab / moose

Multiphysics Object Oriented Simulation Environment
https://www.mooseframework.org
GNU Lesser General Public License v2.1
1.75k stars 1.05k forks source link

Test: mortar/2d.equalgradient #10924

Open fdkong opened 6 years ago

fdkong commented 6 years ago

Rationale

The test is sensitive to the change of floating point operations because the system is singular (nearly).

Description

../../../moose_test-opt -i equalgradient.i -pc_type lu -ksp_error_if_not_converged

Time Step  0, time = 0
                dt = 0

Time Step  1, time = 0.05
                dt = 0.05
 0 Nonlinear |R| = 1.118034e-01
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Zero pivot in LU factorization: http://www.mcs.anl.gov/petsc/documentation/faq.html#zeropivot
[0]PETSC ERROR: Zero pivot row 36 value 2.10942e-15 tolerance 2.22045e-14
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.8.3, unknown 
[0]PETSC ERROR: ../../../moose_test-opt on a arch-darwin-c-debug-maint named FN604208 by kongf Wed Feb 28 09:52:32 2018
[0]PETSC ERROR: Configure options --with-debugging=yes --with-shared-libraries=1 --with-mpi=1 --download-fblaslapack=1 --download-metis=1 --download-parmetis=1 --download-superlu_dist=1 --download-hypre=1 --download-mumps=1 --download-scalapack=1 PETSC_ARCH=arch-darwin-c-debug-maint
[0]PETSC ERROR: #1 MatPivotCheck_none() line 713 in /Users/kongf/projects/petsc/include/petsc/private/matimpl.h
[0]PETSC ERROR: #2 MatPivotCheck() line 730 in /Users/kongf/projects/petsc/include/petsc/private/matimpl.h
[0]PETSC ERROR: #3 MatLUFactorNumeric_SeqAIJ() line 558 in /Users/kongf/projects/petsc/src/mat/impls/aij/seq/aijfact.c
[0]PETSC ERROR: #4 MatLUFactorNumeric() line 3075 in /Users/kongf/projects/petsc/src/mat/interface/matrix.c
[0]PETSC ERROR: #5 PCSetUp_LU() line 131 in /Users/kongf/projects/petsc/src/ksp/pc/impls/factor/lu/lu.c
[0]PETSC ERROR: #6 PCSetUp() line 924 in /Users/kongf/projects/petsc/src/ksp/pc/interface/precon.c
[0]PETSC ERROR: #7 KSPSetUp() line 381 in /Users/kongf/projects/petsc/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: #8 KSPSolve() line 612 in /Users/kongf/projects/petsc/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: #9 SNESSolve_NEWTONLS() line 224 in /Users/kongf/projects/petsc/src/snes/impls/ls/ls.c
[0]PETSC ERROR: #10 SNESSolve() line 4131 in /Users/kongf/projects/petsc/src/snes/interface/snes.c
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1

../../../moose_test-opt -i equalgradient.i -pc_type svd -pc_svd_monitor

Time Step  0, time = 0
                dt = 0

Time Step  1, time = 0.05
                dt = 0.05
 0 Nonlinear |R| = 1.118034e-01
      SVD: condition number 9.392908993773e+16, 1 of 40 singular values are (nearly) zero
      SVD: smallest singular values: 3.996083156132e-17 2.179395859968e-01 2.844113774607e-01 3.296969146968e-01 3.670726968780e-01
      SVD: largest singular values : 3.345536541528e+00 3.619271475718e+00 3.643751373744e+00 3.659984648281e+00 3.753484541710e+00
 1 Nonlinear |R| = 2.847296e-16
 Solve Converged!

Outlier Variable Residual Norms:
  u: 2.771173e-16

Time Step  2, time = 0.1
                dt = 0.05
 0 Nonlinear |R| = 1.203254e-01
      SVD: condition number 9.392908993773e+16, 1 of 40 singular values are (nearly) zero
      SVD: smallest singular values: 3.996083156132e-17 2.179395859968e-01 2.844113774607e-01 3.296969146968e-01 3.670726968780e-01
      SVD: largest singular values : 3.345536541528e+00 3.619271475718e+00 3.643751373744e+00 3.659984648281e+00 3.753484541710e+00
 1 Nonlinear |R| = 3.164265e-16
 Solve Converged!

Outlier Variable Residual Norms:
  u: 3.106631e-16

Time Step  3, time = 0.15
                dt = 0.05
 0 Nonlinear |R| = 1.294912e-01
      SVD: condition number 9.392908993773e+16, 1 of 40 singular values are (nearly) zero
      SVD: smallest singular values: 3.996083156132e-17 2.179395859968e-01 2.844113774607e-01 3.296969146968e-01 3.670726968780e-01
      SVD: largest singular values : 3.345536541528e+00 3.619271475718e+00 3.643751373744e+00 3.659984648281e+00 3.753484541710e+00
 1 Nonlinear |R| = 3.154318e-16
 Solve Converged!

Outlier Variable Residual Norms:
  u: 3.096645e-16

Impact

Affect the test only!

fdkong commented 6 years ago

@dschwen, could you look at this?

dschwen commented 6 years ago

Should I try adding a nonzero shift?

On Wed, Feb 28, 2018, 10:05 AM Fande Kong notifications@github.com wrote:

@dschwen https://github.com/dschwen, could you look at this?

— You are receiving this because you were mentioned.

Reply to this email directly, view it on GitHub https://github.com/idaholab/moose/issues/10924#issuecomment-369307159, or mute the thread https://github.com/notifications/unsubscribe-auth/AAMWPgR1LKysRtWT-8acIxsEiC1Qr-n_ks5tZYY4gaJpZM4SW_sb .