[petsc-users] Nonconforming object sizes using TAO (TAOBNTR)

Alexis Marboeuf alexis.marboeuf at hotmail.fr
Fri Jan 13 19:21:10 CST 2023


Hi Matt,

Here is the output from Petsc when I view the TAO solver:


Tao Object: (Damage_) 4 MPI processes

  type: bntr

    Tao Object: (Damage_tao_bnk_cg_) 4 MPI processes

      type: bncg

        CG Type: ssml_bfgs

        Skipped Stepdirection Updates: 0

        Scaled gradient steps: 0

        Pure gradient steps: 0

        Not a descent direction: 0

        Line search fails: 0

        Matrix has not been preallocated yet

      TaoLineSearch Object: (Damage_tao_bnk_cg_) 4 MPI processes

        type: more-thuente

        maximum function evaluations=30

        tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9

        total number of function evaluations=0

        total number of gradient evaluations=0

        total number of function/gradient evaluations=0

        Termination reason: 0

      Active Set subset type: subvec

      convergence tolerances: gatol=1e-08,       steptol=0.,       gttol=0.

      Residual in Function/Gradient:=0.

      Objective value=0.

      total number of iterations=0,                              (max: 2000)

      Solver terminated: 0

    Rejected BFGS updates: 0

    CG steps: 0

    Newton steps: 11

    BFGS steps: 0

    Scaled gradient steps: 0

    Gradient steps: 0

    KSP termination reasons:

      atol: 4

      rtol: 0

      ctol: 7

      negc: 0

      dtol: 0

      iter: 0

      othr: 0

  TaoLineSearch Object: (Damage_) 4 MPI processes

    type: more-thuente

    maximum function evaluations=30

    tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9

    total number of function evaluations=0

    total number of gradient evaluations=0

    total number of function/gradient evaluations=0

    using variable bounds

    Termination reason: 0

  KSP Object: (Damage_tao_bnk_) 4 MPI processes

    type: stcg

    maximum iterations=10000, nonzero initial guess

    tolerances:  relative=1e-08, absolute=1e-08, divergence=1e+10

    left preconditioning

    using UNPRECONDITIONED norm type for convergence test

  PC Object: (Damage_tao_bnk_) 4 MPI processes

    type: lmvm

    Mat Object: (Damage_tao_bnk_pc_lmvm_) 4 MPI processes

      type: lmvmbfgs

      rows=30634, cols=30634

        Scale type: DIAGONAL

        Scale history: 1

        Scale params: alpha=1., beta=0.5, rho=1.

        Convex factors: phi=0., theta=0.125

        Max. storage: 5

        Used storage: 5

        Number of updates: 11

        Number of rejects: 0

        Number of resets: 0

        Mat Object: (Damage_tao_bnk_pc_lmvm_J0_) 4 MPI processes

          type: lmvmdiagbroyden

          rows=30634, cols=30634

            Scale history: 1

            Scale params: alpha=1., beta=0.5, rho=1.

            Convex factor: theta=0.125

            Max. storage: 1

            Used storage: 1

            Number of updates: 11

            Number of rejects: 0

            Number of resets: 0

    linear system matrix = precond matrix:

    Mat Object: 4 MPI processes

      type: mpiaij

      rows=468, cols=468

      total: nonzeros=2932, allocated nonzeros=2932

      total number of mallocs used during MatSetValues calls=0

        not using I-node (on process 0) routines

  total KSP iterations: 103

  Active Set subset type: subvec

  convergence tolerances: gatol=0.0001,   steptol=0.,   gttol=1e-05

  Residual in Function/Gradient:=9.11153e-05

  Objective value=0.00665458

  total number of iterations=11,                          (max: 50)

  total number of function evaluations=17,                  max: -1

  total number of gradient evaluations=13,                  max: -1

  total number of Hessian evaluations=12

  Solution converged:    ||g(X)|| <= gatol

Thanks again for your help!
Alexis

