[petsc-users] difference in Face Sets in latest petsc release

[Matthew Knepley]

On Tue, Dec 12, 2023 at 7:13 PM Sharan Roongta <s.roongta at mpie.de> wrote:

> Thanks for the clarification.
> However, would it be more consistent to differentiate between 0D (vertex
> sets), 1D (edge sets), 2d (faces) and 3D (cell sets)?
> If I want to now apply boundary condition on a face with tag 1, it would
> contain the 4 edges making up that face and an additional edge with the
> same physical tag?
> Basically, I can’t differentiate between the two entities
>

When we do this in PyLith, we only tag the things we want tagged. In this
mesh, absolutely everything is tagged, and with overlapping tag ranges.
That seems counterproductive.

  Thanks,

     Matt


> Thanks,
> Sharan
>
> *Group - Theory & Simulation*
> *Department of Microstructure Physics & Alloy Design*
>
>
>
> * From: * Matthew Knepley <knepley at gmail.com>
> * To: * Sharan Roongta <s.roongta at mpie.de>
> * Cc: * "petsc-users at mcs.anl.gov" <petsc-users at mcs.anl.gov>
> * Sent: * 13/12/2023 12:51 AM
> * Subject: * Re: [petsc-users] difference in Face Sets in latest petsc
> release
>
> On Tue, Dec 12, 2023 at 5:50 PM Sharan Roongta <s.roongta at mpie.de> wrote:
>
> Hello,
>
> I see discrepancy in the size/value of the 'Face Sets' printed in the
> current release v3.20.2 , and v3.18.6
>
> Attached is the .msh file
>
> -dm_view with v3.18.6
> DM Object: Generated Mesh 1 MPI process
>   type: plex
> Generated Mesh in 3 dimensions:
>   Number of 0-cells per rank: 14
>   Number of 1-cells per rank: 49
>   Number of 2-cells per rank: 60
>   Number of 3-cells per rank: 24
> Labels:
>   celltype: 4 strata with value/size (0 (14), 6 (24), 3 (60), 1 (49))
>   depth: 4 strata with value/size (0 (14), 1 (49), 2 (60), 3 (24))
>   Cell Sets: 1 strata with value/size (1 (24))
>   Face Sets: 5 strata with value/size (1 (4), 2 (4), 3 (4), 4 (4), 5 (4))
>
>
> -dm_view with the current release (commit 4b9a870af96)
>
> DM Object: Generated Mesh 1 MPI process
>   type: plex
> Generated Mesh in 3 dimensions:
>   Number of 0-cells per rank: 14
>   Number of 1-cells per rank: 49
>   Number of 2-cells per rank: 60
>   Number of 3-cells per rank: 24
> Labels:
>   celltype: 4 strata with value/size (0 (14), 6 (24), 3 (60), 1 (49))
>   depth: 4 strata with value/size (0 (14), 1 (49), 2 (60), 3 (24))
>   Cell Sets: 1 strata with value/size (1 (24))
>   Face Sets: 12 strata with value/size (1 (5), 2 (5), 3 (5), 4 (5), 5 (5),
> 6 (1), 7 (1), 8 (1), 9 (1), 10 (1), 11 (1), 12 (1))
>
> I believe the older version printed the correct thing? Has something
> changed in the interpretation of Face Sets?
>
>
> Yes. In the older version, I was only labeling cells, faces, and vertices.
> There were complaints, so I put in the edge labels. If you check, all the
> additional labels are on edges, and checking your .msh file, those edges
> clearly have those labels.
>
>   Thanks,
>
>      Matt
>
>
> Thanks,
> Sharan
>
> *Group - Theory & Simulation*
> *Department of Microstructure Physics & Alloy Design*
>
>
> ------------------------------
> -------------------------------------------------
> Stay up to date and follow us on LinkedIn, Twitter and YouTube.
>
> Max-Planck-Institut für Eisenforschung GmbH
> Max-Planck-Straße 1
> D-40237 Düsseldorf
>
> Handelsregister B 2533
> Amtsgericht Düsseldorf
>
> Geschäftsführung
> Prof. Dr. Gerhard Dehm
> Prof. Dr. Jörg Neugebauer
> Prof. Dr. Dierk Raabe
> Dr. Kai de Weldige
>
> Ust.-Id.-Nr.: DE 11 93 58 514
> Steuernummer: 105 5891 1000
>
>
> Please consider that invitations and e-mails of our institute are
> only valid if they end with …@mpie.de.
> If you are not sure of the validity please contact rco at mpie.de
>
> Bitte beachten Sie, dass Einladungen zu Veranstaltungen und E-Mails
> aus unserem Haus nur mit der Endung …@mpie.de gültig sind.
> In Zweifelsfällen wenden Sie sich bitte an rco at mpie.de
> -------------------------------------------------
>
>
>
> --
> 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/>
>
>
>
> ------------------------------
> -------------------------------------------------
> Stay up to date and follow us on LinkedIn, Twitter and YouTube.
>
> Max-Planck-Institut für Eisenforschung GmbH
> Max-Planck-Straße 1
> D-40237 Düsseldorf
>
> Handelsregister B 2533
> Amtsgericht Düsseldorf
>
> Geschäftsführung
> Prof. Dr. Gerhard Dehm
> Prof. Dr. Jörg Neugebauer
> Prof. Dr. Dierk Raabe
> Dr. Kai de Weldige
>
> Ust.-Id.-Nr.: DE 11 93 58 514
> Steuernummer: 105 5891 1000
>
>
> Please consider that invitations and e-mails of our institute are
> only valid if they end with …@mpie.de.
> If you are not sure of the validity please contact rco at mpie.de
>
> Bitte beachten Sie, dass Einladungen zu Veranstaltungen und E-Mails
> aus unserem Haus nur mit der Endung …@mpie.de gültig sind.
> In Zweifelsfällen wenden Sie sich bitte an rco at mpie.de
> -------------------------------------------------
>


-- 
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/20231212/9eacebef/attachment-0001.html>