Changeset 3031


Ignore:
Timestamp:
02/16/10 08:58:07 (15 years ago)
Author:
Mathieu Morlighem
Message:

better comments

Location:
issm/trunk/src/c/Bamgx/objects
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r3022 r3031  
    7777
    7878                //Vertices
    79                 if(verbose>3) printf("      writing Vertices\n");
     79                if(verbose>5) printf("      writing Vertices\n");
    8080                bamggeom->NumVertices=nbv;
    8181                xfree((void**)&bamggeom->Vertices);
    … …  
    9696
    9797                //Edges
    98                 if(verbose>3) printf("      writing Edges\n");
     98                if(verbose>5) printf("      writing Edges\n");
    9999                bamggeom->NumEdges=nbe;
    100100                xfree((void**)&bamggeom->Edges);
    … …  
    121121
    122122                //CrackedEdges
    123                 if(verbose>3) printf("      writing CrackedEdges\n");
     123                if(verbose>5) printf("      writing CrackedEdges\n");
    124124                bamggeom->NumCrackedEdges=nbcracked;
    125125                xfree((void**)&bamggeom->CrackedEdges);
    … …  
    142142
    143143                //RequiredEdges
    144                 if(verbose>3) printf("      writing RequiredEdges\n");
     144                if(verbose>5) printf("      writing RequiredEdges\n");
    145145                bamggeom->NumRequiredEdges=nbreq;
    146146                xfree((void**)&bamggeom->RequiredEdges);
    … …  
    162162
    163163                //RequiredVertices
    164                 if(verbose>3) printf("      writing RequiredVertices\n");
     164                if(verbose>5) printf("      writing RequiredVertices\n");
    165165                bamggeom->NumRequiredVertices=nbreqv;
    166166                xfree((void**)&bamggeom->RequiredVertices);
    … …  
    180180
    181181                //SubDomains
    182                 if(verbose>3) printf("      writing SubDomains\n");
     182                if(verbose>5) printf("      writing SubDomains\n");
    183183                bamggeom->NumSubDomains=NbSubDomains;
    184184                xfree((void**)&bamggeom->SubDomains);
    … …  
    197197
    198198                //TangentAtEdges
    199                 if(verbose>3) printf("      writing TangentAtEdges\n");
     199                if(verbose>5) printf("      writing TangentAtEdges\n");
    200200                bamggeom->NumTangentAtEdges=nbtan;
    201201                xfree((void**)&bamggeom->TangentAtEdges);
    … …  
    254254                //Vertices
    255255                if (bamggeom->Vertices){
    256                         if(verbose>3) printf("      processing Vertices\n");
     256                        if(verbose>5) printf("      processing Vertices\n");
    257257                        vertices = new GeometricalVertex[nbvx];
    258258                        for (i=0;i<nbv;i++) {
    … …  
    293293                        Real4* len=NULL;
    294294
    295                         if(verbose>3) printf("      processing Edges\n");
     295                        if(verbose>5) printf("      processing Edges\n");
    296296                        edges = new GeometricalEdge[nbe];
    297297
    … …  
    346346                //CrackedEdges
    347347                if (bamggeom->CrackedEdges){
    348                         if(verbose>3) printf("      processing CrackedEdges\n");
     348                        if(verbose>5) printf("      processing CrackedEdges\n");
    349349                        for (i=0;i<bamggeom->NumCrackedEdges;i++){
    350350                                throw ErrorException(__FUNCT__,exprintf("Not implemented yet"));
    … …  
    352352                }
    353353                else{
    354                         if(verbose>3) printf("      no CrackedEdges found\n");
     354                        if(verbose>5) printf("      no CrackedEdges found\n");
    355355                }
    356356
    357357                //hVertices
    358358                if(bamggeom->hVertices){
    359                         if(verbose>3) printf("      processing hVertices\n");
     359                        if(verbose>5) printf("      processing hVertices\n");
    360360                        for (i=0;i< nbv;i++){
    361361                                vertices[i].m=Metric((Real4)bamggeom->hVertices[i]);
    … …  
    363363                }
    364364                else{
    365                         if(verbose>3) printf("      no hVertices found\n");
     365                        if(verbose>5) printf("      no hVertices found\n");
    366366                }
    367367
    368368                //MetricVertices
    369369                if(bamggeom->MetricVertices){
    370                         if(verbose>3) printf("      processing MetricVertices\n");
     370                        if(verbose>5) printf("      processing MetricVertices\n");
    371371                        hvertices=1;
    372372                        for (i=0;i< nbv;i++) {
    … …  
    375375                }
    376376                else{
    377                         if(verbose>3) printf("      no MetricVertices found\n");
     377                        if(verbose>5) printf("      no MetricVertices found\n");
    378378                }
    379379
    380380                //h1h2VpVertices
    381381                if(bamggeom->h1h2VpVertices){
    382                         if(verbose>3) printf("      processing h1h2VpVertices\n");
     382                        if(verbose>5) printf("      processing h1h2VpVertices\n");
    383383                        Real4 h1,h2,v1,v2;
    384384                        hvertices =1;
    … …  
    392392                }
    393393                else{
    394                         if(verbose>3) printf("      no h1h2VpVertices found\n");
     394                        if(verbose>5) printf("      no h1h2VpVertices found\n");
    395395                }
    396396
    397397                //MaximalAngleOfCorner
    398398                if (bamgopts->MaximalAngleOfCorner){
    399                         if(verbose>3) printf("      processing MaximalAngleOfCorner\n");
     399                        if(verbose>5) printf("      processing MaximalAngleOfCorner\n");
    400400                        MaximalAngleOfCorner=bamgopts->MaximalAngleOfCorner*Pi/180;
    401401                }
    402402                else{
    403                         if(verbose>3) printf("      no MaximalAngleOfCorner found\n");
     403                        if(verbose>5) printf("      no MaximalAngleOfCorner found\n");
    404404                }
    405405
    406406                //TangentAtEdges
    407407                if (bamggeom->TangentAtEdges){
    408                         if(verbose>3) printf("      processing TangentAtEdges");
     408                        if(verbose>5) printf("      processing TangentAtEdges");
    409409                        int n,i,j,k;
    410410                        R2 tg;
    … …  
    423423                }
    424424                else{
    425                         if(verbose>3) printf("      no TangentAtEdges found\n");
     425                        if(verbose>5) printf("      no TangentAtEdges found\n");
    426426                }
    427427
    428428                //Corners
    429429                if(bamggeom->Corners){
    430                         if(verbose>3) printf("      processing Corners");
     430                        if(verbose>5) printf("      processing Corners");
    431431                        n=bamggeom->NumCorners;
    432432                        for (i=0;i<n;i++) {     
    … …  
    439439                }
    440440                else{
    441                         if(verbose>3) printf("      no Corners found\n");
     441                        if(verbose>5) printf("      no Corners found\n");
    442442                }
    443443
    444444                //RequiredVertices
    445445                if(bamggeom->RequiredVertices){
    446                         if(verbose>3) printf("      processing RequiredVertices");
     446                        if(verbose>5) printf("      processing RequiredVertices");
    447447                        n=bamggeom->NumRequiredVertices;
    448448                        for (i=0;i<n;i++) {     
    … …  
    454454                }
    455455                else{
    456                         if(verbose>3) printf("      no RequiredVertices found\n");
     456                        if(verbose>5) printf("      no RequiredVertices found\n");
    457457                }
    458458
    459459                //RequiredEdges
    460460                if(bamggeom->RequiredEdges){
    461                         if(verbose>3) printf("      processing RequiredEdges");
     461                        if(verbose>5) printf("      processing RequiredEdges");
    462462                        n=bamggeom->NumRequiredEdges;
    463463                        for (i=0;i<n;i++) {     
    … …  
    470470                }
    471471                else{
    472                         if(verbose>3) printf("      no RequiredEdges found\n");
     472                        if(verbose>5) printf("      no RequiredEdges found\n");
    473473                }
    474474
    … …  
    476476                if(bamggeom->SubDomains){
    477477                        Int4 i0,i1,i2,i3;
    478                         if(verbose>3) printf("      processing SubDomains\n");
     478                        if(verbose>5) printf("      processing SubDomains\n");
    479479                        NbSubDomains=bamggeom->NumSubDomains;
    480480                        subdomains = new GeometricalSubDomain[NbSubDomains];
    … …  
    492492                }
    493493                else{
    494                         if(verbose>3) printf("      no SubDomains found\n");
     494                        if(verbose>5) printf("      no SubDomains found\n");
    495495                }
    496496        }

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

    r3022 r3031  
    247247                //Vertices
    248248                if(bamgmesh->Vertices){
    249                         if(verbose>3) printf("      processing Vertices\n");
     249                        if(verbose>5) printf("      processing Vertices\n");
    250250
    251251                        vertices=(Vertex*)xmalloc(nbv*sizeof(Vertex));
    … …  
    263263                }
    264264                else{
    265                         if(verbose>3) throw ErrorException(__FUNCT__,exprintf("no Vertices found in the initial mesh"));
     265                        if(verbose>5) throw ErrorException(__FUNCT__,exprintf("no Vertices found in the initial mesh"));
    266266                }
    267267
    268268                //Triangles
    269269                if(bamgmesh->Triangles){
    270                         if(verbose>3) printf("      processing Triangles\n");
     270                        if(verbose>5) printf("      processing Triangles\n");
    271271                        triangles =new Triangle[nbtx]; //we cannot allocate only nbt triangles since
    272272                        //other triangles will be added for each edge
    … …  
    281281                }
    282282                else{
    283                         if(verbose>3) throw ErrorException(__FUNCT__,exprintf("no Triangles found in the initial mesh"));
     283                        if(verbose>5) throw ErrorException(__FUNCT__,exprintf("no Triangles found in the initial mesh"));
    284284                }
    285285
    286286                //Quadrilaterals
    287287                if(bamgmesh->Quadrilaterals){
    288                         if(verbose>3) printf("      processing Quadrilaterals\n");
     288                        if(verbose>5) printf("      processing Quadrilaterals\n");
    289289                        Int4 i1,i2,i3,i4,iref;
    290290                        triangles =new Triangle[nbt];
    … …  
    306306                }
    307307                else{
    308                         if(verbose>3) printf("      no Quadrilaterals found\n");
     308                        if(verbose>5) printf("      no Quadrilaterals found\n");
    309309                }
    310310
    311311                //hVertices
    312312                if(bamgmesh->hVertices){
    313                         if(verbose>3) printf("      processing hVertices\n");
     313                        if(verbose>5) printf("      processing hVertices\n");
    314314                        hvertices=1;
    315315                        for (i=0;i< nbv;i++){
    … …  
    320320                }
    321321                else{
    322                         if(verbose>3) printf("      no hVertices found\n");
     322                        if(verbose>5) printf("      no hVertices found\n");
    323323                }
    324324
    325325                //VerVerticesOnGeometricEdge
    326326                if(bamgmesh->VerticesOnGeometricEdge){
    327                         if(verbose>3) printf("      processing VerticesOnGeometricEdge\n");
     327                        if(verbose>5) printf("      processing VerticesOnGeometricEdge\n");
    328328                        NbVerticesOnGeomEdge=bamgmesh->NumVerticesOnGeometricEdge;
    329329                        VerticesOnGeomEdge= new  VertexOnGeom[NbVerticesOnGeomEdge] ;
    … …  
    338338                }
    339339                else{
    340                         if(verbose>3) printf("      no VertexOnGeometricEdge found\n");
     340                        if(verbose>5) printf("      no VertexOnGeometricEdge found\n");
    341341                }
    342342
    … …  
    347347                        Real4 *len =0;
    348348
    349                         if(verbose>3) printf("      processing Edges\n");
     349                        if(verbose>5) printf("      processing Edges\n");
    350350                        nbe=bamgmesh->NumEdges;
    351351                        edges = new Edge[nbe];
    … …  
    411411                }
    412412                else{
    413                         if(verbose>3) printf("      no Edges found\n");
     413                        if(verbose>5) printf("      no Edges found\n");
    414414                }
    415415                //CrackedEdges
    416416                if(bamgmesh->CrackedEdges){
    417                         if(verbose>3) printf("      processing CrackedEdges\n");
     417                        if(verbose>5) printf("      processing CrackedEdges\n");
    418418                        NbCrackedEdges=bamgmesh->NumCrackedEdges;
    419419                        CrackedEdges=new CrackedEdge[NbCrackedEdges];
    … …  
    426426                }
    427427                else{
    428                         if(verbose>3) printf("      no CrackedEdges found\n");
     428                        if(verbose>5) printf("      no CrackedEdges found\n");
    429429                }
    430430
    431431                //EdgeOnGeometricEdge
    432432                if(bamgmesh->EdgesOnGeometricEdge){
    433                         if(verbose>3) printf("      processing EdgesOnGeometricEdge\n");
     433                        if(verbose>5) printf("      processing EdgesOnGeometricEdge\n");
    434434                        int i1,i2,i,j;
    435435                        i2=bamgmesh->NumEdgesOnGeometricEdge;
    … …  
    444444                }
    445445                else{
    446                         if(verbose>3) printf("      no EdgesOnGeometricEdge found\n");
     446                        if(verbose>5) printf("      no EdgesOnGeometricEdge found\n");
    447447                }
    448448
    … …  
    450450                if(bamgmesh->SubDomains){
    451451                        Int4 i3,head,sens;
    452                         if(verbose>3) printf("      processing SubDomains\n");
     452                        if(verbose>5) printf("      processing SubDomains\n");
    453453                        NbSubDomains=bamgmesh->NumSubDomains;
    454454                        subdomains = new SubDomain [ NbSubDomains ];
    … …  
    463463                }
    464464                else{
    465                         if(verbose>3) printf("      no SubDomains found\n");
     465                        if(verbose>5) printf("      no SubDomains found\n");
    466466                }
    467467
    … …  
    483483                //renumber
    484484                if (bamgopts->renumber){
    485                         if(verbose>3) printf("      Renumbering...");
     485                        if(verbose>5) printf("      Renumbering...");
    486486                        Renumber(bamgopts);
    487                         if(verbose>3) printf(" done\n");
     487                        if(verbose>5) printf(" done\n");
    488488                }
    489489
    … …  
    493493
    494494                //Vertices
    495                 if(verbose>3) printf("      writing Vertices\n");
     495                if(verbose>5) printf("      writing Vertices\n");
    496496                bamgmesh->NumVertices=nbv;
    497497                xfree((void**)&bamgmesh->Vertices);
    … …  
    509509
    510510                //Edges
    511                 if(verbose>3) printf("      writing Edges\n");
     511                if(verbose>5) printf("      writing Edges\n");
    512512                bamgmesh->NumEdges=nbe;
    513513                xfree((void**)&bamgmesh->Edges);
    … …  
    530530                //Segments
    531531                bamgmesh->NumSegments=0;
    532                 if(verbose>3) printf("      writing Segments\n");
     532                if(verbose>5) printf("      writing Segments\n");
    533533
    534534                //chaining algorithm
    … …  
    599599
    600600                //CrackedEdges
    601                 if(verbose>3) printf("      writing CrackedEdges\n");
     601                if(verbose>5) printf("      writing CrackedEdges\n");
    602602                bamgmesh->NumCrackedEdges=NbCrackedEdges;
    603603                xfree((void**)&bamgmesh->CrackedEdges);
    … …  
    614614
    615615                //Triangles
    616                 if(verbose>3) printf("      writing Triangles\n");
     616                if(verbose>5) printf("      writing Triangles\n");
    617617                k=nbInT-NbOfQuad*2;
    618618                num=0;
    … …  
    638638
    639639                //Quadrilaterals
    640                 if(verbose>3) printf("      writing Quadrilaterals\n");
     640                if(verbose>5) printf("      writing Quadrilaterals\n");
    641641                bamgmesh->NumQuadrilaterals=NbOfQuad;
    642642                xfree((void**)&bamgmesh->Quadrilaterals);
    … …  
    662662
    663663                //SubDomains
    664                 if(verbose>3) printf("      writing SubDomains\n");
     664                if(verbose>5) printf("      writing SubDomains\n");
    665665                bamgmesh->NumSubDomains=NbSubDomains;
    666666                xfree((void**)&bamgmesh->SubDomains);
    … …  
    679679
    680680                //SubDomainsFromGeom
    681                 if(verbose>3) printf("      writing SubDomainsFromGeom\n");
     681                if(verbose>5) printf("      writing SubDomainsFromGeom\n");
    682682                bamgmesh->NumSubDomainsFromGeom=Gh.NbSubDomains;
    683683                xfree((void**)&bamgmesh->SubDomainsFromGeom);
    … …  
    696696
    697697                //VerticesOnGeomVertex
    698                 if(verbose>3) printf("      writing VerticesOnGeometricVertex\n");
     698                if(verbose>5) printf("      writing VerticesOnGeometricVertex\n");
    699699                bamgmesh->NumVerticesOnGeometricVertex=NbVerticesOnGeomVertex;
    700700                xfree((void**)&bamgmesh->VerticesOnGeometricVertex);
    … …  
    715715
    716716                //VertexOnGeometricEdge
    717                 if(verbose>3) printf("      writing VerticesOnGeometricEdge\n");
     717                if(verbose>5) printf("      writing VerticesOnGeometricEdge\n");
    718718                bamgmesh->NumVerticesOnGeometricEdge=NbVerticesOnGeomEdge;
    719719                xfree((void**)&bamgmesh->VerticesOnGeometricEdge);
    … …  
    735735
    736736                //EdgesOnGeometricEdge
    737                 if(verbose>3) printf("      writing EdgesOnGeometricEdge\n");
     737                if(verbose>5) printf("      writing EdgesOnGeometricEdge\n");
    738738                k=0;
    739739                for (i=0;i<nbe;i++){
    … …  
    28062806
    28072807                //initialize this
    2808                 if (verbose>3) printf("      Generating vertices\n");
     2808                if (verbose>3) printf("      Generating Boundary vertices\n");
    28092809                PreInit(inbvx);
    28102810                nbv=0;
    … …  
    30863086
    30873087                //Insert points inside existing triangles
    3088                 if (verbose>3) printf("      Inserting vertices in mesh\n");
     3088                if (verbose>4) printf("      -- current number of vertices = %i\n",nbv);
     3089                if (verbose>3) printf("      Creating initial Constrained Delaunay Triangulation...\n");
     3090                if (verbose>3) printf("         Inserting points\n");
    30893091                Insert();
    30903092
    30913093                //Force the boundary
    3092                 if (verbose>3) printf("      Forcing boundaries\n");
     3094                if (verbose>3) printf("         Forcing boundaries\n");
    30933095                ForceBoundary();
    30943096
    30953097                //Extract SubDomains
    3096                 if (verbose>3) printf("      Extracting subdomains\n");
     3098                if (verbose>3) printf("         Extracting subdomains\n");
    30973099                FindSubDomain();
    30983100
    30993101                // NewPointsOld(*this) ;
    3100                 if (verbose>3) printf("      Generating mesh properties\n");
     3102                if (verbose>3) printf("      Inserting internal points\n");
    31013103                NewPoints(*this,0) ;
     3104                if (verbose>4) printf("      -- current number of vertices = %i\n",nbv);
    31023105        }
    31033106        /*}}}1*/
    … …  
    39583961                        i=Headt;
    39593962                        next_t=-first_np_or_next_t[i];
    3960                         for(t=&triangles[i];i<nbt;t=&triangles[i=next_t],next_t=-first_np_or_next_t[i]) 
    3961                           { // for each triangle  t
     3963                        for(t=&triangles[i];i<nbt;t=&triangles[i=next_t],next_t=-first_np_or_next_t[i]){
     3964                                // for each triangle  t
    39623965                                // we can change first_np_or_next_t[i]
    39633966                                if (i<0 || i>=nbt ){
    … …  
    39653968                                }
    39663969                                first_np_or_next_t[i] = iter;
    3967                                 for(j=0;j<3;j++)
    3968                                   { // for each edge
     3970                                for(j=0;j<3;j++){
     3971                                        // for each edge
    39693972                                        TriangleAdjacent tj(t,j);
    39703973                                        Vertex & vA = * tj.EdgeVertex(0);
    … …  
    39984001                         break;
    39994002
    4000                         for (i=nbtold;i<nbt;i++)
    4001                          first_np_or_next_t[i]=iter;
     4003                        for (i=nbtold;i<nbt;i++) first_np_or_next_t[i]=iter;
    40024004
    40034005                        Headt = nbt; // empty list
    4004                         for (i=nbvold;i<nbv;i++)
    4005                           { // for all the triangle contening the vertex i
    4006                                 Vertex * s  = vertices + i;
     4006
     4007                        // for all the triangle contening the vertex i
     4008                        for (i=nbvold;i<nbv;i++){
     4009                                Vertex*          s  = vertices + i;
    40074010                                TriangleAdjacent ta(s->t, EdgesVertexTriangle[s->vint][1]);
    4008                                 Triangle * tbegin= (Triangle*) ta;
     4011                                Triangle*        tbegin= (Triangle*) ta;
    40094012                                Int4 kt;
    40104013                                do {
    40114014                                        kt = Number((Triangle*) ta);
    4012                                         if (first_np_or_next_t[kt]>0)
    4013                                          first_np_or_next_t[kt]=-Headt,Headt=kt;
     4015                                        if (first_np_or_next_t[kt]>0){
     4016                                                first_np_or_next_t[kt]=-Headt;
     4017                                                Headt=kt;
     4018                                        }
    40144019                                        if (ta.EdgeVertex(0)!=s){
    40154020                                                throw ErrorException(__FUNCT__,exprintf("ta.EdgeVertex(0)!=s"));
    … …  
    40174022                                        ta = Next(Adj(ta));
    40184023                                } while ( (tbegin != (Triangle*) ta));
    4019                           }   
     4024                        }   
    40204025
    40214026                } while (nbv!=nbvold);
Note: See TracChangeset for help on using the changeset viewer.