Hello! Can somebody explain to me what I am doing wrong here? I took tutorial 28 from KSP folder and modified it a little in order to create manual DM interpolations and set them to PC. When i run the application I receive an error:
[0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: Must call DMShellSetGlobalVector() or DMShellSetCreateGlobalVector() [0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting. [0]PETSC ERROR: Petsc Development GIT revision: v3.17.3-770-gceb9926454c GIT Date: 2022-07-15 21:35:25 +0000 [0]PETSC ERROR: ./a.out on a named MacBook-Pro-Oleg.local by olegshatrov Mon Sep 12 17:53:13 2022 [0]PETSC ERROR: Configure options --prefix=/usr/local CFLAGS="-m64 -O0 -g" CXXFLAGS="-m64 -O0 -g" FCFLAGS="-m64 -O0 -g -ffree-line-length-none" --with-mem-debug=1 --with-debugging=1 --with-mpi-dir=/usr/local/ --with-moab=true [0]PETSC ERROR: #1 DMCreateGlobalVector_Shell() at /Users/olegshatrov/Documents/GitHub/petsc/src/dm/impls/shell/dmshell.c:206 [0]PETSC ERROR: #2 DMCreateGlobalVector() at /Users/olegshatrov/Documents/GitHub/petsc/src/dm/interface/dm.c:998 [0]PETSC ERROR: #3 DMGetGlobalVector() at /Users/olegshatrov/Documents/GitHub/petsc/src/dm/interface/dmget.c:161 [0]PETSC ERROR: #4 KSPCreateVecs() at /Users/olegshatrov/Documents/GitHub/petsc/src/ksp/ksp/interface/iterativ.c:1616 [0]PETSC ERROR: #5 PCSetUp_MG() at /Users/olegshatrov/Documents/GitHub/petsc/src/ksp/pc/impls/mg/mg.c:1129 [0]PETSC ERROR: #6 PCSetUp() at /Users/olegshatrov/Documents/GitHub/petsc/src/ksp/pc/interface/precon.c:991 [0]PETSC ERROR: #7 KSPSetUp() at /Users/olegshatrov/Documents/GitHub/petsc/src/ksp/ksp/interface/itfunc.c:401 [0]PETSC ERROR: #8 main() at ex28.c:67 [0]PETSC ERROR: No PETSc Option Table entries [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-ma...@mcs.anl.gov---------- -------------------------------------------------------------------------- And this is strange, since I did not provide any Shell DM. When i comment out lines with setting interpolation matrix to PCMG everything works fine: [image: image.png] I am testing on MacOS 12.3. MPI: Open MPI: 4.1.4 gcc version 12.1.0 Program is compiled with mpicc ex28.c -I /usr/local/include/ -L /usr/local/lib -lpetsc Thanks in advance, Oleg
ex28.c
Description: Binary data