[Dune] a problem with refinement

arya fallahi arya.fallahi at psi.ch
Wed Jun 22 11:43:42 CEST 2011


Dear Dune support team,

I have a problem with refining a very simple grid, introduced by a dgf file. When I try to a refinement using  " alusimplexgrid->globalRefine(1); " The code crashes and returns the following message:

Assertion failed: (fabs( quadraturTetra3D < VolumeCalc > ( LinearMapping ( this->myvertex(0)->Point(), this->myvertex(1)->Point(), this->myvertex(2)->Point(), this->myvertex(3)->Point())).integrate1 (0.0) - _volume ) < 1e-10), function TetraTop, file serial/gitter_tetra_top.h, line 1114.
[pcp042851:89732] *** Process received signal ***
[pcp042851:89732] Signal: Abort trap (6)
[pcp042851:89732] Signal code:  (0)
[pcp042851:89732] [ 0] 2   libSystem.B.dylib                   0x00007fff8250666a _sigtramp + 26
[pcp042851:89732] [ 1] 3   libSystem.B.dylib                   0x00007fff704855c0 __stack_chk_guard + 0
[pcp042851:89732] [ 2] 4   libSystem.B.dylib                   0x00007fff82581c9c __pthread_markcancel + 0
[pcp042851:89732] [ 3] 5   hades3deltd                         0x000000010033c30d _ZN12ALUGridSpace8TetraTopINS_14GitterBasisPll10ObjectsPll13TetraEmptyPllEEC1EiPNS_6Gitter9Geometric6hface3EiS8_iS8_iS8_iPS4_id + 2109
[pcp042851:89732] [ 4] 6   hades3deltd                         0x000000010033d1a0 _ZN12ALUGridSpace8TetraTopINS_14GitterBasisPll10ObjectsPll13TetraEmptyPllEE9splitISO8Ev + 3632
[pcp042851:89732] [ 5] 7   hades3deltd                         0x000000010033eb85 _ZN12ALUGridSpace8TetraTopINS_14GitterBasisPll10ObjectsPll13TetraEmptyPllEE15refineImmediateENS_6Gitter9Geometric9TetraRuleE + 325
[pcp042851:89732] [ 6] 8   hades3deltd                         0x000000010033f308 _ZN12ALUGridSpace8TetraTopINS_14GitterBasisPll10ObjectsPll13TetraEmptyPllEE6refineEv + 488
[pcp042851:89732] [ 7] 9   hades3deltd                         0x000000010028afbc _ZN12ALUGridSpace6Gitter6refineEv + 188
[pcp042851:89732] [ 8] 10  hades3deltd                         0x00000001002d3410 _ZN12ALUGridSpace13GitterDunePll6refineEv + 1440
[pcp042851:89732] [ 9] 11  hades3deltd                         0x00000001002c56a9 _ZN12ALUGridSpace13GitterDunePll25adaptWithoutLoadBalancingEv + 121
[pcp042851:89732] [10] 12  hades3deltd                         0x00000001002a905c _ZN12ALUGridSpace13GitterDunePll9duneAdaptERNS_6Gitter20AdaptRestrictProlongE + 44
[pcp042851:89732] [11] 13  hades3deltd                         0x0000000100023107 T.94442 + 503
[pcp042851:89732] [12] 14  hades3deltd                         0x000000010002540c main + 5580
[pcp042851:89732] [13] 15  hades3deltd                         0x00000001000071b4 start + 52
[pcp042851:89732] [14] 16  ???                                 0x0000000000000002 0x0 + 2
[pcp042851:89732] *** End of error message ***
Abort trap

Since, the file is a ver simple 12 tetrahedra grid and I can not find any errors in it, I think it may be a problem with dune. I will be thankful if you have a look at the problem and let me know what is mainly the source of error.

Best regards
Arya Fallahi, 
Dr. sc. techn.  ETH, 
Computational Accelerator Scientist 
Paul Scherrer Institute (PSI), 
CH-5232 Villigen, Switzerland, 
arya.fallahi at psi.ch, 
+41(0)56 310 50 17

P.S. The dgf file is as follows:

DGF
vertex
  0.000000e+00   0.000000e+00   0.000000e+00
  0.000000e+00   1.000000e+03   0.000000e+00
  1.000000e+03   0.000000e+00   0.000000e+00
  0.000000e+00   0.000000e+00   1.000000e+03
  0.000000e+00   1.000000e+03   1.000000e+03
  1.000000e+03   0.000000e+00   1.000000e+03
  4.997688e+02   5.002312e+02   0.000000e+00
  4.997796e+02   5.002204e+02   1.000000e+03
  0.000000e+00   5.000000e+02   5.000000e+02
  5.000000e+02   0.000000e+00   5.000000e+02
#
simplex               %  tetrahedra
parameters 1
     1      8      0      6          101
     6      8      9      7          101
     8      3      7      4          101
     6      9      0      2          101
     5      9      3      7          101
     9      3      8      0          101
     6      8      0      9          101
     3      9      8      7          101
     7      6      8      4          101
     1      4      8      6          101
     7      2      9      6          101
     5      2      9      7          101
#
boundarysegments
   105            0      6      2
   105            0      1      6
   105            3      4      7
   105            3      7      5
   105            7      4      6
   105            1      6      4
   105            5      7      2
   105            7      6      2
   105            8      4      3
   105            1      4      8
   105            3      0      8
   105            1      8      0
   105            5      3      9
   105            5      9      2
   105            9      3      0
   105            9      0      2
#
BOUNDARYDOMAIN
default 1   % all other boundary segments have id 1
#
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20110622/a522547d/attachment.htm>


More information about the Dune mailing list