Hi,
I just pulled and built the HEAD (0cceba267df9e8f10a9c739a928c0b109039c888)
this morning. The SNES variational inequality solvers are break on the headnode
when assembling the bounds. For example misc/ex7 now breaks with:
(base) node004:/shared/libmesh/Linux-gnu-O-new/examples/miscellaneous/ex7$
mpirun -machinefile $PBS_NODEFILE -np 2 ./example-opt -snes_type vinewtonrsls
[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR: Object is in wrong state
[0]PETSC ERROR: Invalid vector
[0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.13.0, Mar 29, 2020
[0]PETSC ERROR: ./example-opt on a Linux-gnu-O named node004 by mesgarnejad Tue
Apr 21 16:48:29 2020
[0]PETSC ERROR: Configure options --download-fblaslapack=1 --download-hdf5=1
--download-hypre=1 --download-metis=1 --download-ml=1 --download-mumps=1
--download-parmetis=1 --download-scalapack=1 --download-sprng=1
--download-suitesparse=yes --download-superlu=yes
--download-superlu_dist=/shared/petsc-3.13.0/Linux-gnu-O/externalpackages/v6.3.0.tar.gz
--download-triangle=1 --download-yaml=1 --with-boost-dir=/shared/opt/boost
--with-cmake=yes --with-debugging=0
--with-mpi-dir=/shared/opt/openmpi/gcc-5.4.0 --with-pic
--with-shared-libraries=1 --with-valgrind=0 --with-x11=1 PETSC_ARCH=Linux-gnu-O
[0]PETSC ERROR: #1 DMRestoreGlobalVector() line 286 in
/shared/petsc-3.13.0/src/dm/interface/dmget.c
[0]PETSC ERROR: #2 SNESSetWorkVecs() line 758 in
/shared/petsc-3.13.0/src/snes/interface/snesut.c
[0]PETSC ERROR: #3 SNESSetUp_VI() line 371 in
/shared/petsc-3.13.0/src/snes/impls/vi/vi.c
[0]PETSC ERROR: #4 SNESSetUp_VINEWTONRSLS() line 711 in
/shared/petsc-3.13.0/src/snes/impls/vi/rs/virs.c
[0]PETSC ERROR: #5 SNESSetUp() line 3148 in
/shared/petsc-3.13.0/src/snes/interface/snes.c
And here is the trace out file:
[New LWP 17273]
[New LWP 17271]
[Thread debugging using libthread_db enabled]
0x00000031c34ac90d in waitpid () from /lib64/libc.so.6
To enable execution of this file add
add-auto-load-safe-path /shared/opt/gcc/lib64/libstdc++.so.6.0.21-gdb.py
line to your configuration file "/home/mesgarnejad/.gdbinit".
To completely disable this security protection add
set auto-load safe-path /
line to your configuration file "/home/mesgarnejad/.gdbinit".
For more information about this security protection see the
"Auto-loading safe path" section in the GDB manual. E.g., run from the shell:
info "(gdb)Auto-loading safe path"
#0 0x00000031c34ac90d in waitpid () from /lib64/libc.so.6
#1 0x00000031c343e909 in do_system () from /lib64/libc.so.6
#2 0x00000031c343ec40 in system () from /lib64/libc.so.6
#3 0x00007fc748b56089 in libMesh::print_trace(std::basic_ostream<char,
std::char_traits<char> >&) () from
/shared/libmesh/Linux-gnu-O-new/lib/libmesh_opt.so.0
#4 0x00007fc748b56b2c in libMesh::write_traceout() () from
/shared/libmesh/Linux-gnu-O-new/lib/libmesh_opt.so.0
#5 0x00007fc748b35b1a in libMesh::libmesh_terminate_handler() () from
/shared/libmesh/Linux-gnu-O-new/lib/libmesh_opt.so.0
#6 0x00007fc747d03ba6 in __cxxabiv1::__terminate(void (*)()) () at
../../.././libstdc++-v3/libsupc++/eh_terminate.cc:47
#7 0x00007fc747d03bf1 in std::terminate() () at
../../.././libstdc++-v3/libsupc++/eh_terminate.cc:57
#8 0x00007fc747d03e08 in __cxa_throw () at
../../.././libstdc++-v3/libsupc++/eh_throw.cc:87
#9 0x00007fc7493637d5 in
libMesh::PetscNonlinearSolver<double>::solve(libMesh::SparseMatrix<double>&,
libMesh::NumericVector<double>&, libMesh::NumericVector<double>&, double,
unsigned int) () from /shared/libmesh/Linux-gnu-O-new/lib/libmesh_opt.so.0
#10 0x00007fc7493d352f in libMesh::NonlinearImplicitSystem::solve() () from
/shared/libmesh/Linux-gnu-O-new/lib/libmesh_opt.so.0
#11 0x0000000000411073 in Biharmonic::run() ()
#12 0x000000000041e15f in main ()
My config.log for libMesh is here:
https://www.dropbox.com/s/uhl8isnxbhmr59s/config.log?dl=0
I’ll keep digging to see if I can find the issue.
Thanks,
Ata
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
