shesha.init.target_init Namespace Reference

Initialization of a Target object. More...

Functions
def	target_init (carmaWrap_context ctxt, Telescope telescope, list p_targets, conf.Param_atmos p_atmos, conf.Param_tel p_tel, conf.Param_geom p_geom, dm=None, brahma=False)
	Create a cython target from parametres structures. More...

Detailed Description

Initialization of a Target object.

Author

COMPASS Team https://github.com/ANR-COMPASS

Version

5.0.0

Date

2020/05/18

Copyright

GNU Lesser General Public License

This file is part of COMPASS https://anr-compass.github.io/compass/

Copyright (C) 2011-2019 COMPASS Team https://github.com/ANR-COMPASS All rights reserved. Distributed under GNU - LGPL

COMPASS is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or any later version.

COMPASS: End-to-end AO simulation tool using GPU acceleration The COMPASS platform was designed to meet the need of high-performance for the simulation of AO systems.

The final product includes a software package for simulating all the critical subcomponents of AO, particularly in the context of the ELT and a real-time core based on several control approaches, with performances consistent with its integration into an instrument. Taking advantage of the specific hardware architecture of the GPU, the COMPASS tool allows to achieve adequate execution speeds to conduct large simulation campaigns called to the ELT.

The COMPASS platform can be used to carry a wide variety of simulations to both testspecific components of AO of the E-ELT (such as wavefront analysis device with a pyramid or elongated Laser star), and various systems configurations such as multi-conjugate AO.

COMPASS is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public License along with COMPASS. If not, see https://www.gnu.org/licenses/lgpl-3.0.txt.

Function Documentation

target_init()

def shesha.init.target_init.target_init	(	carmaWrap_context	ctxt,
		Telescope	telescope,
		list	p_targets,
		conf.Param_atmos	p_atmos,
		conf.Param_tel	p_tel,
		conf.Param_geom	p_geom,
			dm = None,
			brahma = False
	)

Create a cython target from parametres structures.

:parameters: ctxt (carmaWrap_context) : telescope (Telescope): Telescope object p_targets (lis of Param_target) : target_settings p_atmos (Param_atmos) : atmos settings p_tel (Param_tel) : telescope settings p_geom (Param_geom) : geom settings dm (Param_dm) : (optional) dm settings brahma (bool): (optional) brahma flag :return: tar (Target): Target object

Definition at line 61 of file target_init.py.

    :return:
        tar: (Target): Target object
    """
    type_target = "atmos"
 
    if (p_targets is not None):
        for p_target in p_targets:
            if (p_target.dms_seen is None and dm is not None):
                p_target.dms_seen = np.arange(len(dm))
 
    sizes = np.ones(len(p_targets), dtype=np.int64) * p_geom.pupdiam
 
    ceiled_pupil = np.ceil(p_geom._spupil)
 
    ceiled_pupil[np.where(ceiled_pupil > 1)] = 1
 
    if (p_target.apod == 1):
        Npts = 0
        # TODO apodizer, Npts=nb element of apodizer>0
        ceiled_apodizer = np.ceil(p_geom._apodizer * p_geom._spupil)
        ceiled_apodizer[np.where(ceiled_apodizer > 1)] = 1
        Npts = np.sum(ceiled_apodizer)
    else:
        Npts = np.sum(ceiled_pupil)
 
    xpos = np.array([p_target.xpos for p_target in p_targets], dtype=np.float32)
    ypos = np.array([p_target.ypos for p_target in p_targets], dtype=np.float32)
    Lambda = np.array([p_target.Lambda for p_target in p_targets], dtype=np.float32)
    mag = np.array([p_target.mag for p_target in p_targets], dtype=np.float32)
    zerop = p_targets[0].zerop
 
    if (brahma):
        target = Target_brahma(ctxt, "target_brahma", telescope, 0, len(p_targets), xpos,
                               ypos, Lambda, mag, zerop, sizes, Npts, ctxt.active_device)
    else:
        target = Target(ctxt, telescope, len(p_targets), xpos, ypos, Lambda, mag, zerop,
                        sizes, Npts, ctxt.active_device)
 
    # cc=i
    for i in range(len(p_targets)):
        p_target = p_targets[i]
        if (p_atmos.nscreens > 0):
            for j in range(p_atmos.nscreens):
                xoff = p_target.xpos * CONST.ARCSEC2RAD * \
                    p_atmos.alt[j] / p_atmos.pupixsize
                yoff = p_target.ypos * CONST.ARCSEC2RAD * \
                    p_atmos.alt[j] / p_atmos.pupixsize
                xoff += float((p_atmos.dim_screens[j] - p_geom._n) / 2)
                yoff += float((p_atmos.dim_screens[j] - p_geom._n) / 2)
                pupdiff = (p_geom._n - p_geom.pupdiam) / 2
                xoff += pupdiff
                yoff += pupdiff
                target.d_targets[i].add_layer(type_target, j, xoff, yoff)
 
        # if (y_dm != []) {
        if (dm is not None):
            # j=ddd
            # for (ddd=1;ddd<=numberof(*y_target(cc).dms_seen);ddd++) {
            for j in range(p_target.dms_seen.size):
                # k=dd
                # dd = (*y_target(cc).dms_seen)(ddd)
                k = p_target.dms_seen[j]
                dims = dm[k]._n2 - dm[k]._n1 + 1
                dim = p_geom._mpupil[2].size
                dim_dm = max(dim, dims)
                xoff = p_target.xpos * CONST.ARCSEC2RAD * \
                    dm[k].alt / p_tel.diam * p_geom.pupdiam
                yoff = p_target.ypos * CONST.ARCSEC2RAD * \
                    dm[k].alt / p_tel.diam * p_geom.pupdiam
 
                xoff += float((dim_dm - p_geom._n) / 2)
                yoff += float((dim_dm - p_geom._n) / 2)
 
                pupdiff = (p_geom._n - p_geom.pupdiam) / 2
                xoff += pupdiff
                yoff += pupdiff
 
                # if (dm[k].type == scons.DmType.KL):
                #     xoff -= 2
                #     yoff -= 2
                target.d_targets[i].add_layer(dm[k].type, k, xoff, yoff)
 
        target.d_targets[i].init_strehlmeter()
 
    return target