Changeset 5340
- Timestamp:
- 08/18/10 08:36:30 (15 years ago)
- Location:
- issm/trunk/src/c/objects/Bamg
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Bamg/GeometricalEdge.h
r5151 r5340 16 16 long ReferenceNumber; 17 17 long CurveNumber; 18 R2 tg[2]; // the 2 tangentes (tg[0] =0 => no continuity)18 R2 tg[2]; // the 2 tangentes (tg[0] =0 => no continuity) 19 19 GeometricalEdge *Adj[2]; 20 int DirAdj[2];20 int AdjVertexNumber[2]; // for a given vertex, this gives the index of the vertex in the adjacent edge (0 or 1) 21 21 int type; 22 22 -
issm/trunk/src/c/objects/Bamg/Geometry.cpp
r5339 r5340 151 151 edges[i].tg[0]=zerovector; 152 152 edges[i].tg[1]=zerovector; 153 edges[i]. DirAdj[0] = edges[i].DirAdj[1] = -1;153 edges[i].AdjVertexNumber[0] = edges[i].AdjVertexNumber[1] = -1; 154 154 edges[i].Adj[0] = edges[i].Adj[1] = NULL; 155 155 edges[i].type = 0; … … 602 602 603 603 edges[i1].Adj[j1] = edges + i; 604 edges[i1]. DirAdj[j1] = j;604 edges[i1].AdjVertexNumber[j1] = j; 605 605 } 606 606 } … … 624 624 * previous and next vertices connected to current vertex 625 625 * normed by the edge length*/ 626 tg = edges[i].v[1-j]->r - edges[i].Adj[j]->v[1-edges[i]. DirAdj[j]]->r;626 tg = edges[i].v[1-j]->r - edges[i].Adj[j]->v[1-edges[i].AdjVertexNumber[j]]->r; 627 627 ltg= Norme2(tg); 628 628 tg = tg *(lAB/ltg); … … 680 680 GeometricalVertex *b=(*e)(k1); 681 681 if (a == b || b->Required() ) break; 682 k0 = e-> DirAdj[k1];// vertex in next edge682 k0 = e->AdjVertexNumber[k1];// vertex in next edge 683 683 e = e->Adj[k1]; // next edge 684 684 } … … 865 865 ge[--bge] =eg0 = eg0->Adj[direction0]; 866 866 ISSMASSERT(bge>=0 && bge<=mxe); 867 direction0 = 1-( directionge[bge] = tmpge-> DirAdj[direction0]);867 direction0 = 1-( directionge[bge] = tmpge->AdjVertexNumber[direction0]); 868 868 } 869 869 while (eg1 != (GeometricalEdge*) vg1 && (*eg1)(direction1) != (GeometricalVertex*) vg1) { … … 881 881 GeometricalEdge* tmpge = eg1; 882 882 ge[++tge] =eg1 = eg1->Adj[direction1]; 883 directionge[tge]= direction1 = 1-tmpge-> DirAdj[direction1];883 directionge[tge]= direction1 = 1-tmpge->AdjVertexNumber[direction1]; 884 884 ISSMASSERT(tge>=0 && tge<=mxe); 885 885 } -
issm/trunk/src/c/objects/Bamg/Mesh.cpp
r5208 r5340 5110 5110 if (b->Required() ) break; 5111 5111 int kprev=k; 5112 k = e-> DirAdj[kprev];// next vertices5112 k = e->AdjVertexNumber[kprev];// next vertices 5113 5113 e = e->Adj[kprev]; 5114 5114 ISSMASSERT(e); -
issm/trunk/src/c/objects/Bamg/Mesh.h
r5208 r5340 170 170 } 171 171 inline void Adj(GeometricalEdge * & on,int &i){ 172 int j=i;i=on-> DirAdj[i];on=on->Adj[j];172 int j=i;i=on->AdjVertexNumber[i];on=on->Adj[j]; 173 173 } 174 174 inline double qualite(const BamgVertex &va,const BamgVertex &vb,const BamgVertex &vc){ -
issm/trunk/src/c/objects/Bamg/Triangle.h
r5151 r5340 18 18 private: 19 19 BamgVertex *vertices[3]; // 3 vertices if t is triangle, t[i] allowed by access function, (*t)[i] if pointer 20 Triangle *adj[3]; // 3 pointers toward the adjacent triangles20 Triangle *adj[3]; // 3 pointers toward the adjacent triangles 21 21 short AdjEdgeNumber[3]; // edge id in the adjacent triangles. The edge number 1 is the edge number AdjEdgeNumber[1] in the Adjacent triangle 1 22 22
Note:
See TracChangeset
for help on using the changeset viewer.