Example image comparison between NFWp and NFWm ellipticity parameterizations (top), with resulting residuals from the PEMD+shear fits (bottom). Caustics for both parameterizations are shown in the bottom right with the source position rescaled as described in the text and indicated as a cross. This example image comes from a system with moderately large ellipticity with an input axis ratio of q_κ = 0.64.