[petsc-users] problem with BT linesearcher,

[Matthew Knepley]

On Tue, Oct 20, 2015 at 8:46 AM, Lukasz Kaczmarczyk <
Lukasz.Kaczmarczyk at glasgow.ac.uk> wrote:

> 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


> [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
>
>
>
>
>


-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151020/b8e3a7cc/attachment.html>