syfi team mailing list archive
-
syfi team
-
Mailing list archive
-
Message #01325
Re: [Branch ~syfi-core/fenics-syfi/main] Rev 521: created new mxpoission output and updated some demos
The mac buildbot now fails on the mxpoisson test :)
Johannes
On Tue, Aug 23, 2011 at 8:33 PM, <noreply@xxxxxxxxxxxxx> wrote:
> ------------------------------------------------------------
> revno: 521
> committer: Kent-Andre Mardal <kent-and@xxxxxxxxx>
> branch nick: fenics-syfi
> timestamp: Tue 2011-08-23 08:31:40 +0200
> message:
> created new mxpoission output and updated some demos
> modified:
> demo/Advection-Diffusion/python/demo.py
> demo/LinearElasticity/cpp/main.cpp
> demo/LinearElasticity/python/demo.py
> tests/cpp/mxpoisson_ex.gar.r
>
>
> --
> lp:fenics-syfi
> https://code.launchpad.net/~syfi-core/fenics-syfi/main
>
> You are subscribed to branch lp:fenics-syfi.
> To unsubscribe from this branch go to https://code.launchpad.net/~syfi-core/fenics-syfi/main/+edit-subscription
>
> === modified file 'demo/Advection-Diffusion/python/demo.py'
> --- demo/Advection-Diffusion/python/demo.py 2011-06-02 20:33:50 +0000
> +++ demo/Advection-Diffusion/python/demo.py 2011-08-23 06:31:40 +0000
> @@ -102,7 +102,7 @@
> # Compute sensitivity to c
> u0.assign(u1)
> sb = assemble(sL)
> - solve(A, su1.vector(), sb)
> + solve(A==sb, su1)
> sout_file << su1
> figure(1)
> plot(su1, title="du/dc")
>
> === modified file 'demo/LinearElasticity/cpp/main.cpp'
> --- demo/LinearElasticity/cpp/main.cpp 2011-06-02 20:33:50 +0000
> +++ demo/LinearElasticity/cpp/main.cpp 2011-08-23 06:31:40 +0000
> @@ -123,11 +123,9 @@
> LinearElasticity::LinearForm L(V);
> L.f = f;
>
> - VariationalProblem problem(a, L, bcs);
> -
> // Solve PDE (using direct solver)
> Function u(V);
> - problem.solve(u);
> + solve(a==L, u, bcs);
>
> // Plot solution
> //plot(u, "displacement");
>
> === modified file 'demo/LinearElasticity/python/demo.py'
> --- demo/LinearElasticity/python/demo.py 2011-06-02 20:33:50 +0000
> +++ demo/LinearElasticity/python/demo.py 2011-08-23 06:31:40 +0000
> @@ -93,15 +93,15 @@
> bcl = DirichletBC(V, c, Left())
>
> # Set up boundary condition at right end
> -r = Rotation(element=V.ufl_element())
> +r = Rotation()
> bcr = DirichletBC(V, r, Right())
>
> # Set up boundary conditions
> bcs = [bcl, bcr]
>
> # Set up PDE and solve
> -problem = VariationalProblem(a, L, bcs)
> -u = problem.solve()
> +u = Function(V)
> +solve(a==L, u, bcs)
>
> # Save solution to VTK format
> vtk_file = File("elasticity.pvd")
>
> === modified file 'tests/cpp/mxpoisson_ex.gar.r'
> Binary files tests/cpp/mxpoisson_ex.gar.r 2009-04-02 07:52:45 +0000 and tests/cpp/mxpoisson_ex.gar.r 2011-08-23 06:31:40 +0000 differ
>
>