linear reconstruction no longer needed for odd resolution

This commit is contained in:
Martin Diehl 2013-06-05 14:05:11 +00:00
parent ac870574ab
commit becacb675a
1 changed files with 2 additions and 7 deletions

View File

@ -463,9 +463,7 @@ for filename in args:
res = numpy.array([len(grid[0]),\
len(grid[1]),\
len(grid[2]),],'i')
if (res[0]%2 != 0 or res[1]%2 != 0 or (res[2] != 1 and res[2]%2 !=0)):
print 'using linear reconstruction for uneven resolution'
options.linearreconstruction = True
dim = numpy.ones(3)
for i,r in enumerate(res):
@ -505,12 +503,9 @@ for filename in args:
tablenodalcoords.labels.index(options.labelnodalcoords[2])]
i = 0
while tablenodalcoords.data_read(): # read next data line of ASCII table
mesh[i,:]=float(tablenodalcoords.data[column[:]])
i += 1
mesh=mesh.reshape(res[0]+1,res[1]+1,res[2]+1,3)
fields = {\