The quadrupole and hexadecapole corrections to the point source magnifications are insignificant when approach a fold from the the outside of a caustic. Bozza et.al. 2018 devised a test (Eq. 47 in the paper) to catch when the source is close to a fold and trigger the full calculation at those points. I've implemented this test but the output looks like this:
No idea what these weird curves emerging out of the fold near the top are. Could be a mistake in my implementation.
The quadrupole and hexadecapole corrections to the point source magnifications are insignificant when approach a fold from the the outside of a caustic. Bozza et.al. 2018 devised a test (Eq. 47 in the paper) to catch when the source is close to a fold and trigger the full calculation at those points. I've implemented this test but the output looks like this:
No idea what these weird curves emerging out of the fold near the top are. Could be a mistake in my implementation.