Ignore:
Timestamp:
08/10/10 13:39:21 (15 years ago)
Author:
Mathieu Morlighem
Message:

improved Bamg objects

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Bamg/ListofIntersectionTriangles.cpp

    r5095 r5120  
    4444                                double ba,bb;
    4545                                TriangleAdjacent edge=CloseBoundaryEdge(a,t,ba,bb);
    46                                 MeshVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);
     46                                BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);
    4747                                NewItem(A,Metric(ba,v0,bb,v1));
    4848                                t=edge;
     
    5050                                if (det(v0.i,v1.i,b)>=0) {
    5151                                        TriangleAdjacent edge=CloseBoundaryEdge(a,t,ba,bb);
    52                                         MeshVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);
     52                                        BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);
    5353                                        NewItem(A,Metric(ba,v0,bb,v1));
    5454                                        return;
     
    8080                                        long int verbose=2;
    8181                                        TriangleAdjacent edge=CloseBoundaryEdge(a,t,ba,bb);
    82                                         MeshVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);
     82                                        BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);
    8383                                        NewItem(A,Metric(ba,v0,bb,v1));
    8484                                        return;
     
    232232                        lIntTria[Size].x = x;
    233233                        Metric m0,m1,m2;
    234                         register MeshVertex * v;
     234                        register BamgVertex * v;
    235235                        if ((v=(*tt)(0))) m0    = v->m;
    236236                        if ((v=(*tt)(1))) m1    = v->m;
     
    307307        /*}}}1*/
    308308        /*FUNCTION ListofIntersectionTriangles::NewPoints{{{1*/
    309         long ListofIntersectionTriangles::NewPoints(MeshVertex* vertices,long &nbv,long nbvx){
     309        long ListofIntersectionTriangles::NewPoints(BamgVertex* vertices,long &nbv,long maxnbv){
    310310                /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/NewPoints)*/
    311311
     
    354354
    355355                        si += sint;
    356                         if ( nbv<nbvx) {
     356                        if ( nbv<maxnbv) {
    357357                                vertices[nbv].r = C;
    358358                                vertices[nbv++].m = Metric(cx,lIntTria[ii-1].m,cy,lIntTria[ii].m);
Note: See TracChangeset for help on using the changeset viewer.