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

Last change on this file since 2177 was 2177, checked in by cmv, 23 years ago

+ de fonctionnalites cmv 11/8/02
File size: 8.8 KB
Line 
1#---------------------- Create Maps and Save
2shell rm -f h4*.fits t4*.ppf h8*.fits t8*.ppf clh*.ppf clt*.ppf almh*.ppf almt*.ppf
3delobjs *
4c++import SkyMap
5c++exec \
6 double theta,phi; \
7 SphereHEALPix<r_4> h4(256); \
8 for(int_4 i=0;i<h4.NbPixels();i++) { \
9 h4.PixThetaPhi(i,theta,phi); \
10 h4(i) = 0.; \
11 double v = (M_PI_2-theta)*(M_PI_2-theta)/(0.50*0.50); \
12 v += (phi-M_PI)*(phi-M_PI)/(M_PI*M_PI/2.); \
13 v = 100.*exp(-v/2.); \
14 v += 10.*NorRand(); \
15 if(phi<20.*M_PI/180.||phi>340.*M_PI/180.) continue; \
16 if(phi<190.*M_PI/180.&&phi>170.*M_PI/180.) continue; \
17 h4(i) = v; \
18 } \
19 KeepObj(h4); \
20 SphereThetaPhi<r_4> t4(400); \
21 for(int_4 i=0;i<t4.NbPixels();i++) { \
22 t4.PixThetaPhi(i,theta,phi); \
23 t4(i) = 0.; \
24 double v = (M_PI_2-theta)*(M_PI_2-theta)/(0.50*0.50); \
25 v += (phi-M_PI)*(phi-M_PI)/(M_PI*M_PI/2.); \
26 v = 100.*exp(-v/2.); \
27 v += 10.*NorRand(); \
28 if(phi<20.*M_PI/180.||phi>340.*M_PI/180.) continue; \
29 if(phi<190.*M_PI/180.&&phi>170.*M_PI/180.) continue; \
30 t4(i) = v; \
31 } \
32 KeepObj(t4); \
33 cout<<"End of Job"<<endl;
34
35savefits h4 !h4.fits
36saveppf t4 t4.ppf
37
38zone 1 2
39disp h4 zoom/4
40disp t4 zoom/4
41
42#---------------------- test Map2Double
43delobjs *
44zone
45openfits h4.fits
46
47map2double h4
48rename h4 h8
49savefits h8 !h8.fits
50
51delobjs *
52zone
53openppf t4.ppf
54map2double t4
55rename t4 t8
56saveppf t8 t8.ppf
57
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
72openfits h4.fits
73openfits h8.fits
74openppf t4.ppf
75openppf t8.ppf
76typemap h4
77typemap h8
78typemap t4
79typemap t8
80
81#---------------------- test Map2Map
82delobjs *
83zone
84openfits h8.fits
85openppf t8.ppf
86map2map h8 t
87map2map t8 h
88
89#---------------------- test MapMult
90delobjs *
91zone
92openfits h8.fits
93openppf t8.ppf
94mapmult h8 1000
95mapmult t8 1000
96
97#---------------------- test MapProj
98delobjs *
99zone
100openfits h8.fits
101maproj h8 h8p
102mapop h8p - h8
103disp h8p zoom/4
104
105maproj h8 h8p64 64
106disp h8p64 zoom/4
107savefits h8p64 !h8p64.fits
108
109delobjs *
110zone
111openppf t8.ppf
112maproj t8 t8p
113mapop t8p - t8
114disp t8p zoom/4
115
116maproj t8 t8p100 100
117disp t8p100 zoom/4
118saveppf t8p100 t8p100.ppf
119
120#---------------------- test Map2Cl
121delobjs *
122zone
123openfits h8p64.fits
124map2cl h8p64 clh64
125saveppf clh64 clh64.ppf
126n/plot clh64.n*(n+1)*val%n ! ! crossmarker5
127
128map2cl h8p64 clh64 50
129n/plot clh64.n*(n+1)*val%n ! ! crossmarker5
130
131delobjs *
132zone
133openppf t8p100.ppf
134map2cl t8p100 clt100 191
135saveppf clt100 clt100.ppf
136n/plot clt100.n*(n+1)*val%n ! ! crossmarker5
137
138#---------------------- test Map2Alm
139delobjs *
140zone
141openfits h8p64.fits
142map2alm h8p64 almh64
143saveppf almh64 almh64.ppf
144disp almh64
145
146delobjs *
147zone
148openppf t8p100.ppf
149map2alm t8p100 almt100 191
150saveppf almt100 almt100.ppf
151disp almt100
152
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
177openppf clh64.ppf
178settypemap h
179cl2map clh64 h8p64frcl 64
180savefits h8p64frcl !h8p64frcl.fits
181disp h8p64frcl
182
183delobjs *
184zone
185openppf clt100.ppf
186settypemap t
187cl2map clt100 t8p100frcl 100
188saveppf t8p100frcl t8p100frcl.ppf
189disp t8p100frcl
190
191#---------------------- test Alm2Map et SetTypeMap
192delobjs *
193zone
194openppf almh64.ppf
195settypemap h
196alm2map almh64 h8p64fralm 64
197savefits h8p64fralm !h8p64fralm.fits
198disp h8p64fralm
199
200delobjs *
201zone
202openppf almt100.ppf
203settypemap t
204alm2map almt100 t8p100fralm 100
205saveppf t8p100fralm t8p100fralm.ppf
206disp t8p100fralm
207
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
309settypemap h
310crmapmask h8m 256 -20,20 100,130 0,1
311savefits h8m !h8m.fits
312disp h8m zoom/4
313
314delobjs *
315zone
316settypemap t
317crmapmask t8m 400 -20,20 100,130 0,1
318saveppf t8m t8m.ppf
319disp t8m zoom/4
320
321#---------------------- test CrMaskFrMap
322delobjs *
323zone
324openfits h8.fits
325crmaskfrmap h8fm 256 h8 -1.e-30,1.e-30 0,1
326savefits h8fm !h8fm.fits
327disp h8fm zoom/4
328
329delobjs *
330zone
331openppf t8.ppf
332crmaskfrmap t8fm 400 t8 -1.e-30,1.e-30 0,1
333saveppf t8fm t8fm.ppf
334disp t8fm zoom/4
335
336#---------------------- test MaskMap
337delobjs *
338zone
339openfits h8.fits
340openfits h8m.fits
341maskmap h8 h8m
342disp h8 zoom/4
343
344delobjs *
345zone
346openppf t8.ppf
347openppf t8m.ppf
348maskmap t8 t8m
349disp t8 zoom/4
350
351delobjs *
352zone
353openfits h8.fits
354openppf t8m.ppf
355maskmap h8 t8m
356disp h8 zoom/4
357
358delobjs *
359zone
360openppf t8.ppf
361openfits h8m.fits
362maskmap t8 h8m
363disp t8 zoom/4
364
365#---------------------- test MapCover
366delobjs *
367zone
368openfits h8m.fits
369mapcover h8m 0.9,1 couvh8m
370echo $couvh8m
371
372delobjs *
373zone
374openppf t8m.ppf
375mapcover t8m 0.9,1 couvt8m
376echo $couvt8m
377
378#---------------------- test Map2Local
379delobjs *
380zone
381openfits h8.fits
382map2local h8 h8loc 200,300 20,30 0,90 !
383disp h8loc
384map2local h8 h8loc 200,300 20,30 0,90 ! 45
385disp h8loc
386
387delobjs *
388zone
389openppf t8.ppf
390map2local t8 t8loc 200,300 20,30 0,90 !
391disp t8loc
392map2local t8 t8loc 200,300 20,30 0,90 ! 45
393disp t8loc
394
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
419disp h8 zoom/4
420
421#---------------------- test MapStat
422delobjs *
423zone
424openfits h8.fits
425openfits h8fm.fits
426mapstat h8 ! mean sig
427echo mean=$mean sig=$sig
428mapstat h8 h8fm mean sig
429echo mean=$mean sig=$sig
430
431openppf t8.ppf
432openppf t8fm.ppf
433mapstat t8 ! mean sig
434echo mean=$mean sig=$sig
435mapstat t8 t8fm mean sig
436echo 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 TracBrowser for help on using the repository browser.