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

[Yann Jobic]

I'd like to add that it only occurs in parallel.

Yann

On 17/12/2018 22:24, Yann Jobic via petsc-users wrote:
> 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
>