/[escript]/trunk/finley/src/generateReferenceElementList.py
ViewVC logotype

Diff of /trunk/finley/src/generateReferenceElementList.py

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

revision 4575 by caltinay, Tue Jul 2 01:44:11 2013 UTC revision 4576 by sshaw, Mon Dec 9 23:35:30 2013 UTC
# Line 173  for name in ["Point1", "Line2", "Line3", Line 173  for name in ["Point1", "Line2", "Line3",
173          else:          else:
174              numSubElements=1              numSubElements=1
175              BasisFunctions=Parametrization              BasisFunctions=Parametrization
176              subElementNodes=[ i for i in xrange(numNodes) ]              subElementNodes=[ i for i in range(numNodes) ]
177    
178          linearTypeId="%s%s"%(z,numLinearNodes)          linearTypeId="%s%s"%(z,numLinearNodes)
179          linearNodes=[ i for i in xrange(numLinearNodes) ]          linearNodes=[ i for i in range(numLinearNodes) ]
180          if isFace: linearTypeId+="Face"          if isFace: linearTypeId+="Face"
181                
182          if isContact:          if isContact:
183               linearTypeId+="_Contact"               linearTypeId+="_Contact"
184               linearNodes+=[ n+i for i in xrange(numLinearNodes) ]               linearNodes+=[ n+i for i in range(numLinearNodes) ]
185                    
186          if isFace:          if isFace:
187             relevantGeoNodes=RELEVANTGEO["%s%s"%(z,n)]             relevantGeoNodes=RELEVANTGEO["%s%s"%(z,n)]
188          else:          else:
189             relevantGeoNodes=[ i for i in xrange(n) ]             relevantGeoNodes=[ i for i in range(n) ]
190    
191          if isContact:          if isContact:
192               faceNodes =[-1]               faceNodes =[-1]
# Line 198  for name in ["Point1", "Line2", "Line3", Line 198  for name in ["Point1", "Line2", "Line3",
198                  shiftNodes = SHIFTNODES["%s%sFace"%(z,n)][0]                  shiftNodes = SHIFTNODES["%s%sFace"%(z,n)][0]
199                  reverseNodes = SHIFTNODES["%s%sFace"%(z,n)][1]                  reverseNodes = SHIFTNODES["%s%sFace"%(z,n)][1]
200             else:             else:
201                  faceNodes=[  i for i in xrange(numNodes) ]                  faceNodes=[  i for i in range(numNodes) ]
202                  shiftNodes = SHIFTNODES["%s%s"%(z,n)][0]                  shiftNodes = SHIFTNODES["%s%s"%(z,n)][0]
203                  reverseNodes = SHIFTNODES["%s%s"%(z,n)][1]                  reverseNodes = SHIFTNODES["%s%s"%(z,n)][1]
204    

Legend:
Removed from v.4575  
changed lines
  Added in v.4576

  ViewVC Help
Powered by ViewVC 1.1.26