Mercurial > fem-fenics-eugenio

changeset 107:0e83f8ce9083

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use the function SubSpace to extract thesubspace and apply correctly the BCs.
author gedeone-octave <marcovass89@hotmail.it>
date Tue, 20 Aug 2013 22:31:16 +0200
parents 6b4a77dc734f
children 5cbc7341ded5
files example/Mixed-Poisson/MixedPoisson.m
diffstat 1 files changed, 7 insertions(+), 6 deletions(-) [+]
line wrap: on
line diff
--- a/example/Mixed-Poisson/MixedPoisson.m	Tue Aug 20 22:23:12 2013 +0200
+++ b/example/Mixed-Poisson/MixedPoisson.m	Tue Aug 20 22:31:16 2013 +0200
@@ -1,5 +1,6 @@
 pkg load msh
 pkg load fem-fenics
+
 fem_init_env ();
 
 problem = "MixedPoisson";
@@ -9,9 +10,9 @@
 mshd = Mesh (msho);
 
 V = MixedPoisson_FunctionSpace (mshd);
-
-bc1 = DirichletBC (V, @(x,y) [0; -sin(5.0*x); 0], 1);
-bc2 = DirichletBC (V, @(x,y) [0;  sin(5.0*x); 0], 3);
+V0 = SubSpace (V, 0);
+bc1 = DirichletBC (V0, @(x,y) [0; -sin(5.0*x)], 1);
+bc2 = DirichletBC (V0, @(x,y) [0;  sin(5.0*x)], 3);
 
 f = Expression ('f', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02));
 
@@ -22,9 +23,9 @@
 
 u = A \ b;
 
-func = Function ('u', V, u);
-uu = Function ('u', func, 0);
-sigma = Function ('sigma', func, 1);
+func = Function ('solution', V, u);
+sigma = Function ('sigma', func, 0);
+uu = Function ('u', func, 1);
 
 fem_plot (sigma);
 fem_plot (uu);