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 GridFunctionSpace
s. The FiniteElementMap
s 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?