Changeset 5461
- Timestamp:
- 08/20/10 16:11:32 (15 years ago)
- Location:
- issm/trunk/src/c/objects/Bamg
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Bamg/Mesh.cpp
r5460 r5461 3792 3792 for (i=0;i< NbVerticesOnGeomVertex;i++) 3793 3793 { 3794 BamgVertex *v = VerticesOnGeomVertex[i].m v;3794 BamgVertex *v = VerticesOnGeomVertex[i].meshvertex; 3795 3795 if (v>=vertices && v < ve) 3796 VerticesOnGeomVertex[i].m v=vertices+renu[GetId(v)];3796 VerticesOnGeomVertex[i].meshvertex=vertices+renu[GetId(v)]; 3797 3797 } 3798 3798 … … 3800 3800 for (i=0;i< NbVerticesOnGeomEdge;i++) 3801 3801 { 3802 BamgVertex *v =VerticesOnGeomEdge[i].m v;3802 BamgVertex *v =VerticesOnGeomEdge[i].meshvertex; 3803 3803 if (v>=vertices && v < ve) 3804 VerticesOnGeomEdge[i].m v=vertices+renu[GetId(v)];3804 VerticesOnGeomEdge[i].meshvertex=vertices+renu[GetId(v)]; 3805 3805 } 3806 3806 … … 4036 4036 tstart[i]=0; 4037 4037 for ( j=0;j<NbVerticesOnGeomVertex;j++ ) 4038 tstart[ GetId(VerticesOnGeomVertex[j].m v)]=&vide;4038 tstart[ GetId(VerticesOnGeomVertex[j].meshvertex)]=&vide; 4039 4039 for ( k=0;k<NbVerticesOnGeomEdge;k++ ) 4040 tstart[ GetId(VerticesOnGeomEdge[k].m v)]=&vide;4040 tstart[ GetId(VerticesOnGeomEdge[k].meshvertex)]=&vide; 4041 4041 if(verbose>2) printf(" SmoothingVertex: nb Iteration = %i, Omega=%g\n",nbiter,omega); 4042 4042 for (k=0;k<nbiter;k++) -
issm/trunk/src/c/objects/Bamg/VertexOnGeom.cpp
r5460 r5461 12 12 /*Constructors/Destructors*/ 13 13 /*FUNCTION VertexOnGeom::VertexOnGeom(){{{1*/ 14 VertexOnGeom::VertexOnGeom(): mv(0),abscisse(0){ 14 VertexOnGeom::VertexOnGeom(){ 15 meshvertex=NULL; 16 curvilincoord=0; 15 17 gv=0; 16 18 } 17 19 /*}}}*/ 18 20 /*FUNCTION VertexOnGeom::VertexOnGeom(BamgVertex & m,GeometricalVertex &g){{{1*/ 19 VertexOnGeom::VertexOnGeom(BamgVertex & m,GeometricalVertex &g) : mv(&m),abscisse(-1){ 21 VertexOnGeom::VertexOnGeom(BamgVertex & m,GeometricalVertex &g){ 22 meshvertex=&m; 23 curvilincoord=-1; 20 24 gv=&g; 21 25 } 22 26 /*}}}*/ 23 27 /*FUNCTION VertexOnGeom::VertexOnGeom(BamgVertex & m,GeometricalEdge &g,double s){{{1*/ 24 VertexOnGeom::VertexOnGeom(BamgVertex & m,GeometricalEdge &g,double s) : mv(&m),abscisse(s){ 28 VertexOnGeom::VertexOnGeom(BamgVertex & m,GeometricalEdge &g,double s){ 29 meshvertex=&m; 30 curvilincoord=s; 25 31 ge=&g; 26 32 } … … 31 37 void VertexOnGeom::Set(const VertexOnGeom & rec,const Mesh & Th ,Mesh & ThNew){ 32 38 *this = rec; 33 m v = ThNew.vertices + Th.GetId(mv);39 meshvertex = ThNew.vertices + Th.GetId(meshvertex); 34 40 if (gv) 35 if ( abscisse< 0 )41 if (curvilincoord < 0 ) 36 42 gv = ThNew.Gh.vertices + Th.Gh.GetId(gv); 37 43 else … … 42 48 /*FUNCTION VertexOnGeom::OnGeomVertex{{{1*/ 43 49 int VertexOnGeom::OnGeomVertex()const{ 44 return this? abscisse<0 :0;50 return this? curvilincoord <0 :0; 45 51 } 46 52 /*}}}*/ 47 53 /*FUNCTION VertexOnGeom::OnGeomEdge{{{1*/ 48 54 int VertexOnGeom::OnGeomEdge() const{ 49 return this? abscisse>=0 :0;55 return this? curvilincoord >=0 :0; 50 56 } 51 57 /*}}}*/ 52 58 /*FUNCTION VertexOnGeom::IsRequiredVertex{{{1*/ 53 59 int VertexOnGeom::IsRequiredVertex() { 54 return this? (( abscisse<0 ? (gv?gv->Required():0):0 )) : 0;60 return this? ((curvilincoord<0 ? (gv?gv->Required():0):0 )) : 0; 55 61 } 56 62 /*}}}*/ 57 63 /*FUNCTION VertexOnGeom::SetOn{{{1*/ 58 64 void VertexOnGeom::SetOn(){ 59 m v->GeometricalEdgeHook=this;60 m v->IndexInTriangle=IsVertexOnGeom;65 meshvertex->GeometricalEdgeHook=this; 66 meshvertex->IndexInTriangle=IsVertexOnGeom; 61 67 } 62 68 /*}}}*/ -
issm/trunk/src/c/objects/Bamg/VertexOnGeom.h
r5150 r5461 16 16 public: 17 17 18 BamgVertex* m v;19 double abscisse;18 BamgVertex* meshvertex; 19 double curvilincoord; 20 20 union{ 21 GeometricalVertex* gv; // if abscisse<0;22 GeometricalEdge* ge; // if abscissein [0..1]21 GeometricalVertex* gv; // if curvilincoord <0; 22 GeometricalEdge* ge; // if curvilincoord in [0..1] 23 23 }; 24 24 … … 29 29 30 30 //Operators 31 operator BamgVertex*() const {return m v;}31 operator BamgVertex*() const {return meshvertex;} 32 32 operator GeometricalVertex * () const {return gv;} 33 33 operator GeometricalEdge * () const {return ge;} 34 operator const double & () const {return abscisse;}34 operator const double & () const {return curvilincoord;} 35 35 36 36 //Methods
Note:
See TracChangeset
for help on using the changeset viewer.