shesha.widgets.testNoise Namespace Reference

Functions
def	testNoise (noise=0.9, niter=100, close_loop=True)

Variables
list	errscreenList = []
list	imagesErrList = []
list	imagesHRErrList = []
list	slopesErrList = []
list	subErrList = []
	errscreen
	imagesErr
	imagesHRErr
	slopesErr
	subErr

Function Documentation

testNoise()

def shesha.widgets.testNoise.testNoise	(		noise = 0.9,
			niter = 100,
			close_loop = True
	)

Definition at line 6 of file testNoise.py.

def testNoise(noise=0.9, niter=100, close_loop=True):
 
    s1 = np.zeros(((niter, int(wao.sim.config.p_controller0.nvalid * 2))))
    s2 = np.zeros(((niter, int(wao.sim.config.p_controller0.nvalid * 2))))
    sub1 = np.zeros(((niter, int(wao.sim.config.p_controller0.nvalid))))
    sub2 = np.zeros(((niter, int(wao.sim.config.p_controller0.nvalid))))
    dimim = wao.sim.wfs.get_pyrimg(0).shape[0]
    dimimHR = wao.sim.wfs.get_pyrimghr(0).shape[0]
    dimscreen = wao.sim.atm.get_screen(0).shape[0]
    im1 = np.zeros(((niter, dimim, dimim)))
    im2 = np.zeros(((niter, dimim, dimim)))
    #im1HR = np.zeros(((niter, dimimHR,dimimHR)))
    #im2HR = np.zeros(((niter, dimimHR,dimimHR)))
    screen1 = np.zeros(((niter, dimscreen, dimscreen)))
    screen2 = np.zeros(((niter, dimscreen, dimscreen)))
 
    wao.sim.rtc.set_open_loop(0, 1)
    wao.sim.atm.set_seed(0, 1234)
    wao.sim.atm.refresh_screen(0)
    wao.sim.wfs.set_noise(0, noise)
    wao.sim.dms.resetdm(b"pzt", 0.)
    wao.sim.dms.resetdm(b"tt", 0.)
    if close_loop:
        wao.sim.rtc.set_open_loop(0, 0)
    #s1, v1, t1, aiData1, psfLE1 = wao.record_ao_circular_buffer(niter)
 
    for k in trange(niter):
        wao.sim.next()
        s1[k, :] = wao.sim.rtc.get_centroids(0)
        sub1[k, :] = wao.sim.rtc.get_intensities(0)
 
        im1[k, :, :] = wao.sim.wfs.get_pyrimg(0)
        # im1HR[k,:,:] = wao.sim.wfs.get_pyrimghr(0)
        screen1[k, :, :] = wao.sim.atm.get_screen(0)
    wao.sim.rtc.set_open_loop(0, 1)
    wao.sim.atm.set_seed(0, 1234)
    wao.sim.atm.refresh_screen(0)
    wao.sim.wfs.set_noise(0, noise)
    wao.sim.dms.resetdm(b"pzt", 0.)
    wao.sim.dms.resetdm(b"tt", 0.)
    if close_loop:
        wao.sim.rtc.set_open_loop(0, 0)
    #s2, v2, t2, aiData2, psfLE2 = wao.record_ao_circular_buffer(niter)
 
    for k in trange(niter):
        wao.sim.next()
        s2[k, :] = wao.sim.rtc.get_centroids(0)
        sub2[k, :] = wao.sim.rtc.get_intensities(0)
        im2[k, :, :] = wao.sim.wfs.get_pyrimg(0)
        screen2[k, :, :] = wao.sim.atm.get_screen(0)
        # im2HR[k,:,:] = wao.sim.wfs.get_pyrimghr(0)
 
    screenErr = np.abs(screen2 - screen1).max() / screen1.max()
    imagesErr = np.abs(im2 - im1).max() / im1.max()
    subErr = np.abs(sub2 - sub1).max() / sub1.max()
    # imagesHRErr = np.abs(im2HR-im1HR).max()/im1HR.max()
    imagesHRErr = 0.
    slopesErr = np.abs(s2 - s1).max() / s1.max()
    print("Screen Error =" + str(screenErr))
    # print("Images HR Error ="+str(imagesHRErr))
    print("Images Error =" + str(imagesErr))
    print("Slopes Error =" + str(slopesErr))
    print("Subsum Error =" + str(subErr))
 
    #print(np.abs(psfLE2-psfLE1).max())
    #print(psfLE1.max())
    #print(psfLE2.max())
    return screenErr, imagesErr, imagesHRErr, slopesErr, subErr
    #plt.figure(1)
    #plt.clf()
    #plt.matshow(s2-s1, aspect="auto", fignum=1)
 
 
#plt.matshow(im2[-1,:,:]-im1[-1,:,:], aspect="auto", fignum=2)
 

Variable Documentation

errscreen

shesha.widgets.testNoise.errscreen

Definition at line 87 of file testNoise.py.

errscreenList

list shesha.widgets.testNoise.errscreenList = []

Definition at line 81 of file testNoise.py.

imagesErr

shesha.widgets.testNoise.imagesErr

Definition at line 87 of file testNoise.py.

imagesErrList

list shesha.widgets.testNoise.imagesErrList = []

Definition at line 82 of file testNoise.py.

imagesHRErr

shesha.widgets.testNoise.imagesHRErr

Definition at line 87 of file testNoise.py.

imagesHRErrList

list shesha.widgets.testNoise.imagesHRErrList = []

Definition at line 83 of file testNoise.py.

slopesErr

shesha.widgets.testNoise.slopesErr

Definition at line 87 of file testNoise.py.

slopesErrList

list shesha.widgets.testNoise.slopesErrList = []

Definition at line 84 of file testNoise.py.

subErr

shesha.widgets.testNoise.subErr

Definition at line 87 of file testNoise.py.

subErrList

list shesha.widgets.testNoise.subErrList = []

Definition at line 85 of file testNoise.py.