[petsc-users] strange error by changing the number of cells.

[Yann Jobic]

Dear petsc users,

I'm using petsc FEM framework, with p4est. I'm getting a strange error, 
when changing the number of point in DMPlexCreateBoxMesh() (2D, 20x20 is 
ok, 22x22 is not) :

[0]PETSC ERROR: Nonconforming object sizes
[0]PETSC ERROR: The section point closure size 20 != dual space dimension 18
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html 
for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.10.2, Oct, 09, 2018
[0]PETSC ERROR: ./gauss on a  named crabe by yann Mon Dec 17 22:17:38 2018
[0]PETSC ERROR: Configure options 
--prefix=/local/lib/petsc/3.10/p2/gcc/openmpi_gcc_all 
--with-single-library=0 --with-debugging=1 --download-scalapack=1 
--download-metis=1 --download-parmetis=1 --download-ptscotch=1 
--download-mumps=1 --download-hypre=1 --download-superlu=1 
--download-superlu_dist=1 --download-fblaslapack=1 --download-metis=1 
--download-ml=1 --download-p4est=1 --download-netcdf=1 
--download-pragmatic=1 --with-cxx-dialect=C++11 --download-eigen=1 
--download-parms=1 --download-triangle=1 --download-hdf5=1 --with-zlib=1 
--download-szlib=1 --download-suitesparse=1 --with-shared-libraries=0 
PETSC_ARCH=openmpi_gcc_all
[0]PETSC ERROR: #1 DMProjectLocal_Generic_Plex() line 400 in 
/local/yann/petsc-3.10.2/src/dm/impls/plex/plexproject.c
[0]PETSC ERROR: #2 DMProjectFunctionLocal_Plex() line 530 in 
/local/yann/petsc-3.10.2/src/dm/impls/plex/plexproject.c
[0]PETSC ERROR: #3 DMProjectFunctionLocal() line 6320 in 
/local/yann/petsc-3.10.2/src/dm/interface/dm.c
[0]PETSC ERROR: #4 DMProjectFunctionLocal_p4est() line 4393 in 
/local/yann/petsc-3.10.2/src/dm/impls/forest/p4est/pforest.c
[0]PETSC ERROR: #5 DMProjectFunctionLocal() line 6320 in 
/local/yann/petsc-3.10.2/src/dm/interface/dm.c
[0]PETSC ERROR: #6 SetupVelocity() line 332 in 
/local/yann/fe-utils/petsc/3.10/gauss.c
[0]PETSC ERROR: #7 SetupVelBC() line 359 in 
/local/yann/fe-utils/petsc/3.10/gauss.c
[0]PETSC ERROR: #8 main() line 652 in 
/local/yann/fe-utils/petsc/3.10/gauss.c

I'm using petscspace_degree == 2 for the primary field, and also for the 
auxiliary one (velocity). The error occurs when i project the velocity 
function into the auxiliary field.

Do you see where the problem could come from ?

Best regards,

Yann