yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #10596
[Branch ~yade-pkg/yade/git-trunk] Rev 3850: fix a bug in PFV code (fictious vetices not defined correctly)
------------------------------------------------------------
revno: 3850
committer: Bruno Chareyre <bruno.chareyre@xxxxxxxxxxxxxxx>
timestamp: Sun 2014-03-23 03:00:05 +0100
message:
fix a bug in PFV code (fictious vetices not defined correctly)
modified:
lib/triangulation/Network.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/Network.ipp'
--- lib/triangulation/Network.ipp 2014-03-21 18:47:45 +0000
+++ lib/triangulation/Network.ipp 2014-03-23 02:00:05 +0000
@@ -39,15 +39,16 @@
int Network<Tesselation>::detectFacetFictiousVertices (CellHandle& cell, int& j)
{
facetNFictious = 0;
+ int nRealVtx=0;
for (int kk=0; kk<3; kk++) {
if (cell->vertex(facetVertices[j][kk])->info().isFictious) {
if (facetNFictious==0) facetF1=kk; else facetF2=kk;
facetNFictious +=1;
} else {
- if (real_vertex==0) facetRe1=kk;
- else if (real_vertex==1) facetRe2=kk;
- else if (real_vertex==2) facetRe3=kk;
- real_vertex+=1;}}
+ if (nRealVtx==0) facetRe1=kk;
+ else if (nRealVtx==1) facetRe2=kk;
+ else if (nRealVtx==2) facetRe3=kk;
+ nRealVtx+=1;}}
return facetNFictious;
}