<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Feb 20, 2018 at 12:30 PM, Danyang Su <span dir="ltr"><<a href="mailto:danyang.su@gmail.com" target="_blank">danyang.su@gmail.com</a>></span> wrote:<br><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">
<p>Hi All,</p>
<p>I tried to compile the DMPlexCreateSection code but got error
information as shown below.<br>
</p>
<p>Error: Symbol 'petsc_null_is' at (1) has no IMPLICIT type</p>
<p>I tried to use PETSC_NULL_OBJECT instead of PETSC_NULL_IS, then
the code can be compiled but run into Segmentation Violation error
in DMPlexCreateSection.<br></p></div></blockquote><div>From the webpage</div><div><br></div><div> <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexCreateSection.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexCreateSection.html</a> <br></div><div><br></div><div>The F90 version is DMPlexCreateSectionF90. Doing this with F77 arrays would have been too painful.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</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"><p>
</p>
<p>dmda_flow%da is distributed dm object that works fine. </p>
<p>The fortran example I follow is
<a class="gmail-m_8962138310016109093moz-txt-link-freetext" href="http://www.mcs.anl.gov/petsc/petsc-dev/src/dm/impls/plex/examples/tutorials/ex1f90.F90" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>petsc-dev/src/dm/impls/plex/<wbr>examples/tutorials/ex1f90.F90</a>.
<br>
</p>
<p>What parameters should I use if passing null to bcField, bcComps,
bcPoints and perm.<br>
</p>
<pre style="color:rgb(0,0,0);font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Vec/PetscErrorCode.html#PetscErrorCode" target="_blank">PetscErrorCode</a> <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/DMPLEX/DMPlexCreateSection.html#DMPlexCreateSection" target="_blank">DMPlexCreateSection</a>(<a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/DM/DM.html#DM" target="_blank">DM</a> dm, <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Vec/PetscInt.html#PetscInt" target="_blank">PetscInt</a> dim, <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Vec/PetscInt.html#PetscInt" target="_blank">PetscInt</a> numFields,const <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Vec/PetscInt.html#PetscInt" target="_blank">PetscInt</a> numComp[],const <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Vec/PetscInt.html#PetscInt" target="_blank">PetscInt</a> numDof[], <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Vec/PetscInt.html#PetscInt" target="_blank">PetscInt</a> numBC,const <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Vec/PetscInt.html#PetscInt" target="_blank">PetscInt</a> bcField[],
const <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/IS/IS.html#IS" target="_blank">IS</a> bcComps[], const <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/IS/IS.html#IS" target="_blank">IS</a> bcPoints[], <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/IS/IS.html#IS" target="_blank">IS</a> perm, <a href="https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/IS/PetscSection.html#PetscSection" target="_blank">PetscSection</a> *section)
</pre>
<p>#include <petsc/finclude/petscis.h><br>
#include <petsc/finclude/petscis.h90><br>
#include <petsc/finclude/petscdmplex.h></p>
<p>...<br>
</p>
<p>#ifdef USG<br>
numFields = 1<br>
numComp(1) = 1<br>
pNumComp => numComp<br>
<br>
do i = 1, numFields*(dmda_flow%dim+1)<br>
numDof(i) = 0<br>
end do<br>
numDof(0*(dmda_flow%dim+1)+1) = dmda_flow%dof<br>
pNumDof => numDof<br>
<br>
numBC = 0<br>
<br>
call
DMPlexCreateSection(dmda_flow%<wbr>da,dmda_flow%dim, &<br>
<wbr>
numFields,pNumComp,pNumDof, <wbr> &<br>
<wbr>
numBC,PETSC_NULL_INTEGER, <wbr> &<br>
<wbr>
PETSC_NULL_IS,PETSC_NULL_IS, <wbr>
& !Error here<br>
<wbr> PETSC_NULL_IS,section,ierr)<br>
CHKERRQ(ierr)<br>
<br>
call PetscSectionSetFieldName(<wbr>section,0,'flow',ierr)<br>
CHKERRQ(ierr)<br>
<br>
call DMSetDefaultSection(dmda_flow%<wbr>da,section,ierr)<br>
CHKERRQ(ierr)<br>
<br>
call PetscSectionDestroy(section,<wbr>ierr)<br>
CHKERRQ(ierr)<br>
#endif</p>
<p>Thanks,</p>
<p>Danyang<br>
</p>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><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.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>