Something is wrong with sub-functions in Python. For the mixed-poisson demo, doing
(sigma, u) = problem.solve().split() plot(u) the plot looks as expected. Doing U = problem.solve() plot(U[1]) the plot is wrong. Also, it's not possible to do f = File("u.pvd") f << U[1] I get the message TypeError: in method 'File___lshift__', argument 2 of type 'dolfin::FunctionPlotData const &' Garth _______________________________________________ DOLFIN-dev mailing list DOLFIN-dev@fenics.org http://www.fenics.org/mailman/listinfo/dolfin-dev