-
Notifications
You must be signed in to change notification settings - Fork 5
/
changelog.txt
6996 lines (4808 loc) · 210 KB
/
changelog.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
From 7e1a8a2 (master):
=====================
In README.rst,
(*) Update with citation and acknowledgement information
From daa77da (master):
=====================
In prisim/__init__.py,
(*) Update PRISim version to 2.2.1
From 7a06295 (master):
=====================
In README.rst,
(*) Update PRISim description
(*) Update links for obtaining PRISim data
From ceaf7a2 (mwa), 7eb5233 (master):
====================================
(*) Merge 'mwa' branch into 'master'
From 76a01b6 (master):
=====================
In scripts/setup_prisim_data.py,
(*) Update outdated and inaccurate comment about the purpose of the
script
From 38132bc (master):
=====================
In prisim/__init__.py,
(*) Update PRISim version to 2.2
From d3a9bed (mwa):
==================
In prisim/examples/ioparms/data_setup_parms.yaml,
(*) Update Google Drive link for PRISim data
In README.rst,
(*) Update Google Drive link for gzipped PRISim data
From dacf5b7 (mwa):
==================
In prisim/interferometry.py,
(*) In function getBaselineInfo(), deprecate support for 'MWA-128T'
layout which is now replaced with 'MWA-I-128T'
From 38132bc (master):
=====================
In prisim/interferometry.py,
(*) In function getBaselineInfo(), add support for updated MWA
layouts using default antenna locations. Update member function
docstring accordingly
In prisim/examples/simparms/defaultparms.yaml,
(*) Add support for updated MWA layouts using default antenna
locations. Update information about latitude and longitude of MWA
From 1e4267d (closure-adv):
==========================
In prisim/examples/codes/BispectrumPhase/,
(*) Update *.ipynb codes corresponding to the bispectrum phase math
paper
From 40c4d19 (master):
=====================
In prisim/interferometry.py,
(*) Import astroutils
(*) In member function save() of class InterferometerArray, save the
githash of the AstroUtils and PRISim versions in the HDF5 output file
(*) In member function __init__() of class InterferometerArray, read
or assign the prisim_githash and astroutils_githash attributes. Update
class and member function docstring accordingly
From 04f7905 (master):
=====================
In prisim/examples/simparms/defaultparms.yaml,
(*) Change default value of 'wait_after_run' from false to true
In prisim/interferometry.py,
(*) In member function observe() of class InterferometerArray, remove
a redundant set of statements
From 9eb2a00 (master):
=====================
In prisim/__init__.py,
(*) Update PRISim version
From fa5a832 (master):
=====================
(*) `git rm -r main` folder from the repository
From d3592b0 (closure-adv):
==========================
In prisim/examples/codes/BispectrumPhase/,
(*) Update *.ipynb codes corresponding to the bispectrum phase data
analysis demo paper
From b3cbcec (closure-adv):
==========================
In prisim/examples/codes/BispectrumPhase/,
(*) Update *.ipynb codes corresponding to the bispectrum phase math
paper
From cf6f5df (parabolic-cylinder):
=================================
(*) In member function duplicate_measurements() of class
InterferometerArray, simplify the process to add a missing baseline
label key
From eba5016 (parabolic-cylinder):
=================================
In prisim/interferometry.py,
(*) In __init__() member function of class InterferometerArray, accept
input labels as a numpy array, and convert if it is a list to a numpy
recArray using datatype of antenna labels in attribute layout
(*) In member function duplicate_measurements() of class
InterferometerArray, do more rigorous type checking of baseline labels
to avoid bugs due to slight mismatch in datatypes, avoid type
promotion errors, and append without errors
(*) In member function save() of class InterferometerArray, determine
baseline label datatypes from antenna label datatypes instead of
finding max string length because the latter may be a smaller string
if some baselines have been skipped
From 01ca9c4 (parabolic-cylinder):
=================================
In prisim/interferometry.py,
(*) Remove some commented lines in member function observe() of class
InterferometerArray
In scripts/run_prisim.py,
(*) Pass whole skymodel instead of only a subset in call to observe()
member function of class InterferometerArray when the MPI chunking is
based on frequency
From 11698dd (parabolic-cylinder):
=================================
In prisim/interferometry.py,
(*) In function antenna_power(), generate sky model spectrum using
explicit call to 'pchip' interpolation method
(*) In member function observe() of class SkyModel, generate sky
model spectrum using explicit call to 'pchip' interpolation method
In scripts/run_prisim.py,
(*) Replace the use of 'fgparm' key with new key 'skyparm' which was
missed in the merge of master into parabolic-cylinder branch
(*) Move the creation of organized directory structure to an earlier
location in the script
(*) Include 'gsm2008' and 'gsm2016' in checks for accepted sky model
(*) Compute sky model and the sky locations only for the rank=0
process and broadcast it to other processes
(*) Always store sky model computed under rank=0
(*) Make 'pyGSM' for the sky model work
(*) Make creation of sky model subsets explicit in calls to
observe() member function of class InterferometerArray
(*) Replace chans_chunk_indices with chans_chunk and explicitly use
'pchip' interpolation method in calls to observe() member function
of class InterferometerArray
From c3f2978 (parabolic-cylinder):
=================================
In prisim/interferometry.py,
(*) Call function rectangle_generator() in function
getBaselineInfo() when array_layout is set to 'CHIME' which was
missed in the merge of master with the parabolic-cylinder branch.
(*) Add comments in observe() member function of class
InterferometerArray
(*) Add comments in function baseline_generator()
From 1a2b107 (parabolic-cylinder):
=================================
In scripts/run_prisim.py,
(*) Remove some leftover code from the merge with the branch
parabolic-cylinder that used old timestamps_JD
(*) Check for presence of PyGSM module and start a set_trace() if
absent
In prisim/interferometry.py,
(*) In member function observe() of class InterferometerArray, if
the previous skymodel has to be updated, delete the contents of
the previous iteration and then store the current contents
From 5378176 (master), 3f7d3bc (parabolic-cylinder):
===================================================
(*) Merge master branch onto parabolic-cylinder
From 6cd37a7 (master):
=====================
(*) Add 21cmforest codes under prisim/examples/codes/21cmforest/
From 6418179 (master):
=====================
In scripts/run_prisim.py,
(*) Add support for processing NVSS catalog as sky model
From 7daa921 (gmrt):
===================
In prisim/examples/simparms/defaultparms.yaml,
(*) Add support for 'ugmrt' as a telescope id
In prisim/examples/simparms/defaultparms_dev.yaml,
(*) Add support for 'ugmrt' as a telescope id
In scripts/run_prisim.py,
(*) Add support for 'ugmrt' as a telescope id
In prisim/interferometry.py,
(*) Update docstrings to indicate support for 'ugmrt' as a telescope id
In prisim/primary_beams.py,
(*) In function primary_beam_generator(), add support for 'ugmrt' as a
telescope id
(*) In function GMRT_primary_beam(), add input keyword instrument that
can support 'gmrt' (default) or 'ugmrt' and provide corresponding
polynomial coefficients for the primary beam. Update function docstring
accordingly
(*) In function GMRT_primary_beam(), check for NaN values and raise a
ValueError if NaN values found
From 8061bbd (master):
=====================
In prisim/primary_beams.py,
(*) In functions VLA_primary_beam() and GMRT_primary_beam(), check if
the returned primary beam exceeds unity and raise ValueError if it
does due to the invalidity of the polynomial equations.
From a2d544f (master):
=====================
In prisim/interferometry.py,
(*) In member function observe() of class InterferometerArray, pass
the right shape of inputs in calls to altaz2hadec()
From bb083bc (closure-adv):
==========================
(*) Add folder prisim/examples/codes/ and its contents to repo
In setup.py,
(*) Include folder prisim/examples/codes/ and its contents during
installation
From acdf710 (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) Modify function save_CPhase_cross_power_spectrum() to save also the
output of function incoherent_cross_power_spectrum_average()
(*) Modify function read_CPhase_cross_power_spectrum() to also read the
what is saved using function save_CPhase_cross_power_spectrum() that
came as the output of function incoherent_cross_power_spectrum_average()
From c8b09de (master):
=====================
In prisim/interferometry.py,
(*) In member function append-settings() of class ROI_parameters, there
was no option to handle achromatic beams when the beams were to be
determined analytically. Now it accepts beam chromaticity parameters in
input roi_info and computes analystical beam at specified frequency and
propagates it to all frequencies
In scripts/run_prisim.py,
(*) Beam chromaticity information in case of analytical beams was not
being passed on to the ROI parameters in call to roi.append_settings()
which is now passed
(*) Minor changes to chunk completion time by process
In prisim/delay_spectrum.py,
(*) In member function beam3Dvol() of class DelayPowerSpectrum, update
code block to determine achromatic analytical beams.
In prisim/primary_beams.py,
(*) In function primary_beam_generator(), ensure the input frequency
is converted to a numpy array of non-zero length even if input as a
scalar. Update function docstring accordingly
From ba17ed5 (closure-adv):
==========================
(*) Add prisim/scriptUtils/write_PRISim_bispectrum_phase_to_npz_util.py
to repo
In prisim/scriptUtils/replicatesim_util.py,
(*) Minor cosmetic changes
In scripts/write_PRISim_bispectrum_phase_to_npz.py,
(*) Move some imported modules and main code block to reside in
prisim/scriptUtils/write_PRISim_bispectrum_phase_to_npz_util.py
From a203161 (replicatesim):
===========================
(*) Add folder prisim/scriptUtils/ to repo
(*) Add prisim/scriptUtils/replicatesim_util.py to repo
In scripts/replicatesim.py,
(*) Move some imported modules and function replicate() to reside in
prisim/scriptUtils/replicatesim_util.py
From 234d167 (master):
=====================
In scripts/replicatesim.py,
(*) Remove unused import of h5py
(*) Restructure __main__ block to call a function replicate() by
passing the parameters
(*) Insert an option to go into debug mode before calling replicate()
In examples/simparms/replicatesim.yaml,
(*) Add option to enter debug mode before calling replicate()
From bcad114 (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In member function smooth_in_tbins() of class ClosurePhase, update
handling of LST bins and weights when specified LST bin size is
smaller than the LST resolution in the data
(*) In member function subsample_differencing() of class ClosurePhase,
update handling of LST bins and weights when specified LST bin size is
smaller than the LST resolution in the data
From 683604f (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In member function smooth_in_tbins() of class ClosurePhase, update
handling of the case when the specified LST bin size is smaller than
the LST resolution in the data, namely, do no binning and pass the
existing bin size, and other keys
From f11e09b (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In member function subsample_differencing() of class ClosurePhase,
handle the computing of eicp_diff when specified LST bin size is
smaller than the LST resolution in the data, namely, do no binning
and pass the existing bin size, and other keys
From bc5041b (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In member function smooth_in_tbins() of class ClosurePhase, handle
the case when the specified LST bin size is smaller than the LST
resolution in the data, namely, do no binning and pass the existing
bin size, and other keys
(*) In member function subsample_differencing() of class ClosurePhase,
handle the case when the specified LST bin size is smaller than the
LST resolution in the data, namely, do no binning and pass the
existing bin size, and other keys
From a3b9a25 (closure-adv):
==========================
In scripts/write_PRISim_bispectrum_phase_to_npz.py,
(*) Import prisim and set default path
(*) Correct the script description and default path of input YAML file
In prisim/examples/ioparms/model_bispectrum_phase_to_npz_parms.yaml,
(*) Update input parameters with updated defaults and updated
descriptions
From 47b70df (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) Remove commented lines in function
incoherent_cross_power_spectrum_average()
From dd0d710 (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In function incoherent_cross_power_spectrum_average(), fix bug to
handle when diagweights in the input dictionaries could either be a
dictionary or a numpy array
From 73f7a47 (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In function incoherent_kbin_averaging(), add a progressbar
From f783389 (master):
=====================
In prisim/bispectrum_phase.py,
(*) In function incoherent_cross_power_spectrum_average(), fix bug in
estimating axes to average in case of subsample differences
From aadeca3 (noise):
====================
In prisim/bispectrum_phase.py,
(*) In member function subsample_differencing() of class ClosurePhase,
fix bug in creation of the outshape
From 8700a47 (noise):
====================
(*) Add prisim/examples/simparms/noise_update_parms.yaml that serves
as a template for scripts/update_PRISim_noise.py
(*) Update .gitignore
In scripts/update_PRISim_noise.py,
(*) Update default location of noise parameters file to
prisim/examples/simparms/noise_update_parms.yaml
From 518ca3c (pyuvdata):
=======================
In scripts/make_redundant_visibilities.py,
(*) Fix bug in reading and saving of attribute projected_baselines
In scripts/update_PRISim_noise.py,
(*) Fix bug in reading and saving of attribute projected_baselines
In scripts/write_PRISim_visibilities.py,
(*) Fix bug in reading and saving of attribute projected_baselines
From fb71114 (noise):
====================
(*) Add scripts/update_PRISim_noise.py to repo, which updates noise
in a PRISim output
In prisim/interferometry.py,
(*) In member function save() of class InterferometerArray, fix bug
in saving Tsysinfo parameters
(*) In member function __init__() of class InterferometerArray, fix
bug in reading Tsysinfo parameters
In scripts/write_PRISim_visibilities.py,
(*) Minor change in description string
From 08d74fd (pyuvdata):
=======================
In scripts/make_redundant_visibilities.py,
(*) Import script scripts/write_PRISim_visibilities.py as a module
(*) Update the __main__() function to set up parameters to be passed
to save() function in imported scripts/write_PRISim_visibilities.py
module to save the outputs in the requested formats. This removes the
duplication of saving output in this script and instead calls the
imported script
(*) Fix minor typo in a printed message of a raised exception
In scripts/write_PRISim_visibilities.py,
(*) Fix minor typo in a printed message of a raised exception
From 7355c0a (pyuvdata):
=======================
(*) Add scripts/write_PRISim_visibilities.py to repo
In scripts/make_redundant_visibilities.py,
(*) Introduce a conditional statement to fix previous erroneous
outputs of attribute projected_baselines in the wrong shape.
From 6a21515 (ephemeris):
========================
In scripts/run_prisim.py,
(*) Add hooks to improve ephemeris in conversion of ROI center alt-az
to RA-Dec, but still commented just to act as a placeholder
From d2d33ae (ephemeris):
========================
In prisim/interferometry.py,
(*) In member function append_settings() of class ROI_parameters,
remove commented lines in the previous commit
From c46e608 (ephemeris):
========================
In prisim/interferometry.py,
(*) In member function append_settings() of class ROI_parameters,
initialize sky coordinates in FK5 frame instead of ICRS to fix bug in
ephemeris and to be consistent with ephemeris calculations in
scripts/run_prisim.py
(*) In member function append_settings() of class ROI_parameters,
comment out the previous initialization of sky coords in ICRS frame
In scripts/run_prisim.py,
(*) Add a comment to clarify the intention behind computing of the
variable roi_center_hadec
From 07e1117 (master):
=====================
In scripts/run_prisim.py,
(*) Modify DSM to be read from the default file independent of the
'freq' parameter and simply use spectral index to extrapolate to
any frequency. Make this update for skymodel when set to 'dsm', 'gsm'
or 'usm'
From dade112 (pyuvdata):
=======================
In prisim/interferometry.py,
(*) In member function save() of class InterferometerArray, fix bug
in writing the correct attribute for projected_baselines, wherein
earlier it only wrote baselines attribute wrongly in place of the
projected_baselines attribute
In scripts/make_redundant_visibilities.py,
(*) Insert a comment that can temporarily fix the bug that was saving
incorrect projected_baselines (see above), which will be a
placeholder until all simulations have got this bug fixed in them
From d9ddcd8 (pyuvdata):
=======================
In scripts/make_redundant_visibilities.py,
(*) Update to save redundant visibilities in UVFITS format
(*) Start work on saving redundant visibilities in uvh5 format
From e64b513 (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In member function FT() of class ClosurePhaseDelaySpectrum, add
ability to handle if the reference visibility is provided only at a
single LST, in which case this scaling applies to all the LST
From d6090ea (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) Import primary_beams module
(*) In member function beam3Dvol() of class ClosurePhaseDelaySpectrum,
fix bug in computation of analytic power pattern
In prisim/primary_beams.py,
(*) In function primary_beam_generator(), pass telecope orientation
as a numpy array in call to function airy_disk_pattern()
From 892760b (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) In member function smooth_in_tbins() of class ClosurePhase, add
ability to handle when LST binning is not requested or if there is
only one LST.
(*) In member function subsample_differencing() of class ClosurePhase,
add ability to handle when LST binning is not requested or if there is
only one LST.
From 6d211aa (noise):
====================
In examples/simparms/replicatesim.yaml,
(*) Add input parameters eff_aprtr and eff_Q
(*) Change default value of Trx to 162 K from 50 K
In scripts/replicate_sim.py,
(*) Accept input parameters eff_aprtr and eff_Q from the YAML file
with template in examples/simparms/replicatesim.yaml
From 58b3903 (pyuvsim):
======================
In scripts/altsim_interface.py,
(*) In function simparms_from_pyuvsim_to_prisim(), change some
defaults, and fix some parameters
From 0d3c19c (noise):
====================
In examples/simparms/defaultparms.yaml,
(*) Add parameters eff_aprtr (for aperture efficiency), and eff_Q
(quantization efficiency) for enabling better estimates of thermal
noise RMS
In scripts/run_prisim.py,
(*) Read and process eff_aprtr, and eff_Q parameters from the input
YAML file and pass appropriate values while creating instances of
class InterferometerArray.
From fd82446 (master):
=====================
In prisim/interferometry.py,
(*) In function thermalNoiseRMS(), revise RMS calculation by removing
the factor 2 inside square root to account for a complex-valued rms
that is made of both real and imaginary parts. Earlier the presence
of the factor 2 inside the square root corresponds to only the real
or the imaginary rms. Update docstring accordingly with a reference
for the equation used from SIRA II
(*) In function generateNoise(), update docstring to make it explicit
that the input RMS corresponds to a complex value that comprises of
both real and imaginary values and therefore, only 1/sqrt(2) of that
goes into the real and imaginary values. Include a reference to the
equations in SIRA II
(*) In member function generate_noise() of class InterferometerArray,
revise RMS calculation by removing the factor 2 inside square root
to account for a complex-valued rms that is made of both real and
imaginary parts. Earlier the presence of the factor 2 inside the
square root corresponds to only the real or the imaginary rms. This
is consistent with other changes (see above) in this commit.
Update docstring accordingly with a reference for the equation used
from SIRA II.
From 9f5487a (master):
=====================
(*) Add script altsim_interface.py to scripts/ folder. This
script converts YAML used in pyuvsim to a YAML that can be used for
PRISim
From 1d10033 (closure-adv):
==========================
In prisim/bispectrum_phase.py,
(*) Add input keyword lst_format to function npz2hdf5() to allow
different unit/format inputs for LST. Update function docstring
accordingly
(*) Add input keyword lst_format to function loadnpz() to allow
different unit/format inputs for LST. Update function docstring
accordingly
From f84aef6 (ephemeris):
========================
In scripts/run_prisim.py,
(*) Remove recently commented lines related to recent changes in
handling ephemeris
In prisim/interferometry.py,
(*) Remove recently commented lines related to recent changes in
handling ephemeris in member function observe() of class
InterferometerArray
From baa9f04 (ephemeris):
========================
In scripts/run_prisim.py,
(*) Minor changes in the way earth location of the osberver is
initialized
(*) Sky model is precessed to the observing epoch
(*) Comment out some redundant sky coordinate conversions
In prisim/interferometry.py,
(*) In member function observe() of class InterferometerArray, take
precessed version of sky model to observing epoch before calculating
altaz coordinates.
From e88db86 (master):
=====================
In scripts/run_prisim.py,
(*) Use updated ephemeris calculations from
astroutils.ephemeris_timing.py
From 4e8154a (emptysim):
=======================
In prisim/interferometry.py,
(*) In member function append_settings() of class ROI_parameters,
initialize telescope and freq attributes if they are found to be None.
In scripts/run_prisim.py,
(*) If no sources are found to be simulated for the entire simulation,
one source is assumed to be present just for purposes of memory
estimates.
(*) Pass chans, telescope, and freq_scale inputs in call to member
function append_settings() of class ROI_parameters even in case of
no sources (empty simulation).
(*) Update code to handle a completely empty simulation.
From 2c562ed (ephemeris):
========================
In prisim/interferometry.py,
(*) Import Galactic from astropy.coordinates
(*) In member function append_settings() of class ROI_parameters,
add input keyword time_jd to pass time in JD for accurate ephemeris
calculation. Update and adjust member function docstring accordingly
(*) In member function append_settings() of class ROI_parameters,
fix bug in inaccurate calculation of Alt-Az and use the astropy
coordinate transformation instead to convert from RA-Dec to Alt-Az
using the location and time information.
In scripts/run_prisim.py,
(*) Pass time in JD as input in call to member function
append_settings() of class RI.ROI_parameters based on the changes
above
From d277b30 (BPS_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In function read_CPhase_cross_power_spectrum(), ensure reading of
the group name (which is a string) is converted to int for certain
cases as appropriate and to ensure comptability.
From fda8c3f (BPS_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) Add function incoherent_kbin_averaging() along with docstring, to
perform incoherent averaging of power spectra in k-bins and also to
compute power spectrum in Delta^2 units.
(*) Update function save_CPhase_cross_power_spectrum() to save cross
power spectrum info of closure phases in a dictionary to an external
HDF5 file.
(*) Update function read_CPhase_cross_power_spectrum() to read cross
power spectrum info of closure phases from an external HDF5 file into
a dictionary
From 55df5dc (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) Add function save_CPhase_cross_power_spectrum() to save cross
power spectrum info of closure phases in a dictionary to an external
HDF5 file.
(*) Add function read_CPhase_cross_power_spectrum() to read cross
power spectrum info of closure phases from an external HDF5 file into
a dictionary
From 7783d74 (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) Add docstring for function
incoherent_cross_power_spectrum_average()
From 72be647 (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In member function subsample_differencing() of class ClosurePhase,
enforce masked values to data difference of the two terms in eicp_diff
values
(*) Update incoherent averaging of subsample differences in function
incoherent_cross_power_spectrum_average()
From fd86c01 (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In save() member function of class InterferometerArray, fix bug
in writing of baseline labels.
From e6a7ebb (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In member function smooth_in_tbins() and subsample_differencing()
of class ClosurePhase, if the LST cross the 24-hour boundary, unwrap
them to remove the discontinuity.
From d59fce7 (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) Remove commented lines due to change made earlier in the way
diagonal offsets are stored in member functions
compute_power_spectrum() and compute_power_spectrum_uncertainty() of
class ClosurePhaseDelaySpectrum
From 5547b9e (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) Develop work on producing incoherent average of power spectrum
through function incoherent_cross_power_spectrum_average()
From 0172ece (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In member functions compute_power_spectrum() and
compute_power_spectrum_uncertainty() of class
ClosurePhaseDelaySpectrum, update diagoffsets to be a dictionary
instead of a list. Update respective docstrings accordingly
From 08df256 (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) Start work on producing incoherent average of power spectrum
through function incoherent_cross_power_spectrum_average()
From 6a887cb (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In member function compute_power_spectrum() of class
ClosurePhaseDelaySpectrum, fix bug in calculation of diagonal weights
along the LST axis while operating on collapse_axes based on the
number of NaN along that axis
(*) In member function compute_power_spectrum_uncertainty() of class
ClosurePhaseDelaySpectrum, fix bug in calculation of diagonal weights
along the LST axis while operating on collapse_axes based on the
number of NaN along that axis
From d2e76c3 (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In member function FT() of class ClosurePhaseDelaySpectrum,
allow visscaleinfo to contain visibilities in the form of a numpy
array or an instance of class InterferometerArray. Update computing
of matching baselines and LST accordingly. Update member function
docstring accordingly.
(*) In member function FT() of class ClosurePhaseDelaySpectrum,
pass key 'bltriplet' in input visscaleinfo. Use this information to
determine the matching baselines in the instance of class
ClosurePhaseDelaySpectrum if the latter has to be used for scaling
visibilities. Update member function docstring accordingly
(*) If LST range in the input instance of class InterferometerArray
does not cover the data LST range, fill those visibilities with NaN
and propagate it to the computed delay spectrum. Update member
function docstring accordingly
From ccae7ac (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) Remove deprecated member function compute_power_spectrum_old()
under class ClosurePhaseDelaySpectrum. Update class docstring
accordingly
From fe688ef (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In member function compute_power_spectrum() of class
ClosurePhaseDelaySpectrum, compute and return diagonal weights in the
output. Update member function docstring accordingly
(*) In member function compute_power_spectrum() of class
ClosurePhaseDelaySpectrum, instead of returning traces of diagonals
which represent the sum, return the means along the traces along with
the diagonal weights.
(*) In member function compute_power_spectrum_uncertainty() of class
ClosurePhaseDelaySpectrcompute_power_spectrum_uncertainty diagonal
weights in the output. Update member function docstring accordingly
(*) In member function compute_power_spectrum_uncertainty() of class
ClosurePhaseDelaySpectrum, instead of returning traces of diagonals
which represent the sum, return the means along the traces along with
the diagonal weights.
From 5d0cbaa (BSP_PS_errorbar):
==============================
In prisim/bispectrum_phase.py,
(*) In member function smooth_in_tbins() of class ClosurePhase,
adjust mean and median of eicp quantities to unit amplitudes.
(*) In member function subsample_differencing() of class ClosurePhase,
adjust mean and median of eicp quantities to unit amplitudes.
(*) In member function subsample_differencing() of class ClosurePhase,
weight the subsample differences by 0.5 to get noise properties with
the right scaling
(*) In member function subsample_differencing() of class ClosurePhase,
change returned list_of_pair_of_pairs attribute to be a numpy array
instead of a list of pairs of sets
(*) In member function FT() of class ClosurePhaseDelaySpectrum,