Note that there appears to be a systematic astrometry difference between MARX and mkpsf--the centroid of the two PSFs (before blurring) are generally offset from each other.
I find this offset is near zero on-axis and increases linearly with off-axis angle at a rate of
0.12 sky pixels per arc minute.
I do not, of course, know which PSF has the more correct astrometry.
One implication of this offset is that one should avoid the situation where a source's position was determined using a PSF from one
system (e.g. via AE's ``correlation position'' estimate or via image reconstruction), and the source's extraction region was constructed using a PSF from the other system.
In a multi-observation project, one should not extract some observations with one flavor of PSF and some with the other
unless you have considered very carefully how that would affect your analysis.