[petsc-users] problem with BT linesearcher,

[Lukasz Kaczmarczyk]

Hello,

I have problem with BT line-searcher. I can do calculations with L2, CP line-searchers, however I get error for BT. Error is at first step. This is plasticity which is initially linear after some time nonlinearity kicks in. Could you advise where to search fro problem,

For L2  or CP I get,
Step 1 Time 0.05 final time   3
  0 SNES Function norm 2.282177322938e-03 
      Line search: lambdas = [1, 0.5, 0], fnorms = [8.50089e-17, 0.00114109, 0.00228218]
      Line search terminated: lambda = 1, fnorms = 8.23981e-17
  1 SNES Function norm 8.239809253872e-17 
number of Newton iterations = 1
reduction delta_time = 1.0000e+00 delta_time = 5.0000e-02
out file out_1.h5m
Step 2 Time 0.1 final time   3
  0 SNES Function norm 2.282177322938e-03 
      Line search: lambdas = [1, 0.5, 0], fnorms = [9.15228e-17, 0.00114109, 0.00228218]
      Line search terminated: lambda = 1, fnorms = 9.15228e-17
  1 SNES Function norm 9.152279314292e-17 
number of Newton iterations = 1
reduction delta_time = 1.0000e+00 delta_time = 5.0000e-02

Fro BT I get,
Step 1 Time 0.05 final time   3
  0 SNES Function norm 2.282177322938e-03 
      Line search: gnorm after quadratic fit 1.141088661469e-03
      Line search: Quadratically determined step, lambda=4.9999999999999706e-01
  1 SNES Function norm 1.141088661469e-03 
[0]PETSC ERROR: --------------------- MoFEM Error Message---------------------------------------------------------------------------
[0]PETSC ERROR: MoFEM version 0.3.0
[0]PETSC ERROR: MoFEM git commit id a0b356adadb599d17aa11ae2e6d354355111ee7d
[0]PETSC ERROR: See http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/guidelines_bug_reporting.html for bug reporting.
[0]PETSC ERROR: See http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/faq_and_bugs.html for trouble shooting.
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Object is in wrong state
[0]PETSC ERROR: Not for unassembled vector
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.6-41-g8f9ce63  GIT Date: 2015-06-20 12:33:17 -0500
[0]PETSC ERROR: ./small_strain_plasticity_problem on a arch-darwin-c-opt named dh23-177.eng.gla.ac.uk by likask Tue Oct 20 14:21:15 2015
[0]PETSC ERROR: Configure options --download-blacs=1 --download-fblaslapack=1 --download-hdf5=1 --download-hypre=1 --download-metis=1 --download-moab=1 --download-mumps=1 --download-netcdf=netcdf-4.3.3.1.tar.gz --download-openmpi=1 --download-parmetis=1 --download-scalapack=1 --download-superlu_dist=1 --with-debugging=0 PETSC_ARCH=arch-darwin-c-opt
[0]PETSC ERROR: #1 VecCopy() line 1630 in /opt/petsc/src/vec/vec/interface/vector.c
[0]PETSC ERROR: #2 KSPSolve_FGMRES() line 290 in /opt/petsc/src/ksp/ksp/impls/gmres/fgmres/fgmres.c
[0]PETSC ERROR: #3 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: #4 SNESSolve_NEWTONLS() line 233 in /opt/petsc/src/snes/impls/ls/ls.c
[0]PETSC ERROR: #5 SNESSolve() line 3894 in /opt/petsc/src/snes/interface/snes.c
[0]PETSC ERROR: #6 main() line 421 in /Users/likask/MyBuild/mofem-cephas/mofem/users_modules/small_strain_plasticity/small_strain_plasticity_problem.cpp
[0]PETSC ERROR: PETSc Option Table entries:
[0]PETSC ERROR: -ksp_type fgmres
[0]PETSC ERROR: -my_H 0.1
[0]PETSC ERROR: -my_delta_time 0.05
[0]PETSC ERROR: -my_file plate_with_hole.h5m
[0]PETSC ERROR: -my_final_time 3
[0]PETSC ERROR: -my_load_history plate_with_hole_history.txt
[0]PETSC ERROR: -my_order 1
[0]PETSC ERROR: -my_poisson_ratio 0.25

Kind regards,
Lukasz Kaczmarczyk

"Prediction is very difficult, especially about the future." -- Niels Bohr
“… there are dark corners in the Bourne shell, and people use all of them." --Chet Ramey

HomePage:
http://www.gla.ac.uk/schools/engineering/staff/lukaszkaczmarczyk/

MoFEM @mofemjoseph 
http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/index.html