Changeset 21629


Ignore:
Timestamp:
03/23/17 16:17:46 (8 years ago)
Author:
Mathieu Morlighem
Message:

CHG: minor

Location:
issm/trunk-jpl/src/c/bamg
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/bamg/BamgQuadtree.cpp

    r21625 r21629  
    77
    88namespace bamg {
     9
    910#define ABS(i) ((i)<0 ?-(i) :(i))
    10 #define MAXDEPTH 30
     11#define MAXDEPTH  30
     12#define MAXISIZE  1073741824 //2^30
     13#define MAXICOORD 1073741823 //2^30 - 1 = =111...111 (29 times one)
     14
    1115        /*DOCUMENTATION What is a BamgQuadtree? {{{
    1216         * A Quadtree is a very simple way to group vertices according
    … …  
    6165
    6266                /*Initialize fields*/
    63                 this->MaxISize      = 1073741824; //2^30
    64                 this->MaxICoord     = 1073741823; //2^30 - 1
    6567                this->NbQuadtreeBox = 0;
    6668                this->NbVertices    = 0;
    … …  
    7678
    7779                /*Initialize fields*/
    78                 this->MaxISize      = 1073741824; //2^30
    79                 this->MaxICoord     = 1073741823; //2^30 - 1
    8080                this->NbQuadtreeBox = 0;
    8181                this->NbVertices    = 0;
    … …  
    8888
    8989                /*Check Sizes*/
    90                 _assert_(this->MaxISize>this->MaxICoord);
     90                _assert_(MAXISIZE>MAXICOORD);
    9191
    9292                /*Add all vertices of the mesh*/
    … …  
    114114
    115115                /*Initialize level*/
    116                 long level=MaxISize;
     116                long level=MAXISIZE;
    117117
    118118                /*Get inital box (the largest)*/
    … …  
    236236                if(!this->root->nbitems) return nearest_v;
    237237
    238                 /*Project coordinates (xi,yi) onto [0,MaxICoord-1] x [0,MaxICoord-1]*/
     238                /*Project coordinates (xi,yi) onto [0,MAXICOORD-1] x [0,MAXICOORD-1]*/
    239239                int xi2 = xi;
    240240                int yi2 = yi;
    241241                if(xi<0)        xi2 = 0;
    242                 if(xi>MaxISize) xi2 = MaxICoord;
     242                if(xi>MAXISIZE) xi2 = MAXICOORD;
    243243                if(yi<0)        yi2 = 0;
    244                 if(yi>MaxISize) yi2 = MaxICoord;
     244                if(yi>MAXISIZE) yi2 = MAXICOORD;
    245245
    246246                /*Get inital box (the largest)*/
    … …  
    249249                /*Initialize level and sizes for largest box*/
    250250                int levelbin = (1L<<MAXDEPTH);// = 2^30
    251                 int        h = this->MaxISize;
    252                 int       hb = this->MaxISize;
     251                int        h = MAXISIZE;
     252                int       hb = MAXISIZE;
    253253                int       i0 = 0;
    254254                int       j0 = 0;
    … …  
    439439                int l=0; // level
    440440                BamgQuadtreeBox* b;
    441                 int hb =  MaxISize;
     441                int hb =  MAXISIZE;
    442442                int i0=0,j0=0;
    443443

  • issm/trunk-jpl/src/c/bamg/BamgQuadtree.h

    r21625 r21629  
    3737
    3838                        /*BamgQuadtree public Fields*/
    39                         int              MaxICoord;       // maximum integer coordinate 2^MaxDepth -1
    40                         int              MaxISize;        // maximum integer coordinate 2^MaxDepth
    4139                        BamgQuadtreeBox *root;            // main box
    4240                        long             NbQuadtreeBox;   // total number of boxes

  • issm/trunk-jpl/src/c/bamg/Geometry.cpp

    r21623 r21629  
    103103                         * coefIcoor = (2^30 -1)/D
    104104                         */
    105                         coefIcoor=(MaxICoor)/(Max(pmax.x-pmin.x,pmax.y-pmin.y));
     105                        int MaxICoord = 1073741823;
     106                        coefIcoor=(MaxICoord)/(Max(pmax.x-pmin.x,pmax.y-pmin.y));
    106107                        if(coefIcoor<=0) _error_("coefIcoor should be positive");
    107108                }
    … …  
    909910                 * coefIcoor = (2^30 -1)/D
    910911                 */
    911                 return  I2( (Icoor1) (coefIcoor*(P.x-pmin.x)) ,(Icoor1) (coefIcoor*(P.y-pmin.y)) );
     912                return  I2( (int) (coefIcoor*(P.x-pmin.x)) ,(int) (coefIcoor*(P.y-pmin.y)) );
    912913        }/*}}}*/
    913914        void Geometry::UnMarkEdges() {/*{{{*/

  • issm/trunk-jpl/src/c/bamg/Geometry.h

    r16237 r21629  
    3131                        Curve         *curves;
    3232                        R2             pmin,pmax;             // domain extrema coordinates
    33                         double         coefIcoor;             // coef to integer Icoor1;
     33                        double         coefIcoor;             // coef to integer coordinates;
    3434                        double         MaxCornerAngle;
    3535

  • issm/trunk-jpl/src/c/bamg/Mesh.cpp

    r21623 r21629  
    15751575
    15761576                //Build Gh.coefIcoor
    1577                 Gh.coefIcoor= (MaxICoor)/(Max(Gh.pmax.x-Gh.pmin.x,Gh.pmax.y-Gh.pmin.y));
     1577                int MaxICoord = 1073741823; //2^30 - 1 = =111...111 (29 times one)
     1578                Gh.coefIcoor= (MaxICoord)/(Max(Gh.pmax.x-Gh.pmin.x,Gh.pmax.y-Gh.pmin.y));
    15781579                if (Gh.coefIcoor<=0){
    15791580                        delete [] colorV;
    … …  
    16651666                //unset adj
    16661667                for (i=0;i<nbt;i++){
    1667                         for ( j=0;j<3;j++){
     1668                        for (j=0;j<3;j++){
    16681669                                triangles[i].SetAdj2(j,0,triangles[i].GetAllflag(j));
    16691670                        }
    … …  
    28632864                        double hx,hy;
    28642865                        vi.m.Box(hx,hy);
    2865                         Icoor1 hi=(Icoor1) (hx*coefIcoor),hj=(Icoor1) (hy*coefIcoor);
     2866                        int hi=(int) (hx*coefIcoor),hj=(int) (hy*coefIcoor);
    28662867                        if(!quadtree->TooClose(&vi,seuil,hi,hj)){
    28672868                                // a good new point
    … …  
    30363037        }
    30373038        /*}}}*/
    3038         BamgVertex* Mesh::NearestVertex(Icoor1 i,Icoor1 j) {/*{{{*/
     3039        BamgVertex* Mesh::NearestVertex(int i,int j) {/*{{{*/
    30393040                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/NearestVertex)*/
    30403041                return  quadtree->NearestVertex(i,j);
    … …  
    36953696
    36963697                //Compute coefIcoor
    3697                 coefIcoor= (MaxICoor)/(Max(pmax.x-pmin.x,pmax.y-pmin.y));
     3698                int MaxICoord = 1073741823; //2^30 - 1 = =111...111 (29 times one)
     3699                coefIcoor= (MaxICoord)/(Max(pmax.x-pmin.x,pmax.y-pmin.y));
    36983700                if (coefIcoor<=0){
    36993701                        _error_("coefIcoor should be positive, a problem in the geometry is likely");
    … …  
    39253927        /*}}}*/
    39263928        I2 Mesh::R2ToI2(const R2 & P) const {/*{{{*/
    3927                 return  I2( (Icoor1) (coefIcoor*(P.x-pmin.x)),(Icoor1) (coefIcoor*(P.y-pmin.y)) );
     3929                return  I2( (int) (coefIcoor*(P.x-pmin.x)),(int) (coefIcoor*(P.y-pmin.y)) );
    39283930        }
    39293931        /*}}}*/

  • issm/trunk-jpl/src/c/bamg/Mesh.h

    r21623 r21629  
    102102                        long GetId(const Edge & t) const;
    103103                        long GetId(const Edge * t) const;
    104                         BamgVertex* NearestVertex(Icoor1 i,Icoor1 j) ;
     104                        BamgVertex* NearestVertex(int i,int j) ;
    105105                        Triangle* TriangleFindFromCoord(const I2 & ,Icoor2 [3],Triangle *tstart=0) const;
    106106                        void ReadMesh(int* index,double* x,double* y,int nods,int nels);

  • issm/trunk-jpl/src/c/bamg/macros.h

    r12934 r21629  
    1919        static const short NextVertex[3] = {1,2,0};
    2020        static const short PreviousVertex[3] = {2,0,1};
    21         const  Icoor1 MaxICoor   = 1073741823; // 2^30-1 =111...111 (29 times one)
    2221}
    2322
Note: See TracChangeset for help on using the changeset viewer.