def test_ne():
""" Check that inequality works as expected."""
spec1 = lambda x: x/1000
spec2 = galsim.LookupTable([400, 550], [0.4, 0.55], interpolant='linear')
spec3 = '3'
# These should all compare unequal.
seds = [galsim.SED(spec1, wave_type='nm', flux_type='flambda'),
galsim.SED(spec1, wave_type='A', flux_type='flambda'),
galsim.SED(spec1, wave_type='nm', flux_type='fnu'),
galsim.SED(spec1, 'nm', 'flambda', redshift=1.0),
galsim.SED(spec2, 'nm', 'flambda'),
galsim.SED(spec3, 'nm', 'flambda')]
all_obj_diff(seds)
def test_ne():
""" Check that inequality works as expected."""
spec1 = lambda x: x/1000
spec2 = galsim.LookupTable([400, 550], [0.4, 0.55], interpolant='linear')
spec3 = '3'
# These should all compare unequal.
seds = [galsim.SED(spec1, wave_type='nm', flux_type='flambda'),
galsim.SED(spec1, wave_type='A', flux_type='flambda'),
galsim.SED(spec1, wave_type='nm', flux_type='fnu'),
galsim.SED(spec1, 'nm', 'flambda', redshift=1.0),
galsim.SED(spec2, 'nm', 'flambda'),
galsim.SED(spec3, 'nm', 'flambda')]
all_obj_diff(seds)
def test_ne():
"""Test Apertures, PhaseScreens, PhaseScreenLists, and PhaseScreenPSFs for not-equals."""
pupil_plane_im = galsim.fits.read(os.path.join(imgdir, pp_file))
# Test galsim.Aperture __ne__
objs = [
galsim.Aperture(diam=1.0),
galsim.Aperture(diam=1.1),
galsim.Aperture(diam=1.0, oversampling=1.5),
galsim.Aperture(diam=1.0, pad_factor=1.5),
galsim.Aperture(diam=1.0, circular_pupil=False),
galsim.Aperture(diam=1.0, obscuration=0.3),
galsim.Aperture(diam=1.0, nstruts=3),
galsim.Aperture(diam=1.0, nstruts=3, strut_thick=0.2),
galsim.Aperture(diam=1.0, nstruts=3, strut_angle=15 * galsim.degrees),
galsim.Aperture(diam=1.0, pupil_plane_im=pupil_plane_im),
galsim.Aperture(diam=1.0,
pupil_plane_im=pupil_plane_im,
pupil_angle=10.0 * galsim.degrees)
]
all_obj_diff(objs)
# Test AtmosphericScreen __ne__
rng = galsim.BaseDeviate(1)
objs = [
galsim.AtmosphericScreen(10.0, rng=rng),
galsim.AtmosphericScreen(10.0, rng=rng, vx=1.0),
galsim.AtmosphericScreen(10.0, rng=rng, vy=1.0),
galsim.AtmosphericScreen(10.0, rng=rng, alpha=0.999, time_step=0.01),
galsim.AtmosphericScreen(10.0, rng=rng, altitude=1.0),
galsim.AtmosphericScreen(10.0, rng=rng, alpha=0.999, time_step=0.02),
galsim.AtmosphericScreen(10.0, rng=rng, alpha=0.998, time_step=0.02),
galsim.AtmosphericScreen(10.0, rng=rng, r0_500=0.1),
galsim.AtmosphericScreen(10.0, rng=rng, L0=10.0),
galsim.AtmosphericScreen(10.0, rng=rng, vx=10.0),
]
all_obj_diff(objs)
# Test OpticalScreen __ne__
objs = [
galsim.OpticalScreen(diam=1.0, ),
galsim.OpticalScreen(diam=1.0, tip=1.0),
galsim.OpticalScreen(diam=1.0, tilt=1.0),
galsim.OpticalScreen(diam=1.0, defocus=1.0),
galsim.OpticalScreen(diam=1.0, astig1=1.0),
galsim.OpticalScreen(diam=1.0, astig2=1.0),
galsim.OpticalScreen(diam=1.0, coma1=1.0),
galsim.OpticalScreen(diam=1.0, coma2=1.0),
galsim.OpticalScreen(diam=1.0, trefoil1=1.0),
galsim.OpticalScreen(diam=1.0, trefoil2=1.0),
galsim.OpticalScreen(diam=1.0, spher=1.0),
galsim.OpticalScreen(diam=1.0, spher=1.0, lam_0=100.0),
galsim.OpticalScreen(diam=1.0, aberrations=[0, 0, 1.1]), # tip=1.1
]
all_obj_diff(objs)
# Test PhaseScreenList __ne__
atm = galsim.Atmosphere(10.0, vx=1.0)
objs = [
galsim.PhaseScreenList(atm),
galsim.PhaseScreenList(objs), # Reuse list of OpticalScreens above
galsim.PhaseScreenList(objs[0:2])
]
all_obj_diff(objs)
# Test PhaseScreenPSF __ne__
psl = galsim.PhaseScreenList(atm)
objs = [galsim.PhaseScreenPSF(psl, 500.0, exptime=0.03, diam=1.0)]
objs += [galsim.PhaseScreenPSF(psl, 700.0, exptime=0.03, diam=1.0)]
objs += [galsim.PhaseScreenPSF(psl, 700.0, exptime=0.03, diam=1.1)]
objs += [
galsim.PhaseScreenPSF(psl, 700.0, exptime=0.03, diam=1.0, flux=1.1)
]
objs += [
galsim.PhaseScreenPSF(psl,
700.0,
exptime=0.03,
diam=1.0,
interpolant='linear')
]
stepk = objs[0].stepk
maxk = objs[0].maxk
objs += [
galsim.PhaseScreenPSF(psl,
700.0,
exptime=0.03,
diam=1.0,
_force_stepk=stepk / 1.5)
]
objs += [
galsim.PhaseScreenPSF(psl,
700.0,
exptime=0.03,
diam=1.0,
_force_maxk=maxk * 2.0)
]
all_obj_diff(objs)