Allow kwds in vector linsolve
second order neumann conditions