Closer to a working Neumann test
vec: Add axpy
Neumann Laplace works