void Foam::multiDirRefinement::addCells ( const primitiveMesh& mesh, const Map<label>& splitMap ) { // Construct inverse addressing: from new to original cell. labelList origCell(mesh.nCells(), -1); forAll(addedCells_, cellI) { const labelList& added = addedCells_[cellI]; forAll(added, i) { label slave = added[i]; if (origCell[slave] == -1) { origCell[slave] = cellI; } else if (origCell[slave] != cellI) { FatalErrorIn ( "multiDirRefinement::addCells(const primitiveMesh&" ", const Map<label>&" ) << "Added cell " << slave << " has two different masters:" << origCell[slave] << " , " << cellI << abort(FatalError); } } }