August 2021 - petsc-dev - mailman3.cels.anl.gov

Re: [petsc-dev] building on Spock
by Paul Mullowney 13 Aug '21

13 Aug '21

Here's a script I use for building Hypre that worked a few days ago. It may help with PETSc. https://github.com/hypre-space/hypre/issues/401 -Paul On Fri, Aug 13, 2021 at 12:02 PM Matthew Knepley <knepley(a)gmail.com> wrote: > On Fri, Aug 13, 2021 at 1:49 PM Mark Adams <mfadams(a)lbl.gov> wrote: > >> I was building on Spock a few weeks ago, but am getting these errors now. >> I have a setup with this environment and get this error. >> >> Any ideas? >> > > Here is the error: > > Executing: cc -o /tmp/petsc-6izfgpkr/config.setCompilers/conftest > -L/opt/rocm-4.2.0/lib -lhsa-runtime64 -L${ROCM_PATH}/lib -lamdhip64 > -lhsa-runtime64 /tmp/petsc-6izfgpkr/config.setCompilers/conftest.o > Possible ERROR while running linker: exit code 1 > stderr: > ld.lld: error: unable to find library -lmpi_gtl_hsa > clang-12: error: linker command failed with exit code 1 (use -v to see > invocation) > Linker output before filtering: > > ld.lld: error: unable to find library -lmpi_gtl_hsa > clang-12: error: linker command failed with exit code 1 (use -v to see > invocation) > : > Linker output after filtering: > > ld.lld: error: unable to find library -lmpi_gtl_hsa > clang-12: error: linker command failed with exit code 1 (use -v to see > invocation): > Error testing C compiler: Cannot compile/link C with cc. > Deleting "CC" > > That -l look like it is put in directly by 'cc'. How does anything work on > that machine? > > Matt > > >> Thanks, >> Mark >> >> module load craype-accel-amd-gfx908 >> module load rocm >> module load emacs >> module load zlib >> module load autoconf automake libtool >> >> ## These must be set before compiling so the executable picks up GTL >> export >> PE_MPICH_GTL_DIR_amd_gfx908="-L/opt/cray/pe/mpich/8.1.4/gtl/lib" >> export PE_MPICH_GTL_LIBS_amd_gfx908="-lmpi_gtl_hsa" >> >> ## These must be set before running >> export MPIR_CVAR_GPU_EAGER_DEVICE_MEM=0 >> export MPICH_GPU_SUPPORT_ENABLED=1 >> export MPICH_SMP_SINGLE_COPY_MODE=CMA >> > > > -- > 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/> >

1 0

Re: [petsc-dev] building on Spock
by Matthew Knepley 13 Aug '21

13 Aug '21

On Fri, Aug 13, 2021 at 1:49 PM Mark Adams <mfadams(a)lbl.gov> wrote: > I was building on Spock a few weeks ago, but am getting these errors now. > I have a setup with this environment and get this error. > > Any ideas? > Here is the error: Executing: cc -o /tmp/petsc-6izfgpkr/config.setCompilers/conftest -L/opt/rocm-4.2.0/lib -lhsa-runtime64 -L${ROCM_PATH}/lib -lamdhip64 -lhsa-runtime64 /tmp/petsc-6izfgpkr/config.setCompilers/conftest.o Possible ERROR while running linker: exit code 1 stderr: ld.lld: error: unable to find library -lmpi_gtl_hsa clang-12: error: linker command failed with exit code 1 (use -v to see invocation) Linker output before filtering: ld.lld: error: unable to find library -lmpi_gtl_hsa clang-12: error: linker command failed with exit code 1 (use -v to see invocation) : Linker output after filtering: ld.lld: error: unable to find library -lmpi_gtl_hsa clang-12: error: linker command failed with exit code 1 (use -v to see invocation): Error testing C compiler: Cannot compile/link C with cc. Deleting "CC" That -l look like it is put in directly by 'cc'. How does anything work on that machine? Matt > Thanks, > Mark > > module load craype-accel-amd-gfx908 > module load rocm > module load emacs > module load zlib > module load autoconf automake libtool > > ## These must be set before compiling so the executable picks up GTL > export PE_MPICH_GTL_DIR_amd_gfx908="-L/opt/cray/pe/mpich/8.1.4/gtl/lib" > export PE_MPICH_GTL_LIBS_amd_gfx908="-lmpi_gtl_hsa" > > ## These must be set before running > export MPIR_CVAR_GPU_EAGER_DEVICE_MEM=0 > export MPICH_GPU_SUPPORT_ENABLED=1 > export MPICH_SMP_SINGLE_COPY_MODE=CMA > -- 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/>

