COMPASS

il y a un léger problème a la ligne 662 du programme yoga_ao.i

En effet, le roll que tu fais a pour but de convoluer les spots obtenus par make_lgs_prod1d avec la PSF turbulente, approximée a une gaussienne de largeur a mi hauteur 2*nrebin, mais en fait comme tmp est un tableau (30,30,32) et tmp2 (30,30) il fait des trucs chelou. Au lieu de simplement faire un roll sur chaque sous pupille il les désorganise ... Or il n'y a rien qui réorganise après.

Pour vérifier ca, j'ai lancé une script_system avec une fonction de centroide de type "model" (même pas besoin d'attendre la fin il faut juste que ca initialise les poids) et affichés les poids et le y_wfs(1)._lgskern directement. On voit que les spots ne sont pas dans la meme direction (apres j'ai fait quelques magouilles pour afficher les pupilles mais ca affiche les même au même endroit).

Ma solution pour changer ca : au lieu de faire un roll global il faut faire une boucle sur le nombre de sous pupilles :

tmp3=array(float,y_wfs(nwfs).npix,y_wfs(nwfs).npix,y_wfs(nwfs).nsspup)
for (i=1;i<=nsspup;i++) {
tmp3(,,i) = roll( (fft(fft(tmp(,,i))*(fft(tmp2)),-1) ).re);
}

et la suite reste inchangée.