ratt-ru / CubiCal

A fast radio interferometric calibration suite.
GNU General Public License v2.0
18 stars 13 forks source link

Amplitude weirdness in tf-plane solver (diag-diag 0) #228

Open bennahugo opened 6 years ago

bennahugo commented 6 years ago

Attempting a tf-plane calibration. The gainphase table that gets dumped has 0 amplitude although it does not look like it gets enforced. The convergence rate is very low as well but it is probably just an snr issue. The images do improve slightly though.

image

  4  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [data] Visibility data options
  5  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - ms ................................................ = msdir/COMBINED.1GC.ms
  6  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - column ............................................ = DATA
  7  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - time-chunk ........................................ = 100
  8  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - freq-chunk ........................................ = 0
  9  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - rebin-time ........................................ = 1
 10  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - rebin-freq ........................................ = 1
 11  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - chunk-by .......................................... = None
 12  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - chunk-by-jump ..................................... = 1.0
 13  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - single-chunk ...................................... =
 14  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [sel] Data selection options
 15  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - field ............................................. = 0
 16  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - ddid .............................................. = None
 17  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - taql .............................................. =
 18  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - chan .............................................. =
 19  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [model] Calibration model options
 20  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - list .............................................. = MODEL_DATA
 21  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - ddes .............................................. = auto
 22  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - beam-pattern ...................................... = meerkat_pb_jones_cube_95channels_$(CORR)_$(REIM).fits
 23  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - beam-l-axis ....................................... = None
 24  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - beam-m-axis ....................................... = None
 25  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [montblanc] Montblanc simulation options
 26  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - device-type ....................................... = CPU
 27  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - dtype ............................................. = float
 28  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - feed-type ......................................... = linear
 29  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - mem-budget ........................................ = 1024
 30  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - verbosity ......................................... = WARNING
 31  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - threads ........................................... = 0
 32  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [weight] Weighting options
 33  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - column ............................................ = WEIGHT
 34  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [flags] General flagging options
 35  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - apply ............................................. = -cubical
 36  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - auto-init ......................................... = legacy
 37  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - save .............................................. = cubical
 38  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - save-legacy ....................................... = auto
 39  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - reinit-bitflags ................................... = False
 40  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - warn-thr .......................................... = 0.3
 41  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - see-no-evil ....................................... = False
 42  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [postmortem] Options for "postmortem" flagging based on solution statistics
 43  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - enable ............................................ = False
 44  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - tf-chisq-median ................................... = 1.2
 45  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - tf-np-median ...................................... = 0.5
 46  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - time-density ...................................... = 0.5
 47  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - chan-density ...................................... = 0.5
 48  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - ddid-density ...................................... = 0.5
 49  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [madmax] Options for the "Mad Max" flagger
 50  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - enable ............................................ = 0
 51  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - estimate .......................................... = corr
 52  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - diag .............................................. = True
 53  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - offdiag ........................................... = True
 54  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - threshold ......................................... = [0, 10]
 55  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - global-threshold .................................. = [0, 12]
 56  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - plot .............................................. = 1
 57  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - plot-frac-above ................................... = 0.01
 58  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - plot-bl ........................................... =
 59  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - flag-ant .......................................... = 0
 60  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - flag-ant-thr ...................................... = 5
 61  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [sol] Solution options which apply at the solver level
 62  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - jones ............................................. = G
 63  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - precision ......................................... = 32