1 0

building on Spock
by Mark Adams 13 Aug '21

13 Aug '21

I was building on Spock a few weeks ago, but am getting these errors now. I have a setup with this environment and get this error. Any ideas? Thanks, Mark module load craype-accel-amd-gfx908 module load rocm module load emacs module load zlib module load autoconf automake libtool ## These must be set before compiling so the executable picks up GTL export PE_MPICH_GTL_DIR_amd_gfx908="-L/opt/cray/pe/mpich/8.1.4/gtl/lib" export PE_MPICH_GTL_LIBS_amd_gfx908="-lmpi_gtl_hsa" ## These must be set before running export MPIR_CVAR_GPU_EAGER_DEVICE_MEM=0 export MPICH_GPU_SUPPORT_ENABLED=1 export MPICH_SMP_SINGLE_COPY_MODE=CMA

1 0

Re: [petsc-dev] PETSC_COMM_WORLD not a PETSc communicator
by Barry Smith 12 Aug '21

12 Aug '21

User visible communicators generally do not have a keyval attached. Rather the keyval is attached to the inner communicator; because we don't want both PETSc and the user doing MPI operations on the same communicator (to prevent mixing up tags). I think PetscShmCommGet() is wrong. I think it should not call MPI_Comm_get_attr(globcomm,Petsc_Counter_keyval,&counter,&flg); but should call PetscCommDuplicate() and use that communicator to stash the pshmcomm; then you would not have the problem you are having. Barry > On Aug 12, 2021, at 11:05 AM, Pierre Jolivet <pierre(a)joliv.et> wrote: > > Hello, > Is there a specific reason why PETSC_COMM_WORLD is not a PETSc communicator, i.e., has no Petsc_Counter_keyval attached? > ierr = PetscOmpCtrlCreate(PETSC_COMM_WORLD,nthreads,&ctrl);CHKERRQ(ierr); > yields > [0]PETSC ERROR: Bad MPI communicator supplied must be a PETSc communicator > [0]PETSC ERROR: #1 PetscShmCommGet() at src/sys/utils/mpishm.c:60 > [0]PETSC ERROR: #2 PetscOmpCtrlCreate() at src/sys/utils/mpishm.c:340 > > Thanks, > Pierre

1 0

PETSC_COMM_WORLD not a PETSc communicator
by Pierre Jolivet 12 Aug '21

12 Aug '21

Hello, Is there a specific reason why PETSC_COMM_WORLD is not a PETSc communicator, i.e., has no Petsc_Counter_keyval attached? ierr = PetscOmpCtrlCreate(PETSC_COMM_WORLD,nthreads,&ctrl);CHKERRQ(ierr); yields [0]PETSC ERROR: Bad MPI communicator supplied must be a PETSc communicator [0]PETSC ERROR: #1 PetscShmCommGet() at src/sys/utils/mpishm.c:60 [0]PETSC ERROR: #2 PetscOmpCtrlCreate() at src/sys/utils/mpishm.c:340 Thanks, Pierre

1 0

Re: [petsc-dev] -ksp_error_if_not_converged
by Barry Smith 08 Aug '21

08 Aug '21

