Photon shooting with OpticalPSF/InterpolatedImage gives different results on Mac ARM and Linux

I noticed that I found different results depending if I run some code on Linux or on my laptop with Mac ARM. I tracedback the issue to photon shooting through OpticalPSF but at some point it gets transformed to an InterpolatedImage. I figure that the difference is at this line. I manage to figure that _xInterp is galsim._galsim.Quintic but I wasn't able to find the shoot method.
Here is simple code that returns different results on Mac ARM and Linux (galsim version: 2.9)

import galsim
import galsim.roman as roman

bp = roman.getBandpasses()["J129"]
psf = galsim.OpticalPSF(
    diam=roman.diameter,
    lam=bp.effective_wavelength,
)

psf.shoot(5, galsim.BaseDeviate(1234)).x

On my Mac I get:

array([ 0.00834068, -0.01800069, -0.10280517,  0.0007287 ,  0.05277417])

On Linux I get:

array([-0.00932278, -0.05793788, -0.10280517, -0.00044896,  0.02334157])

When drawing galaxy profile, the number of photons drawn is the same and the flux of the photons is the same but they end up at different locations.

Start with the reproducer and inspect src/SBInterpolatedImage.cpp around line 1306, following the OpticalPSF-to-InterpolatedImage path. Investigate how _xInterp, identified as galsim._galsim.Quintic, participates in photon shooting, then run the supplied snippet on Mac ARM and Linux. Done means the platform-dependent photon locations are explained and the discrepancy is resolved or documented.