Ignore:
Timestamp:
Aug 11, 2002, 4:53:54 PM (23 years ago)
Author:
cmv
Message:

+ de fonctionnalites cmv 11/8/02

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaPI/ProgPI/skymapmodule.pic

    r2175 r2177  
    1 #----------------------
     1#---------------------- Create Maps and Save
    22shell rm -f h4*.fits t4*.ppf h8*.fits t8*.ppf clh*.ppf clt*.ppf almh*.ppf almt*.ppf
    33delobjs *
     
    3939disp h4 zoom/4
    4040disp t4 zoom/4
    41 zone
    42 
    43 #----------------------
    44 delobjs *
     41
     42#---------------------- test Map2Double
     43delobjs *
     44zone
    4545openfits h4.fits
    4646
     
    5050
    5151delobjs *
     52zone
    5253openppf t4.ppf
    5354map2double t4
     
    5556saveppf t8 t8.ppf
    5657
    57 #----------------------
    58 delobjs *
     58#---------------------- test Map2Float
     59delobjs *
     60zone
     61openfits h8.fits
     62openppf t8.ppf
     63
     64map2float h8
     65map2float t8
     66typemap h8
     67typemap t8
     68
     69#---------------------- test TypeMap
     70delobjs *
     71zone
    5972openfits h4.fits
    6073openfits h8.fits
     
    6679typemap t8
    6780
    68 #----------------------
    69 delobjs *
     81#---------------------- test Map2Map
     82delobjs *
     83zone
    7084openfits h8.fits
    7185openppf t8.ppf
     
    7387map2map t8 h
    7488
    75 #----------------------
    76 delobjs *
     89#---------------------- test MapMult
     90delobjs *
     91zone
    7792openfits h8.fits
    7893openppf t8.ppf
     
    8095mapmult t8 1000
    8196
    82 #----------------------
    83 delobjs *
     97#---------------------- test MapProj
     98delobjs *
     99zone
    84100openfits h8.fits
    85101maproj h8 h8p
     
    92108
    93109delobjs *
     110zone
    94111openppf t8.ppf
    95112maproj t8 t8p
     
    101118saveppf t8p100 t8p100.ppf
    102119
    103 #----------------------
    104 delobjs *
     120#---------------------- test Map2Cl
     121delobjs *
     122zone
    105123openfits h8p64.fits
    106124map2cl h8p64 clh64
     
    112130
    113131delobjs *
     132zone
    114133openppf t8p100.ppf
    115134map2cl t8p100 clt100 191
     
    117136n/plot clt100.n*(n+1)*val%n ! ! crossmarker5
    118137
    119 #----------------------
    120 delobjs *
     138#---------------------- test Map2Alm
     139delobjs *
     140zone
    121141openfits h8p64.fits
    122142map2alm h8p64 almh64
     
    125145
    126146delobjs *
     147zone
    127148openppf t8p100.ppf
    128149map2alm t8p100 almt100 191
     
    130151disp almt100
    131152
    132 #----------------------
    133 delobjs *
     153#---------------------- test Alm2Cl
     154delobjs *
     155zone
     156openppf clh64.ppf
     157openppf almh64.ppf
     158alm2cl almh64 clh64fralm
     159n/plot clh64.n*(n+1)*val%n ! ! crossmarker5
     160n/plot clh64fralm.n*(n+1)*val%n ! ! "red same circlemarker5"
     161c++exec clh64fralm -= clh64;
     162n/plot clh64fralm.n*(n+1)*val%n ! ! crossmarker5
     163
     164delobjs *
     165zone
     166openppf clt100.ppf
     167openppf almt100.ppf
     168alm2cl almt100 clt100fralm
     169n/plot clt100.n*(n+1)*val%n ! ! crossmarker5
     170n/plot clt100fralm.n*(n+1)*val%n ! ! "red same circlemarker5"
     171c++exec clt100fralm -= clt100;
     172n/plot clt100fralm.n*(n+1)*val%n ! ! crossmarker5
     173
     174#---------------------- test Cl2Map et SetTypeMap
     175delobjs *
     176zone
    134177openppf clh64.ppf
    135178settypemap h
     
    139182
    140183delobjs *
     184zone
    141185openppf clt100.ppf
    142186settypemap t
     
    145189disp t8p100frcl
    146190
    147 #----------------------
    148 delobjs *
     191#---------------------- test Alm2Map et SetTypeMap
     192delobjs *
     193zone
    149194openppf almh64.ppf
    150195settypemap h
     
    154199
    155200delobjs *
     201zone
    156202openppf almt100.ppf
    157203settypemap t
     
    160206disp t8p100fralm
    161207
    162 #----------------------
    163 delobjs *
     208#---------------------- test Cl2llCl
     209delobjs *
     210zone
     211openppf clh64.ppf
     212cl2llcl clh64 llclh64
     213n/plot clh64.n*(n+1)*val%n ! ! crossmarker5
     214n/plot llclh64.val%n ! ! "red same circlemarker5"
     215
     216delobjs *
     217zone
     218openppf clt100.ppf
     219cl2llcl clt100 llclt100
     220n/plot clt100.n*(n+1)*val%n ! ! crossmarker5
     221n/plot llclt100.val%n ! ! "red same circlemarker5"
     222
     223#---------------------- test ClMean
     224delobjs *
     225zone
     226openppf clh64.ppf
     227cl2llcl clh64 llclh64
     228n/plot llclh64.val%n ! ! circlemarker5
     229clmean llclh64
     230clmean llclh64 !
     231clmean llclh64 100,9999 mean
     232echo $mean
     233
     234delobjs *
     235zone
     236openppf clt100.ppf
     237cl2llcl clt100 llclt100
     238n/plot llclt100.val%n ! ! circlemarker5
     239clmean llclt100
     240clmean llclt100 !
     241clmean llclt100 100,9999 mean
     242echo $mean
     243
     244#---------------------- test ClMult
     245delobjs *
     246openppf clh64.ppf
     247zone 1 2
     248n/plot clh64.val%n ! ! circlemarker5
     249clmult clh64 1000.
     250n/plot clh64.val%n ! ! circlemarker5
     251
     252delobjs *
     253openppf clt100.ppf
     254zone 1 2
     255n/plot clt100.val%n ! ! circlemarker5
     256clmult clt100 1000.
     257n/plot clt100.val%n ! ! circlemarker5
     258
     259#---------------------- test ClOper
     260delobjs *
     261zone
     262openppf clh64.ppf
     263cp clh64 dum
     264clop dum + clh64 - clh64 * clh64 / clh64
     265n/plot clh64.val%n ! ! crossmarker5
     266n/plot dum.val%n ! ! "same red circlemarker5"
     267
     268delobjs *
     269zone
     270openppf clt100.ppf
     271cp clt100 dum
     272clop dum + clt100 - clt100 * clt100 / clt100
     273n/plot clt100.val%n ! ! crossmarker5
     274n/plot dum.val%n ! ! "same red circlemarker5"
     275
     276delobjs *
     277zone
     278openppf clh64.ppf
     279openppf clt100.ppf
     280clop clh64 - clt100
     281n/plot clh64.val%n ! ! crossmarker5
     282
     283#---------------------- test ClRebin
     284delobjs *
     285zone
     286openppf clh64.ppf
     287cl2llcl clh64 llclh64
     288clrebin llclh64 clntu 10,0
     289
     290n/plot clntu.n%l ! ! crossmarker5
     291
     292n/plot llclh64.val%n  ! ! "crossmarker3"
     293n/plot clntu.clmean%l ! ! "same circlemarker5 red"
     294n/plot clntu.cllin%l  ! ! "same boxmarker5 blue"
     295n/plot clntu.clpar%l  ! ! "same trianglemarker5 violet"
     296
     297n/plot llclh64.val%n  ! ! "crossmarker3"
     298nt2d clntu l clmean 0 sclmean 1 " " "same circlemarker5 red thinline"
     299
     300n/plot llclh64.val%n  ! ! "crossmarker3"
     301nt2d clntu l cllin 0 scllin 1 " " "same boxmarker5 blue thinline"
     302
     303n/plot llclh64.val%n  ! ! "crossmarker3"
     304nt2d clntu l clpar 0 sclpar 1 " " "same trianglemarker5 violet thinline"
     305
     306#---------------------- test CrMaskMap et SetTypeMap
     307delobjs *
     308zone
    164309settypemap h
    165310crmapmask h8m 256 -20,20 100,130 0,1
     
    168313
    169314delobjs *
     315zone
    170316settypemap t
    171317crmapmask t8m 400 -20,20 100,130 0,1
     
    173319disp t8m zoom/4
    174320
    175 #----------------------
    176 delobjs *
     321#---------------------- test CrMaskFrMap
     322delobjs *
     323zone
    177324openfits h8.fits
    178325crmaskfrmap h8fm 256 h8 -1.e-30,1.e-30 0,1
     
    181328
    182329delobjs *
     330zone
    183331openppf t8.ppf
    184332crmaskfrmap t8fm 400 t8 -1.e-30,1.e-30 0,1
     
    186334disp t8fm zoom/4
    187335
    188 #----------------------
    189 delobjs *
     336#---------------------- test MaskMap
     337delobjs *
     338zone
    190339openfits h8.fits
    191340openfits h8m.fits
     
    194343
    195344delobjs *
     345zone
    196346openppf t8.ppf
    197347openppf t8m.ppf
     
    200350
    201351delobjs *
     352zone
    202353openfits h8.fits
    203354openppf t8m.ppf
     
    206357
    207358delobjs *
     359zone
    208360openppf t8.ppf
    209361openfits h8m.fits
     
    211363disp t8 zoom/4
    212364
    213 #----------------------
    214 delobjs *
     365#---------------------- test MapCover
     366delobjs *
     367zone
    215368openfits h8m.fits
    216369mapcover h8m 0.9,1 couvh8m
     
    218371
    219372delobjs *
     373zone
    220374openppf t8m.ppf
    221375mapcover t8m 0.9,1 couvt8m
    222376echo $couvt8m
    223377
    224 #----------------------
    225 delobjs *
     378#---------------------- test Map2Local
     379delobjs *
     380zone
    226381openfits h8.fits
    227382map2local h8 h8loc 200,300 20,30 0,90 !
     
    231386
    232387delobjs *
     388zone
    233389openppf t8.ppf
    234390map2local t8 t8loc 200,300 20,30 0,90 !
     
    237393disp t8loc
    238394
    239 #----------------------
    240 delobjs *
    241 openfits h8.fits
    242 rename h8 h8save
    243 openfits h8.fits
    244 mapop h8 + h8save - h8save * h8save / h8save
     395#---------------------- test MapOper
     396delobjs *
     397zone
     398openfits h8.fits
     399cp h8 h8save
     400mapop h8save + h8 - h8 * h8 / h8
     401disp h8save zoom/4
     402mapop h8save - h8
     403disp h8save zoom/4
     404
     405delobjs *
     406zone
     407openppf t8.ppf
     408cp t8 t8save
     409mapop t8save + t8 - t8 * t8 / t8
     410disp t8save zoom/4
     411mapop t8save - t8
     412disp t8save zoom/4
     413
     414delobjs *
     415zone
     416openfits h8.fits
     417openppf t8.ppf
     418mapop h8 - t8
    245419disp h8 zoom/4
    246420
    247 openppf t8.ppf
    248 rename t8 t8save
    249 openppf t8.ppf
    250 mapop t8 + t8save - t8save * t8save / t8save
    251 disp t8 zoom/4
    252 
    253 mapop h8 - t8save
    254 disp h8 zoom/4
    255 mapop t8 - h8save
    256 disp t8 zoom/4
    257 
    258 #----------------------
    259 delobjs *
     421#---------------------- test MapStat
     422delobjs *
     423zone
    260424openfits h8.fits
    261425openfits h8fm.fits
     
    271435mapstat t8 t8fm mean sig
    272436echo mean=$mean sig=$sig
     437
     438#---------------------- test VarOper
     439varop 1000 + 100 result
     440  echo result is $result
     441varop 1000 - 100 result
     442  echo result is $result
     443varop 1000 * 100 result
     444  echo result is $result
     445varop 1000 / 100 result
     446  echo result is $result
     447varop 1000 / 0 result
     448  echo result is $result
     449
     450varop 1000 sqrt ! result
     451  echo result is $result
     452varop -1000 sqrt ! result
     453  echo result is $result
     454varop 10 exp ! result
     455  echo result is $result
     456varop 1000 log ! result
     457  echo result is $result
     458varop -1000 log ! result
     459  echo result is $result
     460varop 1000 log10 ! result
     461  echo result is $result
     462varop -1000 log10 ! result
     463  echo result is $result
     464
     465varop 1000 pow 2 result
     466  echo result is $result
     467varop 1000 pow 1000 result
     468  echo result is $result
Note: See TracChangeset for help on using the changeset viewer.