Mercurial > fem-fenics-eugenio

changeset 77:ff95326e6f13

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make the example compliant with the fem-fenics pkg
author gedeone-octave <marcovass89@hotmail.it>
date Fri, 02 Aug 2013 23:18:01 +0200
parents 6cb2a3e3be68
children 670a5d91c397
files example/MixedPoisson.m example/NavierStokes/NS.m
diffstat 2 files changed, 8 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/example/MixedPoisson.m	Fri Aug 02 23:01:39 2013 +0200
+++ b/example/MixedPoisson.m	Fri Aug 02 23:18:01 2013 +0200
@@ -1,7 +1,9 @@
 pkg load msh
-
+pkg load fem-fenics
 fem_init_env ();
 
+problem = "MixedPoisson";
+fem_create_all (problem);
 x = y = linspace (0, 1, 32);
 msho = msh2m_structured_mesh (x, y, 1, 1:4);
 mshd = fem_init_mesh (msho);
@@ -11,7 +13,6 @@
 bc1 = fem_bc (V, @(x,y) [0; -sin(5.0*x); 0], 1);
 bc2 = fem_bc (V, @(x,y) [0;  sin(5.0*x); 0], 3);
 
-
 f = fem_coeff ('f', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02));
 
 A = fem_rhs_MixedPoisson (V, bc1, bc2);
@@ -19,7 +20,6 @@
 
 u = A \ b;
 
-
 func = fem_func ('u', V, u);
 
 fem_plot (func);
--- a/example/NavierStokes/NS.m	Fri Aug 02 23:01:39 2013 +0200
+++ b/example/NavierStokes/NS.m	Fri Aug 02 23:18:01 2013 +0200
@@ -1,7 +1,12 @@
 pkg load msh
+pkg load fem-fenics
 
 fem_init_env ();
 
+fem_create_all ("TentativeVelocity");
+fem_create_all ("VelocityUpdate");
+fem_create_all ("PressureUpdate");
+
 run lshape-domain.m;
 mshd = fem_init_mesh (msho);