def test_02():
"""
Particle above surface.
"""
mp = msPSF.m_params
rv = numpy.arange(0.0, 1.01, 0.1)
zv = numpy.arange(-1.0, 1.01, 0.2)
fast_rz = msPSF.gLZRFocalScan(mp, rv, zv, pz=0.1)
slow_rz = msPSF.gLZRFocalScanSlow(mp, rv, zv, pz=0.1)
assert (numpy.allclose(fast_rz, slow_rz, atol=1.0e-4, rtol=1.0e-4))
def test_01():
"""
Particle on surface.
"""
mp = msPSF.m_params
rv = numpy.arange(0.0, 1.01, 0.1)
zv = numpy.arange(-1.0, 1.01, 0.2)
fast_rz = msPSF.gLZRFocalScan(mp, rv, zv)
slow_rz = msPSF.gLZRFocalScanSlow(mp, rv, zv)
assert (numpy.allclose(fast_rz, slow_rz))
def test_03():
"""
Detector offset.
"""
mp = msPSF.m_params
rv = numpy.arange(0.0, 1.01, 0.1)
zv = numpy.arange(-1.0, 1.01, 0.2)
zd = mp["zd0"] + 1000
fast_rz = msPSF.gLZRFocalScan(mp, rv, zv, zd=zd)
slow_rz = msPSF.gLZRFocalScanSlow(mp, rv, zv, zd=zd)
assert (numpy.allclose(fast_rz, slow_rz, atol=1.0e-4, rtol=1.0e-4))