Changeset 3245
- Timestamp:
- 03/10/10 10:01:20 (15 years ago)
- Location:
- issm/trunk/src/c
- Files:
-
- 1 added
- 1 deleted
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/Bamgx/BamgObjects.h
r3243 r3245 10 10 11 11 #include "objects/Metric.h" 12 #include "objects/DoubleAndInt 4.h"12 #include "objects/DoubleAndInt.h" 13 13 #include "objects/Direction.h" 14 14 #include "objects/Vertex.h" … … 104 104 105 105 inline void Triangles::SetVertexFieldOn(){ 106 for ( longi=0;i<nbv;i++)106 for (int i=0;i<nbv;i++) 107 107 vertices[i].onGeometry=0; 108 for ( longj=0;j<NbVerticesOnGeomVertex;j++ )108 for (int j=0;j<NbVerticesOnGeomVertex;j++ ) 109 109 VerticesOnGeomVertex[j].SetOn(); 110 for ( longk=0;k<NbVerticesOnGeomEdge;k++ )110 for (int k=0;k<NbVerticesOnGeomEdge;k++ ) 111 111 VerticesOnGeomEdge[k].SetOn(); 112 112 } 113 113 inline void Triangles::SetVertexFieldOnBTh(){ 114 for ( longi=0;i<nbv;i++)114 for (int i=0;i<nbv;i++) 115 115 vertices[i].onGeometry=0; 116 for ( longj=0;j<NbVertexOnBThVertex;j++ )116 for (int j=0;j<NbVertexOnBThVertex;j++ ) 117 117 VertexOnBThVertex[j].SetOnBTh(); 118 for ( longk=0;k<NbVertexOnBThEdge;k++ )118 for (int k=0;k<NbVertexOnBThEdge;k++ ) 119 119 VertexOnBThEdge[k].SetOnBTh(); 120 120 -
issm/trunk/src/c/Bamgx/objects/Geometry.h
r3243 r3245 69 69 long Number(const GeometricalEdge * t) const { return t - edges;} 70 70 long Number(const Curve * c) const { return c - curves;} 71 void UnMarkEdges() {for ( longi=0;i<nbe;i++) edges[i].SetUnMark();}71 void UnMarkEdges() {for (int i=0;i<nbe;i++) edges[i].SetUnMark();} 72 72 GeometricalEdge * ProjectOnCurve(const Edge & ,double,Vertex &,VertexOnGeom &) const ; 73 73 GeometricalEdge * Contening(const R2 P, GeometricalEdge * start) const; -
issm/trunk/src/c/Bamgx/objects/QuadTree.cpp
r3243 r3245 94 94 throw ErrorException(__FUNCT__,exprintf("MaxISize <= MaxICoor")); 95 95 } 96 for ( longi=0;i<nbv;i++)96 for (int i=0;i<nbv;i++) 97 97 Add(t->vertices[i]); 98 98 } -
issm/trunk/src/c/Bamgx/objects/SetOfE4.cpp
r3243 r3245 20 20 NbOfEdges=0; 21 21 head = new long [nx]; 22 Edges= new Int 4Edge[nbax];22 Edges= new IntEdge[nbax]; 23 23 24 24 //initialize head (-1 everywhere) -
issm/trunk/src/c/Bamgx/objects/SetOfE4.h
r3243 r3245 6 6 class SetOfEdges4; 7 7 8 class Int 4Edge{8 class IntEdge{ 9 9 friend class SetOfEdges4; 10 10 public: … … 18 18 long nx,nbax,NbOfEdges; 19 19 long* head; 20 Int 4Edge* Edges;20 IntEdge* Edges; 21 21 22 22 public: … … 33 33 34 34 //operators 35 Int 4Edge & operator[](long k){return Edges[k];}35 IntEdge & operator[](long k){return Edges[k];} 36 36 }; 37 37 } -
issm/trunk/src/c/Bamgx/objects/Triangles.cpp
r3243 r3245 884 884 885 885 //loop over all the vertices on edges 886 for ( longi=0;i<nbe;i++){886 for (int i=0;i<nbe;i++){ 887 887 for (int j=0;j<2;j++){ 888 888 … … 1085 1085 1086 1086 //loop over all vertices 1087 for ( longi=0;i<nbv;i++){1087 for (int i=0;i<nbv;i++){ 1088 1088 MatVVP2x2 Vp(vertices[i]); 1089 1089 double lmax=Vp.lmax(); … … 1655 1655 1656 1656 //for all Solutions 1657 for ( longnusol=0;nusol<nbsol;nusol++) {1657 for (int nusol=0;nusol<nbsol;nusol++) { 1658 1658 double smin=ss[nusol],smax=ss[nusol]; 1659 1659 … … 1861 1861 1862 1862 //for all Solution 1863 for ( longnusol=0;nusol<nbsol;nusol++) {1863 for (int nusol=0;nusol<nbsol;nusol++) { 1864 1864 1865 1865 double smin=ss[nusol],smax=ss[nusol]; … … 2208 2208 2209 2209 //check that there is no triangle with 0 determinant 2210 for ( longt = 0; t < nbt; t++){2210 for (int t = 0; t < nbt; t++){ 2211 2211 if (!triangles[t].det) k++; 2212 2212 } … … 2217 2217 //Force Edges 2218 2218 TriangleAdjacent ta(0,0); 2219 for ( longi = 0; i < nbe; i++){2219 for (int i = 0; i < nbe; i++){ 2220 2220 2221 2221 //Force edge i … … 2226 2226 if (nbswp) nbfe++; 2227 2227 if ( nbswp < 0 && k < 5){ 2228 for ( longj = 0; j < nbe; j++){2228 for (int j = 0; j < nbe; j++){ 2229 2229 printf("Edge %i: %i %i\n",j,Number(edges[j][0]),Number(edges[j][1])); 2230 2230 } … … 2238 2238 throw ErrorException(__FUNCT__,exprintf("There are %i lost edges, the boundary might be crossing",k)); 2239 2239 } 2240 for ( longj=0;j<nbv;j++){2240 for (int j=0;j<nbv;j++){ 2241 2241 Nbswap += vertices[j].Optim(1,0); 2242 2242 } … … 2257 2257 long k,it; 2258 2258 2259 for ( longitt=0;itt<nbt;itt++)2259 for (int itt=0;itt<nbt;itt++) 2260 2260 triangles[itt].link=0; // par defaut pas de couleur 2261 2261 … … 2415 2415 mark[it]=triangles[it].link ? -1 : -2; 2416 2416 long inew =0; 2417 for ( longi=0;i<NbSubDomains;i++) {2417 for (int i=0;i<NbSubDomains;i++) { 2418 2418 GeometricalEdge &eg = *Gh.subdomains[i].edge; 2419 2419 subdomains[i].ref = Gh.subdomains[i].ref; … … 2745 2745 // We add the vertices one by one 2746 2746 long NbSwap=0; 2747 for ( longicount=2; icount<nbvb; icount++) {2747 for (int icount=2; icount<nbvb; icount++) { 2748 2748 Vertex *vi = ordre[icount]; 2749 2749 Icoor2 dete[3]; … … 3633 3633 if (verbosity>3) printf(" Begining of insertion process...\n"); 3634 3634 3635 for ( longicount=2; icount<nbv; icount++) {3635 for (int icount=2; icount<nbv; icount++) { 3636 3636 3637 3637 //Get new vertex … … 3702 3702 long k3 = rand()%nbvnew; 3703 3703 //loop over the new points 3704 for ( longis3=0; is3<nbvnew; is3++){3704 for (int is3=0; is3<nbvnew; is3++){ 3705 3705 register long j=nbvold +(k3 = (k3+PrimeNumber)%nbvnew); 3706 3706 register long i=nbvold+is3; … … 3808 3808 3809 3809 long nbqq = (nbt*3)/2; 3810 DoubleAndInt 4 *qq = new DoubleAndInt4[nbqq];3810 DoubleAndInt *qq = new DoubleAndInt[nbqq]; 3811 3811 3812 3812 long i,ij; … … 5276 5276 long NbSwap = 0; 5277 5277 Icoor2 dete[3]; 5278 for ( longi=nbvold;i<nbv;i++) {// for all the new point5278 for (int i=nbvold;i<nbv;i++) {// for all the new point 5279 5279 Vertex & vi = vertices[i]; 5280 5280 vi.i = toI2(vi.r); … … 5317 5317 5318 5318 //initialize all triangles as -1 (outside) 5319 for ( longit=0;it<nbt;it++) reft[it]=-1;5319 for (int it=0;it<nbt;it++) reft[it]=-1; 5320 5320 5321 5321 //loop over all subdomains 5322 for ( longi=0;i<NbSubDomains;i++){5322 for (int i=0;i<NbSubDomains;i++){ 5323 5323 5324 5324 //first triangle of the subdomain i -
issm/trunk/src/c/Bamgx/objects/Triangles.h
r3243 r3245 10 10 #include "Metric.h" 11 11 #include "SetOfE4.h" 12 #include "DoubleAndInt 4.h"12 #include "DoubleAndInt.h" 13 13 #include "Direction.h" 14 14 #include "Vertex.h" -
issm/trunk/src/c/Makefile.am
r3240 r3245 336 336 ./Bamgx/objects/Curve.h\ 337 337 ./Bamgx/objects/Direction.h\ 338 ./Bamgx/objects/DoubleAndInt 4.h\338 ./Bamgx/objects/DoubleAndInt.h\ 339 339 ./Bamgx/objects/Edge.h\ 340 340 ./Bamgx/objects/GeometricalEdge.cpp\ … … 702 702 ./Bamgx/objects/Curve.h\ 703 703 ./Bamgx/objects/Direction.h\ 704 ./Bamgx/objects/DoubleAndInt 4.h\704 ./Bamgx/objects/DoubleAndInt.h\ 705 705 ./Bamgx/objects/Edge.h\ 706 706 ./Bamgx/objects/GeometricalEdge.cpp\
Note:
See TracChangeset
for help on using the changeset viewer.