YaspGrid not working with fem.globalRefine.
The following code does not work, i.e. the grid is not refined, and we should at least raise an Exception that this does not work. Or is it simply a bug.
from dune.grid import cartesianDomain
from dune.grid import yaspGrid
from dune.fem.view import adaptiveLeafGridView
from dune.fem.space import finiteVolume
from dune.fem import globalRefine as globalRefine
domain = cartesianDomain([0,0], [1,1], [4, 4])
view = adaptiveLeafGridView( yaspGrid( domain ) )
spc = finiteVolume(view)
uh = spc.interpolate(1, name="uh")
globalRefine(1, uh)
uh.plot()