Advection and rendering velocity lines now seem to work. Also adding boundary velocity and the linear interpolation is moved to a inlined function. Next is the build solution kernel (divergence).