- if(surface!=0.) {
-
- matrix[targetCell].insert(std::make_pair(cellSrcIdx, surface));
-
- bool isSrcFaceColinearWithFaceOfTetraTargetCell = false;
- std::set<TriangleFaceKey>::iterator iter;
- for (iter = listOfTetraFacesColinear.begin(); iter != listOfTetraFacesColinear.end(); ++iter)
- {
- if (listOfTetraFacesTreated.count(*iter) != 1)
- {
- isSrcFaceColinearWithFaceOfTetraTargetCell = false;
- break;
- }
- else
- {
- isSrcFaceColinearWithFaceOfTetraTargetCell = true;
- }
- }
-
- if (isSrcFaceColinearWithFaceOfTetraTargetCell)
- {
- DuplicateFacesType::iterator intersectFacesIter = _intersect_faces.find(cellSrcIdx);
- if (intersectFacesIter != _intersect_faces.end())
- {
- intersectFacesIter->second.insert(targetCell);
- }
- else
- {
- std::set<int> targetCellSet;
- targetCellSet.insert(targetCell);
- _intersect_faces.insert(std::make_pair(cellSrcIdx, targetCellSet));
- }
-
- }
-
- }
-
- delete[] polyNodes;
- delete[] polyCoords;
-
+ if(surface!=0.)
+ {
+
+ matrix[targetCell].insert(std::make_pair(cellSrcIdx, surface));
+
+ bool isSrcFaceColinearWithFaceOfTetraTargetCell = false;
+ std::set<TriangleFaceKey>::iterator iter;
+ for (iter = listOfTetraFacesColinear.begin(); iter != listOfTetraFacesColinear.end(); ++iter)
+ {
+ if (listOfTetraFacesTreated.count(*iter) != 1)
+ {
+ isSrcFaceColinearWithFaceOfTetraTargetCell = false;
+ break;
+ }
+ else
+ {
+ isSrcFaceColinearWithFaceOfTetraTargetCell = true;
+ }
+ }
+
+ if (isSrcFaceColinearWithFaceOfTetraTargetCell)
+ {
+ DuplicateFacesType::iterator intersectFacesIter = _intersect_faces.find(cellSrcIdx);
+ if (intersectFacesIter != _intersect_faces.end())
+ {
+ intersectFacesIter->second.insert(targetCell);
+ }
+ else
+ {
+ std::set<int> targetCellSet;
+ targetCellSet.insert(targetCell);
+ _intersect_faces.insert(std::make_pair(cellSrcIdx, targetCellSet));
+ }
+ }
+ }
+ delete [] polyNodes;
+ delete [] polyCoords;