yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #13002
[Branch ~yade-pkg/yade/git-trunk] Rev 4005: return a consistent set of constrictions in the periodic case (some constrictions were previsousl...
------------------------------------------------------------
revno: 4005
committer: bchareyre <bruno.chareyre@xxxxxxxxxxxxxxx>
timestamp: Thu 2017-02-23 15:14:02 +0100
message:
return a consistent set of constrictions in the periodic case (some constrictions were previsously skipped inconsistenty)
modified:
lib/triangulation/FlowBoundingSphere.ipp
--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk
Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription
=== modified file 'lib/triangulation/FlowBoundingSphere.ipp'
--- lib/triangulation/FlowBoundingSphere.ipp 2016-11-17 15:32:26 +0000
+++ lib/triangulation/FlowBoundingSphere.ipp 2017-02-23 14:14:02 +0000
@@ -663,10 +663,9 @@
RTriangulation& Tri = T[currentTes].Triangulation();
vector<double> constrictions;
for (FiniteFacetsIterator f_it=Tri.finite_facets_begin(); f_it != Tri.finite_facets_end();f_it++){
- //in the periodic case, we skip facets with lowest id out of the base period
- if ( ((f_it->first->info().index <= f_it->first->neighbor(f_it->second)->info().index) && f_it->first->info().isGhost)
- || ((f_it->first->info().index >= f_it->first->neighbor(f_it->second)->info().index) && f_it->first->neighbor(f_it->second)->info().isGhost)
- || f_it->first->info().index == 0 || f_it->first->neighbor(f_it->second)->info().index == 0) continue;
+ //in the periodic case, we retain only the facets incident to at least one non-ghost cell
+ if ( (f_it->first->neighbor(f_it->second)->info().isGhost && f_it->first->info().isGhost)
+ || f_it->first->info().index == 0 || f_it->first->neighbor(f_it->second)->info().index == 0) continue;
constrictions.push_back(computeEffectiveRadius(f_it->first, f_it->second));
}
return constrictions;