kronbichler/ceed_benchmarks_dealii

ceed benchmark evaluation with deal.II library

ceed_benchmarks_dealii's Issues

Move re-numeration algorithm into deal.II

ping @arosch

Renumbering algorithm throws exception

If I enable the renumbering e.g. in bp3/bench.cc

ceed_benchmarks_dealii/bp3/bench.cc

Lines 83 to 86 in 33e5ca1

	// renumber Dofs to minimize the number of partitions in import indices of
	// partitioner
	//Renumber<dim, double> renum(0, 1, 2);
	//renum.renumber(dof_handler, constraints, mf_data);

and run, I get the following error:

$ ./bench 4 5
terminate called after throwing an instance of 'dealii::StandardExceptions::ExcMessage'
  what():  
--------------------------------------------------------
An error occurred in line <132> of file </home/kronbichler/Work/deal_tests/ceed_benchmarks/bp3/../common_code/renumber_dofs_for_mf.h> in function
    void Renumber<dim, Number, VectorizedArrayType>::renumber(dealii::DoFHandler<dim>&, const dealii::AffineConstraints<double>&, const typename dealii::MatrixFree<dim, Number, VectorizedArrayType>::AdditionalData&) const [with int dim = 3; Number = double; VectorizedArrayType = dealii::VectorizedArray<double, 4>; typename dealii::MatrixFree<dim, Number, VectorizedArrayType>::AdditionalData = dealii::MatrixFree<3, double, dealii::VectorizedArray<double, 4> >::AdditionalData]
The violated condition was: 
    new_numbers.size() == dof_handler.locally_owned_dofs().n_elements()
Additional information: 
    Dimension mismatch 1575 vs 2601
--------------------------------------------------------

It seems that the Renumber::grouping function misses to add some of the unknowns for this setting.

What version of Deal.II is needed?

I tried with current 'master' (182e974ca9e2051625c4c19a160c3ad28285a3db) which is described as 9.0.0-pre but have template errors.