source: trunk/external-alma/components/SpectralComponents/test/tSpectralFit.out @ 2980

Last change on this file since 2980 was 2980, checked in by Malte Marquarding, 10 years ago

Add a copy of casacore/components/SpectralComponents to external-alma directory to prepare for its removal from casacore-trunk. DOn;t activate in SConscript yet.

File size: 10.6 KB
Line 
1Test RecordInterface
2---------------------------------------------------
3toRecord
4fromRecord
5
6Test SpectralEstimate
7---------------------------------------------------
8Made spectrum with components:
9GAUSSIAN element:
10  Amplitude: 5, 0
11  Center:    45, 0
12  Sigma:     1.1, 0
13
14GAUSSIAN element:
15  Amplitude: 4, 0
16  Center:    50, 0
17  Sigma:     1.1, 0
18
19GAUSSIAN element:
20  Amplitude: 6, 0
21  Center:    55, 0
22  Sigma:     1.1, 0
23
24*** n 100
25*** y [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2.09099e-40, 3.34739e-35, 2.34499e-30, 7.18876e-26, 9.64375e-22, 5.66131e-18, 1.45434e-14, 1.63492e-11, 8.04274e-09, 1.73137e-06, 0.000163101, 0.00672359, 0.12129, 0.957476, 3.30757, 5.00013, 3.31295, 1.05451, 0.887271, 2.65278, 4.00036, 2.65413, 0.911529, 1.246, 3.97447, 6.00013, 3.96909, 1.14897, 0.145548, 0.00806831, 0.000195721, 2.07765e-06, 9.65129e-09, 1.9619e-11, 1.74521e-14, 6.79357e-18, 1.15725e-21, 8.62651e-26, 2.81398e-30, 4.01687e-35, 2.50918e-40, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
26Found (using rms, cutoff, minsig, q): (0.5, 1, 0.5, 2) 3 estimates
27Estimate 0: GAUSSIAN element:
28  Amplitude: 6.00013, 0
29  Center:    55.0715, 0
30  Sigma:     1.33398, 0
31
32Estimate 1: GAUSSIAN element:
33  Amplitude: 5.00013, 0
34  Center:    44.9152, 0
35  Sigma:     1.32309, 0
36
37Estimate 2: GAUSSIAN element:
38  Amplitude: 4.00036, 0
39  Center:    50.0043, 0
40  Sigma:     1.26335, 0
41
42*** res [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2.09099e-40, 3.34739e-35, 2.34499e-30, 7.18876e-26, 9.64375e-22, 5.66131e-18, 1.45434e-14, 1.63492e-11, 8.04274e-09, 1.73137e-06, 0.000163101, 0.00672359, 0.12129, 0.957476, 3.30757, 5.00013, 3.31295, 1.05451, 0.887271, 2.65278, 4.00036, 2.65413, 0.911529, 1.246, 3.97447, 6.00013, 3.96909, 1.14897, 0.145548, 0.00806831, 0.000195721, 2.07765e-06, 9.65129e-09, 1.9619e-11, 1.74521e-14, 6.79357e-18, 1.15725e-21, 8.62651e-26, 2.81398e-30, 4.01687e-35, 2.50918e-40, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
43Minimum, Maximum residuals: -0.961305, 0.0086807
44Window (using rms, cutoff, minsig, q): (0.5, 1, 0.5, 2) 0 estimates
45Found (using rms, cutoff, minsig, q): (0.5, 1, 0.5, 1) 3 estimates
46Estimate 0: GAUSSIAN element:
47  Amplitude: 6.00013, 0
48  Center:    54.9829, 0
49  Sigma:     1.17129, 0
50
51Estimate 1: GAUSSIAN element:
52  Amplitude: 5.00013, 0
53  Center:    45.0202, 0
54  Sigma:     1.17355, 0
55
56Estimate 2: GAUSSIAN element:
57  Amplitude: 4.00036, 0
58  Center:    50.0182, 0
59  Sigma:     1.38534, 0
60
61Found (using rms, cutoff, minsig, q): (0, 0, 0, 1) 3 estimates
62Estimate 0: GAUSSIAN element:
63  Amplitude: 6.00013, 0
64  Center:    54.9829, 0
65  Sigma:     1.17129, 0
66
67Estimate 1: GAUSSIAN element:
68  Amplitude: 5.00013, 0
69  Center:    45.0202, 0
70  Sigma:     1.17355, 0
71
72Estimate 2: GAUSSIAN element:
73  Amplitude: 4.00036, 0
74  Center:    50.0182, 0
75  Sigma:     1.38534, 0
76
77Found (using rms, cutoff, minsig, q): (0, 0, 0, 2) 3 estimates
78Estimate 0: GAUSSIAN element:
79  Amplitude: 6.00013, 0
80  Center:    55.0715, 0
81  Sigma:     1.33398, 0
82
83Estimate 1: GAUSSIAN element:
84  Amplitude: 5.00013, 0
85  Center:    44.9152, 0
86  Sigma:     1.32309, 0
87
88Estimate 2: GAUSSIAN element:
89  Amplitude: 4.00036, 0
90  Center:    50.0043, 0
91  Sigma:     1.26335, 0
92
93Test SpectralFit
94---------------------------------------------------
95Spectral elements:
96GAUSSIAN element:
97  Amplitude: 1, 0
98  Center:    1402.5, 0
99  Sigma:     0.0195312, 0
100
101GAUSSIAN element:
102  Amplitude: 2.5, 0
103  Center:    1405, 0
104  Sigma:     0.0390625, 0
105
106GAUSSIAN element:
107  Amplitude: 0.5, 0
108  Center:    1407.5, 0
109  Sigma:     0.078125, 0
110
111GAUSSIAN element:
112  Amplitude: 5, 0
113  Center:    1403, 0
114  Sigma:     0.0292969, 0
115
116---------------------------------------------------
117Data for frequencies 1402.93 - 1403.08
1181402.93: 0.219133
1191402.94: -0.480645
1201402.95: -0.281389
1211402.96: 0.0423511
1221402.97: 0.0837382
1231402.98: 1.05514
1241402.99: 3.14053
1251403: 4.93017
1261403.01: 4.53475
1271403.02: 1.54477
1281403.03: 0.247699
1291403.04: -0.50962
1301403.05: 0.0446798
1311403.06: -0.0518962
1321403.07: -0.750935
1331403.08: 0.499586
134---------------------------------------------------
135Specify fitter on the 4 gaussians:
136GAUSSIAN element:
137  Amplitude: 1, 0
138  Center:    1402.5, 0
139  Sigma:     0.0195312, 0
140
141GAUSSIAN element:
142  Amplitude: 2.5, 0
143  Center:    1405, 0
144  Sigma:     0.0390625, 0
145
146GAUSSIAN element:
147  Amplitude: 0.5, 0
148  Center:    1407.5, 0
149  Sigma:     0.078125, 0
150
151GAUSSIAN element:
152  Amplitude: 5, 0
153  Center:    1403, 0
154  Sigma:     0.0292969, 0
155
156---------------------------------------------------
157Execute the fitter, can do (1:Y): 1
158The results:
159GAUSSIAN element:
160  Amplitude: 0.917797, 0.310674
161  Center:    1402.5, 0.00347673
162  Sigma:     0.00882989, 0.00343021
163Parameters: [0.917797, 1402.5, 0.00882989]
164Errors:     [0.310674, 0.00347673, 0.00343021]
165
166GAUSSIAN element:
167  Amplitude: 2.76117, 0.238196
168  Center:    1405, 0.00151172
169  Sigma:     0.0151761, 0.00151172
170Parameters: [2.76117, 1405, 0.0151761]
171Errors:     [0.238196, 0.00151172, 0.00151172]
172
173GAUSSIAN element:
174  Amplitude: 0.612907, 0.18968
175  Center:    1407.48, 0.0085527
176  Sigma:     0.0239303, 0.0085527
177Parameters: [0.612907, 1407.48, 0.0239303]
178Errors:     [0.18968, 0.0085527, 0.0085527]
179
180GAUSSIAN element:
181  Amplitude: 5.24752, 0.272235
182  Center:    1403, 0.000695977
183  Sigma:     0.0116179, 0.000695945
184Parameters: [5.24752, 1403, 0.0116179]
185Errors:     [0.272235, 0.000695977, 0.000695945]
186
187>>>
188Number of iterations: 29
189<<<
190---------------------------------------------------
191GAUSSIAN element:
192  Amplitude: 1, 0
193  Center:    1402.5, 0
194  Sigma:     0.0195312, 0
195
196GAUSSIAN element:
197  Amplitude: 2.5, 0
198  Center:    1405, 0
199  Sigma:     0.0390625, 0
200
201GAUSSIAN element:
202  Amplitude: 0.5, 0
203  Center:    1407.5, 0
204  Sigma:     0.078125, 0
205
206GAUSSIAN element:
207  Amplitude: 5, 0 (fixed)
208  Center:    1403, 0
209  Sigma:     0.0292969, 0
210
211---------------------------------------------------
212Execute the fitter again with a fixed ampl, can do (1:Y): 1
213The results:
214GAUSSIAN element:
215  Amplitude: 0.917772, 0.310662
216  Center:    1402.5, 0.00347618
217  Sigma:     0.00883038, 0.00342961
218Parameters: [0.917772, 1402.5, 0.00883038]
219Errors:     [0.310662, 0.00347618, 0.00342961]
220
221GAUSSIAN element:
222  Amplitude: 2.76117, 0.23818
223  Center:    1405, 0.00151162
224  Sigma:     0.0151761, 0.00151162
225Parameters: [2.76117, 1405, 0.0151761]
226Errors:     [0.23818, 0.00151162, 0.00151162]
227
228GAUSSIAN element:
229  Amplitude: 0.61293, 0.18968
230  Center:    1407.48, 0.00855088
231  Sigma:     0.0239286, 0.00855088
232Parameters: [0.61293, 1407.48, 0.0239286]
233Errors:     [0.18968, 0.00855088, 0.00855088]
234
235GAUSSIAN element:
236  Amplitude: 5, 0 (fixed)
237  Center:    1403, 0.000741387
238  Sigma:     0.0119707, 0.00060533
239Parameters: [5, 1403, 0.0119707]
240Errors:     [0, 0.000741387, 0.00060533]
241
242>>>
243Number of iterations: 26
244<<<
245---------------------------------------------------
246Different start values:
247GAUSSIAN element:
248  Amplitude: 0.9, 0
249  Center:    1402.5, 0
250  Sigma:     0.0185547, 0
251
252GAUSSIAN element:
253  Amplitude: 2.25, 0
254  Center:    1405, 0
255  Sigma:     0.0371094, 0
256
257GAUSSIAN element:
258  Amplitude: 0.45, 0
259  Center:    1407.5, 0
260  Sigma:     0.0742188, 0
261
262GAUSSIAN element:
263  Amplitude: 4.5, 0
264  Center:    1403, 0
265  Sigma:     0.027832, 0
266
267Execute the fitter, can do (1:Y): 1
268The results:
269GAUSSIAN element:
270  Amplitude: 0.917785, 0.310679
271  Center:    1402.5, 0.00347653
272  Sigma:     0.00883013, 0.00342999
273
274GAUSSIAN element:
275  Amplitude: 2.76117, 0.238196
276  Center:    1405, 0.00151172
277  Sigma:     0.0151761, 0.00151172
278
279GAUSSIAN element:
280  Amplitude: 0.61293, 0.189693
281  Center:    1407.48, 0.00855143
282  Sigma:     0.0239285, 0.00855143
283
284GAUSSIAN element:
285  Amplitude: 5.24752, 0.272235
286  Center:    1403, 0.000695977
287  Sigma:     0.0116179, 0.000695945
288
289>>>
290Number of iterations: 32
291<<<
292---------------------------------------------------
293Differences:
294Min difference: -1.09994, max: 1.01271, average: 0.00978336, sigma: 0.0100325
295---------------------------------------------------
296---------------------------------------------------
297Estimates:
298Found 20 estimates
299Estimate 0: GAUSSIAN element:
300  Amplitude: 4.93017, 0
301  Center:    307.195, 0
302  Sigma:     3.05174, 0
303
304Estimate 1: GAUSSIAN element:
305  Amplitude: 2.79525, 0
306  Center:    511.819, 0
307  Sigma:     3.29148, 0
308
309Estimate 2: GAUSSIAN element:
310  Amplitude: 0.906691, 0
311  Center:    766.653, 0
312  Sigma:     2.31423, 0
313
314Estimate 3: GAUSSIAN element:
315  Amplitude: 0.90563, 0
316  Center:    501.372, 0
317  Sigma:     3.60214, 0
318
319Estimate 4: GAUSSIAN element:
320  Amplitude: 0.861827, 0
321  Center:    839.949, 0
322  Sigma:     3.43708, 0
323
324Estimate 5: GAUSSIAN element:
325  Amplitude: 0.846625, 0
326  Center:    112.878, 0
327  Sigma:     2.28534, 0
328
329Estimate 6: GAUSSIAN element:
330  Amplitude: 0.751603, 0
331  Center:    871.157, 0
332  Sigma:     2.78939, 0
333
334Estimate 7: GAUSSIAN element:
335  Amplitude: 0.699624, 0
336  Center:    36.4464, 0
337  Sigma:     3.65743, 0
338
339Estimate 8: GAUSSIAN element:
340  Amplitude: 0.682235, 0
341  Center:    395.917, 0
342  Sigma:     3.5893, 0
343
344Estimate 9: GAUSSIAN element:
345  Amplitude: 0.663379, 0
346  Center:    450.12, 0
347  Sigma:     0.548263, 0
348
349Estimate 10: GAUSSIAN element:
350  Amplitude: 0.640783, 0
351  Center:    827.254, 0
352  Sigma:     3.35359, 0
353
354Estimate 11: GAUSSIAN element:
355  Amplitude: 0.616448, 0
356  Center:    245.555, 0
357  Sigma:     3.28022, 0
358
359Estimate 12: GAUSSIAN element:
360  Amplitude: 0.578568, 0
361  Center:    120.043, 0
362  Sigma:     1.45528, 0
363
364Estimate 13: GAUSSIAN element:
365  Amplitude: 0.569486, 0
366  Center:    663.971, 0
367  Sigma:     1.81295, 0
368
369Estimate 14: GAUSSIAN element:
370  Amplitude: 0.56377, 0
371  Center:    611.848, 0
372  Sigma:     2.99822, 0
373
374Estimate 15: GAUSSIAN element:
375  Amplitude: 0.561451, 0
376  Center:    937.91, 0
377  Sigma:     2.32834, 0
378
379Estimate 16: GAUSSIAN element:
380  Amplitude: 0.54666, 0
381  Center:    336.559, 0
382  Sigma:     1.8761, 0
383
384Estimate 17: GAUSSIAN element:
385  Amplitude: 0.524891, 0
386  Center:    443.644, 0
387  Sigma:     0.809077, 0
388
389Estimate 18: GAUSSIAN element:
390  Amplitude: 0.523377, 0
391  Center:    645.266, 0
392  Sigma:     3.01063, 0
393
394Estimate 19: GAUSSIAN element:
395  Amplitude: 0.515985, 0
396  Center:    446.61, 0
397  Sigma:     1.18317, 0
398
399---------------------------------------------------
400Limited number of estimates:
401Found 4 estimates
402Estimate 0: GAUSSIAN element:
403  Amplitude: 4.93017, 0
404  Center:    307.195, 0
405  Sigma:     3.05174, 0
406
407Estimate 1: GAUSSIAN element:
408  Amplitude: 2.79525, 0
409  Center:    511.819, 0
410  Sigma:     3.29148, 0
411
412Estimate 2: GAUSSIAN element:
413  Amplitude: 0.906691, 0
414  Center:    766.653, 0
415  Sigma:     2.31423, 0
416
417Estimate 3: GAUSSIAN element:
418  Amplitude: 0.90563, 0
419  Center:    501.372, 0
420  Sigma:     3.60214, 0
421
422---------------------------------------------------
Note: See TracBrowser for help on using the repository browser.