dune/alugrid/3d/indexsets.hh · d4912038d11030f0f3835fc03ac9f24374338443 · extensions / dune-alugrid

id type = int becomes quickly too small, e.g., assert in · d4912038

Andreas Dedner authored Nov 14, 2019

dune/alugrid/3d/indexsets.hh:504 fails
alugrid_assert ( creatorNumber < nEntitiesFactor[creatorCodim] * nElements );
with for example nEntitiesFactor[creatorCodim]=2 and
const typename IdType::IntegerType nElements = std::pow(nChildren,level);
                                             = 5^13 = 1220703125
The compare in the assert is then done after converting to int - which causes problems

The question is if unsigned int is big enough or if we want unsigned long?

d4912038