sm/smbinterp
1
0
Fork 0
Code Issues Pull Requests Activity

minor, updated the driver to handle linear case

Browse Source
This commit is contained in:
Stephen Mardson McQuay 2010-01-31 21:15:47 -07:00
parent 0f4daa9a9c
commit f9973e7580
1 changed files with 4 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
bin/driver.py
View File

@ -100,14 +100,13 @@ if __name__ == '__main__':
S = [mesh_source.points[i] for i in indicies[3:] ] S = [mesh_source.points[i] for i in indicies[3:] ]
Sq = [mesh_source.q[i] for i in indicies[3:] ] Sq = [mesh_source.q[i] for i in indicies[3:] ]
s_mesh = grid.grid(S, Sq) s_mesh = grid.grid(S, Sq)
print len(S)
answer = baker.run_baker(X, r_mesh, s_mesh, options.verbose) answer = baker.run_baker(X, r_mesh, s_mesh, options.verbose)
print "]\n[".join([str(i) for i in [answer, r_mesh, s_mesh]]) if answer['a'] == None:
print errors.append(0)
print continue
print
exact = exact_func(X[0], X[1]) exact = exact_func(X[0], X[1])
if np.abs(exact - answer['final']) < np.abs(exact - answer['qlin']): if np.abs(exact - answer['final']) < np.abs(exact - answer['qlin']):
success += 1 success += 1