[petsc-users] Problems with PCMGSetInterpolation

[Oleg Shatrov]

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-maint at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220912/f4986cce/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 88252 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220912/f4986cce/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex28.c
Type: application/octet-stream
Size: 4630 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220912/f4986cce/attachment-0001.obj>