[petsc-users] non-manifold DMPLEX

Wed Dec 15 23:47:09 CST 2021

Dear Nicolas :

Sorry for cutting in!

I'd like to indicate that I have also encountered the same problem. Pls
refer to the follows

https://www.mail-archive.com/petsc-users@mcs.anl.gov/msg42462.html

I have read those kinds of mesh successfully  but my program cracked when
doing DMPlexInterpolate. It would be great that you could solve this
problem.

Yuan

2021年12月13日(月) 6:36 TARDIEU Nicolas via petsc-users <petsc-users at mcs.anl.gov
>:

> Dear Patrick and Matthew,
>
> Thank you very much for your answers. I am gonna try to set up such a test
> by assigning cell types.
> Shall I open a MR in order to contribute this example ?
>
> Regards,
> Nicolas
>
> ------------------------------
> *De :* knepley at gmail.com <knepley at gmail.com>
> *Envoyé :* dimanche 12 décembre 2021 12:17
> *À :* Patrick Sanan <patrick.sanan at gmail.com>
> *Cc :* TARDIEU Nicolas <nicolas.tardieu at edf.fr>; petsc-users at mcs.anl.gov <
> petsc-users at mcs.anl.gov>
> *Objet :* Re: [petsc-users] non-manifold DMPLEX
>
> On Sun, Dec 12, 2021 at 6:11 AM Patrick Sanan <patrick.sanan at gmail.com>
> wrote:
>
> 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
> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fpetsc.org%2Frelease%2Fdocs%2Fmanual%2Fdmplex%2F%23representing-unstructured-grids&data=04%7C01%7Cnicolas.tardieu%40edf.fr%7Cb7faa53e924149df02bb08d9bd610c10%7Ce242425b70fc44dc9ddfc21e304e6c80%7C1%7C0%7C637749046784571371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=NX0gPHyjCX3kU%2BZIAZwZQ4951FpJCdri36OzDfRoLbk%3D&reserved=0>
>
> 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.
>
>   Thanks,
>
>     Matt
>
>
> 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/
> <https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.cse.buffalo.edu%2F~knepley%2F&data=04%7C01%7Cnicolas.tardieu%40edf.fr%7Cb7faa53e924149df02bb08d9bd610c10%7Ce242425b70fc44dc9ddfc21e304e6c80%7C1%7C0%7C637749046784571371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=4MLthiWkfhugmu%2BQJ2%2BHxTUR8JbsgEZvlK1b7fUehvo%3D&reserved=0>
>
>
> 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.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20211216/10fabd0a/attachment.html>