LevelGridView examples

[1]:

from dune.grid import cartesianDomain
from dune.alugrid import aluConformGrid, aluSimplexGrid
domain = cartesianDomain([0,0],[1,1],[1,1])

# first construct a grid using quartering as refinement strategy
aluView = aluSimplexGrid(domain)
hGrid = aluView.hierarchicalGrid
hGrid.globalRefine(5)
for level in range(hGrid.maxLevel):
    print("level:",level, "number of elements:",hGrid.levelView(level).size(0))

# now construct a grid using bisection as refinement strategy
aluView = aluConformGrid(domain)
hGrid = aluView.hierarchicalGrid
hGrid.globalRefine(10)
for level in range(hGrid.maxLevel):
    print("level:",level, "number of elements:",hGrid.levelView(level).size(0))

level: 0 number of elements: 2
level: 1 number of elements: 8
level: 2 number of elements: 32
level: 3 number of elements: 128
level: 4 number of elements: 512
level: 0 number of elements: 2
level: 1 number of elements: 4
level: 2 number of elements: 8
level: 3 number of elements: 16
level: 4 number of elements: 32
level: 5 number of elements: 64
level: 6 number of elements: 128
level: 7 number of elements: 256
level: 8 number of elements: 512
level: 9 number of elements: 1024

This page was generated from the notebook levelgridview_nb.ipynb and is part of the tutorial for the dune-fem python bindings DOI