[petsc-users] problem with BT linesearcher,

Matthew Knepley

2015-10-20 14:46:28 UTC

On Tue, Oct 20, 2015 at 8:46 AM, Lukasz Kaczmarczyk <

Post by 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

Right here it says that the vector F, your residual, is unassembled. Do you
assemble it at the end
of your routine?

Thanks,

Matt

Post by Lukasz Kaczmarczyk
[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: -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
http://www.gla.ac.uk/schools/engineering/staff/lukaszkaczmarczyk/
http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/index.html

--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener