[petsc-users] "Exact-ish" preallocation for (mostly) unstructured mesh with periodic boundary conditions?

[Ramsey, James J CIV (US)]

________________________________________

for ( ... loop over elements ...) {

   wrapNodeIds(*mesh, nIdVec);  // Accounts for periodic boundary conditions
   mesh->nodes2globalDOF(nIdVec, globalDOFsPerElem);

________________________________________

Sorry, in cutting and pasting snippets of code, I forgot to indicate that "nIdVec" here is the vector of node labels of the current element.