Mercurial > fem-fenics-eugenio

changeset 103:2d90592f73ab

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
MAint: move the file to the new specific folder
author gedeone-octave <marcovass89@hotmail.it>
date Wed, 14 Aug 2013 11:56:58 +0200
parents d880c19f691b
children 68262f409ac4
files example/Mixed-Poisson/MixedPoisson.m example/Mixed-Poisson/MixedPoisson.ufl example/MixedPoisson.m example/MixedPoisson.ufl
diffstat 4 files changed, 42 insertions(+), 42 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/example/Mixed-Poisson/MixedPoisson.m	Wed Aug 14 11:56:58 2013 +0200
@@ -0,0 +1,30 @@
+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 = 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);
+
+f = Expression ('f', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02));
+
+a = MixedPoisson_BilinearForm (V);
+L = MixedPoisson_LinearForm (V, f);
+
+[A, b] = assemble_system (a, L, bc1, bc2);
+
+u = A \ b;
+
+func = Function ('u', V, u);
+uu = Function ('u', func, 0);
+sigma = Function ('sigma', func, 1);
+
+fem_plot (sigma);
+fem_plot (uu);
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/example/Mixed-Poisson/MixedPoisson.ufl	Wed Aug 14 11:56:58 2013 +0200
@@ -0,0 +1,12 @@
+BDM = FiniteElement("BDM", triangle, 1)
+DG  = FiniteElement("DG", triangle, 0)
+W = BDM * DG
+
+(sigma, u) = TrialFunctions(W)
+(tau, v)   = TestFunctions(W)
+
+CG = FiniteElement("CG", triangle, 1)
+f = Coefficient(CG)
+
+a = (dot(sigma, tau) + div(tau)*u + div(sigma)*v)*dx
+L = - f*v*dx
--- a/example/MixedPoisson.m	Wed Aug 14 11:56:36 2013 +0200
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,30 +0,0 @@
-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 = 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);
-
-f = Expression ('f', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02));
-
-a = MixedPoisson_BilinearForm (V);
-L = MixedPoisson_LinearForm (V, f);
-
-[A, b] = assemble_system (a, L, bc1, bc2);
-
-u = A \ b;
-
-func = Function ('u', V, u);
-uu = Function ('u', func, 0);
-sigma = Function ('sigma', func, 1);
-
-fem_plot (sigma);
-fem_plot (uu);
--- a/example/MixedPoisson.ufl	Wed Aug 14 11:56:36 2013 +0200
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,12 +0,0 @@
-BDM = FiniteElement("BDM", triangle, 1)
-DG  = FiniteElement("DG", triangle, 0)
-W = BDM * DG
-
-(sigma, u) = TrialFunctions(W)
-(tau, v)   = TestFunctions(W)
-
-CG = FiniteElement("CG", triangle, 1)
-f = Coefficient(CG)
-
-a = (dot(sigma, tau) + div(tau)*u + div(sigma)*v)*dx
-L = - f*v*dx