[Dune] reading multiple tags with GmshReader

Benedikt Oswald benedikt.oswald at lspr.ch
Fri Jun 19 09:48:56 CEST 2015


Hi Oliver, yes. Just ask Aleksejs, I am currently on the move. Greetings, Benedikt 

dr. benedikt oswald - Suisse

Am 19.06.2015 um 09:04 schrieb Oliver Sander <sander at igpm.rwth-aachen.de>:

> Gruezi Benedikt,
> I was going to ask you about your GmshReader anyway, because I am currently interested
> in higher-order support.  Is the code somewhere I can look at?
> Thanks,
> Oliver
> 
> Am 19.06.2015 um 08:39 schrieb Benedikt Oswald:
>> Dear Ganesh,
>> 
>> We have implemented a parallel gumsh reader that reads all the tags. Why not try this? Ask Aleksejs, you'll find his coordinates on the mailing list. Also, why not use our dune-curvilineargrid? It is parallel, reads gmsh files including higher order elements.
>> Over and out, Benedikt
>> 
>> 
>> dr. benedikt oswald - Suisse
>> 
>> Am 19.06.2015 um 07:40 schrieb Oliver Sander <sander at igpm.rwth-aachen.de>:
>> 
>>> Am 18.06.2015 um 20:56 schrieb Christian Engwer:
>>>> Dear Ganesh,
>>>> 
>>>> sadly in the current state none of this works. Oli recently removed
>>>> the half-implement reader for more than one tag, as it was never
>>>> finished. Reading other things than elements or boundary elements is
>>>> currently not possible.
>>>> 
>>>> In the issue you reference, nobody ever offered a patch and further
>>>> more the author never replied to my questions.
>>>> 
>>>> But perhaps you can answer some of them:
>>>> * Which information can be added as a tag. Is this restricted to int,
>>>> or can it be arbitrary data, e.g. double values?
>>> 
>>> As I understand the spec, and element can have an arbitrary positive number of
>>> integer tags.
>>> 
>>>> * Is it only possible to use additional tags for cells, or also for
>>>> boundary entries? Reading your comment below I guess there is also
>>>> some way to add tags to other entities!?
>>> 
>>> Yes, that is possible.  The boundarySegmentToPhysicalEntity field actually
>>> returns the first tag of each codim-1 entity.
>>> 
>>>> 
>>>> and which program can I use to generate such files? My Gmsh
>>>> installation only uses the physical entity tag.
>>>> 
>>>> If you are interested in providing patches, I'll happily guide you
>>>> through the code to implement the missing features. Reagrding the
>>>> cells should be fairly simple once we know what type of data can be
>>>> attached. We basically have to revert Olis recent patch and add an
>>>> additional constructor which takes more than one vector.
>>> 
>>> Maybe we can think a little before adding a new constructor.  Just adding
>>> another vector argument still does not cover all the generality of the
>>> file format.
>>> 
>>>>> partition info tags as well. Additionally, is there a way to extract the
>>>>> tags (physical, geometry and partition) for the edge-element (denoted by a
>>>>> type 1 element in Gmsh, eg. a 2-noded edge of a 4-noded quadrilateral) as
>>>>> well?
>>> 
>>> The 'physical' for edge elements are returned in the boundarySegmentToPhysicalEntity
>>> field.
>>> 
>>> --
>>> Oliver
>>> 
>>>>> 
>>>>> I see that this issue was considered before (FS#935) but I do not know if
>>>>> it was resolved and a patch was applied.
>>>>> 
>>>>> 
>>>>> Ganesh
>>>> 
>>>>> _______________________________________________
>>>>> Dune mailing list
>>>>> Dune at dune-project.org
>>>>> http://lists.dune-project.org/mailman/listinfo/dune
>>> 
>>> 
>>> _______________________________________________
>>> Dune mailing list
>>> Dune at dune-project.org
>>> http://lists.dune-project.org/mailman/listinfo/dune
> 
> 




More information about the Dune mailing list