My original test was to have the successive `applyShear` calls in the python script and the shear addition in the yaml.  But it failed.  So I had to switch it to this version.

The problem is that two shears do not really add to become a single net shear.  There is also a slight rotation.  The Miralda Escude formula, for example, is just the shear part, ignoring the rotation.

As I mentioned here, it is normally fine, since we cannot observe the original galaxy orientation, so we don't care about the rotation part.  But it means that unless you start with a circular source galaxy, the resulting image is different in the two cases.