Sorry it’s me again.
I want to ask some questions regarding the class Function. It seems like this class is quite different between the old DILFIN, and many functions are deleted. As before, the users can do operation on the DOFs of a function easily, for example, we can do :
f->vector()->set_local(f_vec, dofmap,size_dofmap)
or
f->vector()->axpy(a,g)
or
f->vector() = g->vector()
Now, it seems to be more difficult, especially in CPP interface. We can access the DOFs of a function using f->x(), which gives la::vector. If I want to do some operatons on whole vectors of functions, say if I have function f and g, I want to add them up or assign g to f, I have tried:
f->x()->mutable_array() = g->x()->mutable_array(),
but it doesn’t work. Is there any suggestion? Thanks!