Is there any problem in adding a function Vector& Function::vector() ? It would be useful for the nonlinear solvers. Garth