# HG changeset patch # User gedeone-octave # Date 1373575572 -7200 # Node ID c6ede799884c4c49edaf5202c53f02197a7b6e5a # Parent fe2264f1e789ec22d66a789e53f190bb0f045d0e Now include also call to fem_fs diff -r fe2264f1e789 -r c6ede799884c test/laplace.m --- a/test/laplace.m Thu Jul 11 22:44:38 2013 +0200 +++ b/test/laplace.m Thu Jul 11 22:46:12 2013 +0200 @@ -3,18 +3,18 @@ # Create mesh fem_init_env (); +pkg load msh; -x = y = linspace (0, 1, 32) ; -msho = msh2m_structured_mesh (x, y, 1, [ 1 1 2 2 ]); +msho = msh2m_structured_mesh (0:0.05:1, 0:0.05:1, 1, 1:4); mshd = fem_init_mesh (msho); #Compile the problem defined in the .ufl file -fem_ffc ("Laplace.ufl"); +#fem_ffc ("Laplace.ufl"); # #Import the problem inside Octave -# V = FunctionSpace (mshd, "Laplace"); -# sd = SubDomain (mshd, @(x, y) abs (x - 1.0) < eps, "Laplace"); -# bc = DirichletBC (V, @(x, y) 0.0, sd, "Laplace"); +V = fem_fs (mshd); +#Set BC equal to zero on sides with label 2 and 4 +#bc = fem_bc (V, @(x, y) 0.0, [2,4]); # #Get the matrices # A = BilinearForm (V, V, bc, "Laplace");