FS#33 Box models do not run with ALUGrid
Metadata
Property |
Value |
Project |
dumux |
Category |
Implicit models |
Reported by |
Bernd Flemisch (bernd@iws.uni-stuttgart.de) |
Reported at |
Dec 16, 2009 16:32 |
Type |
Bug Report |
Version |
Git |
Last edited by |
Bernd Flemisch (bernd@iws.uni-stuttgart.de) |
Last edited at |
May 25, 2010 16:32 |
Closed by |
Bernd Flemisch (bernd@iws.uni-stuttgart.de) |
Closed at |
May 25, 2010 16:32 |
Closed in version |
unknown (Id=0) |
Resolution |
Fixed |
Description
If used with ALUGrid, the box models drive an assertion like
test_1p: /temp/bernd/DUMUX/dune-grid/dune/grid/alugrid/2d/entity_imp.cc:182: const typename Dune::ALU2dGridEntity<0, dim, GridImp>::Geometry& Dune::ALU2dGridEntity<0, dim, GridImp>::geometry() const [with int dim = 2, GridImp = const Dune::ALU2dGrid<2, 2>]: Assertion `item_ != 0' failed.
This happens basically in PDELab's constraints.hh in the ConstraintsVisitNodeMetaProgram. There, it tries to construct a FVElementGeometry on an empty element.
Hard to say whether it is a PDELab or a Dumux bug. Maybe one should try to test PDELab with ALUGrid first to be sure.
The traceback for test_1p is
[tobo:28300] *** Process received signal ***
[tobo:28300] Signal: Aborted (6)
[tobo:28300] Signal code: (-6)
[tobo:28300] [ 0] /lib64/libpthread.so.0 [0x7f6f59ce7c00]
[tobo:28300] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7f6f599b04e5]
[tobo:28300] [ 2] /lib64/libc.so.6(abort+0x180) [0x7f6f599b19b0]
[tobo:28300] [ 3] /lib64/libc.so.6(__assert_fail+0xea) [0x7f6f599a924a]
[tobo:28300] [ 4] ./test_1p(_ZNK4Dune15ALU2dGridEntityILi0ELi2EKNS_9ALU2dGridILi2ELi2EEEE8geometryEv+0x32) [0x74030e]
[tobo:28300] [ 5] ./test_1p(_ZNK4Dune6EntityILi0ELi2EKNS_9ALU2dGridILi2ELi2EEENS_15ALU2dGridEntityEE8geometryEv+0x15) [0x74038d]
[tobo:28300] [ 6] ./test_1p(_ZN4Dune17FVElementGeometryINS_14ALUSimplexGridILi2ELi2EEENS_6PDELab23P1LocalFiniteElementMapIddLi2EEEE6updateERKNS_6EntityILi0ELi2EKNS_9ALU2dGridILi2ELi2EEENS_15ALU2dGridEntityEEE+0x28) [0x7c5e70]
[tobo:28300] [ 7] ./test_1p(_ZNK4Dune23BoxDirichletConstraintsINS_10Properties4TTag15OnePTestProblemEE8boundaryIN25BoundaryIndexHelperPDELabIS3_E5ChildILi0EE4TypeENS_12IntersectionIKNS_9ALU2dGridILi2ELi2EEENS_23LeafIntersectionWrapperEEENS_6PDELab22LocalFunctionSpaceNodeINSH_17GridFunctionSpaceINS_8GridViewINS_25DefaultLeafGridViewTraitsIKNS_14ALUSimplexGridILi2ELi2EEELNS_21PartitionIteratorTypeE4EEEEENSH_23P1LocalFiniteElementMapIddLi2EEENS_37NonoverlappingBoxDirichletConstraintsIS3_EENSH_17ISTLVectorBackendILi1EEENSH_25GridFunctionGeneralMapperEEEEENSH_25ConstraintsTransformationImdEEEEvRKT_RKNSH_20IntersectionGeometryIT0_EERKT1_RT2_+0x73) [0x7dacdb]
[tobo:28300] [ 8] ./test_1p(_ZN4Dune6PDELab23ConstraintsCallBoundaryINS_37NonoverlappingBoxDirichletConstraintsINS_10Properties4TTag15OnePTestProblemEEELb1EE8boundaryIN25BoundaryIndexHelperPDELabIS5_E5ChildILi0EE4TypeENS_12IntersectionIKNS_9ALU2dGridILi2ELi2EEENS_23LeafIntersectionWrapperEEENS0_22LocalFunctionSpaceNodeINS0_17GridFunctionSpaceINS_8GridViewINS_25DefaultLeafGridViewTraitsIKNS_14ALUSimplexGridILi2ELi2EEELNS_21PartitionIteratorTypeE4EEEEENS0_23P1LocalFiniteElementMapIddLi2EEES6_NS0_17ISTLVectorBackendILi1EEENS0_25GridFunctionGeneralMapperEEEEENS0_25ConstraintsTransformationImdEEEEvRKS6_RKT_RKNS0_20IntersectionGeometryIT0_EERKT1_RT2_+0x38) [0x7db064]
[tobo:28300] [ 9] ./test_1p(_ZN4Dune6PDELab31ConstraintsVisitNodeMetaProgramIN25BoundaryIndexHelperPDELabINS_10Properties4TTag15OnePTestProblemEE5ChildILi0EE4TypeELb1ENS0_22LocalFunctionSpaceNodeINS0_17GridFunctionSpaceINS_8GridViewINS_25DefaultLeafGridViewTraitsIKNS_14ALUSimplexGridILi2ELi2EEELNS_21PartitionIteratorTypeE4EEEEENS0_23P1LocalFiniteElementMapIddLi2EEENS_37NonoverlappingBoxDirichletConstraintsIS5_EENS0_17ISTLVectorBackendILi1EEENS0_25GridFunctionGeneralMapperEEEEELb1EE8boundaryINS0_25ConstraintsTransformationImdEENS_12IntersectionIKNS_9ALU2dGridILi2ELi2EEENS_23LeafIntersectionWrapperEEEEEvRKS9_RKSS_RT_RKNS0_20IntersectionGeometryIT0_EE+0x48) [0x7dbd8f]
[tobo:28300] [10] ./test_1p(_ZN4Dune6PDELab32ConstraintsVisitChildMetaProgramI25BoundaryIndexHelperPDELabINS_10Properties4TTag15OnePTestProblemEENS0_23PowerLocalFunctionSpaceINS0_26PowerGridFunctionSpaceBaseINS0_17GridFunctionSpaceINS_8GridViewINS_25DefaultLeafGridViewTraitsIKNS_14ALUSimplexGridILi2ELi2EEELNS_21PartitionIteratorTypeE4EEEEENS0_23P1LocalFiniteElementMapIddLi2EEENS_37NonoverlappingBoxDirichletConstraintsIS5_EENS0_17ISTLVectorBackendILi1EEENS0_25GridFunctionGeneralMapperEEELi1ENS0_32GridFunctionSpaceBlockwiseMapperEEEEELi1ELi0EE8boundaryINS0_25ConstraintsTransformationImdEENS_12IntersectionIKNS_9ALU2dGridILi2ELi2EEENS_23LeafIntersectionWrapperEEEEEvRKS6_RKSS_RT_RKNS0_20IntersectionGeometryIT0_EE+0x47) [0x7dbe1d]
[tobo:28300] [11] ./test_1p(_ZN4Dune6PDELab31ConstraintsVisitNodeMetaProgramI25BoundaryIndexHelperPDELabINS_10Properties4TTag15OnePTestProblemEELb0ENS0_23PowerLocalFunctionSpaceINS0_26PowerGridFunctionSpaceBaseINS0_17GridFunctionSpaceINS_8GridViewINS_25DefaultLeafGridViewTraitsIKNS_14ALUSimplexGridILi2ELi2EEELNS_21PartitionIteratorTypeE4EEEEENS0_23P1LocalFiniteElementMapIddLi2EEENS_37NonoverlappingBoxDirichletConstraintsIS5_EENS0_17ISTLVectorBackendILi1EEENS0_25GridFunctionGeneralMapperEEELi1ENS0_32GridFunctionSpaceBlockwiseMapperEEEEELb0EE8boundaryINS0_25ConstraintsTransformationImdEENS_12IntersectionIKNS_9ALU2dGridILi2ELi2EEENS_23LeafIntersectionWrapperEEEEEvRKS6_RKSS_RT_RKNS0_20IntersectionGeometryIT0_EE+0x2d) [0x7dbe61]
[tobo:28300] [12] ./test_1p(_ZN4Dune6PDELab11constraintsI25BoundaryIndexHelperPDELabINS_10Properties4TTag15OnePTestProblemEENS0_22PowerGridFunctionSpaceINS0_17GridFunctionSpaceINS_8GridViewINS_25DefaultLeafGridViewTraitsIKNS_14ALUSimplexGridILi2ELi2EEELNS_21PartitionIteratorTypeE4EEEEENS0_23P1LocalFiniteElementMapIddLi2EEENS_37NonoverlappingBoxDirichletConstraintsIS5_EENS0_17ISTLVectorBackendILi1EEENS0_25GridFunctionGeneralMapperEEELi1ENS0_32GridFunctionSpaceBlockwiseMapperEEENS0_25ConstraintsTransformationImdEEEEvRT_RKT0_RT1_b+0x2e8) [0x7dc14b]
[tobo:28300] [13] ./test_1p(_ZN15AssemblerPDELabIN4Dune10Properties4TTag15OnePTestProblemEEC1ERNS0_12OnePBoxModelIS3_EERNS0_15OnePTestProblemIS3_EE+0x20f) [0x7de7eb]
[tobo:28300] [14] ./test_1p(_ZN4Dune9BoxSchemeINS_10Properties4TTag15OnePTestProblemENS_12OnePBoxModelIS3_EEE14allocateStuff_Ev+0x3c) [0x7de93c]
[tobo:28300] [15] ./test_1p(_ZN4Dune9BoxSchemeINS_10Properties4TTag15OnePTestProblemENS_12OnePBoxModelIS3_EEE7initialEv+0x2a) [0x7deae8]
[tobo:28300] [16] ./test_1p(_ZN4Dune10BoxProblemINS_10Properties4TTag15OnePTestProblemENS_15OnePTestProblemIS3_EEE4initEv+0x1d) [0x7eae39]
[tobo:28300] [17] ./test_1p(_ZN4Dune11TimeManagerINS_10BoxProblemINS_10Properties4TTag15OnePTestProblemENS_15OnePTestProblemIS4_EEE7EpisodeEE13runSimulationIS6_EEvRT_+0x5e) [0x7ee6a8]
[tobo:28300] [18] ./test_1p(_ZN4Dune10BoxProblemINS_10Properties4TTag15OnePTestProblemENS_15OnePTestProblemIS3_EEE8simulateEdd+0x53) [0x7ee96f]
[tobo:28300] [19] ./test_1p(main+0x7ea) [0x6fc0de]
[tobo:28300] [20] /lib64/libc.so.6(__libc_start_main+0xfd) [0x7f6f5999ca7d]
[tobo:28300] [21] ./test_1p [0x6fb469]