Changeset 17507
- Timestamp:
- 03/20/14 16:34:36 (11 years ago)
- Location:
- issm/trunk-jpl/src/c/bamg
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/bamg/BamgQuadtree.cpp
r16521 r17507 150 150 151 151 /*Get integer coodinate of current point w*/ 152 registerlong i=w.i.x, j=w.i.y;152 long i=w.i.x, j=w.i.y; 153 153 154 154 /*Initialize level*/ 155 registerlong level=MaxISize;155 long level=MaxISize; 156 156 157 157 /*Get inital box (the largest)*/ … … 426 426 int pi[ MaxDepth ]; 427 427 Icoor1 ii[ MaxDepth ], jj [ MaxDepth]; 428 registerint l=0; // level429 registerBamgQuadtreeBox * b;428 int l=0; // level 429 BamgQuadtreeBox * b; 430 430 Icoor1 hb = MaxISize; 431 431 Icoor1 i0=0,j0=0; … … 444 444 b= pb[l]; 445 445 while (pi[l]--){ 446 registerint k = pi[l];446 int k = pi[l]; 447 447 448 448 if (b->nbitems>0){ // BamgVertex BamgQuadtreeBox none empty … … 457 457 } 458 458 else{ // Pointer BamgQuadtreeBox 459 registerBamgQuadtreeBox *b0=b;459 BamgQuadtreeBox *b0=b; 460 460 if ((b=b->b[k])){ 461 461 hb >>=1 ; // div by 2 462 registerlong iii = ii[l]+I_IJ(k,hb);463 registerlong jjj = jj[l]+J_IJ(k,hb);462 long iii = ii[l]+I_IJ(k,hb); 463 long jjj = jj[l]+J_IJ(k,hb); 464 464 465 465 if (INTER_SEG(iii,iii+hb,i-hx,i+hx) && INTER_SEG(jjj,jjj+hb,j-hy,j+hy)){ -
issm/trunk-jpl/src/c/bamg/BamgVertex.cpp
r15104 r17507 125 125 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/Smoothing)*/ 126 126 127 registerBamgVertex* s=this;127 BamgVertex* s=this; 128 128 BamgVertex &vP = *s,vPsave=vP; 129 129 130 registerTriangle* tbegin= t , *tria = t , *ttc;131 132 registerint k=0,kk=0,j = EdgesVertexTriangle[IndexInTriangle][0],jc;130 Triangle* tbegin= t , *tria = t , *ttc; 131 132 int k=0,kk=0,j = EdgesVertexTriangle[IndexInTriangle][0],jc; 133 133 R2 P(s->r),PNew(0,0); 134 134 do { -
issm/trunk-jpl/src/c/bamg/Geometry.cpp
r16237 r17507 598 598 long *po=p; 599 599 long n=*p; 600 registerfloat angleold=-1000 ; // angle = - infinity600 float angleold=-1000 ; // angle = - infinity 601 601 ord=0; exch=0; 602 602 … … 878 878 double sg; 879 879 if (eg0 == eg1) { 880 registerdouble s0=vg0,s1=vg1;880 double s0=vg0,s1=vg1; 881 881 sg = s0*(1.0-s) + s*s1; 882 882 on=eg0; -
issm/trunk-jpl/src/c/bamg/ListofIntersectionTriangles.cpp
r16233 r17507 83 83 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/NewItem)*/ 84 84 85 registerint n;85 int n; 86 86 R2 x(0,0); 87 87 if ( d0) x = (*tt)[0].r * d0; … … 97 97 lIntTria[Size].x = x; 98 98 Metric m0,m1,m2; 99 registerBamgVertex * v;99 BamgVertex * v; 100 100 if ((v=(*tt)(0))) m0 = v->m; 101 101 if ((v=(*tt)(1))) m1 = v->m; … … 111 111 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/NewItem)*/ 112 112 113 registerint n;113 int n; 114 114 if(!Size || Norme2_2(lIntTria[Size-1].x-A)) { 115 115 if (Size==MaxSize) ReShape(); … … 183 183 void ListofIntersectionTriangles::ReShape(){ 184 184 185 registerint newsize = MaxSize*2;185 int newsize = MaxSize*2; 186 186 IntersectionTriangles* nw = new IntersectionTriangles[newsize]; 187 187 _assert_(nw); -
issm/trunk-jpl/src/c/bamg/Mesh.cpp
r16989 r17507 2933 2933 //loop over the new points 2934 2934 for (int is3=0; is3<nbvnew; is3++){ 2935 registerlong j=nbvold +(k3 = (k3+PrimeNumber)%nbvnew);2936 registerlong i=nbvold+is3;2935 long j=nbvold +(k3 = (k3+PrimeNumber)%nbvnew); 2936 long i=nbvold+is3; 2937 2937 orderedvertices[i]= vertices + j; 2938 2938 orderedvertices[i]->ReferenceNumber=i; … … 3722 3722 3723 3723 long *renu= new long[nbt]; 3724 registerTriangle *t0,*t,*te=triangles+nbt;3725 registerlong k=0,it,i,j;3724 Triangle *t0,*t,*te=triangles+nbt; 3725 long k=0,it,i,j; 3726 3726 3727 3727 for ( it=0;it<nbt;it++) … … 3914 3914 for(i=Head0;i>=0;i=first_np_or_next_t0[ip=i],first_np_or_next_t0[ip]=-1) { 3915 3915 // pour tous les triangles autour du sommet s 3916 registerTriangle* t= vertices[i].t;3916 Triangle* t= vertices[i].t; 3917 3917 if (!t){ 3918 3918 _error_("!t"); … … 4154 4154 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/ConsRefTriangle)*/ 4155 4155 4156 registerTriangle *t0,*t;4157 registerlong k=0, num;4156 Triangle *t0,*t; 4157 long k=0, num; 4158 4158 4159 4159 //initialize all triangles as -1 (outside) … … 4707 4707 L=0;// length of the curve 4708 4708 long i=0;// index of new points on the curve 4709 registerGeomVertex * GA0 = *(*peequi)[k0equi].GeomEdgeHook;4709 GeomVertex * GA0 = *(*peequi)[k0equi].GeomEdgeHook; 4710 4710 BamgVertex *A0; 4711 4711 A0 = GA0->MeshVertexHook; // the vertex in new mesh … … 4781 4781 if (ee[k1].GeomEdgeHook->IsRequiredVertex()) { 4782 4782 _assert_(eeequi[k1equi].GeomEdgeHook->IsRequiredVertex()); 4783 registerGeomVertex * GA1 = *eeequi[k1equi].GeomEdgeHook;4783 GeomVertex * GA1 = *eeequi[k1equi].GeomEdgeHook; 4784 4784 A1=GA1->MeshVertexHook;// the vertex in new mesh 4785 4785 _assert_(A1-vertices>=0 && A1-vertices<nbv); -
issm/trunk-jpl/src/c/bamg/Triangle.cpp
r16237 r17507 231 231 // Copy all the mark 232 232 a &= 3; 233 registerTriangle *tt=adj[a];233 Triangle *tt=adj[a]; 234 234 AdjEdgeIndex [a] &= 55; // remove MarkUnSwap 235 registershort aatt = AdjEdgeIndex[a] & 3;235 short aatt = AdjEdgeIndex[a] & 3; 236 236 if(tt){ 237 237 tt->adj[aatt]=this; … … 255 255 void Triangle::SetHidden(int a){ 256 256 //Get Adjacent Triangle number a 257 registerTriangle* t = adj[a];257 Triangle* t = adj[a]; 258 258 //if it exist 259 259 //C|=D -> C=(C|D) bitwise inclusive OR … … 264 264 void Triangle::SetLocked(int a){ 265 265 //mark the edge as on Boundary 266 registerTriangle * t = adj[a];266 Triangle * t = adj[a]; 267 267 t->AdjEdgeIndex[AdjEdgeIndex[a] & 3] |=4; 268 268 AdjEdgeIndex[a] |= 4; … … 270 270 /*FUNCTION Triangle::SetMarkUnSwap{{{*/ 271 271 void Triangle::SetMarkUnSwap(int a){ 272 registerTriangle * t = adj[a];272 Triangle * t = adj[a]; 273 273 t->AdjEdgeIndex[AdjEdgeIndex[a] & 3] |=8; 274 274 AdjEdgeIndex[a] |=8 ; … … 282 282 /*FUNCTION Triangle::SetUnMarkUnSwap{{{*/ 283 283 void Triangle::SetUnMarkUnSwap(int a){ 284 registerTriangle * t = adj[a];284 Triangle * t = adj[a]; 285 285 t->AdjEdgeIndex[AdjEdgeIndex[a] & 3] &=55; // 23 + 32 286 286 AdjEdgeIndex[a] &=55 ; … … 292 292 if(a/4 !=0) return 0;// arete lock or MarkUnSwap 293 293 294 registerTriangle *t1=this,*t2=adj[a];// les 2 triangles adjacent295 registershort a1=a,a2=AdjEdgeIndex[a];// les 2 numero de l arete dans les 2 triangles294 Triangle *t1=this,*t2=adj[a];// les 2 triangles adjacent 295 short a1=a,a2=AdjEdgeIndex[a];// les 2 numero de l arete dans les 2 triangles 296 296 if(a2/4 !=0) return 0; // arete lock or MarkUnSwap 297 297 298 registerBamgVertex *sa=t1->vertices[VerticesOfTriangularEdge[a1][0]];299 registerBamgVertex *sb=t1->vertices[VerticesOfTriangularEdge[a1][1]];300 registerBamgVertex *s1=t1->vertices[OppositeVertex[a1]];301 registerBamgVertex *s2=t2->vertices[OppositeVertex[a2]];298 BamgVertex *sa=t1->vertices[VerticesOfTriangularEdge[a1][0]]; 299 BamgVertex *sb=t1->vertices[VerticesOfTriangularEdge[a1][1]]; 300 BamgVertex *s1=t1->vertices[OppositeVertex[a1]]; 301 BamgVertex *s2=t2->vertices[OppositeVertex[a2]]; 302 302 303 303 Icoor2 det1=t1->det , det2=t2->det ; … … 328 328 if(kopt) { 329 329 // critere de Delaunay pure isotrope 330 registerIcoor2 xb1 = sb->i.x - s1->i.x,330 Icoor2 xb1 = sb->i.x - s1->i.x, 331 331 x21 = s2->i.x - s1->i.x, 332 332 yb1 = sb->i.y - s1->i.y, … … 336 336 yba = sb->i.y - sa->i.y, 337 337 y2a = s2->i.y - sa->i.y; 338 registerdouble338 double 339 339 cosb12 = double(xb1*x21 + yb1*y21), 340 340 cosba2 = double(xba*x2a + yba*y2a) ,
Note:
See TracChangeset
for help on using the changeset viewer.