No boundary projection for 3d, simplex, conforming grids.
Using the most recent master branch (updated today), it appears that the boundary projection block from the attached dgf
file is ignored.
I've also attached test code test.cc
which computes the volume of the unit ball by refining a coarse triangulation specified in ball.dgf
using both alugrid and alberta.
The output I get is:
$ ./test
Dune :: ALUGrid< 3, 3, Dune::simplex, Dune::conforming >
Loading macro bulk grid: ball.dgf
Created serial ALUGrid<3,3,simplex,conforming> from macro grid file 'ball.dgf'.
level: 0 elements: 48 volume: 1.5396 error: 2.64919
level: 1 elements: 384 volume: 1.5396 error: 2.64919
level: 2 elements: 3072 volume: 1.5396 error: 2.64919
level: 3 elements: 24576 volume: 1.5396 error: 2.64919
level: 4 elements: 196608 volume: 1.5396 error: 2.64919
level: 5 elements: 1572864 volume: 1.5396 error: 2.64919
Dune :: AlbertaGrid< 3 >
Loading macro bulk grid: ball.dgf
macro_test_3d: WARNING in ./../3d/macro_3d.c, line 368
not implemented for 3d yet: no check is performed
check_mesh: checking mesh 'DUNE AlbertaGrid'
checking done; no error detected
level: 0 elements: 48 volume: 3.26599 error: 0.922804
level: 1 elements: 384 volume: 3.91975 error: 0.269036
level: 2 elements: 3072 volume: 4.11879 error: 0.0699986
level: 3 elements: 24576 volume: 4.17111 error: 0.0176772
level: 4 elements: 196608 volume: 4.18436 error: 0.0044305
level: 5 elements: 1572864 volume: 4.18768 error: 0.00110833