It's a bit cumbersome to pass around mesh, dofmap, and form to be able to create functions in finite element spaces. Maybe a FunctionSpace class could collect the necessary info?
# Instead of: uvec = Vector() u = cpp_Function(mesh, uvec, dofmap, form, 0) vvec = Vector() v = cpp_Function(mesh, vvec, dofmap, form, 0) # It would be nice to write: U = FunctionSpace(mesh, dofmapset, form, 0) # or something similar u = Function(U) v = Function(U) This concept is probably usable in C++ as well. (Note that I'm not volunteering to implement this myself now, don't have the time). -- Martin _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
