closed_loop_mpi Namespace Reference

Functions
def	loop (n)

Variables
	rank = int(os.environ['OMPI_COMM_WORLD_RANK'])
	c = ch.carmaWrap_context()
	comm = MPI.COMM_WORLD
	comm_size = comm.Get_size()
string	error = 'command line should be:"python test.py parameters_filename"\n with "parameters_filename" the path to the parameters file'
	param_file = sys.argv[1]
	filename = param_file.split('/')[-1]
	param_path = param_file.split(filename)[0]
string	simul_name = ""
dictionary	matricesToLoad = {}
int	clean = 1
	param_dict = h5u.params_dictionary(config)
	wfs
	tel
	atm
	dms = ao.dm_init(config.p_dms, config.p_wfss, wfs, config.p_geom, config.p_tel)
	tar
	rtc
int	mimg = 0.

Function Documentation

loop()

def closed_loop_mpi.loop

(

n

)

Definition at line 120 of file closed_loop_mpi.py.

def loop(n):
    # if(rank==0):
    #fig,((turbu,image),(shak,defMir))=pl.subplots(2,2, figsize=(15,15))
    # pl.ion()
    # pl.show()
 
    t0 = time.time()
    for i in range(n):
        if (rank == 0):
            atm.move_atmos()
            for t in range(config.p_target.ntargets):
                tar.atmos_trace(t, atm, tel)
                tar.dmtrace(t, dms)
            for w in range(len(config.p_wfss)):
                wfs.sensors_trace(w, "all", tel, atm, dms)
        wfs.Bcast_dscreen()
        for w in range(len(config.p_wfss)):
            wfs.sensors_compimg(w)
            wfs.gather_bincube(w)
        if (rank == 0):
            rtc.docentroids(0)
            rtc.docontrol(0)
            rtc.applycontrol(0, dms)
 
        if ((i + 1) % 50 == 0):
            # s=rtc.get_centroids(0)
            if (rank == 0):
                """ FOR DEBUG PURPOSE
                turbu.clear()
                image.clear()
                shak.clear()
                defMir.clear()
 
                screen=atm.get_screen(0.)
 
                im=tar.get_image(0,"se")
                im=np.roll(im,im.shape[0]/2,axis=0)
                im=np.roll(im,im.shape[1]/2,axis=1)
 
                #sh=wfs.get_binimg(0)
 
                dm=dms.get_dm("pzt",0.)
 
                f1=turbu.matshow(screen,cmap='Blues_r')
                f2=image.matshow(im,cmap='Blues_r')
                #f3=shak.matshow(sh,cmap='Blues_r')
                f4=defMir.matshow(dm)
                pl.draw()
 
 
                c=rtc.get_command(0)
                v=rtc.get_voltages(0)
 
                sh_file="dbg/shak_"+str(i)+"_np_"+str(comm.Get_size())+".npy"
                im_file="dbg/imag_"+str(i)+"_np_"+str(comm.Get_size())+".npy"
                dm_file="dbg/DM_"+str(i)+"_np_"+str(comm.Get_size())+".npy"
                s_file="dbg/cent_"+str(i)+"_np_"+str(comm.Get_size())+".npy"
                c_file="dbg/comm_"+str(i)+"_np_"+str(comm.Get_size())+".npy"
                v_file="dbg/volt_"+str(i)+"_np_"+str(comm.Get_size())+".npy"
 
                np.save(sh_file,sh)
                np.save(im_file,im)
                np.save(dm_file,dm)
                np.save(s_file,s)
                np.save(c_file,c)
                np.save(v_file,v)
                """
 
                strehltmp = tar.get_strehl(0)
                print("%5d" % (i + 1), "  %1.5f" % strehltmp[0],
                      "  %1.5f" % strehltmp[1])
 
    t1 = time.time()
    print(rank, "| loop execution time:", t1 - t0, "  (", n, "iterations), ",
          (t1 - t0) / n, "(mean)  ", n / (t1 - t0), "Hz")
 
 
loop(config.p_loop.niter)

Variable Documentation

atm

closed_loop_mpi.atm

Initial value:

=  ao.atmos_init(c, config.p_atmos, config.p_tel, config.p_geom, config.p_loop,
                    rank=rank, load=matricesToLoad)

Definition at line 76 of file closed_loop_mpi.py.

c

closed_loop_mpi.c = ch.carmaWrap_context()

Definition at line 13 of file closed_loop_mpi.py.

clean

int closed_loop_mpi.clean = 1

Definition at line 61 of file closed_loop_mpi.py.

comm

closed_loop_mpi.comm = MPI.COMM_WORLD

Definition at line 22 of file closed_loop_mpi.py.

comm_size

closed_loop_mpi.comm_size = comm.Get_size()

Definition at line 23 of file closed_loop_mpi.py.

dms

closed_loop_mpi.dms = ao.dm_init(config.p_dms, config.p_wfss, wfs, config.p_geom, config.p_tel)

Definition at line 81 of file closed_loop_mpi.py.

error

string closed_loop_mpi.error = 'command line should be:"python test.py parameters_filename"\n with "parameters_filename" the path to the parameters file'

Definition at line 29 of file closed_loop_mpi.py.

filename

closed_loop_mpi.filename = param_file.split('/')[-1]

Definition at line 35 of file closed_loop_mpi.py.

matricesToLoad

closed_loop_mpi.matricesToLoad = {}

Definition at line 59 of file closed_loop_mpi.py.

mimg

int closed_loop_mpi.mimg = 0.

Definition at line 115 of file closed_loop_mpi.py.

param_dict

closed_loop_mpi.param_dict = h5u.params_dictionary(config)

Definition at line 64 of file closed_loop_mpi.py.

param_file

closed_loop_mpi.param_file = sys.argv[1]

Definition at line 33 of file closed_loop_mpi.py.

param_path

closed_loop_mpi.param_path = param_file.split(filename)[0]

Definition at line 36 of file closed_loop_mpi.py.

rank

closed_loop_mpi.rank = int(os.environ['OMPI_COMM_WORLD_RANK'])

Definition at line 12 of file closed_loop_mpi.py.

rtc

closed_loop_mpi.rtc

Initial value:

=  ao.rtc_init(tel, wfs, config.p_wfss, dms, config.p_dms, config.p_geom,
                  config.p_rtc, config.p_atmos, atm, config.p_tel, config.p_loop,
                  clean=clean, simul_name=simul_name, load=matricesToLoad)

Definition at line 90 of file closed_loop_mpi.py.

simul_name

string closed_loop_mpi.simul_name = ""

Definition at line 52 of file closed_loop_mpi.py.

tar

closed_loop_mpi.tar

Initial value:

=  ao.target_init(c, tel, config.p_target, config.p_atmos, config.p_geom,
                     config.p_tel, config.p_dms)

Definition at line 85 of file closed_loop_mpi.py.

tel

closed_loop_mpi.tel

Definition at line 71 of file closed_loop_mpi.py.

wfs

closed_loop_mpi.wfs

Definition at line 71 of file closed_loop_mpi.py.