NoConstOrderingSize in 2.4 release
I've come across an error in our code (UDG) concerning the ordering/blocking of dofs in the releases/2.4 branch. I have a PowerGridFunctionSpace with GridFunctionSpaces. The FiniteElementMaps of these child GFSs have 0 or k codim 0 dofs (i.e. fixedSize returns false). Using the ISTLBackend, I want to block the dofs of each child GFS together, resulting in a (wrapped) BCRSVector<FieldVector<T,k>>.
Up to commit PDELab/dune-pdelab@f341d78 I used the Chunked decorator and set the Params of the LeafOrderingTag to
struct MyLeafOrderingParams
: public Dune::PDELab::NoConstOrderingSize<true>
, public Dune::PDELab::AllPartitionSelector
{};
This used to worked fine. PDELab/dune-pdelab@f341d78 removed the usage of NoConstOrderingSize. Without this, the ordering somehow assumes a fixed size which results in a wrong size of the GFS. If I revert the two hunks
dune/pdelab/ordering/directleaflocalordering.hh
- _fixed_size_possible = !OrderingTag::no_const_ordering_size;
+ _fixed_size_possible = true;
and
dune/pdelab/ordering/leaflocalordering.hh
- this->_fixed_size_possible = !OrderingTag::no_const_ordering_size;
+ this->_fixed_size_possible = true;
everything works fine again.
Is there a replacement for the NoConstOrderingSize in the 2.4 release or a different way to force the ordering not to assume a fixed size?