<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<br>
<div class="moz-forward-container">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<div class="moz-cite-prefix">Le 18/12/2018 à 14:26, Matthew
Knepley a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:CAMYG4GkvV4Mx7UTEg_PDPwzGMfRjf=S9qCpoBK78NmNMM3tUYw@mail.gmail.com">
<div dir="ltr">
<div class="gmail_quote">
<div dir="ltr">On Tue, Dec 18, 2018 at 5:23 AM Yann JOBIC
via petsc-users <<a
href="mailto:petsc-users@mcs.anl.gov"
moz-do-not-send="true">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">
<div bgcolor="#FFFFFF">
<div class="gmail-m_-2741456036267599212moz-cite-prefix">Le
18/12/2018 à 03:01, Matthew Knepley a écrit :<br>
</div>
<blockquote type="cite">
<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"
target="_blank" moz-do-not-send="true">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>
</div>
</blockquote>
The problem arise in parallel, and for a 6x6 grid size.
I'm doing something wrong, but i don't know where. It's
working fine in older petsc version.<br>
You will find in attachment the log files.<br>
<br>
The code for the section is the following :<br>
ierr =
DMGetSection(dmAux,§ion);CHKERRQ(ierr);<br>
PetscSectionView(section,
PETSC_VIEWER_STDOUT_WORLD);<br>
<br>
Thanks for the help !<br>
</div>
</blockquote>
<div><br>
</div>
<div>On 1 proc, you have 120 cells, but on 2 procs you have
126? Can you also send the output with -dm_view?</div>
<div>I can't see what is going on with the parallelization.
Are you specfiying some overlap in the partition?</div>
</div>
</div>
</blockquote>
I don't have an overlap. My DM is created by :<br>
<br>
ierr = DMPlexCreateBoxMesh(comm, dim, user->simplex,
user->cells, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE,
DM_BOUNDARY_NONE, interpolate, dm);CHKERRQ(ierr);<br>
<br>
I then distribute the DM :<br>
<br>
ierr = DMPlexDistribute(*dm, 0, NULL,
&distributedMesh);CHKERRQ(ierr);<br>
<br>
And finally i convert to DMP4EST the resulting DM.<br>
I copy this mesh to create the auxiliary one, for the auxiliary
velocity field.<br>
<br>
Is that right ?<br>
<br>
Thanks,<br>
<br>
Yann<br>
<br>
<blockquote type="cite"
cite="mid:CAMYG4GkvV4Mx7UTEg_PDPwzGMfRjf=S9qCpoBK78NmNMM3tUYw@mail.gmail.com">
<div dir="ltr">
<div class="gmail_quote">
<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">
<div bgcolor="#FFFFFF"> Yann<br>
<br>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_quote">
<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"
moz-do-not-send="true">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>
</div>
</blockquote>
<br>
</div>
</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/%7Eknepley/"
target="_blank" moz-do-not-send="true">https://www.cse.buffalo.edu/~knepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<p><br>
</p>
<pre class="moz-signature" cols="72">--
___________________________
Yann JOBIC
HPC engineer
Polytech Marseille DME
IUSTI-CNRS UMR 6595
Technopôle de Château Gombert
5 rue Enrico Fermi
13453 Marseille cedex 13
Tel : (33) 4 91 10 69 39
ou (33) 4 91 10 69 43
Fax : (33) 4 91 10 69 69
</pre>
</div>
</body>
</html>