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))