/[escript]/trunk/pycad/py_src/Triangle.py
ViewVC logotype

Diff of /trunk/pycad/py_src/Triangle.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3773 by caltinay, Fri Dec 9 04:48:08 2011 UTC revision 3774 by jfenwick, Wed Jan 18 06:29:34 2012 UTC
# Line 133  class Design(design.Design): Line 133  class Design(design.Design):
133              see U{http://www.cs.cmu.edu/~quake/triangle.switch.html}              see U{http://www.cs.cmu.edu/~quake/triangle.switch.html}
134          """          """
135          if self.__cmdLineArgs == "":          if self.__cmdLineArgs == "":
136              print "warning: using default command line arguments for Triangle"              print("warning: using default command line arguments for Triangle")
137          exe="triangle %s %%s"%self.__cmdLineArgs          exe="triangle %s %%s"%self.__cmdLineArgs
138          return exe          return exe
139    
# Line 152  class Design(design.Design): Line 152  class Design(design.Design):
152              ret=0              ret=0
153          ret=getMPIWorldMax(ret)          ret=getMPIWorldMax(ret)
154          if ret > 0:          if ret > 0:
155            raise RuntimeError, "Could not build mesh: %s"%" ".join(args)            raise RuntimeError("Could not build mesh: %s"%" ".join(args))
156          else:          else:
157              # <hack> so that users can set the mesh filename they want.              # <hack> so that users can set the mesh filename they want.
158              name=self.getScriptFileName()              name=self.getScriptFileName()
# Line 197  class Design(design.Design): Line 197  class Design(design.Design):
197                         pts=list(PSp.getControlPoints())                         pts=list(PSp.getControlPoints())
198                         for pt in pts:                         for pt in pts:
199                             c=pt.getCoordinates()                             c=pt.getCoordinates()
200                             if pt not in ctrlPts.keys():                             if pt not in list(ctrlPts.keys()):
201                                 vertCnt+=1                                 vertCnt+=1
202                                 vertices+="%s %s %s %s\n"%(vertCnt,c[0],c[1],p.getID())                                 vertices+="%s %s %s %s\n"%(vertCnt,c[0],c[1],p.getID())
203                                 ctrlPts[pt]=vertCnt                                 ctrlPts[pt]=vertCnt
# Line 222  class Design(design.Design): Line 222  class Design(design.Design):
222                             pts=list(curve.getControlPoints())                             pts=list(curve.getControlPoints())
223                             for pt in pts:                             for pt in pts:
224                                 c=pt.getCoordinates()                                 c=pt.getCoordinates()
225                                 if pt not in ctrlPts.keys():                                 if pt not in list(ctrlPts.keys()):
226                                     vertCnt+=1                                     vertCnt+=1
227                                     vertices+="%s %s %s %s\n"%(vertCnt,c[0],c[1],p.getID())                                     vertices+="%s %s %s %s\n"%(vertCnt,c[0],c[1],p.getID())
228                                     ctrlPts[pt]=vertCnt                                     ctrlPts[pt]=vertCnt
# Line 246  class Design(design.Design): Line 246  class Design(design.Design):
246                                 pts=list(curve.getControlPoints())                                 pts=list(curve.getControlPoints())
247                                 for pt in pts:                                 for pt in pts:
248                                     c=pt.getCoordinates()                                     c=pt.getCoordinates()
249                                     if pt not in ctrlPts.keys():                                     if pt not in list(ctrlPts.keys()):
250                                         vertCnt+=1                                         vertCnt+=1
251                                         vertices+="%s %s %s %s\n"%(vertCnt,c[0],c[1],p.getID())                                         vertices+="%s %s %s %s\n"%(vertCnt,c[0],c[1],p.getID())
252                                         ctrlPts[pt]=vertCnt                                         ctrlPts[pt]=vertCnt
# Line 276  class Design(design.Design): Line 276  class Design(design.Design):
276                                 vectors[i].append(self.__getVector(A,B))                                 vectors[i].append(self.__getVector(A,B))
277                                 vectors[i].append(self.__getVector(A,C))                                 vectors[i].append(self.__getVector(A,C))
278                             # get angle between vectors at each vertex                             # get angle between vectors at each vertex
279                             for i in vectors.keys():                             for i in list(vectors.keys()):
280                                 angle=self.__getAngle(vectors[i][0],vectors[i][1])                                 angle=self.__getAngle(vectors[i][0],vectors[i][1])
281                                 vectors[i].append(angle)                                 vectors[i].append(angle)
282                             # find the vertex with the smallest angle                             # find the vertex with the smallest angle
283                             minAngle=360.                             minAngle=360.
284                             indx=0                             indx=0
285                             for i in vectors.keys():                             for i in list(vectors.keys()):
286                                 if vectors[i][2] < minAngle:                                 if vectors[i][2] < minAngle:
287                                     indx=i                                     indx=i
288                                     minAngle=vectors[i][2]                                     minAngle=vectors[i][2]

Legend:
Removed from v.3773  
changed lines
  Added in v.3774

  ViewVC Help
Powered by ViewVC 1.1.26