finite element demo with ALUCubeGrid is not working
Using ALUCubeGrid
in finiteelements.py
with
grid_m = grid.get("ALUCubeGrid", dimgrid=2, dimworld=2) # Grid module
gridView = grid.leafGrid("../data/unitcube-2d.dgf", grid_m)
causes the following error:
[misg@manjaro build-cmake]$ python demo/finiteelements.pyc
WARNING (ignored): Could not open file 'alugrid.cfg', using default values 0 < [balance] < 1.2, partitioning method 'ALUGRID_SpaceFillingCurve(9)'.
You are using DUNE-ALUGrid, please don't forget to cite the paper:
Alkaemper, Dedner, Kloefkorn, Nolte. The DUNE-ALUGrid Module, 2016.
Created parallel ALUGrid<2,2,cube,nonconforming> from macro grid file '../data/unitcube-2d.dgf'.
----------------------------------
Number of unknowns: 1089
Assembling...
Traceback (most recent call last):
File "/home/misg/gsoc2016/dune-dev/dune-corepy/demo/finiteelements.py", line 148, in <module>
p1.assemble()
File "/home/misg/gsoc2016/dune-dev/dune-corepy/demo/finiteelements.py", line 66, in assemble
jacInvTra.dot(basis[j].evaluateGradient(p.position), grad2)
File "/home/misg/gsoc2016/dune-dev/dune-corepy/demo/shapefunctions.py", line 45, in __getitem__
return self.f1[i-1]
IndexError: list index out of range
You are using DUNE-ALUGrid, please don't forget to cite the paper:
Alkaemper, Dedner, Kloefkorn, Nolte. The DUNE-ALUGrid Module, 2016.
[misg@manjaro build-cmake]$