[Dune] Gmsh reader crashes when using the grid factory read for ALUSimpleGRid<3, 3> with boundary id's
Oliver Sander
sander at mi.fu-berlin.de
Mon Oct 12 09:29:16 CEST 2009
Hi Benedikt!
Please file this as a bug, even though it may not be one. I'll have a look.
--
Oliver
Andreas Dedner schrieb:
> I do not know anything about the gmsh reader...
>
> Oswald Benedikt wrote:
>
>> Hi Andreas, thanks. Actually, importing boundary Id's is essential
>> when we read meshes in the GMSH format.
>>
>> Would it be a lot work to implement it ?
>>
>> Greetings, Benedikt
>>
>>
>>
>> ------------------------------------------------------------------------------------------------------------------
>> Benedikt Oswald, Dr. sc. techn., dipl. El. Ing. ETH, www.psi.ch, Computational Accelerator Scientist
>> Paul Scherrer Institute (PSI), CH-5232 Villigen, Suisse, benedikt.oswald at psi.ch, +41(0)56 310 32 12
>> "Passion is required for any great work, and for the Revolution passion and audacity are required in big doses.",
>> Ernesto 'Che' Guevara, Letter to his parents.
>> http://maxwell.psi.ch/amaswiki/index.php/User:BenediktOswald
>> ------------------------------------------------------------------------------------------------------------------
>>
>>
>>
>>
>> -----Original Message-----
>> From: dune-bounces+benedikt.oswald=psi.ch at dune-project.org on behalf of Andreas Dedner
>> Sent: Sun 11.10.2009 15:13
>> To: Dune
>> Subject: Re: [Dune] Gmsh reader crashes when using the grid factory read for ALUSimpleGRid<3, 3> with boundary id's
>>
>> Hi Benedict,
>> that is because Gmesh uses a method on the Factory which is non standard
>> - and
>> not implemented in this form in ALU.
>> As I see it, the gmsh reader can not be used with ALU if boundary ids
>> are required.
>> Andreas
>>
>>
>>
>> Oswald Benedikt wrote:
>>
>>
>>> Hi Dune developers, I just tried to use the gmsh reader to read gmsh formatted file, including
>>> boundary and element id, in this way (using the Dune head version):
>>>
>>> typedef Dune::ALUSimplexGrid< 3, 3 > ALUSimplexGridType;
>>>
>>>
>>> /** \brief provide a grid factory object for a grid of the ALUGSimplexGrid<3,3> type */
>>> Dune::GridFactory<ALUSimplexGrid< 3, 3 > > factory;
>>>
>>>
>>> /** \brief open the GMSH formatted tetrahedral mesh file */
>>>
>>> /** \brief open the GMSH formatted tetrahedral mesh file into a grid factory */
>>> std::vector<int> boundaryId2physicalEntity;
>>> boundaryId2physicalEntity.clear();
>>>
>>> std::vector<int> elementIndex2PhysicalEntity;
>>> elementIndex2PhysicalEntity.clear();
>>>
>>> Dune::GmshReader< ALUSimplexGridType >::read<int>(factory,
>>> GMSH_FILE_NAME,
>>> boundaryId2physicalEntity,
>>> elementIndex2PhysicalEntity,
>>> true,
>>> true);
>>>
>>>
>>> When I specify that the boundary id's of the GMSH file should be read, the code crashes as follows:
>>>
>>> version 2 Gmsh file detected
>>> file contains 2645 nodes
>>> file contains 15119 elements
>>> number of real vertices = 2645
>>> number of boundary elements = 328
>>> number of elements = 14791
>>> terminate called after throwing an instance of 'Dune::GridError'
>>> [benedikt-oswalds-macbook-pro:83658] *** Process received signal ***
>>> [benedikt-oswalds-macbook-pro:83658] Signal: Abort trap (6)
>>> [benedikt-oswalds-macbook-pro:83658] Signal code: (0)
>>> [benedikt-oswalds-macbook-pro:83658] [ 0] 2 libSystem.B.dylib 0x912bc2bb _sigtramp + 43
>>> [benedikt-oswalds-macbook-pro:83658] [ 1] 3 ??? 0xffffffff 0x0 + 4294967295
>>> [benedikt-oswalds-macbook-pro:83658] [ 2] 4 libSystem.B.dylib 0x9133023a raise + 26
>>> [benedikt-oswalds-macbook-pro:83658] [ 3] 5 libSystem.B.dylib 0x9133c679 abort + 73
>>> [benedikt-oswalds-macbook-pro:83658] [ 4] 6 libstdc++.6.dylib 0x96366005 _ZN9__gnu_cxx27__verbose_terminate_handlerEv + 445
>>> [benedikt-oswalds-macbook-pro:83658] [ 5] 7 libstdc++.6.dylib 0x9636410c __gxx_personality_v0 + 1108
>>> [benedikt-oswalds-macbook-pro:83658] [ 6] 8 libstdc++.6.dylib 0x9636414b _ZSt9terminatev + 29
>>> [benedikt-oswalds-macbook-pro:83658] [ 7] 9 libstdc++.6.dylib 0x96364261 __cxa_throw + 101
>>> [benedikt-oswalds-macbook-pro:83658] [ 8] 10 gmshreadertest 0x0005f99b _ZN4Dune20GridFactoryInterfaceINS_14ALUSimplexGridILi3ELi3EEEE21insertBoundarySegmentESt6vectorIjSaIjEEPKNS_15BoundarySegmentILi3ELi3EEE + 277
>>> [benedikt-oswalds-macbook-pro:83658] [ 9] 11 gmshreadertest 0x000c0950 _ZN4Dune13GmshReaderImpINS_14ALUSimplexGridILi3ELi3EEELi3EE4readIiEEvRNS_11GridFactoryIS2_EERKSsbbRSt6vectorIT_SaISB_EESE_ + 11326
>>> [benedikt-oswalds-macbook-pro:83658] [10] 12 gmshreadertest 0x000c1523 _ZN4Dune10GmshReaderINS_14ALUSimplexGridILi3ELi3EEEE4readERNS_11GridFactoryIS2_EERKSsbb + 147
>>> [benedikt-oswalds-macbook-pro:83658] [11] 13 gmshreadertest 0x00003802 main + 3248
>>> [benedikt-oswalds-macbook-pro:83658] [12] 14 gmshreadertest 0x00001cc6 start + 54
>>> [benedikt-oswalds-macbook-pro:83658] *** End of error message ***
>>>
>>>
>>>
>>> Any idea, on what I do wrong ?
>>>
>>> I had a look at the files alu3dgridfactory.hh and *.cc but could not find a method
>>> called insertBoundarySegment there.
>>>
>>> Thanks for any suggestions! Benedikt
>>>
>>>
>>>
>>>
>>> ------------------------------------------------------------------------------------------------------------------
>>> Benedikt Oswald, Dr. sc. techn., dipl. El. Ing. ETH, www.psi.ch, Computational Accelerator Scientist
>>> Paul Scherrer Institute (PSI), CH-5232 Villigen, Suisse, benedikt.oswald at psi.ch, +41(0)56 310 32 12
>>> "Passion is required for any great work, and for the Revolution passion and audacity are required in big doses.",
>>> Ernesto 'Che' Guevara, Letter to his parents.
>>> http://maxwell.psi.ch/amaswiki/index.php/User:BenediktOswald
>>> ------------------------------------------------------------------------------------------------------------------
>>>
>>>
>>>
>>> _______________________________________________
>>> 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
>>
>>
>>
>
>
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune
>
--
************************************************************************
* Oliver Sander ** email: sander at mi.fu-berlin.de *
* Freie Universität Berlin ** phone: + 49 (30) 838 75348 *
* Institut für Mathematik ** URL : page.mi.fu-berlin.de/~sander *
* Arnimallee 6 ** -------------------------------------*
* 14195 Berlin, Germany ** Member of MATHEON (www.matheon.de) *
************************************************************************
More information about the Dune
mailing list