Now loads b/-rho (self.b)
Added dot_product function
Changed subtraction in some kernels to addition
Now runs to completion, but output does not appear to be correct.
Debugging needed.