Changeset 22213


Ignore:
Timestamp:
11/02/17 14:28:36 (7 years ago)
Author:
Mathieu Morlighem
Message:

CHG: cosmetics

File:
1 edited

Legend:

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

    r21931 r22213  
    11341134
    11351135        }/*}}}*/
    1136         void  Mesh::BoundAnisotropy(double anisomax,double hminaniso) {/*{{{*/
     1136        void Mesh::BoundAnisotropy(double anisomax,double hminaniso) {/*{{{*/
    11371137                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/BoundAnisotropy)*/
    11381138
     
    23352335                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/FindSubDomain)*/
    23362336
    2337                 long int verbose=0;
     2337                int verbose=0;
    23382338
    23392339                if (verbose >2){
    23402340                        if (OutSide) _printf_("   Find all external sub-domain\n");
    23412341                        else _printf_("   Find all internal sub-domain\n");
    2342                   }
     2342                }
    23432343                short * HeapArete = new short[nbt];
    23442344                Triangle  **  HeapTriangle = new Triangle*  [nbt];
     
    23462346                long k,it;
    23472347
    2348                 for (int itt=0;itt<nbt;itt++)
    2349                  triangles[itt].link=0; // par defaut pas de couleur
     2348                /*No color by default*/
     2349                for(int itt=0;itt<nbt;itt++) triangles[itt].link=0;
    23502350
    23512351                long  NbSubDomTot =0;
    2352                 for ( it=0;it<nbt;it++)  {
    2353                         if ( ! triangles[it].link  ) {
     2352                for(it=0;it<nbt;it++)  {
     2353                        if( !triangles[it].link ){
    23542354                                t = triangles + it;
    23552355                                NbSubDomTot++;; // new composante connexe
     
    23602360                                HeapArete[i] = 3;
    23612361
    2362                                 while (i >= 0) // boucle sur la pile
    2363                                   { while ( HeapArete[i]--) // boucle sur les 3 aretes
     2362                                while(i >= 0) // boucle sur la pile
     2363                                  { while (HeapArete[i]--) // boucle sur les 3 aretes
    23642364                                          {
    23652365                                                int na =  HeapArete[i];
     
    24452445                                                        _error_("k!=nbsubdomains");
    24462446                                                }
    2447                                                 if(OutSide)
    2448                                                   {
     2447                                                if(OutSide){
    24492448                                                        //  to remove all the sub domain by parity adjacents
    24502449                                                        //  because in this case we have only the true boundary edge
    2451                                                         // so teh boundary is manifold
     2450                                                        //  so the boundary is manifold
    24522451                                                        long nbk = nbsubdomains;
    24532452                                                        while (nbk)
     
    30843083                                        Triangle* ta=tadjj;
    30853084
    3086                                         //if the adjacent triangle is a boundary triangle, continur
     3085                                        //if the adjacent triangle is a boundary triangle, continue
    30873086                                        if (ta->det<0) continue;         
    30883087
Note: See TracChangeset for help on using the changeset viewer.