Implement discrete differentiation on images.
Started to evaluate Jacobi. Needs some cleanup first...