problem with importing and load-balancing simplex, non-conforming grid
I try to set up a load-balanced simplex grid for an overlapping computation like this:
using GridType = Dune::ALUGrid<dim,dim,Dune::ALUGridElementType::simplex,Dune::ALUGridRefinementType::nonconforming,Dune::ALUGridMPIComm>;
Dune::GridFactory<GridType> factory;
if(helper.rank() == 0){
Dune::GmshReader<GridType>::read(factory,meshfilename,true,true);
}
std::shared_ptr<GridType> grid(factory.createGrid());
grid->loadBalance();
grid->globalRefine(level);
I set up a GridFunctionSpace with an OverlappingEntitySet and P0ParallelConstraints. Apparently, the load balancing does not work, as calling Dune::PDELab::constraints()
on the GFS constraints a different number of DOFs on every process:
constraints:
40 constrained degrees of freedom
constraints:
39 constrained degrees of freedom
Afterwards, the GridOperators are assembled and throw a segfault.