016 - Added the gradient calculation

Now the qpm algorithm returns also the xy gradient of the calculated
phase