Changeset 3407 for issm/trunk/src/c/Bamgx/objects/Geometry.cpp
- Timestamp:
- 04/06/10 11:19:09 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/Bamgx/objects/Geometry.cpp
r3335 r3407 120 120 } 121 121 else{ 122 ISSMERROR(exprintf("No Vertex provided"));122 ISSMERROR(exprintf("No MeshVertex provided")); 123 123 } 124 124 … … 438 438 double eps=1e-20; 439 439 QuadTree quadtree; // build quadtree to find duplicates 440 Vertex* v0=vertices;440 MeshVertex* v0=vertices; 441 441 GeometricalVertex* v0g=(GeometricalVertex*) (void*)v0; 442 442 … … 457 457 458 458 //find nearest vertex already present in the quadtree (NULL if empty) 459 Vertex* v=quadtree.NearestVertex(vertices[i].i.x,vertices[i].i.y);459 MeshVertex* v=quadtree.NearestVertex(vertices[i].i.x,vertices[i].i.y); 460 460 461 461 //if there is a vertex found that is to close to vertices[i] -> error … … 465 465 j=vg-v0g; 466 466 //check that the clostest vertex is not itself... 467 if ( v != &( Vertex &) vertices[j]){468 ISSMERROR(exprintf(" v != &( Vertex &) vertices[j]"));467 if ( v != &(MeshVertex &) vertices[j]){ 468 ISSMERROR(exprintf(" v != &(MeshVertex &) vertices[j]")); 469 469 } 470 470 vertices[i].link = vertices + j; … … 830 830 /*}}}1*/ 831 831 /*FUNCTION Geometry::ProjectOnCurve {{{1*/ 832 GeometricalEdge* Geometry::ProjectOnCurve(const Edge &e,double s, Vertex &V,VertexOnGeom &GV) const {832 GeometricalEdge* Geometry::ProjectOnCurve(const Edge &e,double s,MeshVertex &V,VertexOnGeom &GV) const { 833 833 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, MeshGeom.cpp/ProjectOnCurve)*/ 834 834 /*Add a vertex on an existing geometrical edge according to the metrics of the two vertices constituting the edge*/ … … 849 849 850 850 //Get the two vertices of the edge 851 const Vertex &v0=e[0];852 const Vertex &v1=e[1];851 const MeshVertex &v0=e[0]; 852 const MeshVertex &v1=e[1]; 853 853 854 854 //Get position of V0, V1 and vector v0->v1 … … 916 916 917 917 if ((*eg0)(sens0)==(GeometricalVertex*)vg0) 918 vg0=VertexOnGeom(*( Vertex*) vg0,*eg0,sens0); //vg0 = absisce918 vg0=VertexOnGeom(*(MeshVertex*) vg0,*eg0,sens0); //vg0 = absisce 919 919 920 920 if ((*eg1)(sens1)==(GeometricalVertex*)vg1) 921 vg1=VertexOnGeom(*( Vertex*) vg1,*eg1,sens1);921 vg1=VertexOnGeom(*(MeshVertex*) vg1,*eg1,sens1); 922 922 923 923 double sg;
Note:
See TracChangeset
for help on using the changeset viewer.