← Back to team overview

yade-dev team mailing list archive

[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;
 }