https://gitlab.com/petsc/petsc/-/merge_requests/4219 > On Aug 8, 2021, at 7:42 AM, Pierre Jolivet <pierre(a)joliv.et> wrote: > > Hello, > I’d like to trigger an error if a linear solve does not converge due to DIVERGED_ITS. > Since this MR https://gitlab.com/petsc/petsc/-/merge_requests/1294 <https://gitlab.com/petsc/petsc/-/merge_requests/1294>, no error is triggered anymore if ksp->reason == KSP_DIVERGED_ITS. > This seems to be inconsistent with -snes_error_if_not_converged. > Is there a workaround? Shouldn’t this be clarified in the manual page? > > Thank you, > Pierre > > $ src/ksp/ksp/tutorials/ex1 -ksp_max_it 1 -ksp_converged_reason -ksp_error_if_not_converged && echo "No error" > Linear solve did not converge due to DIVERGED_ITS iterations 1 > KSP Object: 1 MPI processes > type: gmres > restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement > happy breakdown tolerance 1e-30 > maximum iterations=1, initial guess is zero > tolerances: relative=1e-05, absolute=1e-50, divergence=10000. > left preconditioning > using PRECONDITIONED norm type for convergence test > PC Object: 1 MPI processes > type: jacobi > type DIAGONAL > linear system matrix = precond matrix: > Mat Object: 1 MPI processes > type: seqaij > rows=10, cols=10 > total: nonzeros=28, allocated nonzeros=50 > total number of mallocs used during MatSetValues calls=0 > not using I-node routines > Norm of error 2.95296, Iterations 1 > No error

1 0

-ksp_error_if_not_converged
by Pierre Jolivet 08 Aug '21

08 Aug '21

Hello, I’d like to trigger an error if a linear solve does not converge due to DIVERGED_ITS. Since this MR https://gitlab.com/petsc/petsc/-/merge_requests/1294 <https://gitlab.com/petsc/petsc/-/merge_requests/1294>, no error is triggered anymore if ksp->reason == KSP_DIVERGED_ITS. This seems to be inconsistent with -snes_error_if_not_converged. Is there a workaround? Shouldn’t this be clarified in the manual page? Thank you, Pierre $ src/ksp/ksp/tutorials/ex1 -ksp_max_it 1 -ksp_converged_reason -ksp_error_if_not_converged && echo "No error" Linear solve did not converge due to DIVERGED_ITS iterations 1 KSP Object: 1 MPI processes type: gmres restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement happy breakdown tolerance 1e-30 maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using PRECONDITIONED norm type for convergence test PC Object: 1 MPI processes type: jacobi type DIAGONAL linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=10, cols=10 total: nonzeros=28, allocated nonzeros=50 total number of mallocs used during MatSetValues calls=0 not using I-node routines Norm of error 2.95296, Iterations 1 No error

1 0

petsc-3.15.3 now available
by Satish Balay 06 Aug '21

06 Aug '21

Dear PETSc users, The patch release petsc-3.15.3 is now available for download. https://petsc.org/release/download/ Satish

1 0

Re: [petsc-dev] Broken snes_tutorials-ex20_2d_p1_gmg_vcycle_cr and CRSetup_Private()
by Pierre Jolivet 04 Aug '21

04 Aug '21

