Ignore:
Timestamp:
09/09/11 14:58:09 (14 years ago)
Author:
seroussi
Message:

keep building mesh

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified issm/trunk/src/m/model/mesh/findsegments.m

    r9729 r9733  
    1919                error(' ''mesh.elementconnectivity'' option does not have thge right size.');
    2020        else
    21                 mesh.elementconnectivity=ElementConnectivity(md.elements,md.mesh.vertexconnectivity);
     21                mesh.elementconnectivity=ElementConnectivity(md.mesh.elements,md.mesh.vertexconnectivity);
    2222        end
    2323end
     
    4343
    4444                %get nodes of el1
    45                 nods1=md.elements(el1,:);
     45                nods1=md.mesh.elements(el1,:);
    4646
    4747                %find the common vertices to the two elements connected to el1 (1 or 2)
    48                 flag=intersect(md.elements(els2(1),:),md.elements(els2(2),:));
     48                flag=intersect(md.mesh.elements(els2(1),:),md.mesh.elements(els2(2),:));
    4949
    5050                %get the vertices on the boundary and build segment
     
    5353
    5454                %swap segment nodes if necessary
    55                 ord1=find(nods1(1)==md.elements(el1,:));
    56                 ord2=find(nods1(2)==md.elements(el1,:));
     55                ord1=find(nods1(1)==md.mesh.elements(el1,:));
     56                ord2=find(nods1(2)==md.mesh.elements(el1,:));
    5757                if ( (ord1==1 & ord2==2) | (ord1==2 & ord2==3) | (ord1==3 & ord2==1) ),
    5858                        temp=segments(count,1);
     
    6666        else
    6767                %get nodes of el1
    68                 nods1=md.elements(el1,:);
     68                nods1=md.mesh.elements(el1,:);
    6969
    7070                %find the vertex  the el1 to not share with els2
    71                 flag=setdiff(nods1,md.elements(els2,:));
     71                flag=setdiff(nods1,md.mesh.elements(els2,:));
    7272
    7373                for j=1:3,
     
    7878
    7979                                %swap segment nodes if necessary
    80                                 ord1=find(nods(1)==md.elements(el1,:));
    81                                 ord2=find(nods(2)==md.elements(el1,:));
     80                                ord1=find(nods(1)==md.mesh.elements(el1,:));
     81                                ord2=find(nods(2)==md.mesh.elements(el1,:));
    8282                                if ( (ord1==1 & ord2==2) | (ord1==2 & ord2==3) | (ord1==3 & ord2==1) ),
    8383                                        temp=segments(count,1);
Note: See TracChangeset for help on using the changeset viewer.