<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Mon, Dec 17, 2018 at 4:24 PM Yann Jobic via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear petsc users,<br>
<br>
I'm using petsc FEM framework, with p4est. I'm getting a strange error, <br>
when changing the number of point in DMPlexCreateBoxMesh() (2D, 20x20 is <br>
ok, 22x22 is not) :<br>
<br>
[0]PETSC ERROR: Nonconforming object sizes<br>
[0]PETSC ERROR: The section point closure size 20 != dual space dimension 18<br></blockquote><div><br></div><div>Okay, here is what the nomenclature means. The closure size is the size I get from the Section</div><div>when I take the closure of the given point. So for example, if we have quads and Q2 elements,</div><div>that would have a closure size of</div><div><br></div><div>  3 x 3 = 9</div><div><br></div><div>per component. So a 2-component Q2 element would have closure size 18. The dual space</div><div>dimension is the number of dual space basis vectors. If we use a Lagrange element, there</div><div>would be 1 dof per vertex, 1 per edge and 1 per cell for quadratic fields, so on a quad that is</div><div>9 dofs, and for 2 component fields that is 18. This is why the two numbers are designed to match.</div><div>Could you give</div><div><br></div><div>  -petscds_view</div><div><br></div><div>so I can see the element you are using, and also PetscSectionView() on the DM section (or use</div><div>-petscsection_view although that catches all sections, so it hard to find the right one).</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> <br>
for trouble shooting.<br>
[0]PETSC ERROR: Petsc Release Version 3.10.2, Oct, 09, 2018<br>
[0]PETSC ERROR: ./gauss on a  named crabe by yann Mon Dec 17 22:17:38 2018<br>
[0]PETSC ERROR: Configure options <br>
--prefix=/local/lib/petsc/3.10/p2/gcc/openmpi_gcc_all <br>
--with-single-library=0 --with-debugging=1 --download-scalapack=1 <br>
--download-metis=1 --download-parmetis=1 --download-ptscotch=1 <br>
--download-mumps=1 --download-hypre=1 --download-superlu=1 <br>
--download-superlu_dist=1 --download-fblaslapack=1 --download-metis=1 <br>
--download-ml=1 --download-p4est=1 --download-netcdf=1 <br>
--download-pragmatic=1 --with-cxx-dialect=C++11 --download-eigen=1 <br>
--download-parms=1 --download-triangle=1 --download-hdf5=1 --with-zlib=1 <br>
--download-szlib=1 --download-suitesparse=1 --with-shared-libraries=0 <br>
PETSC_ARCH=openmpi_gcc_all<br>
[0]PETSC ERROR: #1 DMProjectLocal_Generic_Plex() line 400 in <br>
/local/yann/petsc-3.10.2/src/dm/impls/plex/plexproject.c<br>
[0]PETSC ERROR: #2 DMProjectFunctionLocal_Plex() line 530 in <br>
/local/yann/petsc-3.10.2/src/dm/impls/plex/plexproject.c<br>
[0]PETSC ERROR: #3 DMProjectFunctionLocal() line 6320 in <br>
/local/yann/petsc-3.10.2/src/dm/interface/dm.c<br>
[0]PETSC ERROR: #4 DMProjectFunctionLocal_p4est() line 4393 in <br>
/local/yann/petsc-3.10.2/src/dm/impls/forest/p4est/pforest.c<br>
[0]PETSC ERROR: #5 DMProjectFunctionLocal() line 6320 in <br>
/local/yann/petsc-3.10.2/src/dm/interface/dm.c<br>
[0]PETSC ERROR: #6 SetupVelocity() line 332 in <br>
/local/yann/fe-utils/petsc/3.10/gauss.c<br>
[0]PETSC ERROR: #7 SetupVelBC() line 359 in <br>
/local/yann/fe-utils/petsc/3.10/gauss.c<br>
[0]PETSC ERROR: #8 main() line 652 in <br>
/local/yann/fe-utils/petsc/3.10/gauss.c<br>
<br>
I'm using petscspace_degree == 2 for the primary field, and also for the <br>
auxiliary one (velocity). The error occurs when i project the velocity <br>
function into the auxiliary field.<br>
<br>
Do you see where the problem could come from ?<br>
<br>
Best regards,<br>
<br>
Yann<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>