[petsc-users] non-manifold DMPLEX

Matthew Knepley knepley at gmail.com
Sun Dec 12 05:17:41 CST 2021

On Sun, Dec 12, 2021 at 6:11 AM Patrick Sanan <patrick.sanan at gmail.com>

> Here you have the following "points":
> - 1 3-cell (the cube volume)
> - 7 2-cells (the 6 faces of the cube plus the extra one)
> - 16 1-cells  (the 12 edges of the cube, plus 3 extra ones from the extra
> face, plus the extra edge)
> - 11 0-cells (the 8 vertices of the cube, pus 2 extra ones from the extra
> face, plus the extra vertex)
> You could encode your mesh as here, by directly specifying relationships
> between these points in the Hasse diagram:
> https://petsc.org/release/docs/manual/dmplex/#representing-unstructured-grids
> Then, maybe the special relation is captured because you've defined the
> "cone" or "support" for each "point", which tells you about the local
> topology everywhere. E.g. to take the simpler case, three of the faces have
> the yellow edge in their "cone", or equivalently the yellow edge has those
> three faces in its "support".

This is correct. I can help you make this if you want. I think if you
assign cell types, you can even get Plex to automatically interpolate.

Note that with this kind of mesh, algorithms which assume a uniform cell
dimension will break, but I am guessing you would not
be interested in those anyway.



> Am Fr., 10. Dez. 2021 um 17:04 Uhr schrieb TARDIEU Nicolas via petsc-users
> <petsc-users at mcs.anl.gov>:
>> Dear PETSc Team,
>> Following a previous discussion on the mailing list, I'd like to
>> experiment with DMPLEX with a very simple non-manifold mesh as shown in the
>> attached picture : a cube connected to a square by an edge and to an edge
>> by a point.
>> I have read some of the papers that Matthew et al. have written, but I
>> must admit that I do not see how to start...
>> I see how the define the different elements but I do not see how to
>> specify the special relationship between the cube and the square and
>> between the cube and the edge.
>> Once it will have been set correctly, what I am hoping is to be able to
>> use all the nice features of the DMPLEX object.
>> Best regards,
>> Nicolas
>> Ce message et toutes les pièces jointes (ci-après le 'Message') sont
>> établis à l'intention exclusive des destinataires et les informations qui y
>> figurent sont strictement confidentielles. Toute utilisation de ce Message
>> non conforme à sa destination, toute diffusion ou toute publication totale
>> ou partielle, est interdite sauf autorisation expresse.
>> Si vous n'êtes pas le destinataire de ce Message, il vous est interdit de
>> le copier, de le faire suivre, de le divulguer ou d'en utiliser tout ou
>> partie. Si vous avez reçu ce Message par erreur, merci de le supprimer de
>> votre système, ainsi que toutes ses copies, et de n'en garder aucune trace
>> sur quelque support que ce soit. Nous vous remercions également d'en
>> avertir immédiatement l'expéditeur par retour du message.
>> Il est impossible de garantir que les communications par messagerie
>> électronique arrivent en temps utile, sont sécurisées ou dénuées de toute
>> erreur ou virus.
>> ____________________________________________________
>> This message and any attachments (the 'Message') are intended solely for
>> the addressees. The information contained in this Message is confidential.
>> Any use of information contained in this Message not in accord with its
>> purpose, any dissemination or disclosure, either whole or partial, is
>> prohibited except formal approval.
>> If you are not the addressee, you may not copy, forward, disclose or use
>> any part of it. If you have received this message in error, please delete
>> it and all copies from your system and notify the sender immediately by
>> return message.
>> E-mail communication cannot be guaranteed to be timely secure, error or
>> virus-free.

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-users/attachments/20211212/95cbc3f8/attachment.html>

More information about the petsc-users mailing list