Isn't the code: ...................... Matrix A; Vector b,x; assemble(A,a,mesh); assemble(b,L,mesh); solver.solve(A,x,b); Function up(mesh,x,a); //up=(u,p) velocities and pressure Function u=up[0]; Function p=up[1]; ........................... supposed to work? this gives an error: ---------------------- conversion from 'dolfin::SubFunction' to non-scalar type 'dolfin::Function' requested --------------------
On the other hand I know that this works: Function u(up[0]); Function p(up[1]); -- Nuno David Lopes _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
