Changeset 13496


Ignore:
Timestamp:
10/01/12 12:09:25 (12 years ago)
Author:
jschierm
Message:

CHG: cosmetic.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/mesh/bamg.py

    r13495 r13496  
    159159#                                                       x=det([det([x1 y1; x2 y2])  x1-x2;det([x3 y3; x4 y4])  x3-x4])/det([x1-x2 y1-y2;x3-x4 y3-y4]);
    160160#                                                       y=det([det([x1 y1; x2 y2])  y1-y2;det([x3 y3; x4 y4])  y3-y4])/det([x1-x2 y1-y2;x3-x4 y3-y4]);
    161                                                         x=numpy.linalg.det(numpy.array([[numpy.linalg.det(numpy.array([[x1,y1],[x2,y2]])),x1-x2],[numpy.linalg.det(numpy.array([[x3,y3],[x4,y4]])),x3-x4]]))/numpy.linalg.det(numpy.array([[x1-x2,y1-y2],[x3-x4,y3-y4]]));
    162                                                         y=numpy.linalg.det(numpy.array([[numpy.linalg.det(numpy.array([[x1,y1],[x2,y2]])),y1-y2],[numpy.linalg.det(numpy.array([[x3,y3],[x4,y4]])),y3-y4]]))/numpy.linalg.det(numpy.array([[x1-x2,y1-y2],[x3-x4,y3-y4]]));
     161                                                        x=numpy.linalg.det(numpy.array([[numpy.linalg.det(numpy.array([[x1,y1],[x2,y2]])),x1-x2],[numpy.linalg.det(numpy.array([[x3,y3],[x4,y4]])),x3-x4]]))/numpy.linalg.det(numpy.array([[x1-x2,y1-y2],[x3-x4,y3-y4]]))
     162                                                        y=numpy.linalg.det(numpy.array([[numpy.linalg.det(numpy.array([[x1,y1],[x2,y2]])),y1-y2],[numpy.linalg.det(numpy.array([[x3,y3],[x4,y4]])),y3-y4]]))/numpy.linalg.det(numpy.array([[x1-x2,y1-y2],[x3-x4,y3-y4]]))
    163163
    164164                                                        segdis= sqrt((x4-x3)**2+(y4-y3)**2)
     
    205205                                                                        numpy.array([[count,count+1,2]]),\
    206206                                                                        numpy.hstack((numpy.arange(count+1,count+nods).reshape(-1,1),numpy.arange(count+2,count+nods+1).reshape(-1,1),(1+i)*numpy.ones((nods-1,1))))))
    207                                                                 count=count+nods;
     207                                                                count+=nods
    208208
    209209                                                                break
     
    327327        md.mesh.dimension=2
    328328        md.mesh.numberofelements=numpy.size(md.mesh.elements,axis=0)
    329         md.mesh.numberofvertices=numpy.size(md.mesh.x);
     329        md.mesh.numberofvertices=numpy.size(md.mesh.x)
    330330        md.mesh.numberofedges=numpy.size(md.mesh.edges,axis=0)
    331331        md.mesh.z=numpy.zeros(md.mesh.numberofvertices)
Note: See TracChangeset for help on using the changeset viewer.