Mercurial > fem-fenics-eugenio

changeset 100:5c5ef0f7498d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
save the solution instead of plotting it
author gedeone-octave <marcovass89@hotmail.it>
date Mon, 12 Aug 2013 22:13:37 +0200
parents 9818e6302e7b
children 6d05e6c85724
files example/NavierStokes/NS.m
diffstat 1 files changed, 5 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/example/NavierStokes/NS.m	Mon Aug 12 20:47:50 2013 +0200
+++ b/example/NavierStokes/NS.m	Mon Aug 12 22:13:37 2013 +0200
@@ -28,10 +28,11 @@
 a1 = TentativeVelocity_BilinearForm (V, k);
 A1 = assemble (a1, bcnoslip);
 a3 = VelocityUpdate_BilinearForm (V);
-A3 = assemble (a3, bcnoslip)
+A3 = assemble (a3, bcnoslip);
 
 # solve the problem for each time step
 # We need to update only the lhs
+i = 0;
 while t < T
   t += dt;
   bcin = DirichletBC (Q, @(x,y) sin(3.0*t), 1);
@@ -56,6 +57,6 @@
   ut = A3 \ b3;
   u0 = Function ('u0', V, ut);
 
-  #press Q to make the plot continue
-  fem_plot (u0);
-end  
+  name = sprintf ("u_%3.3d", ++i);
+  fem_save (u0, name);
+end