Untested version of ImageAlloc::applyCD

This function takes four matrices with the shift parameters, (2dmax+1)^2 dimensional
if we go out to a maximum separation of dmax pixels, applies the charge shifts to the
image and returns the result.
Merge branch '#524' of github.com:GalSim-developers/GalSim into #524