I’ve flagged the test as broken on top of my feature branch (https://gitlab.com/petsc/petsc/-/commit/b02fd907789ed54b7deef6e9395f48df009… <https://gitlab.com/petsc/petsc/-/commit/b02fd907789ed54b7deef6e9395f48df009…>). Satish will have to make the executive decision of merging this or not, as the pipeline is now successful. Thanks, Pierre > On 4 Aug 2021, at 2:50 PM, Stefano Zampini <stefano.zampini(a)gmail.com> wrote: > > Matt, > > I’m fine flagging the test as broken for now > Actually, this test has always been broken…. > As Pierre said, if we make these matrices inherit from MatShell, we get shift and scale for free > > >> On Aug 4, 2021, at 10:13 AM, Pierre Jolivet <pierre(a)joliv.et <mailto:[email protected]>> wrote: >> >> Hello, >> Up until very recently, MatGetOwnershipRange() was returning only zeros for any MATNORMAL (thanks Stefano for the help fixing that https://gitlab.com/petsc/petsc/-/merge_requests/4182#note_640848528 <https://gitlab.com/petsc/petsc/-/merge_requests/4182#note_640848528>). >> In CRSetup_Private(), there is a call to MatCreateNormal() followed by MatShift(). >> Since ops->shift = NULL for MATNORMAL, code goes into MatShift_Basic(). >> There, a for loop goes from “start” to “end”, but the values were previously only 0. >> Now that MatGetOwnershipRange() has been fixed for MATNORMAL, the bounds are nonzero, and snes_tutorials-ex20_2d_p1_gmg_vcycle_cr fails because MATNORMAL has no ops->setvalues >> What is the proper fix or how can I bypass the issue? Would it be OK to (temporarily) mark this test as broken? >> Previous code was obviously not doing what was intended in theory, but yet I’m guessing numerically everything was OK? >> I’m thinking https://gitlab.com/petsc/petsc/-/issues/972 <https://gitlab.com/petsc/petsc/-/issues/972> should do the trick, but this is somehow orthogonal to my developments, which are now blocked because of failing pipelines in snes_tutorials-ex20_2d_p1_gmg_vcycle_cr (https://gitlab.com/petsc/petsc/-/jobs/1475595231 <https://gitlab.com/petsc/petsc/-/jobs/1475595231>). >> >> Thanks, >> Pierre >> >> $ make -f gmakefile test query='name' queryval="*snes_tutorials-ex20_2d_p1_gmg_vcycle_cr*" >> Using MAKEFLAGS: queryval=*snes_tutorials-ex20_2d_p1_gmg_vcycle_cr* query=name >> TEST arch-darwin-c-debug-real/tests/counts/snes_tutorials-ex20_2d_p1_gmg_vcycle_cr.counts >> not ok snes_tutorials-ex20_2d_p1_gmg_vcycle_cr # Error code: 56 >> # [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- >> # [0]PETSC ERROR: No support for this operation for this object type >> # [0]PETSC ERROR: Mat type normal >> # [0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html <https://www.mcs.anl.gov/petsc/documentation/faq.html> for trouble shooting. >> # [0]PETSC ERROR: Petsc Development GIT revision: v3.15.2-633-gdc9bbe6e6a GIT Date: 2021-08-03 20:53:25 +0200 >> # [0]PETSC ERROR: ../ex20 on a arch-darwin-c-debug-real named aardvark.local by jolivet Wed Aug 4 10:04:59 2021 >> # [0]PETSC ERROR: Configure options --download-c2html --download-elemental --download-hpddm --download-hpddm-commit=origin/omp --download-htool --download-htool-commit=origin/main --download-hypre --download-metis --download-suitesparse --download-mmg --download-mumps --download-hwloc --download-p4est --download-parmetis --download-parmmg --download-ptscotch --download-slepc --download-slepc-commit=origin/main --download-slepc-configure-arguments="--download-arpack=https://github.com/prj-/arpack-ng/archive/b64dccb.tar.gz <https://github.com/prj-/arpack-ng/archive/b64dccb.tar.gz>" --download-tetgen --download-triangle --download-zlib --with-blaslapack-include=/opt/intel/oneapi/mkl/latest/include --with-blaslapack-lib="[/opt/intel/oneapi/compiler/latest/mac/compiler/lib/libiomp5.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_intel_lp64.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_intel_thread.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_core.dylib]" --with-fortran-bindings=1 --with-mkl_cpardiso-include=/opt/intel/oneapi/mkl/latest/include --with-mkl_cpardiso=1 --with-mkl_pardiso-include=/opt/intel/oneapi/mkl/latest/include --with-mkl_pardiso=1 --with-opengl=1 --with-petsc4py --with-scalapack-include=/opt/intel/oneapi/mkl/latest/include --with-scalapack-lib="[/opt/intel/oneapi/mkl/latest/lib/libmkl_scalapack_lp64.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_blacs_mpich_lp64.dylib]" --with-scalar-type=real --with-shared-libraries=1 CFLAGS="-Wall -g" CXXFLAGS="-Wall -g" FFLAGS="-Wall -ffree-line-length-512 -g" PETSC_ARCH=arch-darwin-c-debug-real >> # [0]PETSC ERROR: #1 MatSetValues() at /Users/jolivet/Documents/repositories/petsc/src/mat/interface/matrix.c:1376 >> # [0]PETSC ERROR: #2 MatShift_Basic() at /Users/jolivet/Documents/repositories/petsc/src/mat/utils/gcreate.c:26 >> # [0]PETSC ERROR: #3 MatShift() at /Users/jolivet/Documents/repositories/petsc/src/mat/utils/axpy.c:309 >> # [0]PETSC ERROR: #4 CRSetup_Private() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/mg/mg.c:300 >> # [0]PETSC ERROR: #5 PCSetUp_Shell() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/shell/shellpc.c:101 >> # [0]PETSC ERROR: #6 PCSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:1016 >> # [0]PETSC ERROR: #7 PCApply() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:441 >> # [0]PETSC ERROR: #8 PCApply_Composite_Multiplicative() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/composite/composite.c:50 >> # [0]PETSC ERROR: #9 PCApply() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:445 >> # [0]PETSC ERROR: #10 KSP_PCApply() at /Users/jolivet/repo/petsc/include/petsc/private/kspimpl.h:365 >> # [0]PETSC ERROR: #11 KSPInitialResidual() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itres.c:65 >> # [0]PETSC ERROR: #12 KSPSolve_GMRES() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/impls/gmres/gmres.c:246 >> # [0]PETSC ERROR: #13 KSPSolve_Private() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:912 >> # [0]PETSC ERROR: #14 KSPSolve() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:1084 >> # [0]PETSC ERROR: #15 KSPSetUp_Chebyshev() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/impls/cheby/cheby.c:87 >> # [0]PETSC ERROR: #16 KSPSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:365 >> # [0]PETSC ERROR: #17 PCSetUp_MG() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/mg/mg.c:1247 >> # [0]PETSC ERROR: #18 PCSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:1016 >> # [0]PETSC ERROR: #19 KSPSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:406 >> # [0]PETSC ERROR: #20 KSPSolve_Private() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:850 >> # [0]PETSC ERROR: #21 KSPSolve() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:1084 >> # [0]PETSC ERROR: #22 SNESSolve_NEWTONLS() at /Users/jolivet/Documents/repositories/petsc/src/snes/impls/ls/ls.c:225 >> # [0]PETSC ERROR: #23 SNESSolve() at /Users/jolivet/Documents/repositories/petsc/src/snes/interface/snes.c:4769 >> # [0]PETSC ERROR: #24 main() at /Users/jolivet/Documents/repositories/petsc/src/snes/tutorials/ex20.c:175 >> # [0]PETSC ERROR: PETSc Option Table entries: >> # [0]PETSC ERROR: -dm_plex_box_faces 2,2 >> # [0]PETSC ERROR: -dm_refine_hierarchy 3 >> # [0]PETSC ERROR: -ksp_rtol 5e-10 >> # [0]PETSC ERROR: -mg_levels_cr_ksp_converged_rate >> # [0]PETSC ERROR: -mg_levels_cr_ksp_converged_rate_type error >> # [0]PETSC ERROR: -mg_levels_cr_ksp_max_it 5 >> # [0]PETSC ERROR: -mg_levels_esteig_ksp_max_it 10 >> # [0]PETSC ERROR: -mg_levels_esteig_ksp_type cg >> # [0]PETSC ERROR: -mg_levels_ksp_chebyshev_esteig 0,0.05,0,1.05 >> # [0]PETSC ERROR: -mg_levels_ksp_max_it 5 >> # [0]PETSC ERROR: -mg_levels_ksp_norm_type preconditioned >> # [0]PETSC ERROR: -pc_mg_adapt_cr >> # [0]PETSC ERROR: -pc_type mg >> # [0]PETSC ERROR: -potential_petscspace_degree 1 >> # [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint(a)mcs.anl.gov <mailto:[email protected]>---------- >> # Abort(56) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0 >> ok snes_tutorials-ex20_2d_p1_gmg_vcycle_cr # SKIP Command failed so no diff >> >> >> # FAILED snes_tutorials-ex20_2d_p1_gmg_vcycle_cr >> # >> # To rerun failed tests: >> # /usr/bin/make -f gmakefile test test-fail=1 >

1 0

Re: [petsc-dev] Broken snes_tutorials-ex20_2d_p1_gmg_vcycle_cr and CRSetup_Private()
by Stefano Zampini 04 Aug '21

04 Aug '21

Matt, I’m fine flagging the test as broken for now Actually, this test has always been broken…. As Pierre said, if we make these matrices inherit from MatShell, we get shift and scale for free > On Aug 4, 2021, at 10:13 AM, Pierre Jolivet <pierre(a)joliv.et> wrote: > > Hello, > Up until very recently, MatGetOwnershipRange() was returning only zeros for any MATNORMAL (thanks Stefano for the help fixing that https://gitlab.com/petsc/petsc/-/merge_requests/4182#note_640848528 <https://gitlab.com/petsc/petsc/-/merge_requests/4182#note_640848528>). > In CRSetup_Private(), there is a call to MatCreateNormal() followed by MatShift(). > Since ops->shift = NULL for MATNORMAL, code goes into MatShift_Basic(). > There, a for loop goes from “start” to “end”, but the values were previously only 0. > Now that MatGetOwnershipRange() has been fixed for MATNORMAL, the bounds are nonzero, and snes_tutorials-ex20_2d_p1_gmg_vcycle_cr fails because MATNORMAL has no ops->setvalues > What is the proper fix or how can I bypass the issue? Would it be OK to (temporarily) mark this test as broken? > Previous code was obviously not doing what was intended in theory, but yet I’m guessing numerically everything was OK? > I’m thinking https://gitlab.com/petsc/petsc/-/issues/972 <https://gitlab.com/petsc/petsc/-/issues/972> should do the trick, but this is somehow orthogonal to my developments, which are now blocked because of failing pipelines in snes_tutorials-ex20_2d_p1_gmg_vcycle_cr (https://gitlab.com/petsc/petsc/-/jobs/1475595231 <https://gitlab.com/petsc/petsc/-/jobs/1475595231>). > > Thanks, > Pierre > > $ make -f gmakefile test query='name' queryval="*snes_tutorials-ex20_2d_p1_gmg_vcycle_cr*" > Using MAKEFLAGS: queryval=*snes_tutorials-ex20_2d_p1_gmg_vcycle_cr* query=name > TEST arch-darwin-c-debug-real/tests/counts/snes_tutorials-ex20_2d_p1_gmg_vcycle_cr.counts > not ok snes_tutorials-ex20_2d_p1_gmg_vcycle_cr # Error code: 56 > # [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- > # [0]PETSC ERROR: No support for this operation for this object type > # [0]PETSC ERROR: Mat type normal > # [0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html <https://www.mcs.anl.gov/petsc/documentation/faq.html> for trouble shooting. > # [0]PETSC ERROR: Petsc Development GIT revision: v3.15.2-633-gdc9bbe6e6a GIT Date: 2021-08-03 20:53:25 +0200 > # [0]PETSC ERROR: ../ex20 on a arch-darwin-c-debug-real named aardvark.local by jolivet Wed Aug 4 10:04:59 2021 > # [0]PETSC ERROR: Configure options --download-c2html --download-elemental --download-hpddm --download-hpddm-commit=origin/omp --download-htool --download-htool-commit=origin/main --download-hypre --download-metis --download-suitesparse --download-mmg --download-mumps --download-hwloc --download-p4est --download-parmetis --download-parmmg --download-ptscotch --download-slepc --download-slepc-commit=origin/main --download-slepc-configure-arguments="--download-arpack=https://github.com/prj-/arpack-ng/archive/b64dccb.tar.gz <https://github.com/prj-/arpack-ng/archive/b64dccb.tar.gz>" --download-tetgen --download-triangle --download-zlib --with-blaslapack-include=/opt/intel/oneapi/mkl/latest/include --with-blaslapack-lib="[/opt/intel/oneapi/compiler/latest/mac/compiler/lib/libiomp5.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_intel_lp64.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_intel_thread.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_core.dylib]" --with-fortran-bindings=1 --with-mkl_cpardiso-include=/opt/intel/oneapi/mkl/latest/include --with-mkl_cpardiso=1 --with-mkl_pardiso-include=/opt/intel/oneapi/mkl/latest/include --with-mkl_pardiso=1 --with-opengl=1 --with-petsc4py --with-scalapack-include=/opt/intel/oneapi/mkl/latest/include --with-scalapack-lib="[/opt/intel/oneapi/mkl/latest/lib/libmkl_scalapack_lp64.dylib,/opt/intel/oneapi/mkl/latest/lib/libmkl_blacs_mpich_lp64.dylib]" --with-scalar-type=real --with-shared-libraries=1 CFLAGS="-Wall -g" CXXFLAGS="-Wall -g" FFLAGS="-Wall -ffree-line-length-512 -g" PETSC_ARCH=arch-darwin-c-debug-real > # [0]PETSC ERROR: #1 MatSetValues() at /Users/jolivet/Documents/repositories/petsc/src/mat/interface/matrix.c:1376 > # [0]PETSC ERROR: #2 MatShift_Basic() at /Users/jolivet/Documents/repositories/petsc/src/mat/utils/gcreate.c:26 > # [0]PETSC ERROR: #3 MatShift() at /Users/jolivet/Documents/repositories/petsc/src/mat/utils/axpy.c:309 > # [0]PETSC ERROR: #4 CRSetup_Private() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/mg/mg.c:300 > # [0]PETSC ERROR: #5 PCSetUp_Shell() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/shell/shellpc.c:101 > # [0]PETSC ERROR: #6 PCSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:1016 > # [0]PETSC ERROR: #7 PCApply() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:441 > # [0]PETSC ERROR: #8 PCApply_Composite_Multiplicative() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/composite/composite.c:50 > # [0]PETSC ERROR: #9 PCApply() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:445 > # [0]PETSC ERROR: #10 KSP_PCApply() at /Users/jolivet/repo/petsc/include/petsc/private/kspimpl.h:365 > # [0]PETSC ERROR: #11 KSPInitialResidual() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itres.c:65 > # [0]PETSC ERROR: #12 KSPSolve_GMRES() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/impls/gmres/gmres.c:246 > # [0]PETSC ERROR: #13 KSPSolve_Private() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:912 > # [0]PETSC ERROR: #14 KSPSolve() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:1084 > # [0]PETSC ERROR: #15 KSPSetUp_Chebyshev() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/impls/cheby/cheby.c:87 > # [0]PETSC ERROR: #16 KSPSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:365 > # [0]PETSC ERROR: #17 PCSetUp_MG() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/impls/mg/mg.c:1247 > # [0]PETSC ERROR: #18 PCSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/pc/interface/precon.c:1016 > # [0]PETSC ERROR: #19 KSPSetUp() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:406 > # [0]PETSC ERROR: #20 KSPSolve_Private() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:850 > # [0]PETSC ERROR: #21 KSPSolve() at /Users/jolivet/Documents/repositories/petsc/src/ksp/ksp/interface/itfunc.c:1084 > # [0]PETSC ERROR: #22 SNESSolve_NEWTONLS() at /Users/jolivet/Documents/repositories/petsc/src/snes/impls/ls/ls.c:225 > # [0]PETSC ERROR: #23 SNESSolve() at /Users/jolivet/Documents/repositories/petsc/src/snes/interface/snes.c:4769 > # [0]PETSC ERROR: #24 main() at /Users/jolivet/Documents/repositories/petsc/src/snes/tutorials/ex20.c:175 > # [0]PETSC ERROR: PETSc Option Table entries: > # [0]PETSC ERROR: -dm_plex_box_faces 2,2 > # [0]PETSC ERROR: -dm_refine_hierarchy 3 > # [0]PETSC ERROR: -ksp_rtol 5e-10 > # [0]PETSC ERROR: -mg_levels_cr_ksp_converged_rate > # [0]PETSC ERROR: -mg_levels_cr_ksp_converged_rate_type error > # [0]PETSC ERROR: -mg_levels_cr_ksp_max_it 5 > # [0]PETSC ERROR: -mg_levels_esteig_ksp_max_it 10 > # [0]PETSC ERROR: -mg_levels_esteig_ksp_type cg > # [0]PETSC ERROR: -mg_levels_ksp_chebyshev_esteig 0,0.05,0,1.05 > # [0]PETSC ERROR: -mg_levels_ksp_max_it 5 > # [0]PETSC ERROR: -mg_levels_ksp_norm_type preconditioned > # [0]PETSC ERROR: -pc_mg_adapt_cr > # [0]PETSC ERROR: -pc_type mg > # [0]PETSC ERROR: -potential_petscspace_degree 1 > # [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint(a)mcs.anl.gov <mailto:[email protected]>---------- > # Abort(56) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0 > ok snes_tutorials-ex20_2d_p1_gmg_vcycle_cr # SKIP Command failed so no diff > > > # FAILED snes_tutorials-ex20_2d_p1_gmg_vcycle_cr > # > # To rerun failed tests: > # /usr/bin/make -f gmakefile test test-fail=1

1 0