source: Sophya/trunk/SophyaPI/ProgPI/skymapmodule.pic@ 3891

#---------------------- Create Maps and Save
shell rm -f h4*.fits t4*.ppf h8*.fits t8*.ppf clh*.ppf clt*.ppf almh*.ppf almt*.ppf
delobjs *
c++import SkyMap
c++exec \
  double theta,phi; \
  SphereHEALPix<r_4>  h4(256); \
  for(int_4 i=0;i<h4.NbPixels();i++) { \
    h4.PixThetaPhi(i,theta,phi); \
    h4(i) = 0.; \
    double v = (M_PI_2-theta)*(M_PI_2-theta)/(0.50*0.50); \
    v += (phi-M_PI)*(phi-M_PI)/(M_PI*M_PI/2.); \
    v = 100.*exp(-v/2.); \
    v += 10.*NorRand(); \
    if(phi<20.*M_PI/180.||phi>340.*M_PI/180.) continue; \
    if(phi<190.*M_PI/180.&&phi>170.*M_PI/180.) continue; \
    h4(i) = v; \
  } \
  KeepObj(h4); \
  SphereThetaPhi<r_4> t4(400); \
  for(int_4 i=0;i<t4.NbPixels();i++) { \
    t4.PixThetaPhi(i,theta,phi); \
    t4(i) = 0.; \
    double v = (M_PI_2-theta)*(M_PI_2-theta)/(0.50*0.50); \
    v += (phi-M_PI)*(phi-M_PI)/(M_PI*M_PI/2.); \
    v = 100.*exp(-v/2.); \
    v += 10.*NorRand(); \
    if(phi<20.*M_PI/180.||phi>340.*M_PI/180.) continue; \
    if(phi<190.*M_PI/180.&&phi>170.*M_PI/180.) continue; \
    t4(i) = v; \
  } \
  KeepObj(t4); \
  cout<<"End of Job"<<endl;

savefits h4 !h4.fits
saveppf t4 t4.ppf

zone 1 2
disp h4 zoom/4
disp t4 zoom/4

#---------------------- test Map2Double
delobjs *
zone
openfits h4.fits

map2double h4
rename h4 h8
savefits h8 !h8.fits

delobjs *
zone
openppf t4.ppf
map2double t4
rename t4 t8
saveppf t8 t8.ppf

#---------------------- test Map2Float
delobjs *
zone
openfits h8.fits
openppf t8.ppf

map2float h8
map2float t8
typemap h8
typemap t8

#---------------------- test TypeMap
delobjs *
zone
openfits h4.fits
openfits h8.fits
openppf t4.ppf
openppf t8.ppf
typemap h4
typemap h8
typemap t4
typemap t8

#---------------------- test Map2Map
delobjs *
zone
openfits h8.fits
openppf t8.ppf
map2map h8 t
map2map t8 h

#---------------------- test MapMult
delobjs *
zone
openfits h8.fits
openppf t8.ppf
mapmult h8 1000
mapmult t8 1000

#---------------------- test MapProj
delobjs *
zone
openfits h8.fits
maproj h8 h8p
mapop h8p - h8
disp h8p zoom/4

maproj h8 h8p64 64
disp h8p64 zoom/4
savefits h8p64 !h8p64.fits

delobjs *
zone
openppf t8.ppf
maproj t8 t8p
mapop t8p - t8
disp t8p zoom/4

maproj t8 t8p100 100
disp t8p100 zoom/4
saveppf t8p100 t8p100.ppf

#---------------------- test Map2Cl
delobjs *
zone
openfits h8p64.fits
map2cl h8p64 clh64
saveppf clh64 clh64.ppf
n/plot clh64.n*(n+1)*val%n ! ! crossmarker5

map2cl h8p64 clh64 50
n/plot clh64.n*(n+1)*val%n ! ! crossmarker5

delobjs *
zone
openppf t8p100.ppf
map2cl t8p100 clt100 191
saveppf clt100 clt100.ppf
n/plot clt100.n*(n+1)*val%n ! ! crossmarker5

#---------------------- test Map2Alm
delobjs *
zone
openfits h8p64.fits
map2alm h8p64 almh64
saveppf almh64 almh64.ppf
disp almh64

