### Script pour comparer et verifier la nouvelle version de calcul
## des transforme en Ylm (optimise ds LambdaLMBuilder) et l'ancienne
##  Reza, Mai 2006

traceon
timingon
#  Creation d'une carte a partir d'un c(l) et calcul transforme inverse
#  Choix de type de cartes -> ThetaPhi
settypemap T
LMAX = 511
# lmax+1
CLSIZE = $LMAX+1
# Param. pixelisation Theta,Phi 
M_TP = 400
# Une courbe de C(l)
newvec vcl $CLSIZE 1./(x+10.)*(1+sin(x*M_PI/50)) 
# SynFast : C(l) -> map 
cl2map vcl map $M_TP
# AnaFast : map -> alm
echo "---> map2alm map mtxalm $LMAX 1"
map2alm map mtxalm $LMAX 1
echo "---> alm2map mtxalm mapfalm $M_TP"
alm2map mtxalm mapfalm $M_TP
print map 
disp map

#  On sauve la carte et alm ds un fichier
#  A faire sur la version de reference du programme/SOPHYA (ancien)
savelist mapfalm mtxalm map vcl cksamba.ppf 

###  VERIFICATION 
delobjs *   
# Lecture du fichier cree ci-dessus (avec la version de reference)
openppf cksamba.ppf 
#  Choix de type de cartes -> ThetaPhi 
settypemap T
#  Fabrication carte a partir des alm du fichier
echo "---> alm2map mtxalm mapfalm2 $M_TP"
alm2map mtxalm mapfalm2 $M_TP
print map 
print mapfalm2 
#  Calcul des alm sur la carte du fichier
map2alm map mtxalm2 $LMAX 1

#  Ce bout de code ne peut etre execute qu'avec les dernieres modifs du 2 juin 2006
c++import Samba SkyMap
c++exec \ 
  TMatrix< complex<r_8> > diffalm = mtxalm-mtxalm2; \
  KeepObj(diffalm); \
  SphereThetaPhi<r_8> diffmap = mapfalm-mapfalm2; \
  KeepObj(diffmap); 

zone 1 2 
disp diffmap 'lut=lin,-0.001,0.001'
disp diffalm 'lut=lin,0.,0.05'

expmeansig diffmap val
expmeansig diffalm mod

saveobjs * cksambaA.ppf 


#### Pour verifier en version r_8 et r_4
## Modifier le type dans le fichier cksphylm.cc
c++import Samba SkyMap
delobjs *
c++execfrf cksphylm.cc 
zone 1 
disp clin
disp clfm 'same red'
sleep 4
disp mapdiff
expmeansig mapdiff val