<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Oct 20, 2015 at 8:46 AM, Lukasz Kaczmarczyk <span dir="ltr"><<a href="mailto:Lukasz.Kaczmarczyk@glasgow.ac.uk" target="_blank">Lukasz.Kaczmarczyk@glasgow.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
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,<br>
<br>
For L2 or CP I get,<br>
Step 1 Time 0.05 final time 3<br>
0 SNES Function norm 2.282177322938e-03<br>
Line search: lambdas = [1, 0.5, 0], fnorms = [8.50089e-17, 0.00114109, 0.00228218]<br>
Line search terminated: lambda = 1, fnorms = 8.23981e-17<br>
1 SNES Function norm 8.239809253872e-17<br>
number of Newton iterations = 1<br>
reduction delta_time = 1.0000e+00 delta_time = 5.0000e-02<br>
out file out_1.h5m<br>
Step 2 Time 0.1 final time 3<br>
0 SNES Function norm 2.282177322938e-03<br>
Line search: lambdas = [1, 0.5, 0], fnorms = [9.15228e-17, 0.00114109, 0.00228218]<br>
Line search terminated: lambda = 1, fnorms = 9.15228e-17<br>
1 SNES Function norm 9.152279314292e-17<br>
number of Newton iterations = 1<br>
reduction delta_time = 1.0000e+00 delta_time = 5.0000e-02<br>
<br>
Fro BT I get,<br>
Step 1 Time 0.05 final time 3<br>
0 SNES Function norm 2.282177322938e-03<br>
Line search: gnorm after quadratic fit 1.141088661469e-03<br>
Line search: Quadratically determined step, lambda=4.9999999999999706e-01<br>
1 SNES Function norm 1.141088661469e-03<br>
[0]PETSC ERROR: --------------------- MoFEM Error Message---------------------------------------------------------------------------<br>
[0]PETSC ERROR: MoFEM version 0.3.0<br>
[0]PETSC ERROR: MoFEM git commit id a0b356adadb599d17aa11ae2e6d354355111ee7d<br>
[0]PETSC ERROR: See <a href="http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/guidelines_bug_reporting.html" rel="noreferrer" target="_blank">http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/guidelines_bug_reporting.html</a> for bug reporting.<br>
[0]PETSC ERROR: See <a href="http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/faq_and_bugs.html" rel="noreferrer" target="_blank">http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/faq_and_bugs.html</a> for trouble shooting.<br>
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>
[0]PETSC ERROR: Object is in wrong state<br>
[0]PETSC ERROR: Not for unassembled vector<br>
[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>
[0]PETSC ERROR: Petsc Development GIT revision: v3.6-41-g8f9ce63 GIT Date: 2015-06-20 12:33:17 -0500<br>
[0]PETSC ERROR: ./small_strain_plasticity_problem on a arch-darwin-c-opt named <a href="http://dh23-177.eng.gla.ac.uk" rel="noreferrer" target="_blank">dh23-177.eng.gla.ac.uk</a> by likask Tue Oct 20 14:21:15 2015<br>
[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<br>
[0]PETSC ERROR: #1 VecCopy() line 1630 in /opt/petsc/src/vec/vec/interface/vector.c<br>
[0]PETSC ERROR: #2 KSPSolve_FGMRES() line 290 in /opt/petsc/src/ksp/ksp/impls/gmres/fgmres/fgmres.c<br>
[0]PETSC ERROR: #3 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c<br></blockquote><div><br></div><div>Right here it says that the vector F, your residual, is unassembled. Do you assemble it at the end</div><div>of your routine?</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
[0]PETSC ERROR: #4 SNESSolve_NEWTONLS() line 233 in /opt/petsc/src/snes/impls/ls/ls.c<br>
[0]PETSC ERROR: #5 SNESSolve() line 3894 in /opt/petsc/src/snes/interface/snes.c<br>
[0]PETSC ERROR: #6 main() line 421 in /Users/likask/MyBuild/mofem-cephas/mofem/users_modules/small_strain_plasticity/small_strain_plasticity_problem.cpp<br>
[0]PETSC ERROR: PETSc Option Table entries:<br>
[0]PETSC ERROR: -ksp_type fgmres<br>
[0]PETSC ERROR: -my_H 0.1<br>
[0]PETSC ERROR: -my_delta_time 0.05<br>
[0]PETSC ERROR: -my_file plate_with_hole.h5m<br>
[0]PETSC ERROR: -my_final_time 3<br>
[0]PETSC ERROR: -my_load_history plate_with_hole_history.txt<br>
[0]PETSC ERROR: -my_order 1<br>
[0]PETSC ERROR: -my_poisson_ratio 0.25<br>
<br>
Kind regards,<br>
Lukasz Kaczmarczyk<br>
<br>
"Prediction is very difficult, especially about the future." -- Niels Bohr<br>
“… there are dark corners in the Bourne shell, and people use all of them." --Chet Ramey<br>
<br>
HomePage:<br>
<a href="http://www.gla.ac.uk/schools/engineering/staff/lukaszkaczmarczyk/" rel="noreferrer" target="_blank">http://www.gla.ac.uk/schools/engineering/staff/lukaszkaczmarczyk/</a><br>
<br>
MoFEM @mofemjoseph<br>
<a href="http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/index.html" rel="noreferrer" target="_blank">http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/index.html</a><br>
<br>
<br>
<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>