delobjs *
zone
openppf t8p100.ppf
map2alm t8p100 almt100 191
saveppf almt100 almt100.ppf
disp almt100

#---------------------- test Alm2Cl
delobjs *
zone
openppf clh64.ppf
openppf almh64.ppf
alm2cl almh64 clh64fralm
n/plot clh64.n*(n+1)*val%n ! ! crossmarker5
n/plot clh64fralm.n*(n+1)*val%n ! ! "red same circlemarker5"
c++exec clh64fralm -= clh64;
n/plot clh64fralm.n*(n+1)*val%n ! ! crossmarker5

delobjs *
zone
openppf clt100.ppf
openppf almt100.ppf
alm2cl almt100 clt100fralm
n/plot clt100.n*(n+1)*val%n ! ! crossmarker5
n/plot clt100fralm.n*(n+1)*val%n ! ! "red same circlemarker5"
c++exec clt100fralm -= clt100;
n/plot clt100fralm.n*(n+1)*val%n ! ! crossmarker5

#---------------------- test Cl2Map et SetTypeMap
delobjs *
zone
openppf clh64.ppf
settypemap h
cl2map clh64 h8p64frcl 64
savefits h8p64frcl !h8p64frcl.fits
disp h8p64frcl

delobjs *
zone
openppf clt100.ppf
settypemap t
cl2map clt100 t8p100frcl 100
saveppf t8p100frcl t8p100frcl.ppf
disp t8p100frcl

#---------------------- test Alm2Map et SetTypeMap
delobjs *
zone
openppf almh64.ppf
settypemap h
alm2map almh64 h8p64fralm 64
savefits h8p64fralm !h8p64fralm.fits
disp h8p64fralm

delobjs *
zone
openppf almt100.ppf
settypemap t
alm2map almt100 t8p100fralm 100
saveppf t8p100fralm t8p100fralm.ppf
disp t8p100fralm

#---------------------- test Cl2llCl
delobjs *
zone
openppf clh64.ppf
cl2llcl clh64 llclh64
n/plot clh64.n*(n+1)*val%n ! ! crossmarker5
n/plot llclh64.val%n ! ! "red same circlemarker5"

delobjs *
zone
openppf clt100.ppf
cl2llcl clt100 llclt100
n/plot clt100.n*(n+1)*val%n ! ! crossmarker5
n/plot llclt100.val%n ! ! "red same circlemarker5"

#---------------------- test ClMean
delobjs *
zone
openppf clh64.ppf
cl2llcl clh64 llclh64
n/plot llclh64.val%n ! ! circlemarker5
clmean llclh64
clmean llclh64 !
clmean llclh64 100,9999 mean
echo $mean

delobjs *
zone
openppf clt100.ppf
cl2llcl clt100 llclt100
n/plot llclt100.val%n ! ! circlemarker5
clmean llclt100
clmean llclt100 !
clmean llclt100 100,9999 mean
echo $mean

#---------------------- test ClMult
delobjs *
openppf clh64.ppf
zone 1 2
n/plot clh64.val%n ! ! circlemarker5
clmult clh64 1000.
n/plot clh64.val%n ! ! circlemarker5

delobjs *
openppf clt100.ppf
zone 1 2
n/plot clt100.val%n ! ! circlemarker5
clmult clt100 1000.
n/plot clt100.val%n ! ! circlemarker5

#---------------------- test ClOper
delobjs *
zone
openppf clh64.ppf
cp clh64 dum
clop dum + clh64 - clh64 * clh64 / clh64
n/plot clh64.val%n ! ! crossmarker5
n/plot dum.val%n ! ! "same red circlemarker5"

delobjs *
zone
openppf clt100.ppf
cp clt100 dum
clop dum + clt100 - clt100 * clt100 / clt100
n/plot clt100.val%n ! ! crossmarker5
n/plot dum.val%n ! ! "same red circlemarker5"

delobjs *
zone
openppf clh64.ppf
openppf clt100.ppf
clop clh64 - clt100
n/plot clh64.val%n ! ! crossmarker5

#---------------------- test ClRebin
delobjs *
zone
openppf clh64.ppf
cl2llcl clh64 llclh64
clrebin llclh64 clntu 10,0

n/plot clntu.n%l ! ! crossmarker5

