This would be a particularly useful example for testing the PETSc library, as most applications are set on manifolds. There’s a recently-reformatted chapter in the developers docs on how the testing works - this explains the /* TEST */ blocks you see. If parts of these docs are confusing, that’s very useful feedback for us! https://petsc.org/release/developers/testing
Matthew Knepley <[email protected]> schrieb am Mo. 13. Dez. 2021 um 03:12: > On Sun, Dec 12, 2021 at 4:36 PM TARDIEU Nicolas <[email protected]> > wrote: > >> 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 ? >> > > Yes, that would be great. > > Thanks, > > Matt > > >> 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. >> > > > -- > 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/> >