- 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - delta-g ........................................... = 1e-06
 65  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - delta-chi ......................................... = 1e-06
 66  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - chi-int ........................................... = 5
 67  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - last-rites ........................................ = True
 68  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - stall-quorum ...................................... = 0.99
 69  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - diag-diag ......................................... = False
 70  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - term-iters ........................................ =
 71  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - min-bl ............................................ = 50.0
 72  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - max-bl ............................................ = 0.0
 73  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - subset ............................................ =
 74  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [bbc] Options for baseline-based corrections (a.k.a. BBCs, a.k.a. interferometer gains).
 75  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - load-from ......................................... =
 76  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - compute-2x2 ....................................... = False
 77  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - apply-2x2 ......................................... = False
 78  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - save-to ........................................... = {data[ms]}/BBC-field_{sel[field]}-ddid_{sel[ddid]}.parmdb
 79  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - per-chan .......................................... = True
 80  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - plot .............................................. = True
 81  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [dist] Parallelization and distribution options
 82  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - ncpu .............................................. = 16
 83  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - nworker ........................................... = 8
 84  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - nthread ........................................... = 0
 85  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - max-chunks ........................................ = 0
 86  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - min-chunks ........................................ = 0
 87  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - pin ............................................... = 0
 88  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - pin-io ............................................ = False
 89  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - pin-main .......................................... = io
 90  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [out] Options for output products
 91  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - name .............................................. = 2gc_phasecal
 92  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - mode .............................................. = sc
 93  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - column ............................................ = CORRECTED_DATA
 94  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - model-column ...................................... =
 95  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - reinit-column ..................................... = False
 96  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - subtract-model .................................... = 0
 97  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - subtract-dirs ..................................... = :
 98  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - plots ............................................. = 1
 99  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - casa-gaintables ................................... = True