________________________________
De : Matthew Knepley <knepley at gmail.com>
Envoyé : samedi 14 janvier 2023 01:38
À : Alexis Marboeuf <alexis.marboeuf at hotmail.fr>
Cc : petsc-users at mcs.anl.gov <petsc-users at mcs.anl.gov>
Objet : Re: [petsc-users] Nonconforming object sizes using TAO (TAOBNTR)

On Fri, Jan 13, 2023 at 11:22 AM Alexis Marboeuf <alexis.marboeuf at hotmail.fr<mailto:alexis.marboeuf at hotmail.fr>> wrote:
Hi all,

In a variational approach of brittle fracture setting, I try to solve a bound constraint minimization problem using TAO. I checkout on the main branch of Petsc. Minimization with respect to the bounded variable (damage) is achieved through the Bounded Newton Trust Region (TAOBNTR). All other TAO parameters are set by default. On a Linux machine, I get the following error with a 4 processors run:

Can you view the solver?

  Thanks,

    Matt


[3]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[3]PETSC ERROR: Nonconforming object sizes
[3]PETSC ERROR: Preconditioner number of local rows 1122 does not equal input vector size 1161
[3]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.
[3]PETSC ERROR: Petsc Development GIT revision: v3.18.3-342-gdab44c92d91  GIT Date: 2023-01-04 13:37:04 +0000
[3]PETSC ERROR: /home/marboeua/Developpement/mef90/arch-darwin-c/bin/vDefTAO on a arch-darwin-c named bb01 by marboeua Thu Jan 12 16:55:18 2023
[2]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[3]PETSC ERROR: Configure options --FFLAGS=-ffree-line-length-none --COPTFLAGS="-O3 -march=znver3 -g" --CXXOPTFLAGS="-O3 -march=znver3 -g" --FOPTFLAGS="-O3 -march=znver3 -g" --download-fblaslapack=1 --download-mumps=1 --download-chaco=1 --download-exodusii=1 --download-hypre=1 --download-ml=1 --download-triangle --download-scalapack=1 --download-superlu=1 --download-sowing=1 --download-sowing-cc=/opt/rh/devtoolset-9/root/usr/bin/gcc --download-sowing-cxx=/opt/rh/devtoolset-9/root/usr/bin/g++ --download-sowing-cpp=/opt/rh/devtoolset-9/root/usr/bin/cpp --download-sowing-cxxcpp=/opt/rh/devtoolset-9/root/usr/bin/cpp --download-yaml=1 --download-bison=1 --download-hdf5=1 --download-metis=1 --download-parmetis=1 --download-netcdf=1 --download-pnetcdf=1 --download-zlib=1 --with-cmake=1 --with-debugging=0 --with-mpi-dir=/opt/HPC/mvapich2/2.3.7-gcc11.2.1 --with-ranlib=ranlib --with-shared-libraries=1 --with-sieve=1 --download-p4est=1 --with-pic --with-mpiexec=srun --with-x11=0 PETSC_ARCH=arch-darwin-c
[3]PETSC ERROR: #1 PCApply() at /1/home/marboeua/Developpement/petsc/src/ksp/pc/interface/precon.c:434
[3]PETSC ERROR: #2 KSP_PCApply() at /home/marboeua/Developpement/petsc/include/petsc/private/kspimpl.h:380
[3]PETSC ERROR: #3 KSPCGSolve_STCG() at /1/home/marboeua/Developpement/petsc/src/ksp/ksp/impls/cg/stcg/stcg.c:76
[3]PETSC ERROR: #4 KSPSolve_Private() at /1/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:898
[3]PETSC ERROR: #5 KSPSolve() at /1/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:1070
[3]PETSC ERROR: #6 TaoBNKComputeStep() at /1/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bnk.c:459
[3]PETSC ERROR: #7 TaoSolve_BNTR() at /1/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bntr.c:138
[3]PETSC ERROR: #8 TaoSolve() at /1/home/marboeua/Developpement/petsc/src/tao/interface/taosolver.c:177
[2]PETSC ERROR: Nonconforming object sizes
[2]PETSC ERROR: Preconditioner number of local rows 1229 does not equal input vector size 1254
[2]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.
[2]PETSC ERROR: Petsc Development GIT revision: v3.18.3-342-gdab44c92d91  GIT Date: 2023-01-04 13:37:04 +0000
[2]PETSC ERROR: /home/marboeua/Developpement/mef90/arch-darwin-c/bin/vDefTAO on a arch-darwin-c named bb01 by marboeua Thu Jan 12 16:55:18 2023
[2]PETSC ERROR: Configure options --FFLAGS=-ffree-line-length-none --COPTFLAGS="-O3 -march=znver3 -g" --CXXOPTFLAGS="-O3 -march=znver3 -g" --FOPTFLAGS="-O3 -march=znver3 -g" --download-fblaslapack=1 --download-mumps=1 --download-chaco=1 --download-exodusii=1 --download-hypre=1 --download-ml=1 --download-triangle --download-scalapack=1 --download-superlu=1 --download-sowing=1 --download-sowing-cc=/opt/rh/devtoolset-9/root/usr/bin/gcc --download-sowing-cxx=/opt/rh/devtoolset-9/root/usr/bin/g++ --download-sowing-cpp=/opt/rh/devtoolset-9/root/usr/bin/cpp --download-sowing-cxxcpp=/opt/rh/devtoolset-9/root/usr/bin/cpp --download-yaml=1 --download-bison=1 --download-hdf5=1 --download-metis=1 --download-parmetis=1 --download-netcdf=1 --download-pnetcdf=1 --download-zlib=1 --with-cmake=1 --with-debugging=0 --with-mpi-dir=/opt/HPC/mvapich2/2.3.7-gcc11.2.1 --with-ranlib=ranlib --with-shared-libraries=1 --with-sieve=1 --download-p4est=1 --with-pic --with-mpiexec=srun --with-x11=0 PETSC_ARCH=arch-darwin-c
[2]PETSC ERROR: #1 PCApply() at /1/home/marboeua/Developpement/petsc/src/ksp/pc/interface/precon.c:434
[2]PETSC ERROR: #2 KSP_PCApply() at /home/marboeua/Developpement/petsc/include/petsc/private/kspimpl.h:380
[2]PETSC ERROR: #3 KSPCGSolve_STCG() at /1/home/marboeua/Developpement/petsc/src/ksp/ksp/impls/cg/stcg/stcg.c:76
[2]PETSC ERROR: #4 KSPSolve_Private() at /1/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:898
[2]PETSC ERROR: #5 KSPSolve() at /1/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:1070
[2]PETSC ERROR: #6 TaoBNKComputeStep() at /1/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bnk.c:459
[2]PETSC ERROR: #7 TaoSolve_BNTR() at /1/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bntr.c:138
[2]PETSC ERROR: #8 TaoSolve() at /1/home/marboeua/Developpement/petsc/src/tao/interface/taosolver.c:177
[3]PETSC ERROR: #9 /home/marboeua/Developpement/mef90/vDef/vDefTAO.F90:370
application called MPI_Abort(MPI_COMM_SELF, 60) - process 0
[2]PETSC ERROR: #9 /home/marboeua/Developpement/mef90/vDef/vDefTAO.F90:370
application called MPI_Abort(MPI_COMM_SELF, 60) - process 0
slurmstepd: error: *** STEP 5034.0 ON bb01 CANCELLED AT 2023-01-12T17:21:07 ***
srun: Job step aborted: Waiting up to 32 seconds for job step to finish.
srun: error: bb01: tasks 0-1: Killed
srun: error: bb01: tasks 2-3: Exited with exit code 1

The error is raised in the middle of the computation after many successful calls of TAOSolve and TAO iterations. My guess is that TAO computes the preconditioner during its first iteration with all variables in the active set. But the preconditioner is never updated when some variables are moved to the inactive set during the next TAO iterations. Am I right? Can you help me with that?

Thanks a lot for your help and your time.
Regards,
Alexis


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

https://www.cse.buffalo.edu/~knepley/<http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230114/5692a32c/attachment-0001.html>


More information about the petsc-users mailing list