Mercurial > fem-fenics-eugenio

changeset 164:7190852cf57f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update the examples giving two input parameters to the BilinearForm.
author gedeone-octave <marcovass89@hotmail.it>
date Sat, 05 Oct 2013 10:40:06 +0100
parents 5f22232b9fcc
children 3ae8a7398f43
files inst/example/Biharmonic/Biharmonic.m inst/example/Evolution/Evolution.m inst/example/Mixed-Poisson/MixedPoisson.m inst/example/NavierStokes/NS.m inst/example/Poisson/Poisson.m
diffstat 5 files changed, 7 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
--- a/inst/example/Biharmonic/Biharmonic.m	Sat Oct 05 10:30:39 2013 +0100
+++ b/inst/example/Biharmonic/Biharmonic.m	Sat Oct 05 10:40:06 2013 +0100
@@ -14,7 +14,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);
--- a/inst/example/Evolution/Evolution.m	Sat Oct 05 10:30:39 2013 +0100
+++ b/inst/example/Evolution/Evolution.m	Sat Oct 05 10:40:06 2013 +0100
@@ -17,7 +17,7 @@
 
 u0 = Expression ('u0', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02));
 
-a = BilinearForm (problem, V, k);
+a = BilinearForm (problem, V, V, k);
 L = LinearForm (problem, V, k, u0);
 
 A = assemble (a, bc);
--- a/inst/example/Mixed-Poisson/MixedPoisson.m	Sat Oct 05 10:30:39 2013 +0100
+++ b/inst/example/Mixed-Poisson/MixedPoisson.m	Sat Oct 05 10:40:06 2013 +0100
@@ -9,7 +9,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
--- a/inst/example/NavierStokes/NS.m	Sat Oct 05 10:30:39 2013 +0100
+++ b/inst/example/NavierStokes/NS.m	Sat Oct 05 10:40:06 2013 +0100
@@ -46,13 +46,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);
--- a/inst/example/Poisson/Poisson.m	Sat Oct 05 10:30:39 2013 +0100
+++ b/inst/example/Poisson/Poisson.m	Sat Oct 05 10:40:06 2013 +0100
@@ -13,7 +13,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