100  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [log] Options related to logging
101  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - memory ............................................ = True
102  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - boring ............................................ = False
103  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - append ............................................ = False
104  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - verbose ........................................... = 0
105  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - file-verbose ...................................... = None
106  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [debug] Debugging options for the discerning masochist
107  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - pdb ............................................... = False
108  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - panic-amplitude ................................... = 0.0
109  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - stop-before-solver ................................ = False
110  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [misc] Miscellaneous options
111  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - random-seed ....................................... = None
112  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - parset-version .................................... = 0.1
113  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [g] Options for G-Jones term
114  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - solvable .......................................... = True
115  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - type .............................................. = tf-plane
116  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - load-from ......................................... =
117  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - xfer-from ......................................... =
118  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - save-to ........................................... = {data[ms]}/{JONES}-field_{sel[field]}-ddid_{sel[ddid]}.parmdb
119  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - dd-term ........................................... = 0
120  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - fix-dirs .......................................... =
121  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - diag-diag ......................................... = 0
122  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - update-type ....................................... = full
123  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - time-int .......................................... = 20
124  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - freq-int .......................................... = 0
- 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - clip-high ......................................... = 10.0
129  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - clip-after ........................................ = 5
130  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - max-iter .......................................... = 500
131  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - conv-quorum ....................................... = 0.99
132  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - ref-ant ........................................... = None
133  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - prop-flags ........................................ = default
134  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] [de] de
135  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - dd-term ........................................... = 1
136  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - clip-low .......................................... = 0.0
137  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb]  - clip-high ......................................... = 0.0
138  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] Enabling G-Jones
139  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] mode: solve and correct
140  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb] reading MS msdir/COMBINED.1GC.ms
141  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   detected linear (xy) feeds
142  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   fields are *0: J1131-1955
143  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   MS contains 1 spectral windows
144  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   max freq chunk size is -- channels and/or -- MHz
145  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   DDID 0: 1366/1366 channels selected
146  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]     found 1 frequency chunks: 0 1366
147  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   applying TAQL query 'FIELD_ID == 0 && DATA_DESC_ID IN [0]' (86343/86343 rows selected)
148  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   16 antennas, 86343 rows, 1/1 DDIDs, 865 timeslots, 4 corrs
149  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   DDID central frequencies are at 1.28 GHz
150  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   inserting new column BITFLAG
151  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   inserting new column BITFLAG_ROW
152  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   Will auto-fill new BITFLAG 'legacy' (1) from FLAG/FLAG_ROW
153  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]     BITFLAG column defines the following flagsets: legacy:1
154  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]     Excluding flagset cubical
155  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]     flagset 'cubical' not found -- ignoring
156  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   Applying BITFLAG mask 1 to input data
157  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   Will save output flags into BITFLAG 'cubical' (2), and into FLAG/FLAG_ROW
158  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb] Using 1 model(s) for 1 directions(s) (DDEs explicitly disabled)
159  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   model 0 (weight WEIGHT):
160  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]     direction 0: MODEL_DATA
161  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] subtraction directions set to slice(None, None, None)
162  - 02:42:37 - param_db           [0.1/0.1 4.9/4.9 0.2Gb] creating msdir/COMBINED.1GC.ms/G-field_0-ddid_None.parmdb in fragmented mode
163  - 02:42:37 - gain_machine       [0.1/0.1 4.9/4.9 0.2Gb] G solutions will be saved to msdir/COMBINED.1GC.ms/G-field_0-ddid_None.parmdb
162  - 02:42:37 - param_db           [0.1/0.1 4.9/4.9 0.2Gb] creating msdir/COMBINED.1GC.ms/G-field_0-ddid_None.parmdb in fragmented mode
163  - 02:42:37 - gain_machine       [0.1/0.1 4.9/4.9 0.2Gb] G solutions will be saved to msdir/COMBINED.1GC.ms/G-field_0-ddid_None.parmdb
164  - 02:42:37 - param_db           [0.1/0.1 4.9/4.9 0.2Gb] creating msdir/COMBINED.1GC.ms/BBC-field_0-ddid_None.parmdb in fragmented mode
165  - 02:42:37 - gain_machine       [0.1/0.1 4.9/4.9 0.2Gb] will compute & save suggested baseline-based corrections (BBCs) to msdir/COMBINED.1GC.ms/BBC-field_0-ddid_None.parmdb
166  - 02:42:37 - gain_machine       [0.1/0.1 4.9/4.9 0.2Gb]   (these can optionally be applied in a subsequent CubiCal run)
167  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] multi-process mode: --dist-nworker 8 (+1), 1 OMP threads
168  - 02:42:37 - main               [0.1/0.1 4.9/4.9 0.2Gb] defining chunks (time 100, freq 0)
169  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   read indexing columns (86343 total rows)
170  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   built timeslot index (865 unique timestamps)
171  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   max chunk size is 100 timeslots and/or -- seconds
172  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   found 9 time chunks: 0 100 200 300 400 500 600 700 800 865
173  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   generated 9 row chunks based on time and DDID
174  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   row chunks yield 9 potential tiles
175  - 02:42:37 - data_handler       [0.1/0.1 4.9/4.9 0.2Gb]   coarsening this to 1 tiles (max 9 chunks per tile, based on 8/0 requested)
176  - 02:42:37 - main               [0.1/0.1 4.6/5.1 0.2Gb] waiting for I/O on tile #1/1
177  - 02:42:37 - main               [io] [0.1/0.1 4.6/4.6 0.2Gb] loading tile #1/1
178  - 02:42:37 - data_handler       [io] [0.1/0.1 4.6/4.6 0.2Gb] tile #1/1: reading MS rows 0~86342
179  - 02:42:52 - data_handler       [io] [5.9/7.2 11.3/11.7 0.2Gb]   applying solvable baseline cutoff deselects 2794 rows
180  - 02:42:52 - data_handler       [io] [5.9/7.2 11.3/11.7 0.2Gb]   3.24% visibilities deselected via specificed subset and/or baseline cutoffs
181  - 02:42:53 - data_handler       [io] [7.6/7.6 13.0/13.0 0.2Gb]   auto-filling BITFLAG/BITFLAG_ROW of shape (86343, 1366, 4) from FLAG/FLAG_ROW
182  - 02:43:06 - data_handler       [io] [9.3/11.5 13.9/16.1 0.2Gb]   51.81% input visibilities flagged and/or deselected
183  - 02:43:11 - data_handler       [io] [12.8/15.5 21.0/21.0 6.4Gb]   reading MODEL_DATA for model 0 direction 0
184  - 02:43:24 - main               [io] [1.9/19.9 6.5/28.0 9.9Gb] I/O job(s) complete
185  - 02:43:24 - main               [0.1/0.1 4.6/5.1 9.9Gb] submitting solver jobs for tile #1/1
186  - 03:27:52 - solver             [x08] [6.8/8.4 11.7/17.8 9.9Gb] D0T7F0 G: 500 iters, conv 0.33%, d/fl 6.25%, PGE 0.0275, stall 47.77%, chi^2 0.9484 -> 0.9117 (1.103), noise 2.55 -> 2.31
187  - 03:28:00 - solver             [x06] [6.8/8.4 11.7/17.9 10.3Gb] D0T5F0 G: 500 iters, conv 0.88%, d/fl 6.25%, PGE 0.0191, stall 54.00%, chi^2 0.938 -> 0.9165 (1.093), noise 2.4 -> 2.19
188  - 03:28:11 - solver             [x05] [6.8/8.4 11.7/17.8 11.0Gb] D0T4F0 G: 500 iters, conv 1.08%, d/fl 6.25%, PGE 0.0508, stall 57.05%, chi^2 0.962 -> 0.9408 (1.097), noise 2.43 -> 2.24
189  - 03:28:22 - solver             [x04] [6.8/8.4 11.7/17.9 11.6Gb] D0T3F0 G: 500 iters, conv 1.23%, d/fl 6.25%, PGE 0.0137, stall 59.96%, chi^2 0.9866 -> 0.9607 (1.097), noise 2.45 -> 2.28
190  - 03:28:25 - solver             [x07] [6.8/8.4 11.7/17.8 11.6Gb] D0T6F0 G: 500 iters, conv 0.56%, d/fl 6.25%, PGE 0.0466, stall 45.89%, chi^2 0.9232 -> 0.897 (1.092), noise 2.34 -> 2.1
191  - 03:28:30 - solver             [x01] [6.8/8.4 11.7/17.9 12.2Gb] D0T0F0 G: 500 iters, conv 0.88%, d/fl 6.25%, PGE 0.0104, stall 55.35%, chi^2 1.029 -> 0.987 (1.097), noise 2.38 -> 2.25
192  - 03:28:34 - solver             [x03] [6.8/8.4 11.7/17.9 12.2Gb] D0T2F0 G: 500 iters, conv 1.05%, d/fl 6.25%, PGE 0.0134, stall 56.55%, chi^2 1.009 -> 0.966 (1.093), noise 2.38 -> 2.23
193  - 03:28:50 - solver             [x02] [6.8/8.4 11.7/17.8 13.9Gb] D0T1F0 G: 500 iters, conv 0.95%, d/fl 6.25%, PGE 0.0153, stall 55.26%, chi^2 1.014 -> 0.9711 (1.095), noise 2.37 -> 2.22
194  - 03:49:38 - solver             [x08] [5.8/10.0 10.7/18.1 14.5Gb] D0T8F0 G: 500 iters, conv 0.14%, d/fl 6.25%, PGE 0.25, stall 40.89%, chi^2 0.9583 -> 0.9028 (1.103), noise 2.59 -> 2.33
195  - 03:49:43 - main               [0.2/0.2 4.8/5.1 14.9Gb] finished processing tile #1/1
196  - 03:49:43 - main               [io] [1.9/19.9 6.5/28.0 14.9Gb] saving tile #1/1
o-smirnov commented 6 years ago

As just discussed with @bennahugo, this is probably a CASA table export issue.

CubiCal in this instance produces G:phase, which is a real number (the phase, in radians). Sticking the number into a CASA gain table as is makes the latter interpret it as a real amplitude, no phase.

Solution is to export exp{1j*G_phase} into the CASA table.