Skip to content
Snippets Groups Projects
Commit 19827488 authored by Franck Pérignon's avatar Franck Pérignon
Browse files

fix in demo_2D

parent d54fe2d7
No related branches found
No related tags found
No related merge requests found
......@@ -70,14 +70,17 @@ def run(nb=257):
## Setting solver to Problem
pb.setUp(finalTime, timeStep)
scalar_initial = np.copy(scal.discreteField[0].data)
## We copy the first discretisation of scalar
scalar_initial = np.copy(scal.discreteFields.keys()[0].data)
## Solve problem
pb.solve()
scal.discreteField[0].toHost()
print 'Erreur : ', np.max(np.abs(scalar_initial - scal.discreteField[0].data))/np.max(scalar_initial)
print np.linalg.norm(scalar_initial - scal.discreteField[0].data, ord=2)/np.linalg.norm(scalar_initial, ord=2)
scal_disc = scal.discreteFields.keys()[0]
scal_disc.toHost()
print 'Erreur : ', np.max(np.abs(scalar_initial -
scal_disc.data)) / np.max(scalar_initial)
print np.linalg.norm(scalar_initial - scal_disc.data, ord=2)/np.linalg.norm(scalar_initial, ord=2)
pb.finalize()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment