dolfin team mailing list archive

["Martin Sandve Alnæs" <martinal@xxxxxxxxx>]

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