exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Additional output before "single" edge error

Browse source
This commit is contained in:
Vuko Vukcevic 2017-11-27 08:32:17 +01:00 committed by Vuko Vukcevic
parent 91a892854e
commit b829763587
1 changed files with 15 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

15
src/dynamicMesh/dynamicMesh/directTopoChange/directTopoChange/directActions/removeFaces.C
View file

@ -874,8 +874,19 @@ void Foam::removeFaces::setRefinement
} }
else else
{ {
// Get edge and all its cells
const edge& e = mesh_.edges()[edgeI]; const edge& e = mesh_.edges()[edgeI];
const labelListList& ec = mesh_.edgeCells();
// Get point cells
const labelListList& pc = mesh_.pointCells();
// Write mesh before termination
mesh_.write();
// Write all edge point cells
// Additionally write all edge cells for debugging
FatalErrorIn("removeFaces::setRefinement") FatalErrorIn("removeFaces::setRefinement")
<< "Problem : edge has too few face neighbours:" << "Problem : edge has too few face neighbours:"
<< eFaces << endl << eFaces << endl
@ -883,6 +894,10 @@ void Foam::removeFaces::setRefinement
<< " vertices:" << e << " vertices:" << e
<< " coords:" << mesh_.points()[e[0]] << " coords:" << mesh_.points()[e[0]]
<< mesh_.points()[e[1]] << mesh_.points()[e[1]]
<< endl
<< "Edge cells: " << ec[edgeI] << nl
<< "First point cells: " << pc[e[0]] << nl
<< "Second point cells: " << pc[e[1]] << nl
<< abort(FatalError); << abort(FatalError);
} }
} }