Changeset 5091


Ignore:
Timestamp:
08/09/10 14:25:16 (15 years ago)
Author:
Mathieu Morlighem
Message:

renamed some fields

Location:
issm/trunk/src/c
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/io/WriteData.cpp

    r4779 r5091  
    167167        fnames[++i] = "Vertices";                 fsize[i][0]=bm->VerticesSize[0];                  fsize[i][1]=bm->VerticesSize[1];                 fpointer[i]=&bm->Vertices;
    168168        fnames[++i] = "Edges";                    fsize[i][0]=bm->EdgesSize[0];                     fsize[i][1]=bm->EdgesSize[1];                    fpointer[i]=&bm->Edges;
    169         fnames[++i] = "Segments";                 fsize[i][0]=bm->SegmentsSize[0];                  fsize[i][1]=bm->SegmentsSize[1];                 fpointer[i]=&bm->Segments;
    170         fnames[++i] = "ElementEdges";             fsize[i][0]=bm->ElementEdgesSize[0];              fsize[i][1]=bm->ElementEdgesSize[1];             fpointer[i]=&bm->ElementEdges;
     169        fnames[++i] = "IssmSegments";             fsize[i][0]=bm->IssmSegmentsSize[0];              fsize[i][1]=bm->IssmSegmentsSize[1];             fpointer[i]=&bm->IssmSegments;
     170        fnames[++i] = "IssmEdges";                fsize[i][0]=bm->IssmEdgesSize[0];                 fsize[i][1]=bm->IssmEdgesSize[1];                fpointer[i]=&bm->IssmEdges;
    171171        fnames[++i] = "Quadrilaterals";           fsize[i][0]=bm->QuadrilateralsSize[0];            fsize[i][1]=bm->QuadrilateralsSize[1];           fpointer[i]=&bm->Quadrilaterals;
    172172        fnames[++i] = "VerticesOnGeometricVertex";fsize[i][0]=bm->VerticesOnGeometricVertexSize[0]; fsize[i][1]=bm->VerticesOnGeometricVertexSize[1];fpointer[i]=&bm->VerticesOnGeometricVertex;

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

    r3912 r5091  
    1717        double coef;
    1818        int    verbosity;
    19         int    NbSmooth;
     19        int    nbsmooth;
    2020
    2121        /*intermediary*/
    … …  
    2424        double costheta=2;
    2525        double hminaniso=1e-100;
    26 
    2726        Triangles* Thr=NULL;
    2827        Triangles* Thb=NULL;
    2928
    3029        /*Bamg options*/
    31         NbSmooth=bamgopts->NbSmooth;
    32         coef=bamgopts->coef;
    33         maxnbv=bamgopts->maxnbv;
     30        nbsmooth =bamgopts->nbsmooth;
     31        coef     =bamgopts->coeff;
     32        maxnbv   =bamgopts->maxnbv;
    3433        verbosity=bamgopts->verbose;
    3534
    3635        // no metric -> no smoothing
    37         if (bamgopts->metric==NULL){
    38                 NbSmooth=0;
    39         }
     36        if (bamgopts->metric==NULL) nbsmooth=0;
    4037
    4138        /*If no mesh in input, generate one*/
    42 
    4339        if(bamgmesh_in->TrianglesSize[0]==0){
    4440                /*Mesh generation {{{1*/
    … …  
    8783                /*}}}*/
    8884        }
    89 
    9085        else{
    9186                /*Anisotropic mesh adaptation {{{1*/
    … …  
    185180
    186181                //Smooth vertices
    187                 if(NbSmooth>0) Th.SmoothingVertex(NbSmooth,bamgopts->omega);
     182                if(nbsmooth>0) Th.SmoothingVertex(nbsmooth,bamgopts->omega);
    188183
    189184                //display info

  • issm/trunk/src/c/objects/Bamg/BamgGeom.cpp

    r3683 r5091  
    88        bamggeom->hVertices=NULL;
    99        bamggeom->MetricVertices=NULL;
    10         bamggeom->h1h2VpVertices=NULL;
    1110        bamggeom->TangentAtEdgesSize[0]=0,    bamggeom->TangentAtEdgesSize[1]=0;    bamggeom->TangentAtEdges=NULL;
    1211        bamggeom->CornersSize[0]=0,           bamggeom->CornersSize[1]=0;           bamggeom->Corners=NULL;

  • issm/trunk/src/c/objects/Bamg/BamgGeom.h

    r3683 r5091  
    1515        double* hVertices;
    1616        double* MetricVertices;
    17         double* h1h2VpVertices;
    1817
    1918        int     TangentAtEdgesSize[2];

  • issm/trunk/src/c/objects/Bamg/BamgMesh.cpp

    r3683 r5091  
    44void BamgMeshInit(BamgMesh* bamgmesh){
    55
    6         bamgmesh->TrianglesSize[0]=0,     bamgmesh->TrianglesSize[1]=0;      bamgmesh->Triangles=NULL;
    76        bamgmesh->VerticesSize[0]=0,      bamgmesh->VerticesSize[1]=0;       bamgmesh->Vertices=NULL;
    87        bamgmesh->EdgesSize[0]=0,         bamgmesh->EdgesSize[1]=0;          bamgmesh->Edges=NULL;
    9         bamgmesh->ElementEdgesSize[0]=0,  bamgmesh->ElementEdgesSize[1]=0;   bamgmesh->ElementEdges=NULL;
    10         bamgmesh->SegmentsSize[0]=0,      bamgmesh->SegmentsSize[1]=0;       bamgmesh->Segments=NULL;
     8        bamgmesh->TrianglesSize[0]=0,     bamgmesh->TrianglesSize[1]=0;      bamgmesh->Triangles=NULL;
    119        bamgmesh->QuadrilateralsSize[0]=0,bamgmesh->QuadrilateralsSize[1]=0; bamgmesh->Quadrilaterals=NULL;
    1210        bamgmesh->VerticesOnGeometricVertexSize[0]=0, bamgmesh->VerticesOnGeometricVertexSize[1]=0;bamgmesh->VerticesOnGeometricVertex=NULL;
    … …  
    1614        bamgmesh->SubDomainsFromGeomSize[0]=0, bamgmesh->SubDomainsFromGeomSize[1]=0;  bamgmesh->SubDomainsFromGeom=NULL;
    1715        bamgmesh->hVertices=NULL;
     16        bamgmesh->IssmEdgesSize[0]=0,  bamgmesh->IssmEdgesSize[1]=0;   bamgmesh->IssmEdges=NULL;
     17        bamgmesh->IssmSegmentsSize[0]=0,             bamgmesh->IssmSegmentsSize[1]=0;             bamgmesh->IssmSegments=NULL;
    1818        bamgmesh->ElementConnectivitySize[0]=0,      bamgmesh->ElementConnectivitySize[1]=0;      bamgmesh->ElementConnectivity=NULL;
    1919        bamgmesh->NodalConnectivitySize[0]=0,        bamgmesh->NodalConnectivitySize[1]=0;        bamgmesh->NodalConnectivity=NULL;

  • issm/trunk/src/c/objects/Bamg/BamgMesh.h

    r3683 r5091  
    77struct BamgMesh{
    88
    9         int     TrianglesSize[2];
    10         double* Triangles;
    11 
     9        /*Bamg input/output*/
    1210        int     VerticesSize[2];
    1311        double* Vertices;
    … …  
    1614        double* Edges;
    1715
    18         int     ElementEdgesSize[2];
    19         double* ElementEdges;
    20 
    21         int     SegmentsSize[2];
    22         double* Segments;
     16        int     TrianglesSize[2];
     17        double* Triangles;
    2318
    2419        int     QuadrilateralsSize[2];
    … …  
    4035        double* SubDomainsFromGeom;
    4136
     37        int     CrackedVerticesSize[2];
     38        double* CrackedVertices;
     39
     40        int     CrackedEdgesSize[2];
     41        double* CrackedEdges;
     42
    4243        double* hVertices;
     44
     45        /*Output for ISSM*/
     46        int     IssmEdgesSize[2];
     47        double* IssmEdges;
     48
     49        int     IssmSegmentsSize[2];
     50        double* IssmSegments;
    4351
    4452        int     ElementConnectivitySize[2];
    … …  
    5159        double* NodalElementConnectivity;
    5260
    53         int     CrackedVerticesSize[2];
    54         double* CrackedVertices;
    55 
    56         int     CrackedEdgesSize[2];
    57         double* CrackedEdges;
    5861};
    5962

  • issm/trunk/src/c/objects/Bamg/BamgOpts.cpp

    r3775 r5091  
    88        bamgopts->iso=0;
    99        bamgopts->maxnbv=0;
    10         bamgopts->MaximalAngleOfCorner=0;
     10        bamgopts->MaxCornerAngle=0;
    1111        bamgopts->Hessiantype=0;
    1212        bamgopts->Metrictype=0;
    … …  
    1616        bamgopts->power=0;
    1717        bamgopts->anisomax=0;
    18         bamgopts->NbSmooth=0;
     18        bamgopts->nbsmooth=0;
    1919        bamgopts->nbjacobi=0;
    2020        bamgopts->omega=0;
    … …  
    3030        bamgopts->err=NULL;
    3131        bamgopts->errg=0;
    32         bamgopts->coef=0;
     32        bamgopts->coeff=0;
    3333        bamgopts->metric=NULL;
    3434        bamgopts->field=NULL;
    … …  
    4141        int i;
    4242
    43         if (bamgopts->coef==0) ISSMERROR("'coef' should be positive");
     43        if (bamgopts->coeff==0) ISSMERROR("'coeff' should be positive");
    4444        if (bamgopts->maxsubdiv<=1) ISSMERROR("'maxsubdiv' should be >1");
    4545        if (bamgopts->Crack!=0  && bamgopts->Crack!=1) ISSMERROR("'Crack' supported options are 0 and 1");
    … …  
    5050        if (bamgopts->nbjacobi<=0) ISSMERROR("'nbjacobi' option should be >0");
    5151        if (bamgopts->geometricalmetric!=0  && bamgopts->geometricalmetric!=1) ISSMERROR("'geometricalmetric' supported options are 0 and 1");
    52         if (bamgopts->NbSmooth<=0) ISSMERROR("'NbSmooth' option should be >0");
     52        if (bamgopts->nbsmooth<=0) ISSMERROR("'nbsmooth' option should be >0");
    5353        if (bamgopts->maxnbv<3) ISSMERROR("'maxnbv' option should be >3");
    5454        if (bamgopts->hmin<=0) ISSMERROR("'hmin' option should be >0");

  • issm/trunk/src/c/objects/Bamg/BamgOpts.h

    r3683 r5091  
    1010        int     iso;
    1111        int     maxnbv;
    12         double  MaximalAngleOfCorner;
     12        double  MaxCornerAngle;
    1313        int     Crack;
    1414        int     Hessiantype;
    … …  
    1818        double  power;
    1919        double  anisomax;
    20         int     NbSmooth;
     20        int     nbsmooth;
    2121        int     nbjacobi;
    2222        double  omega;
    … …  
    3232        double* err;
    3333        double  errg;
    34         double  coef;
     34        double  coeff;
    3535        double* metric;
    3636        double* field;

  • issm/trunk/src/c/objects/Bamg/Geometry.cpp

    r5016 r5091  
    4343                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, MeshGeom.cpp/~Geometry)*/
    4444
    45                 long int verbosity=0;
     45                long int verbose=0;
    4646
    4747                if (NbRef>0){
    … …  
    203203                }
    204204
    205                 //h1h2VpVertices
    206                 if(bamggeom->h1h2VpVertices){
    207                         if(verbose>5) printf("      processing h1h2VpVertices\n");
    208                         double h1,h2,v1,v2;
    209                         for (i=0;i< nbv;i++) {
    210                                 h1=(double)bamggeom->MetricVertices[i*4+0];
    211                                 h2=(double)bamggeom->MetricVertices[i*4+1];
    212                                 v1=(double)bamggeom->MetricVertices[i*4+2];
    213                                 v2=(double)bamggeom->MetricVertices[i*4+3];
    214                                 vertices[i].m = Metric(MatVVP2x2(1/(h1*h1),1/(h2*h2),D2(v1,v2)));
    215                         }
    216                 }
    217 
    218                 //MaximalAngleOfCorner
    219                 if (bamgopts->MaximalAngleOfCorner){
    220                         if(verbose>5) printf("      processing MaximalAngleOfCorner\n");
    221                         MaximalAngleOfCorner=bamgopts->MaximalAngleOfCorner*Pi/180;
     205                //MaxCornerAngle
     206                if (bamgopts->MaxCornerAngle){
     207                        if(verbose>5) printf("      processing MaxCornerAngle\n");
     208                        MaxCornerAngle=bamgopts->MaxCornerAngle*Pi/180;
    222209                }
    223210
    … …  
    431418                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, MeshGeom.cpp/AfterRead)*/
    432419
    433                 long int verbosity=0;
     420                long int verbose=0;
    434421
    435422                long i,j,k;
    … …  
    621608                                float angle2= !j2 ? OppositeAngle(eangle[i2]) : eangle[i2];
    622609                                float da12 = Abs(angle2-angle1);
    623                                 if (( da12 >= MaximalAngleOfCorner ) && (da12 <= 2*Pi -MaximalAngleOfCorner)) {
     610                                if (( da12 >= MaxCornerAngle ) && (da12 <= 2*Pi -MaxCornerAngle)) {
    624611                                        vertices[i].SetCorner() ;
    625612                                }
    … …  
    805792                printf("   pmax (x,y): (%g %g)\n",pmax.x,pmax.y);
    806793                printf("   coefIcoor: %g\n",coefIcoor);
    807                 printf("   MaximalAngleOfCorner: %g\n",MaximalAngleOfCorner);
     794                printf("   MaxCornerAngle: %g\n",MaxCornerAngle);
    808795
    809796                return;
    … …  
    829816                //  BeginOfCurve=0;
    830817                subdomains=0;
    831                 MaximalAngleOfCorner = 10*Pi/180;
     818                MaxCornerAngle = 10*Pi/180;
    832819        }
    833820        /*}}}1*/

  • issm/trunk/src/c/objects/Bamg/Geometry.h

    r3913 r5091  
    3232                        R2 pmin,pmax; // extrema
    3333                        double coefIcoor;  // coef to integer Icoor1;
    34                         double MaximalAngleOfCorner;
     34                        double MaxCornerAngle;
    3535
    3636                        //Constructor/Destructors

  • issm/trunk/src/c/objects/Bamg/ListofIntersectionTriangles.cpp

    r3913 r5091  
    1515                Triangle *tbegin, *t;
    1616
    17                 long int verbosity=2;
     17                long int verbose=2;
    1818                Icoor2 deta[3], deti,detj;
    1919                double ba[3];
    … …  
    7878                                if (t == tbegin) { //
    7979                                        double ba,bb;
    80                                         long int verbosity=2;
     80                                        long int verbose=2;
    8181                                        TriangleAdjacent edge=CloseBoundaryEdge(a,t,ba,bb);
    8282                                        MeshVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);

  • issm/trunk/src/c/objects/Bamg/MeshVertex.cpp

    r3913 r5091  
    123123                double hmin=(bamgopts->hmin);
    124124                double hmax=(bamgopts->hmax);
    125                 double coef=bamgopts->coef;
     125                double coef=bamgopts->coeff;
    126126                int    Metrictype=(bamgopts->Metrictype);
    127127

  • issm/trunk/src/c/objects/Bamg/Triangles.cpp

    r5018 r5091  
    225225                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/Triangles)*/
    226226
    227                 long int verbosity=2;
     227                long int verbose=2;
    228228                //if(vertices)  delete [] vertices; //TEST  crash if not commented
    229229                if(edges)     delete [] edges;
    … …  
    584584                bamgmesh->EdgesSize[0]=nbe;
    585585                bamgmesh->EdgesSize[1]=3;
    586                 int NumSegments=0;
     586                int NumIssmSegments=0;
    587587                if (nbe){
    588588                        bamgmesh->Edges=(double*)xmalloc(3*nbe*sizeof(double));
    … …  
    592592                                bamgmesh->Edges[i*3+2]=edges[i].ref;
    593593                                if(edges[i].onGeometry){
    594                                         NumSegments++;
     594                                        NumIssmSegments++;
    595595                                }
    596596                        }
    … …  
    624624                        }
    625625                }
    626                 bamgmesh->ElementEdgesSize[0]=edge4->nb();
    627                 bamgmesh->ElementEdgesSize[1]=4;
    628                 bamgmesh->ElementEdges=(double*)xmalloc(4*edge4->nb()*sizeof(double));
     626                bamgmesh->IssmEdgesSize[0]=edge4->nb();
     627                bamgmesh->IssmEdgesSize[1]=4;
     628                bamgmesh->IssmEdges=(double*)xmalloc(4*edge4->nb()*sizeof(double));
    629629                for (i=0;i<edge4->nb();i++){
    630630                        /*Invert first two vertices if necessary*/
    … …  
    634634                                        if (triangles[(int)elemedge[2*i+0]]((j+1)%3)==vertices+edge4->j(i)){
    635635                                                //trigonometric direction
    636                                                 bamgmesh->ElementEdges[i*4+0]=edge4->i(i)+1;// back to M indexing
    637                                                 bamgmesh->ElementEdges[i*4+1]=edge4->j(i)+1;// back to M indexing
     636                                                bamgmesh->IssmEdges[i*4+0]=edge4->i(i)+1;// back to M indexing
     637                                                bamgmesh->IssmEdges[i*4+1]=edge4->j(i)+1;// back to M indexing
    638638                                        }
    639639                                        else{
    640                                                 bamgmesh->ElementEdges[i*4+0]=edge4->j(i)+1;// back to M indexing
    641                                                 bamgmesh->ElementEdges[i*4+1]=edge4->i(i)+1;// back to M indexing
     640                                                bamgmesh->IssmEdges[i*4+0]=edge4->j(i)+1;// back to M indexing
     641                                                bamgmesh->IssmEdges[i*4+1]=edge4->i(i)+1;// back to M indexing
    642642                                        }
    643643                                        found=true;
    … …  
    646646                        }
    647647                        ISSMASSERT(found);
    648                         bamgmesh->ElementEdges[i*4+2]=elemedge[2*i+0]+1; // back to M indexing
    649                         bamgmesh->ElementEdges[i*4+3]=elemedge[2*i+1]+1; // back to M indexing
     648                        bamgmesh->IssmEdges[i*4+2]=elemedge[2*i+0]+1; // back to M indexing
     649                        bamgmesh->IssmEdges[i*4+3]=elemedge[2*i+1]+1; // back to M indexing
    650650                }
    651651                //clean up
    … …  
    653653                xfree((void**)&elemedge);
    654654
    655                 /*Segments*/
    656                 if(verbose>5) printf("      writing Segments\n");
    657                 bamgmesh->SegmentsSize[0]=NumSegments;
    658                 bamgmesh->SegmentsSize[1]=4;
    659                 bamgmesh->Segments=(double*)xmalloc(4*NumSegments*sizeof(double));
     655                /*IssmSegments*/
     656                if(verbose>5) printf("      writing IssmSegments\n");
     657                bamgmesh->IssmSegmentsSize[0]=NumIssmSegments;
     658                bamgmesh->IssmSegmentsSize[1]=4;
     659                bamgmesh->IssmSegments=(double*)xmalloc(4*NumIssmSegments*sizeof(double));
    660660                num=0;
    661661                for (i=0;i<nbe;i++){
    … …  
    669669                                                if (Number(triangles[(int)j/3][k])==i1){
    670670                                                        if (Number(triangles[(int)j/3][(int)((k+1)%3)])==i2){
    671                                                                 bamgmesh->Segments[num*4+0]=Number(edges[i][0])+1; //back to M indexing
    672                                                                 bamgmesh->Segments[num*4+1]=Number(edges[i][1])+1; //back to M indexing
    673                                                                 bamgmesh->Segments[num*4+2]=(int)j/3+1;            //back to M indexing
    674                                                                 bamgmesh->Segments[num*4+3]=edges[i].ref;
     671                                                                bamgmesh->IssmSegments[num*4+0]=Number(edges[i][0])+1; //back to M indexing
     672                                                                bamgmesh->IssmSegments[num*4+1]=Number(edges[i][1])+1; //back to M indexing
     673                                                                bamgmesh->IssmSegments[num*4+2]=(int)j/3+1;            //back to M indexing
     674                                                                bamgmesh->IssmSegments[num*4+3]=edges[i].ref;
    675675                                                                num+=1;
    676676                                                                stop=true;
    … …  
    678678                                                        }
    679679                                                        if (Number(triangles[(int)j/3][(int)((k+2)%3)])==i2){
    680                                                                 bamgmesh->Segments[num*4+0]=Number(edges[i][1])+1; //back to M indexing
    681                                                                 bamgmesh->Segments[num*4+1]=Number(edges[i][0])+1; //back to M indexing
    682                                                                 bamgmesh->Segments[num*4+2]=(int)j/3+1;            //back to M indexing
    683                                                                 bamgmesh->Segments[num*4+3]=edges[i].ref;
     680                                                                bamgmesh->IssmSegments[num*4+0]=Number(edges[i][1])+1; //back to M indexing
     681                                                                bamgmesh->IssmSegments[num*4+1]=Number(edges[i][0])+1; //back to M indexing
     682                                                                bamgmesh->IssmSegments[num*4+2]=(int)j/3+1;            //back to M indexing
     683                                                                bamgmesh->IssmSegments[num*4+3]=edges[i].ref;
    684684                                                                num+=1;
    685685                                                                stop=true;
    … …  
    859859                int* connectivitysize_2=NULL;
    860860                int  connectivitymax_2=0;
    861                 i1=bamgmesh->ElementEdgesSize[0];
    862                 i2=bamgmesh->ElementEdgesSize[1];
     861                i1=bamgmesh->IssmEdgesSize[0];
     862                i2=bamgmesh->IssmEdgesSize[1];
    863863                head_2=(int*)xmalloc(nbv*sizeof(int));
    864864                next_2=(int*)xmalloc(2*i1*sizeof(int));
    … …  
    871871                for (i=0;i<i1;i++) {
    872872                        for (j=0;j<2;j++){
    873                                 int v=(int)bamgmesh->ElementEdges[i*i2+j]-1; //back to C indexing
     873                                int v=(int)bamgmesh->IssmEdges[i*i2+j]-1; //back to C indexing
    874874                                if (k>2*i1-1 || k<0) ISSMERROR("Index exceed matrix dimensions (k=%i not in [0 %i]",k,2*i1-1);
    875875                                next_2[k]=head_2[v];
    … …  
    892892                        for(j=head_2[i];j!=-1;j=next_2[j]){
    893893                                ISSMASSERT(connectivitymax_2*i+k < connectivitymax_2*nbv);
    894                                 num=(int)bamgmesh->ElementEdges[int(j/2)*i2+0];
     894                                num=(int)bamgmesh->IssmEdges[int(j/2)*i2+0];
    895895                                if (i+1==num){ //carefull, ElementEdge is in M indexing
    896896                                        //i is the first vertex of the edge, it is therefore connected to the second vertex
    897                                         bamgmesh->NodalConnectivity[connectivitymax_2*i+k]=bamgmesh->ElementEdges[int(j/2)*i2+1];
     897                                        bamgmesh->NodalConnectivity[connectivitymax_2*i+k]=bamgmesh->IssmEdges[int(j/2)*i2+1];
    898898                                }
    899899                                else{
    … …  
    948948                double hmin = Max(bamgopts->hmin,MinimalHmin());
    949949                double hmax = Min(bamgopts->hmax,MaximalHmax());
    950                 double coef = bamgopts->coef;
     950                double coef = bamgopts->coeff;
    951951
    952952                //for now we only use j==3
    … …  
    996996
    997997                /*Get options*/
    998                 int    verbosity=bamgopts->verbose;
     998                int    verbose=bamgopts->verbose;
    999999                double anisomax =bamgopts->anisomax;
    10001000                double errg     =bamgopts->errg;
    … …  
    12101210                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/BoundAnisotropy)*/
    12111211
    1212                 long int verbosity=0;
     1212                long int verbose=0;
    12131213                double lminaniso = 1/ (Max(hminaniso*hminaniso,1e-100));
    12141214
    12151215                //display info
    1216                 if (verbosity > 1)  printf("   BoundAnisotropy by %g\n",anisomax);
     1216                if (verbose > 1)  printf("   BoundAnisotropy by %g\n",anisomax);
    12171217
    12181218                double h1=1.e30,h2=1e-30;
    … …  
    12291229
    12301230                        //info to be displayed
    1231                         if (verbosity>2){
     1231                        if (verbose>2){
    12321232                                h1 =Min(h1,Vp.lmin());
    12331233                                h2 =Max(h2,Vp.lmax());
    … …  
    12401240
    12411241                //display info
    1242                 if (verbosity>2){
     1242                if (verbose>2){
    12431243                        printf("      input:  Hmin = %g, Hmax = %g, factor of anisotropy max  = %g\n",pow(h2,-0.5),pow(h1,-0.5),pow(rx,0.5));
    12441244                        printf("      output: Hmin = %g, Hmax = %g, factor of anisotropy max  = %g\n",pow(hn2,-0.5),pow(hn1,-0.5),pow(rnx,0.5));
    … …  
    12541254                /*Intermediary*/
    12551255                int i,j,k,kk,it,jt;
    1256                 int    verbosity=0;
     1256                int    verbose=0;
    12571257                double cutoffradian=10*Pi/180;
    12581258
    12591259                /*Recover options*/
    12601260                if (bamgopts){
    1261                         verbosity=bamgopts->verbose;
    1262                         cutoffradian=bamgopts->MaximalAngleOfCorner*Pi/180;
     1261                        verbose=bamgopts->verbose;
     1262                        cutoffradian=bamgopts->MaxCornerAngle*Pi/180;
    12631263                }
    12641264
    12651265                //display info
    1266                 if (verbosity>1) printf("   construction of the geometry from the 2d mesh\n");
     1266                if (verbose>1) printf("   construction of the geometry from the 2d mesh\n");
    12671267
    12681268                //check that the mesh is not empty
    … …  
    12711271                }
    12721272
    1273                 //Gh is the geometry of the mesh (this), initialize MaximalAngleOfCorner
    1274                 if (cutoffradian>=0) Gh.MaximalAngleOfCorner = cutoffradian;
     1273                //Gh is the geometry of the mesh (this), initialize MaxCornerAngle
     1274                if (cutoffradian>=0) Gh.MaxCornerAngle = cutoffradian;
    12751275
    12761276                /*Construction of the edges*/
    … …  
    13371337
    13381338                //display info
    1339                 if(verbosity>5) {
     1339                if(verbose>5) {
    13401340                        printf("         info on Mesh:\n");
    13411341                        printf("            - number of vertices    = %i \n",nbv);
    … …  
    13741374
    13751375                        //display info
    1376                         if(verbosity>4) printf("   Construction of the edges %i\n",nbe);
     1376                        if(verbose>4) printf("   Construction of the edges %i\n",nbe);
    13771377
    13781378                        for (i=0;i<nbedges;i++){
    … …  
    15151515                        }
    15161516                }
    1517                 if (verbosity> 3) printf("      The Number of sub domain = %i\n",NbSubDomains);
     1517                if (verbose> 3) printf("      The Number of sub domain = %i\n",NbSubDomains);
    15181518
    15191519                //build subdomains
    … …  
    15661566                Gh.NbSubDomains = NbSubDomains;
    15671567                Gh.subdomains = new GeometricalSubDomain[NbSubDomains];
    1568                 if (verbosity>3) printf("   number of vertices = %i\n   number of edges = %i\n",Gh.nbv,Gh.nbe);
     1568                if (verbose>3) printf("   number of vertices = %i\n   number of edges = %i\n",Gh.nbv,Gh.nbe);
    15691569                NbVerticesOnGeomVertex = Gh.nbv;
    15701570                VerticesOnGeomVertex = new VertexOnGeom[NbVerticesOnGeomVertex];
    … …  
    17001700                double* s=NULL;
    17011701                long    nbsol;
    1702                 int     verbosity;
     1702                int     verbose;
    17031703
    17041704                int   i,j,k,iA,iB,iC;
    … …  
    17061706
    17071707                /*Recover options*/
    1708                 verbosity=bamgopts->verbose;
     1708                verbose=bamgopts->verbose;
    17091709
    17101710                /*Get and process fields*/
    … …  
    17311731
    17321732                //display infos
    1733                 if(verbosity>1) {
     1733                if(verbose>1) {
    17341734                        printf("   Construction of Metric: number of field: %i (nbt=%i, nbv=%i)\n",nbsol,nbt,nbv);
    17351735                }
    … …  
    18041804
    18051805                        //display info
    1806                         if(verbosity>2) printf("      Solution %i, Min = %g, Max = %g, Delta = %g\n",nusol,smin,smax,sdelta);
     1806                        if(verbose>2) printf("      Solution %i, Min = %g, Max = %g, Delta = %g\n",nusol,smin,smax,sdelta);
    18071807
    18081808                        //skip constant field
    … …  
    19041904                long nbsol;
    19051905                int NbJacobi;
    1906                 int verbosity;
     1906                int verbose;
    19071907
    19081908                /*Recover options*/
    1909                 verbosity=bamgopts->verbose;
     1909                verbose=bamgopts->verbose;
    19101910                NbJacobi=bamgopts->nbjacobi;
    19111911
    … …  
    19301930
    19311931                //display infos
    1932                 if(verbosity>1) {
     1932                if(verbose>1) {
    19331933                        printf("   Construction of Metric: number of field: %i (nbt=%i, nbv=%i)\n",nbsol,nbt,nbv);
    19341934                }
    … …  
    20142014
    20152015                        //display info
    2016                         if(verbosity>2) printf("      Solution %i, Min = %g, Max = %g, Delta = %g, number of fields = %i\n",nusol,smin,smax,sdelta,nbsol);
     2016                        if(verbose>2) printf("      Solution %i, Min = %g, Max = %g, Delta = %g, number of fields = %i\n",nusol,smin,smax,sdelta,nbsol);
    20172017
    20182018                        //skip constant field
    20192019                        if (sdelta < 1.0e-10*Max(absmax,1e-20) ){
    2020                                 if (verbosity>2) printf("      Solution %i is constant, skipping...\n",nusol);
     2020                                if (verbose>2) printf("      Solution %i is constant, skipping...\n",nusol);
    20212021                                continue;
    20222022                        }
    … …  
    23352335                        /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/ForceBoundary)*/
    23362336
    2337                         long int verbosity=2;
     2337                        long int verbose=2;
    23382338                        int k=0;
    23392339                        int nbfe=0,nbswp=0,Nbswap=0;
    23402340
    23412341                        //display
    2342                         if (verbosity > 2) printf("   ForceBoundary  nb of edge: %i\n",nbe);
     2342                        if (verbose > 2) printf("   ForceBoundary  nb of edge: %i\n",nbe);
    23432343
    23442344                        //check that there is no triangle with 0 determinant
    … …  
    23712371                                Nbswap +=  vertices[j].Optim(1,0);
    23722372                        }
    2373                         if (verbosity > 3) printf("      number of inforced edge = %i, number of swap= %i\n",nbfe,Nbswap);
     2373                        if (verbose > 3) printf("      number of inforced edge = %i, number of swap= %i\n",nbfe,Nbswap);
    23742374                }
    23752375        /*}}}1*/
    … …  
    23782378                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/FindSubDomain)*/
    23792379
    2380                 long int verbosity=0;
    2381 
    2382                 if (verbosity >2){
     2380                long int verbose=0;
     2381
     2382                if (verbose >2){
    23832383                        if (OutSide) printf("   Find all external sub-domain\n");
    23842384                        else printf("   Find all internal sub-domain\n");
    … …  
    25252525                                                                 }//while (t)
    25262526                                                                }
    2527                                                         if(verbosity>4) printf("      Number of removes subdomains (OutSideMesh) = %i\n",NbSubDomains-j);
     2527                                                        if(verbose>4) printf("      Number of removes subdomains (OutSideMesh) = %i\n",NbSubDomains-j);
    25282528                                                        NbSubDomains=j;
    25292529                                                  }
    … …  
    26102610
    26112611                                if (inew < NbSubDomains) {
    2612                                         if (verbosity>5) printf("WARNING: %i SubDomains are being removed\n",NbSubDomains-inew);
     2612                                        if (verbose>5) printf("WARNING: %i SubDomains are being removed\n",NbSubDomains-inew);
    26132613                                        NbSubDomains=inew;}
    26142614
    … …  
    33753375
    33763376                /*Get options*/
    3377                 long int verbosity=2;
     3377                long int verbose=2;
    33783378
    33793379                //Display info
    3380                 if (verbosity>2) printf("   Insert initial %i vertices\n",nbv);
     3380                if (verbose>2) printf("   Insert initial %i vertices\n",nbv);
    33813381
    33823382                //Compute integer coordinates and determinants for the existing vertices (from Geometry)
    … …  
    34693469
    34703470                long NbSwap=0;
    3471                 if (verbosity>3) printf("   Begining of insertion process...\n");
     3471                if (verbose>3) printf("   Begining of insertion process...\n");
    34723472
    34733473                for (int icount=2; icount<nbv; icount++) {
    … …  
    34913491
    34923492                //Display info
    3493                 if (verbosity>3) {
     3493                if (verbose>3) {
    34943494                        printf("      NbSwap of insertion: %i\n",NbSwap);
    34953495                        printf("      NbSwap/nbv:          %i\n",NbSwap/nbv);
    … …  
    35063506                        for (int is1=0; is1<nbv; is1++)
    35073507                         NbSwap += ordre[is1]->Optim(0,0);
    3508                         if (verbosity>3) {
     3508                        if (verbose>3) {
    35093509                                printf("      Optim Loop: %i\n",Nbloop);
    35103510                                printf("      NbSwap/nbv:          %i\n",NbSwap/nbv);
    … …  
    35213521                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/InsertNewPoints)*/
    35223522
    3523                 long int verbosity=0;
     3523                long int verbose=0;
    35243524                double seuil= 1.414/2 ;// for two close point
    35253525                long i;
    … …  
    35313531
    35323532                //display info if required
    3533                 if (verbosity>5) printf("      Try to Insert %i new points\n",nbvnew);
     3533                if (verbose>5) printf("      Try to Insert %i new points\n",nbvnew);
    35343534
    35353535                //return if no new points
    … …  
    35803580                        }
    35813581                }
    3582                 if (verbosity>3) {
     3582                if (verbose>3) {
    35833583                        printf("         number of new points: %i\n",iv);
    35843584                        printf("         number of to close (?) points: %i\n",nbv-iv);
    … …  
    35883588
    35893589                for (i=nbvold;i<nbv;i++) NbSwap += vertices[i].Optim(1); 
    3590                 if (verbosity>3) printf("   NbSwap=%i\n",NbSwap);
     3590                if (verbose>3) printf("   NbSwap=%i\n",NbSwap);
    35913591
    35923592                NbTSwap +=  NbSwap ;
    … …  
    36413641                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/MakeQuadrangles)*/
    36423642
    3643                 long int verbosity=0;
    3644 
    3645                 if (verbosity>2) printf("MakeQuadrangles costheta = %g\n",costheta);
     3643                long int verbose=0;
     3644
     3645                if (verbose>2) printf("MakeQuadrangles costheta = %g\n",costheta);
    36463646
    36473647                if (costheta >1) {
    3648                         if (verbosity>5) printf("   do nothing: costheta > 1\n");
     3648                        if (verbose>5) printf("   do nothing: costheta > 1\n");
    36493649                }
    36503650
    … …  
    36723672                          }
    36733673                        NbOfQuad = kk;
    3674                         if (verbosity>2){
     3674                        if (verbose>2){
    36753675                                printf("   number of quadrilaterals    = %i\n",NbOfQuad);
    36763676                                printf("   number of triangles         = %i\n",nbt-NbOutT- NbOfQuad*2);
    … …  
    36843684                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/MakeQuadTree)*/
    36853685
    3686                 long int verbosity=0;
     3686                long int verbose=0;
    36873687                if (  !quadtree )  quadtree = new QuadTree(this);
    36883688
    … …  
    36933693                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/MaxSubDivision)*/
    36943694
    3695                 long int verbosity=0;
     3695                long int verbose=0;
    36963696
    36973697                const  double maxsubdiv2 = maxsubdiv*maxsubdiv;
    3698                 if(verbosity>1) printf("   Limit the subdivision of a edges in the new mesh by %g\n",maxsubdiv);
     3698                if(verbose>1) printf("   Limit the subdivision of a edges in the new mesh by %g\n",maxsubdiv);
    36993699                // for all the edges
    37003700                // if the len of the edge is to long
    … …  
    37383738                        }
    37393739                }
    3740                 if(verbosity>3){
     3740                if(verbose>3){
    37413741                        printf("      number of metric changes = %i, maximum number of subdivision of a edges before change = %g\n",nbchange,pow(lmax,0.5));
    37423742                }
    … …  
    37843784
    37853785                /*Recover options*/
    3786                 int verbosity=bamgopts->verbose;
     3786                int verbose=bamgopts->verbose;
    37873787
    37883788                /*First, insert old points if requested*/
    37893789                if (KeepVertices && (&Bh != this) && (nbv+Bh.nbv< nbvx)){
    3790                         if (verbosity>5) printf("         Inserting initial mesh points\n");
     3790                        if (verbose>5) printf("         Inserting initial mesh points\n");
    37913791                        for (i=0;i<Bh.nbv;i++){
    37923792                                MeshVertex &bv=Bh[i];
    … …  
    38073807                // Big loop (most time consuming)
    38083808                int iter=0;
    3809                 if (verbosity>5) printf("         Big loop\n");
     3809                if (verbose>5) printf("         Big loop\n");
    38103810                do {
    38113811                        /*Update variables*/
    … …  
    38943894        /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/PreInit)*/
    38953895
    3896         long int verbosity=0;
     3896        long int verbose=0;
    38973897
    38983898        srand(19999999);
    … …  
    41024102
    41034103        /*Intermediary*/
    4104         int verbosity=0;
     4104        int verbose=0;
    41054105
    41064106        // generation of the integer coordinate
    … …  
    41084108        // find extrema coordinates of vertices pmin,pmax
    41094109        long i;
    4110         if(verbosity>2) printf("      Reconstruct mesh of %i vertices\n",nbv);
     4110        if(verbose>2) printf("      Reconstruct mesh of %i vertices\n",nbv);
    41114111
    41124112        //initialize ordre
    … …  
    41644164
    41654165        //Display info if required
    4166         if(verbosity>5) {
     4166        if(verbose>5) {
    41674167                printf("         info of Mesh:\n");
    41684168                printf("            - number of vertices    = %i \n",nbv);
    … …  
    44004400                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/ReNumberingTheTriangleBySubDomain)*/
    44014401
    4402                 long int verbosity=0;
     4402                long int verbose=0;
    44034403                long *renu= new long[nbt];
    44044404                register Triangle *t0,*t,*te=triangles+nbt;
    … …  
    47174717        /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/SmoothingVertex)*/
    47184718
    4719         long int verbosity=0;
     4719        long int verbose=0;
    47204720        //  if quatree exist remove it end reconstruct
    47214721        if (quadtree) delete quadtree;
    … …  
    47364736        for ( k=0;k<NbVerticesOnGeomEdge;k++ )
    47374737         tstart[ Number(VerticesOnGeomEdge[k].mv)]=&vide;
    4738         if(verbosity>2) printf("   SmoothingVertex: nb Iteration = %i, Omega=%g\n",nbiter,omega);
     4738        if(verbose>2) printf("   SmoothingVertex: nb Iteration = %i, Omega=%g\n",nbiter,omega);
    47394739        for (k=0;k<nbiter;k++)
    47404740          {
    … …  
    47484748                  if (tstart[i] != &vide) // not a boundary vertex
    47494749                        NbSwap += vertices[i].Optim(1);
    4750                 if (verbosity>3) printf("      move max = %g, iteration = %i, nb of swap = %i\n",pow(delta,0.5),k,NbSwap);
     4750                if (verbose>3) printf("      move max = %g, iteration = %i, nb of swap = %i\n",pow(delta,0.5),k,NbSwap);
    47514751          }
    47524752
    … …  
    47594759        /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/SmoothMetric)*/
    47604760
    4761         long int verbosity=0;
     4761        long int verbose=0;
    47624762
    47634763        if(raisonmax<1.1) return;
    4764         if(verbosity > 1) printf("   Triangles::SmoothMetric raisonmax = %g\n",raisonmax);
     4764        if(verbose > 1) printf("   Triangles::SmoothMetric raisonmax = %g\n",raisonmax);
    47654765        ReMakeTriangleContainingTheVertex();
    47664766        long i,j,kch,kk,ip;
    … …  
    48284828                Exchange(first_np_or_next_t0,first_np_or_next_t1);
    48294829        }
    4830         if(verbosity>2) printf("      number of iterations = %i\n",kch);
     4830        if(verbose>2) printf("      number of iterations = %i\n",kch);
    48314831        delete [] first_np_or_next_t0;
    48324832        delete [] first_np_or_next_t1;
    … …  
    48374837                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, MeshQuad.cpp/SplitElement)*/
    48384838
    4839                 long int verbosity=0;
     4839                long int verbose=0;
    48404840
    48414841                Direction NoDirOfSearch;
    … …  
    53205320                ReconstructExistingMesh();
    53215321
    5322                 if (verbosity>2){
     5322                if (verbose>2){
    53235323                        printf("   number of quadrilaterals    = %i\n",NbOfQuad);
    53245324                        printf("   number of triangles         = %i\n",nbt-NbOutT- NbOfQuad*2);
    … …  
    53515351        long it;
    53525352        long nbvold=nbv;
    5353         long int verbosity=2;
     5353        long int verbose=2;
    53545354        for (it=0;it<nbt;it++){
    53555355                Triangle &t=triangles[it];
    … …  
    54025402                        //      }
    54035403        }
    5404         if (verbosity>3) {
     5404        if (verbose>3) {
    54055405                printf("   number of points: %i\n",iv);
    54065406                printf("   number of swap to  split internal edges with border vertices: %i\n",NbSwap);
    … …  
    54095409}
    54105410if (NbSplitEdge>nbv-nbvold) printf("WARNING: not enough vertices  to split all internal edges, we lost %i edges...\n",NbSplitEdge - ( nbv-nbvold));
    5411 if (verbosity>2) printf("SplitInternalEdgeWithBorderVertices: Number of splited edge %i\n",NbSplitEdge);
     5411if (verbose>2) printf("SplitInternalEdgeWithBorderVertices: Number of splited edge %i\n",NbSplitEdge);
    54125412
    54135413return  NbSplitEdge;
    … …  
    54185418        /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/ConsRefTriangle)*/
    54195419
    5420         long int verbosity=0;
     5420        long int verbose=0;
    54215421        register Triangle *t0,*t;
    54225422        register long k=0, num;   
Note: See TracChangeset for help on using the changeset viewer.