[petsc-users] Ignore points outside the mesh in DMInterpolationSetUP

[Swarnava Ghosh]

Hi Petsc team,

I am trying to setup an interpolator. Is there any way to force
DMInterpolationSetUp  to ignore points which are not inside the mesh. My
code fails with the following error:

[1]PETSC ERROR: --------------------- Error Message
------------------------------------\
--------------------------
[1]PETSC ERROR: Petsc has generated inconsistent data
[1]PETSC ERROR: Point 0: -5.4387 -2.15493 18.465 not located in mesh
[1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shoo\
ting.
[1]PETSC ERROR: Petsc Release Version 3.8.4, Mar, 24, 2018
[1]PETSC ERROR: ../lib/cgdft on a arch-linux2-c-opt named hpc-22-08 by
swarnava Wed Aug \
28 21:52:24 2019
[1]PETSC ERROR: Configure options --prefix=/software/PETSc/3.8.4-intel
--with-cc=mpicc -\
-with-cxx=mpicxx --with-fc=mpif90
--with-blaslapack-dir=/software/Intel/2018.1/compilers\
_and_libraries_2018.1.163/linux/mkl/lib/intel64_lin --with-debugging=no
--with-shared-li\
braries=0 --download-metis --download-parmetis --download-superlu_dist
[1]PETSC ERROR: #1 DMInterpolationSetUp() line 177 in
/groups/hpc-support/install/PETSc/\
petsc-3.8.4_intel_no-debug/src/snes/utils/dmplexsnes.c
[1]PETSC ERROR: #2 SetupInterpolator() line 250 in ./src/cgdft_finemesh.cc

Sincerely,
Swarnava
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190828/42f0023d/attachment.html>