test_pyrhr.py

import shesha.config as conf
 
# simul_name="bench_scao_pyr40_8pix"
 
# loop
p_loop = conf.Param_loop()
 
p_loop.set_niter(1000)
p_loop.set_ittime(0.002)  # =1/500
p_loop.set_devices([0, 1, 2, 3])
# geom
p_geom = conf.Param_geom()
 
p_geom.set_zenithangle(0.)
#p_geom.set_pupdiam(288)
 
# tel
p_tel = conf.Param_tel()
 
p_tel.set_diam(8.0)
p_tel.set_cobs(0.12)
 
# atmos
p_atmos = conf.Param_atmos()
 
p_atmos.set_r0(0.16)
p_atmos.set_nscreens(1)
p_atmos.set_frac([1.0])
p_atmos.set_alt([0.0])
p_atmos.set_windspeed([10.0])
p_atmos.set_winddir([0])
p_atmos.set_L0([100.])
 
# target
p_target = conf.Param_target()
p_targets = [p_target]
p_target.set_xpos(0.)
p_target.set_ypos(0.)
p_target.set_Lambda(1.65)
p_target.set_mag(10.)
 
# wfs
p_wfs0 = conf.Param_wfs()
p_wfss = [p_wfs0]
 
p_wfs0.set_type("pyrhr")
p_wfs0.set_nxsub(16)
p_wfs0.set_fssize(1.5)
p_wfs0.set_fracsub(0.8)
p_wfs0.set_xpos(0.)
p_wfs0.set_ypos(0.)
p_wfs0.set_Lambda(0.5)
p_wfs0.set_gsmag(5.)
p_wfs0.set_optthroughput(0.5)
p_wfs0.set_zerop(1.e11)
p_wfs0.set_noise(-1)
p_wfs0.set_fstop("round")
p_wfs0.set_pyr_npts(16)
p_wfs0.set_pyr_ampl(3)
p_wfs0.set_pyr_pup_sep(p_wfs0.nxsub)
p_wfs0.set_atmos_seen(1)
 
# dm
p_dm0 = conf.Param_dm()
p_dm1 = conf.Param_dm()
p_dms = [p_dm0, p_dm1]
p_dm0.set_type("pzt")
nact = p_wfs0.nxsub + 1
p_dm0.set_nact(nact)
p_dm0.set_alt(0.)
p_dm0.set_thresh(0.3)
p_dm0.set_coupling(0.2)
p_dm0.set_unitpervolt(0.01)
p_dm0.set_push4imat(100.)
 
p_dm1.set_type("tt")
p_dm1.set_alt(0.)
p_dm1.set_unitpervolt(0.0005)
p_dm1.set_push4imat(100)
 
# centroiders
p_centroider0 = conf.Param_centroider()
p_centroiders = [p_centroider0]
 
p_centroider0.set_nwfs(0)
p_centroider0.set_type("maskedpix")
# p_centroider0.set_type("corr")
# p_centroider0.set_type_fct("model")
 
# controllers
p_controller0 = conf.Param_controller()
p_controllers = [p_controller0]
 
p_controller0.set_type("ls")
p_controller0.set_nwfs([0])
p_controller0.set_ndm([0, 1])
p_controller0.set_maxcond(100.)
p_controller0.set_delay(1)
p_controller0.set_gain(0.4)