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/[email protected]/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 <[email protected] >: > 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 :* [email protected] <[email protected]> > *Envoyé :* dimanche 12 décembre 2021 12:17 > *À :* Patrick Sanan <[email protected]> > *Cc :* TARDIEU Nicolas <[email protected]>; [email protected] < > [email protected]> > *Objet :* Re: [petsc-users] non-manifold DMPLEX > > On Sun, Dec 12, 2021 at 6:11 AM Patrick Sanan <[email protected]> > 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 > <[email protected]>: > > 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. >
