Mercurial > fem-fenics-eugenio
changeset 165:3ae8a7398f43
Update the examples for the tutorial.
author | gedeone-octave <marcovass89@hotmail.it> |
---|---|
date | Sat, 05 Oct 2013 22:25:52 +0100 |
parents | 7190852cf57f |
children | 40ae9e5dfb93 |
files | inst/femfenics_examples.m |
diffstat | 1 files changed, 6 insertions(+), 6 deletions(-) [+] |
line wrap: on
line diff
--- a/inst/femfenics_examples.m Sat Oct 05 10:40:06 2013 +0100 +++ b/inst/femfenics_examples.m Sat Oct 05 22:25:52 2013 +0100 @@ -79,7 +79,7 @@ %! f = Expression ('f', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02)); %! g = Expression ('g', @(x,y) sin (5.0 * x)); %! -%! a = BilinearForm ('Poisson', V); +%! a = BilinearForm ('Poisson', V, V); %! L = LinearForm ('Poisson', V, f, g); %! %! # Compute solution @@ -119,7 +119,7 @@ %! %! f = Expression ('f', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02)); %! -%! a = BilinearForm ('MixedPoisson', V); +%! a = BilinearForm ('MixedPoisson', V, V); %! L = LinearForm ('MixedPoisson', V, f); %! %! # Define essential boundary @@ -275,13 +275,13 @@ %! f = Constant ('f', [0; 0]); %! %! # Tentative velocity step. -%! a1 = BilinearForm ('TentativeVelocity', V, k); +%! a1 = BilinearForm ('TentativeVelocity', V, V, k); %! %! # Pressure update. -%! a2 = BilinearForm ('PressureUpdate', Q); +%! a2 = BilinearForm ('PressureUpdate', Q, Q); %! %! # Velocity update -%! a3 = BilinearForm ('VelocityUpdate', V); +%! a3 = BilinearForm ('VelocityUpdate', V, V); %! %! # Assemble matrices %! A1 = assemble (a1, noslip); @@ -358,7 +358,7 @@ %! f = Expression ('f', @(x,y) 4.0*pi^4.*sin(pi.*x).*sin(pi.*y)); %! g = Expression ('alpha', @(x,y) 8); %! -%! a = BilinearForm (problem, V, g); +%! a = BilinearForm (problem, V, V, g); %! L = LinearForm (problem, V, f); %! %! [A, b] = assemble_system (a, L, bc);