n/plot llclh64.val%n  ! ! "crossmarker3"
n/plot clntu.clmean%l ! ! "same circlemarker5 red"
n/plot clntu.cllin%l  ! ! "same boxmarker5 blue"
n/plot clntu.clpar%l  ! ! "same trianglemarker5 darkviolet"

n/plot llclh64.val%n  ! ! "crossmarker3"
nt2d clntu l clmean 0 sclmean 1 " " "same circlemarker5 red thinline"

n/plot llclh64.val%n  ! ! "crossmarker3"
nt2d clntu l cllin 0 scllin 1 " " "same boxmarker5 blue thinline"

n/plot llclh64.val%n  ! ! "crossmarker3"
nt2d clntu l clpar 0 sclpar 1 " " "same trianglemarker5 darkviolet thinline"

#---------------------- test CrMaskMap et SetTypeMap
delobjs *
zone
settypemap h
crmapmask h8m 256 -20,20 100,130 0,1
savefits h8m !h8m.fits
disp h8m zoom/4

delobjs *
zone
settypemap t
crmapmask t8m 400 -20,20 100,130 0,1
saveppf t8m t8m.ppf
disp t8m zoom/4

#---------------------- test CrMaskFrMap
delobjs *
zone
openfits h8.fits
crmaskfrmap h8fm 256 h8 -1.e-30,1.e-30 0,1
savefits h8fm !h8fm.fits
disp h8fm zoom/4

delobjs *
zone
openppf t8.ppf
crmaskfrmap t8fm 400 t8 -1.e-30,1.e-30 0,1
saveppf t8fm t8fm.ppf
disp t8fm zoom/4

#---------------------- test MaskMap
delobjs *
zone
openfits h8.fits
openfits h8m.fits
maskmap h8 h8m
disp h8 zoom/4

delobjs *
zone
openppf t8.ppf
openppf t8m.ppf
maskmap t8 t8m
disp t8 zoom/4

delobjs *
zone
openfits h8.fits
openppf t8m.ppf
maskmap h8 t8m
disp h8 zoom/4

delobjs *
zone
openppf t8.ppf
openfits h8m.fits
maskmap t8 h8m
disp t8 zoom/4

#---------------------- test MapCover
delobjs *
zone
openfits h8m.fits
mapcover h8m 0.9,1 couvh8m
echo $couvh8m

delobjs *
zone
openppf t8m.ppf
mapcover t8m 0.9,1 couvt8m
echo $couvt8m

#---------------------- test Map2Local
delobjs *
zone
openfits h8.fits
map2local h8 h8loc 200,300 20,30 0,90 !
disp h8loc
map2local h8 h8loc 200,300 20,30 0,90 ! 45
disp h8loc

delobjs *
zone
openppf t8.ppf
map2local t8 t8loc 200,300 20,30 0,90 !
disp t8loc
map2local t8 t8loc 200,300 20,30 0,90 ! 45
disp t8loc

#---------------------- test MapOper
delobjs *
zone
openfits h8.fits
cp h8 h8save
mapop h8save + h8 - h8 * h8 / h8
disp h8save zoom/4
mapop h8save - h8
disp h8save zoom/4

delobjs *
zone
openppf t8.ppf
cp t8 t8save
mapop t8save + t8 - t8 * t8 / t8
disp t8save zoom/4
mapop t8save - t8
disp t8save zoom/4

delobjs *
zone
openfits h8.fits
openppf t8.ppf
mapop h8 - t8
disp h8 zoom/4

#---------------------- test MapStat
delobjs *
zone
openfits h8.fits
openfits h8fm.fits
mapstat h8 !    mean sig
echo mean=$mean sig=$sig
mapstat h8 h8fm mean sig
echo mean=$mean sig=$sig

openppf t8.ppf
openppf t8fm.ppf
mapstat t8 !    mean sig
echo mean=$mean sig=$sig
mapstat t8 t8fm mean sig
echo mean=$mean sig=$sig

#---------------------- test ExpVal
set x 1
set y 2
set z 3
eval ! x +y +z
eval ! $x +$y +$z
eval r x +y +z
 echo result is $r

eval r 1e150
 echo result is $r
eval r 1e10
 echo result is $r
eval r 0
 echo result is $r

eval r sqrt(-1.)
 echo $r