Changeset 22252


Ignore:
Timestamp:
11/13/17 07:17:18 (7 years ago)
Author:
tsantos
Message:

CHG: setting bamgopts->verbose in some Mesh.cpp methods, fixing bad 'int' variable (number too big to be int, changed to long).

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

Legend:

Unmodified
Added
Removed

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

    r22251 r22252  
    3737
    3838                /*Fill holes and generate mesh properties*/
    39                 ReconstructExistingMesh();
    40         }
    41         /*}}}*/
    42         Mesh::Mesh(int* index,double* x,double* y,int nods,int nels):Gh(*(new Geometry())),BTh(*this){/*{{{*/
     39                ReconstructExistingMesh(bamgopts);
     40        }
     41        /*}}}*/
     42        Mesh::Mesh(int* index,double* x,double* y,int nods,int nels,BamgOpts* bamgopts):Gh(*(new Geometry())),BTh(*this){/*{{{*/
    4343
    4444                Init(0);
    45                 ReadMesh(index,x,y,nods,nels);
     45                ReadMesh(index,x,y,nods,nels,bamgopts);
    4646                SetIntCoor();
    47                 ReconstructExistingMesh();
    48         }
    49         /*}}}*/
    50         Mesh::Mesh(double* x,double* y,int nods):Gh(*(new Geometry())),BTh(*this){/*{{{*/
    51                 Triangulate(x,y,nods);
     47                ReconstructExistingMesh(bamgopts);
     48        }
     49        /*}}}*/
     50        Mesh::Mesh(double* x,double* y,int nods,BamgOpts* bamgopts):Gh(*(new Geometry())),BTh(*this){/*{{{*/
     51                Triangulate(x,y,nods,bamgopts);
    5252        }
    5353        /*}}}*/
    … …  
    143143                  Gh.PostRead();
    144144                  SetIntCoor();
    145                   ReconstructExistingMesh();
     145                  ReconstructExistingMesh(bamgopts);
    146146
    147147                  /*Final checks*/
    … …  
    247247
    248248        /*IO*/
    249         void Mesh::ReadMesh(int* index,double* x,double* y,int nods,int nels){/*{{{*/
     249        void Mesh::ReadMesh(int* index,double* x,double* y,int nods,int nels,BamgOpts* bamgopts){/*{{{*/
    250250
    251251                long i1,i2,i3;
    … …  
    258258                maxnbv=nbv;
    259259                nbt=nels;
     260                if(bamgopts) verbose=bamgopts->verbose;
    260261
    261262                //Vertices
    … …  
    310311        void Mesh::ReadMesh(BamgMesh* bamgmesh, BamgOpts* bamgopts){/*{{{*/
    311312
    312                 int    verbose;
     313                int    verbose=0;
    313314                double Hmin = HUGE_VAL;    // the infinie value
    314315                long   i1,i2,i3;
    … …  
    316317                Metric M1(1);
    317318
    318                 verbose=bamgopts->verbose;
     319                /*Check needed pointer*/
     320                _assert_(bamgmesh);
     321
     322                if(bamgopts) verbose=bamgopts->verbose;
    319323
    320324                nbv=bamgmesh->VerticesSize[0];
    … …  
    496500                int* connectivitysize_1=NULL;
    497501                int  connectivitymax_1=0;
     502                int verbose=0;
     503       
     504                /*Check needed pointer*/
     505                _assert_(bamgmesh);
    498506
    499507                /*Get options*/
    500                 int verbose=bamgopts->verbose;
     508                if(bamgopts) verbose=bamgopts->verbose;
    501509
    502510                /*Build reft that holds the number the subdomain number of each triangle, and the real numbering of the elements*/
    … …  
    898906                /*Intermediary*/
    899907                int  i,j;
    900 
    901                 if(bamgopts->verbose>3) _printf_("      processing metric\n");
     908                int verbose=0;
     909
     910                /*Check pointer*/
     911                _assert_(bamgopts);
     912
     913                /*Get options*/
     914                verbose=bamgopts->verbose;
     915
     916                if(verbose>3) _printf_("      processing metric\n");
    902917                double hmin = Max(bamgopts->hmin,MinimalHmin());
    903918                double hmax = Min(bamgopts->hmax,MaximalHmax());
    … …  
    933948        void Mesh::WriteMetric(BamgOpts* bamgopts) {/*{{{*/
    934949                int i;
     950                _assert_(bamgopts);
    935951                xDelete<double>(bamgopts->metric);
    936952                bamgopts->metric=xNew<double>(3*nbv);
    … …  
    9841000                //  Hessiantype = 0 =>  H is computed using double L2 projection
    9851001                //  Hessiantype = 1 =>  H is computed with green formula
     1002
     1003                /*Check pointer*/
     1004                _assert_(bamgopts);
    9861005
    9871006                /*Options*/
    … …  
    11341153
    11351154        }/*}}}*/
    1136         void Mesh::BoundAnisotropy(double anisomax,double hminaniso) {/*{{{*/
     1155        void Mesh::BoundAnisotropy(BamgOpts* bamgopts,double anisomax,double hminaniso) {/*{{{*/
    11371156                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/BoundAnisotropy)*/
    11381157
    1139                 long int verbose=0;
    1140                 double lminaniso = 1/ (Max(hminaniso*hminaniso,1e-100));
    1141 
     1158                int verbose=0;
     1159                double lminaniso = 1./ (Max(hminaniso*hminaniso,1e-100));
     1160       
     1161                //Get options
     1162                if(bamgopts) verbose=bamgopts->verbose;
     1163                       
    11421164                //display info
    11431165                if (verbose > 1)  _printf_("   BoundAnisotropy by " << anisomax << "\n");
    … …  
    11831205
    11841206                /*Recover options*/
    1185                 if(bamgopts){
    1186                         verbose=bamgopts->verbose;
    1187                 }
     1207                if(bamgopts) verbose=bamgopts->verbose;
    11881208
    11891209                //display info
    … …  
    15161536
    15171537                //Build Gh.coefIcoor
    1518                 int MaxICoord = 1073741823; //2^30 - 1 = =111...111 (29 times one)
     1538                long MaxICoord = 1073741823; //2^30 - 1 = =111...111 (29 times one)
    15191539                Gh.coefIcoor= (MaxICoord)/(Max(Gh.pmax.x-Gh.pmin.x,Gh.pmax.y-Gh.pmin.y));
    15201540                if (Gh.coefIcoor<=0){
    … …  
    16221642                double* s=NULL;
    16231643                long    nbsol;
    1624                 int     verbose;
     1644                int     verbose=0;
    16251645
    16261646                int   i,j,k,iA,iB,iC;
    16271647                int   iv;
     1648
     1649                /*Check pointer*/
     1650                _assert_(bamgopts);
    16281651
    16291652                /*Recover options*/
    … …  
    18301853                int verbose;
    18311854
     1855                /*Check pointer*/
     1856                _assert_(bamgopts);
     1857               
    18321858                /*Recover options*/
    18331859                verbose=bamgopts->verbose;
    … …  
    21282154                int i1,i2;
    21292155                int j1,j2;
     2156                int verbose=0;
    21302157
    21312158                /*Options*/
    2132                 int verbose=bamgopts->verbose;
     2159                if(bamgopts) verbose=bamgopts->verbose;
    21332160
    21342161                //  computed the number of cracked edge
    … …  
    22902317        }
    22912318        /*}}}*/
    2292         void Mesh::ForceBoundary() {/*{{{*/
     2319        void Mesh::ForceBoundary(BamgOpts* bamgopts) {/*{{{*/
    22932320                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/ForceBoundary)*/
    22942321
    2295                 long int verbose=2;
     2322                int verbose=0;
    22962323                int k=0;
    22972324                int nbfe=0,nbswp=0,Nbswap=0;
     2325
     2326                /*Get options*/
     2327                if(bamgopts) verbose=bamgopts->verbose;
    22982328
    22992329                //display
    … …  
    23322362        }
    23332363        /*}}}*/
    2334         void Mesh::FindSubDomain(int OutSide) {/*{{{*/
     2364        void Mesh::FindSubDomain(BamgOpts* bamgopts,int OutSide) {/*{{{*/
    23352365                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/FindSubDomain)*/
    23362366
    23372367                int verbose=0;
     2368
     2369                /*Get options*/
     2370                if(bamgopts) verbose=bamgopts->verbose;
    23382371
    23392372                if (verbose >2){
    … …  
    26452678        }
    26462679        /*}}}*/
    2647         void Mesh::Insert(){/*{{{*/
     2680        void Mesh::Insert(BamgOpts* bamgopts){/*{{{*/
    26482681                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/Insert)*/
    26492682
    … …  
    26522685                //Intermediary
    26532686                int i;
     2687                int verbose=0;
    26542688
    26552689                /*Get options*/
    2656                 long int verbose=2;
     2690                if(bamgopts) verbose=bamgopts->verbose;
    26572691
    26582692                //Display info
    … …  
    27762810        }
    27772811        /*}}}*/
    2778         long Mesh::InsertNewPoints(long nbvold,long & NbTSwap) {/*{{{*/
     2812        long Mesh::InsertNewPoints(long nbvold,long & NbTSwap,BamgOpts* bamgopts) {/*{{{*/
    27792813                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/InsertNewPoints)*/
    27802814
    2781                 long int verbose=0;
    2782                 double seuil= 1.414/2 ;// for two close point
     2815                int verbose=0;
     2816                double seuil= 1.414/2.;// for two close point
    27832817                long i;
    27842818                long NbSwap=0;
    27852819                long long det3[3];
     2820
     2821                /*Get options*/
     2822                if(bamgopts) verbose=bamgopts->verbose;
    27862823
    27872824                //number of new points
    … …  
    29112948        }
    29122949        /*}}}*/
    2913         void  Mesh::MaxSubDivision(double maxsubdiv) {/*{{{*/
     2950        void  Mesh::MaxSubDivision(BamgOpts* bamgopts,double maxsubdiv) {/*{{{*/
    29142951                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/MaxSubDivision)*/
    2915 
    2916                 long int verbose=0;
     2952               
     2953                /*Intermediaries*/
     2954                int verbose=0;
     2955
     2956                /*Get options*/
     2957                if(bamgopts) verbose=bamgopts->verbose;
    29172958
    29182959                const  double maxsubdiv2 = maxsubdiv*maxsubdiv;
    … …  
    29983039
    29993040                int i,j,k;
     3041                int verbose=0;
    30003042                long NbTSwap=0;
    30013043                long nbtold=nbt;
    … …  
    30073049
    30083050                /*Recover options*/
    3009                 int verbose=bamgopts->verbose;
     3051                if(bamgopts) verbose=bamgopts->verbose;
    30103052
    30113053                /*First, insert old points if requested*/
    … …  
    30403082                        //if(pointsoutside) _printf_("WARNING: One or more points of the initial mesh fall outside of the geometric boundary\n");
    30413083                        Bh.CreateSingleVertexToTriangleConnectivity();     
    3042                         InsertNewPoints(nbvold,NbTSwap);
     3084                        InsertNewPoints(nbvold,NbTSwap,bamgopts);
    30433085                }
    30443086                else Bh.CreateSingleVertexToTriangleConnectivity();     
    … …  
    30983140                        }// for triangle   
    30993141
    3100                         if (!InsertNewPoints(nbvold,NbTSwap)) break;
     3142                        if (!InsertNewPoints(nbvold,NbTSwap,bamgopts)) break;
    31013143                        for (i=nbtold;i<nbt;i++) first_np_or_next_t[i]=iter;
    31023144                        Headt = nbt; // empty list
    … …  
    32553297        }                 
    32563298        /*}}}*/
    3257         void Mesh::ReconstructExistingMesh(){/*{{{*/
     3299        void Mesh::ReconstructExistingMesh(BamgOpts* bamgopts){/*{{{*/
    32583300                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/FillHoleInMesh)*/
    32593301
    … …  
    32683310                /*Intermediary*/
    32693311                int verbose=0;
     3312
     3313                /*Get options*/
     3314                if(bamgopts) verbose=bamgopts->verbose;
    32703315
    32713316                // generation of the integer coordinate
    … …  
    34513496                }
    34523497
    3453                 FindSubDomain(1);
     3498                FindSubDomain(bamgopts,1);
    34543499                // remove all the hole
    34553500                // remove all the good sub domain
    … …  
    35253570                        _error_("number of triangles edges alone = " << k);
    35263571                }
    3527                 FindSubDomain();
     3572                FindSubDomain(bamgopts);
    35283573
    35293574                delete edge4;
    … …  
    36483693
    36493694                //Compute coefIcoor
    3650                 int MaxICoord = 1073741823; //2^30 - 1 = =111...111 (29 times one)
     3695                long MaxICoord = 1073741823; //2^30 - 1 = =111...111 (29 times one)
    36513696                coefIcoor= (MaxICoord)/(Max(pmax.x-pmin.x,pmax.y-pmin.y));
    36523697                if (coefIcoor<=0){
    … …  
    36903735        }
    36913736        /*}}}*/
    3692         void Mesh::SmoothingVertex(int nbiter,double omega ) { /*{{{*/
     3737        void Mesh::SmoothingVertex(BamgOpts* bamgopts,int nbiter,double omega ) { /*{{{*/
    36933738                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/SmoothingVertex)*/
    3694 
    3695                 long int verbose=0;
     3739               
     3740                /*Intermediaries*/
     3741                int verbose=0;
     3742       
     3743                /*Get options*/
     3744                if(bamgopts) verbose=bamgopts->verbose;
     3745
    36963746                //  if quatree exist remove it end reconstruct
    36973747                if (quadtree) delete quadtree;
    … …  
    37303780        }
    37313781        /*}}}*/
    3732         void Mesh::SmoothMetric(double raisonmax) { /*{{{*/
     3782        void Mesh::SmoothMetric(BamgOpts* bamgopts,double raisonmax) { /*{{{*/
    37333783                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/SmoothMetric)*/
    3734 
    3735                 long int verbose=0;
     3784               
     3785                /*Intermediaries*/
     3786                int verbose=0;
     3787
     3788                /*Get options*/
     3789                if(bamgopts) verbose=bamgopts->verbose;
    37363790
    37373791                if(raisonmax<1.1) return;
    … …  
    40254079        }
    40264080        /*}}}*/
    4027         void Mesh::Triangulate(double* x,double* y,int nods){/*{{{*/
     4081        void Mesh::Triangulate(double* x,double* y,int nods,BamgOpts* bamgopts){/*{{{*/
    40284082
    40294083                int verbose=0;
    40304084                int i;
    40314085                Metric M1(1);
     4086       
     4087                /*Get options*/
     4088                if(bamgopts) verbose=bamgopts->verbose;
    40324089
    40334090                /*Initialize mesh*/
    … …  
    40474104
    40484105                /*Insert Vertices*/
    4049                 Insert();
     4106                Insert(bamgopts);
    40504107        }
    40514108        /*}}}*/
    … …  
    40564113                /*Intermediaries*/
    40574114                int                i,k;
     4115                int                                      verbose=0;
    40584116                int                nbcurves    = 0;
    40594117                int                NbNewPoints,NbEdgeCurve;
    … …  
    40704128
    40714129                /*Get options*/
    4072                 int verbose=bamgopts->verbose;
     4130                if(bamgopts) verbose=bamgopts->verbose;
    40734131
    40744132                //build background mesh flag (1 if background, else 0)
    … …  
    43484406                if (verbose>3) _printf_("      Creating initial Constrained Delaunay Triangulation...\n");
    43494407                if (verbose>3) _printf_("         Inserting boundary points\n");
    4350                 Insert();
     4408                Insert(bamgopts);
    43514409
    43524410                //Force the boundary
    43534411                if (verbose>3) _printf_("         Forcing boundaries\n");
    4354                 ForceBoundary();
     4412                ForceBoundary(bamgopts);
    43554413
    43564414                //Extract SubDomains
    43574415                if (verbose>3) _printf_("         Extracting subdomains\n");
    4358                 FindSubDomain();
     4416                FindSubDomain(bamgopts);
    43594417
    43604418                if (verbose>3) _printf_("      Inserting internal points\n");
    … …  
    43664424                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/GeomToTriangles1)*/
    43674425
     4426                /*Intermediaries*/
     4427                int verbose=0;
     4428
    43684429                /*Get options*/
    4369                 int verbose=bamgopts->verbose;
     4430                if(bamgopts) verbose=bamgopts->verbose;
    43704431
    43714432                Gh.NbRef++;// add a ref to Gh
    … …  
    46754736                if (verbose>3) _printf_("      Creating initial Constrained Delaunay Triangulation...\n");
    46764737                if (verbose>3) _printf_("         Inserting boundary points\n");
    4677                 Insert();
     4738                Insert(bamgopts);
    46784739
    46794740                //Force the boundary
    46804741                if (verbose>3) _printf_("         Forcing boundaries\n");
    4681                 ForceBoundary();
     4742                ForceBoundary(bamgopts);
    46824743
    46834744                //Extract SubDomains
    46844745                if (verbose>3) _printf_("         Extracting subdomains\n");
    4685                 FindSubDomain();
     4746                FindSubDomain(bamgopts);
    46864747
    46874748                if (verbose>3) _printf_("      Inserting internal points\n");

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

    r22251 r22252  
    5757                        //Constructors/Destructors
    5858                        Mesh(BamgGeom* bamggeom,BamgMesh* bamgmesh,BamgOpts* bamgopts);
    59                         Mesh(int* index,double* x,double* y,int nods,int nels);/*MeshConvert*/
    60                         Mesh(double* x,double* y,int nods); /*BamgTriangulate*/
     59                        Mesh(int* index,double* x,double* y,int nods,int nels,BamgOpts* bamgopts);/*MeshConvert*/
     60                        Mesh(double* x,double* y,int nods,BamgOpts* bamgopts); /*BamgTriangulate*/
    6161                        Mesh(Mesh &,Geometry * pGh=0,Mesh* pBTh=0,long maxnbv_in=0 ); //copy operator
    6262                        Mesh(const Mesh &,const int *flag,const int *bb,BamgOpts* bamgopts); // truncature
    … …  
    7878                        R2 I2ToR2(const I2 & P) const;
    7979                        void AddVertex(BamgVertex & s,Triangle * t,long long *  =0) ;
    80                         void Insert();
     80                        void Insert(BamgOpts* bamgopts);
    8181                        void Echo(void);
    82                         void ForceBoundary();
    83                         void FindSubDomain(int OutSide=0);
     82                        void ForceBoundary(BamgOpts* bamgopts);
     83                        void FindSubDomain(BamgOpts* bamgopts,int OutSide=0);
    8484                        long TriangleReferenceList(long*) const;
    8585                        void TriangleIntNumbering(long* renumbering);
    8686                        void CrackMesh(BamgOpts* bamgopts);
    87                         void SmoothMetric(double raisonmax) ;
    88                         void BoundAnisotropy(double anisomax,double hminaniso= 1e-100) ;
     87                        void SmoothMetric(BamgOpts* bamgopts,double raisonmax) ;
     88                        void BoundAnisotropy(BamgOpts* bamgopts,double anisomax,double hminaniso= 1e-100) ;
    8989                        Edge** MakeGeomEdgeToEdge();
    9090                        long SplitInternalEdgeWithBorderVertices();
    9191                        void MakeBamgQuadtree();
    92                         void MaxSubDivision(double maxsubdiv);
     92                        void MaxSubDivision(BamgOpts* bamgopts,double maxsubdiv);
    9393                        void NewPoints(Mesh &,BamgOpts* bamgopts,int KeepVertices=1);
    94                         long InsertNewPoints(long nbvold,long & NbTSwap);
     94                        long InsertNewPoints(long nbvold,long & NbTSwap,BamgOpts* bamgopts);
    9595                        void TrianglesRenumberBySubDomain(bool justcompress=false);
    96                         void SmoothingVertex(int =3,double=0.3);
     96                        void SmoothingVertex(BamgOpts* bamgopts,int =3,double=0.3);
    9797                        Metric MetricAt (const R2 &);
    9898                        GeomEdge* ProjectOnCurve( Edge & AB, BamgVertex &  A, BamgVertex & B,double theta, BamgVertex & R,VertexOnEdge & BR,VertexOnGeom & GR);
    … …  
    105105                        BamgVertex* NearestVertex(int i,int j) ;
    106106                        Triangle* TriangleFindFromCoord(const I2 & ,long long [3],Triangle *tstart=0);
    107                         void ReadMesh(int* index,double* x,double* y,int nods,int nels);
     107                        void ReadMesh(int* index,double* x,double* y,int nods,int nels,BamgOpts* bamgopts);
    108108                        void ReadMesh(BamgMesh* bamgmesh, BamgOpts* bamgopts);
    109109                        void WriteMesh(BamgMesh* bamgmesh,BamgOpts* bamgopts);
    … …  
    115115                        void BuildMetric1(BamgOpts* bamgopts);
    116116                        void BuildGeometryFromMesh(BamgOpts* bamgopts=NULL);
    117                         long  RandomNumber(long max);
    118                         void ReconstructExistingMesh();
     117                        long RandomNumber(long max);
     118                        void ReconstructExistingMesh(BamgOpts* bamgopts);
    119119
    120120                        //Inline methods
    … …  
    142142                        void TriangulateFromGeom1(BamgOpts* bamgopts,int KeepVertices=1);// the real constructor mesh adaption
    143143                        void TriangulateFromGeom0(BamgOpts* bamgopts);// the real constructor mesh generator
    144                         void Triangulate(double* x,double* y,int nods);
     144                        void Triangulate(double* x,double* y,int nods,BamgOpts* bamgopts);
    145145                        void Init(long);
    146146                        int ForceEdge(BamgVertex &a, BamgVertex & b,AdjacentTriangle & taret) ;

  • issm/trunk-jpl/src/c/classes/AmrBamg.cpp

    r22241 r22252  
    8383
    8484        /*Read father mesh and create geometry*/
    85         Mesh* Th=new Mesh(elements,x,y,numberofvertices,numberofelements);
     85        Mesh* Th=new Mesh(elements,x,y,numberofvertices,numberofelements,this->options);
    8686
    8787        /*Write geometry*/

  • issm/trunk-jpl/src/c/modules/BamgConvertMeshx/BamgConvertMeshx.cpp

    r15061 r22252  
    1717
    1818        /*read mesh*/
    19         Mesh Th(index,x,y,nods,nels);
     19        Mesh Th(index,x,y,nods,nels,bamgopts);
    2020
    2121        /*write mesh and geometry*/

  • issm/trunk-jpl/src/c/modules/BamgTriangulatex/BamgTriangulatex.cpp

    r15065 r22252  
    1414int BamgTriangulatex(int** pindex,int* pnels,double* x,double* y,int nods){
    1515
    16         Mesh Th(x,y,nods);
     16        BamgOpts* bamgopts=new BamgOpts();//use bamgopts->verbose>5 to debug bamg::Mesh()
     17        Mesh Th(x,y,nods,bamgopts);
    1718        Th.WriteIndex(pindex,pnels);
     19        delete bamgopts;
    1820        //delete &Th;
    1921        return 0;
    20 
    2122}

  • issm/trunk-jpl/src/c/modules/Bamgx/Bamgx.cpp

    r21864 r22252  
    150150
    151151                //Smoothe metric
    152                 BTh.SmoothMetric(bamgopts->gradation);
     152                BTh.SmoothMetric(bamgopts,bamgopts->gradation);
    153153
    154154                //Control element subdivision
    155                 BTh.MaxSubDivision(bamgopts->maxsubdiv);
     155                BTh.MaxSubDivision(bamgopts,bamgopts->maxsubdiv);
    156156
    157157                //Bound anisotropy
    158                 BTh.BoundAnisotropy(bamgopts->anisomax,hminaniso);
     158                BTh.BoundAnisotropy(bamgopts,bamgopts->anisomax,hminaniso);
    159159
    160160                //Build new mesh
    … …  
    171171
    172172                //Smooth vertices
    173                 if(nbsmooth>0) Th.SmoothingVertex(nbsmooth,bamgopts->omega);
     173                if(nbsmooth>0) Th.SmoothingVertex(bamgopts,nbsmooth,bamgopts->omega);
    174174
    175175                //display info

  • issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.cpp

    r21856 r22252  
    3030        double aa,bb;
    3131        long long dete[3];
     32        BamgOpts* bamgopts=new BamgOpts();//if necessary debug bamg::Mesh, set bamgopts->verbose>5
    3233
    3334        /*Checks*/
    … …  
    4950
    5051        /*read background mesh*/
    51         Mesh* Th=new Mesh(index_data,x_data,y_data,nods_data,nels_data);
     52        bamgopts->verbose=0;
     53        Mesh* Th=new Mesh(index_data,x_data,y_data,nods_data,nels_data,bamgopts);
    5254
    5355        /*Get reference number (for subdomains)*/
    … …  
    166168        /*clean-up and return*/
    167169        delete Th;
     170        delete bamgopts;
    168171        xDelete<long>(reft);
    169172        xDelete<int>(connectivity);
Note: See TracChangeset for help on using the changeset viewer.