[petsc-users] non-manifold DMPLEX

Thu Dec 16 02:09:53 CST 2021

Dear Yuan,

I had read with great interest your posts on the mailing-list. I am also a structural mechanics practitioner and I share your need to manage non-maniforld meshes.
I will let you know if we manage to do stuff with this kind of meshes in DMPLEX.

Regards,
Nicolas

________________________________
De : yuanxi at advancesoft.jp <yuanxi at advancesoft.jp>
Envoyé : jeudi 16 décembre 2021 06:47
À : TARDIEU Nicolas <nicolas.tardieu at edf.fr>
Cc : knepley at gmail.com <knepley at gmail.com>; Patrick Sanan <patrick.sanan at gmail.com>; petsc-users at mcs.anl.gov <petsc-users at mcs.anl.gov>
Objet : Re: [petsc-users] non-manifold DMPLEX

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<https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.mail-archive.com%2Fpetsc-users%40mcs.anl.gov%2Fmsg42462.html&data=04%7C01%7Cnicolas.tardieu%40edf.fr%7C06d686d98ff44f43d12508d9c0578bae%7Ce242425b70fc44dc9ddfc21e304e6c80%7C1%7C0%7C637752305255610253%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=qAaDA96DUChiUCYf0VEO96AIHj8tMWv%2BGNe0v2oQrkc%3D&reserved=0>

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<mailto: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<mailto:knepley at gmail.com> <knepley at gmail.com<mailto:knepley at gmail.com>>
Envoyé : dimanche 12 décembre 2021 12:17
À : Patrick Sanan <patrick.sanan at gmail.com<mailto:patrick.sanan at gmail.com>>
Cc : TARDIEU Nicolas <nicolas.tardieu at edf.fr<mailto:nicolas.tardieu at edf.fr>>; petsc-users at mcs.anl.gov<mailto:petsc-users at mcs.anl.gov> <petsc-users at mcs.anl.gov<mailto: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<mailto: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%7C06d686d98ff44f43d12508d9c0578bae%7Ce242425b70fc44dc9ddfc21e304e6c80%7C1%7C0%7C637752305255610253%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=lgHWjIjpkXyRFrz%2BD22nhSCc%2FVLnDbDAi%2BXohB0mTnM%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<mailto: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%7C06d686d98ff44f43d12508d9c0578bae%7Ce242425b70fc44dc9ddfc21e304e6c80%7C1%7C0%7C637752305255610253%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=yg4p5f6yo0PsCT5VaGaX%2B54WIeJvwH%2F1J65tXoDxyKQ%3D&reserved=0>