Mercurial > fem-fenics-eugenio
changeset 238:b96f6b12f8ca
Move ufl to m-file in Poisson example
author | Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
---|---|
date | Tue, 24 Jun 2014 16:43:13 +0200 |
parents | 418a5119047b |
children | cb6f11a09667 |
files | inst/example/Poisson/Poisson.m inst/example/Poisson/Poisson.ufl |
diffstat | 2 files changed, 12 insertions(+), 11 deletions(-) [+] |
line wrap: on
line diff
--- a/inst/example/Poisson/Poisson.m Mon Jun 23 20:28:28 2014 +0200 +++ b/inst/example/Poisson/Poisson.m Tue Jun 24 16:43:13 2014 +0200 @@ -14,7 +14,18 @@ ## this program; if not, see <http://www.gnu.org/licenses/>. pkg load fem-fenics msh -import_ufl_Problem ('Poisson') + +ufl start Poisson +ufl 'element = FiniteElement("Lagrange", triangle, 1)' +ufl +ufl u = TrialFunction(element) +ufl v = TestFunction(element) +ufl f = Coefficient(element) +ufl g = Coefficient(element) +ufl +ufl "a = inner(grad(u), grad(v))*dx" +ufl L = f*v*dx + g*v*ds +ufl end # Create mesh and define function space x = y = linspace (0, 1, 33);
--- a/inst/example/Poisson/Poisson.ufl Mon Jun 23 20:28:28 2014 +0200 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,10 +0,0 @@ -# Copyright (C) 2005-2009 Anders Logg -element = FiniteElement("Lagrange", triangle, 1) - -u = TrialFunction(element) -v = TestFunction(element) -f = Coefficient(element) -g = Coefficient(element) - -a = inner(grad(u), grad(v))*dx -L = f*v*dx + g*v*ds