[dune-pdelab] PDELab's testordering is broken

Christoph GrĂ¼ninger christoph.grueninger at iws.uni-stuttgart.de
Mon Jun 16 09:31:39 CEST 2014


Hi Steffen et alii,
for me PDELab's testordering fails with this assert:


Created serial ALUGrid<2,2,simplex,conforming>.

4
testordering: 
../../../dune/pdelab/ordering/leaforderingbase.hh:176: void 
Dune::PDELab::LeafOrderingBase<Dune::PDELab::DirectLeafLocalOrdering<Dune::PDELab::LeafOrderingTag<Dune::PDELab::DefaultLeafOrderingParams>, 
Dune::PDELab::VariableMonomLocalFiniteElementMap<Dune::SingleCodimSingleGeomTypeMapper<Dune::GridView<Dune::DefaultLeafGridViewTraits<const 
Dune::ALUGrid<2, 2, 0, 0, Dune::No_Comm>, 4> >, 0>, float, 
double, 2, 6>, 
Dune::GridView<Dune::DefaultLeafGridViewTraits<const 
Dune::ALUGrid<2, 2, 0, 0, Dune::No_Comm>, 4> >, 
Dune::PDELab::DOFIndex<unsigned long, 1, 2>, 
Dune::PDELab::MultiIndex<unsigned long, 1> > 
>::map_lfs_indices(const ItIn, const ItIn, ItOut) const 
[LocalOrdering = 
Dune::PDELab::DirectLeafLocalOrdering<Dune::PDELab::LeafOrderingTag<Dune::PDELab::DefaultLeafOrderingParams>, 
Dune::PDELab::VariableMonomLocalFiniteElementMap<Dune::SingleCodimSingleGeomTypeMapper<Dune::GridView<Dune::DefaultLeafGridViewTraits<const 
Dune::ALUGrid<2, 2, 0, 0, Dune::No_Comm>, 4> >, 0>, float, 
double, 2, 6>, 
Dune::GridView<Dune::DefaultLeafGridViewTraits<const 
Dune::ALUGrid<2, 2, 0, 0, Dune::No_Comm>, 4> >, 
Dune::PDELab::DOFIndex<unsigned long, 1, 2>, 
Dune::PDELab::MultiIndex<unsigned long, 1> >, ItIn = 
__gnu_cxx::__normal_iterator<Dune::PDELab::DOFIndex<unsigned 
long, 1, 2> *, std::vector<Dune::PDELab::DOFIndex<unsigned 
long, 1, 2>, std::allocator<Dune::PDELab::DOFIndex<unsigned 
long, 1, 2> > > >, ItOut = 
__gnu_cxx::__normal_iterator<Dune::PDELab::MultiIndex<unsigned 
long, 1> *, std::vector<Dune::PDELab::MultiIndex<unsigned 
long, 1>, std::allocator<Dune::PDELab::MultiIndex<unsigned 
long, 1> > > >]: Assertion 
`localOrdering()._gt_used[geometry_type_index]' failed.


A related remark: In the Makefile.am testordering does not get 
added the
ALUGrid flags, but it compiles nevertheless.

Everthing else in typetree and PDELab looks good to me. I have 
not yet
checked the howto.

Bye
Christoph




More information about the dune-pdelab mailing list