[petsc-dev] DMPlex + Cohesive Cells Setup

[Jacob Faibussowitsch]

I’m getting an argument out of bounds error (attached) but I suspect it’s because I’m feeding in bad input to either calls (I am marking __all__ vertices in the mesh in my input label to DMPlexLabelCohesiveComplete). Are the cohesive grains allowed to intersect? Does a cohesive region have to be closed? Does it have to be entirely enclosed by the mesh, i.e. not touching any of the boundaries?

Best regards,

Jacob Faibussowitsch
(Jacob Fai - booss - oh - vitch)
Cell: (312) 694-3391


> On Sep 2, 2020, at 09:28, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Mon, Aug 31, 2020 at 5:39 PM Jacob Faibussowitsch <jacob.fai at gmail.com <mailto:jacob.fai at gmail.com>> wrote:
> Hello All,
> 
> What is the simplest possible configuration required to insert cohesive cells into a plex? Plex/tests/ex5.c would be the example to ape here but there is a lot there that I don’t (think) I need like a DS or any kind of physics. For now, I am essentially looking to:
> 
> 1. Read the mesh
> 2. Perform some mesh analysis
> 3. Slap some cohesive cells along one/multiple areas identified in the previous step (maybe even everywhere?)
> 4. Save the new mesh
> 
> This doesn’t even have to parallel. 
> 
> You can try it out easily. If it does not give you what you want, we can fix that up.
> 
> https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexConstructCohesiveCells.html <https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexConstructCohesiveCells.html> will give you back a mesh
> with the cohesive cells inserted. For input, you need a suitable label.
> 
> https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexLabelCohesiveComplete.html#DMPlexLabelCohesiveComplete <https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexLabelCohesiveComplete.html#DMPlexLabelCohesiveComplete> starts
> with a label marking vertices on an internal surface and creates the label needed by the call above.
> 
> Let me know if that works.
> 
>   Thanks,
> 
>      Matt
>  
> Best regards,
> 
> Jacob Faibussowitsch
> (Jacob Fai - booss - oh - vitch)
> Cell: (312) 694-3391
> 
> 
> 
> -- 
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener
> 
> https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20200902/bb31c3a9/attachment.html>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: err.txt
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20200902/bb31c3a9/attachment.txt>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20200902/bb31c3a9/attachment-0001.html>