def test_calibrate_arlexecute_empty(self):
amp_errors = {'T': 0.0, 'G': 0.0}
phase_errors = {'T': 1.0, 'G': 0.0}
self.actualSetUp(amp_errors=amp_errors, phase_errors=phase_errors)
for v in self.blockvis_list:
v.data['vis'][...] = 0.0 + 0.0j
controls = create_calibration_controls()
controls['T']['first_selfcal'] = 0
controls['T']['timescale'] = 'auto'
calibrate_list = \
calibrate_list_arlexecute_workflow(self.error_blockvis_list, self.blockvis_list,
calibration_context='T', controls=controls, do_selfcal=True,
global_solution=False)
calibrate_list = arlexecute.compute(calibrate_list, sync=True)
assert len(calibrate_list[1][0]) > 0
def test_calibrate_arlexecute_global(self):
amp_errors = {'T': 0.0, 'G': 0.0}
phase_errors = {'T': 1.0, 'G': 0.0}
self.actualSetUp(amp_errors=amp_errors, phase_errors=phase_errors)
controls = create_calibration_controls()
controls['T']['first_selfcal'] = 0
controls['T']['timescale'] = 'auto'
calibrate_list = \
calibrate_list_arlexecute_workflow(self.error_blockvis_list, self.blockvis_list,
calibration_context='T', controls=controls, do_selfcal=True,
global_solution=True)
calibrate_list = arlexecute.compute(calibrate_list, sync=True)
assert numpy.max(calibrate_list[1]['T'].residual) < 7e-6, numpy.max(
calibrate_list[1]['T'].residual)
assert numpy.max(
numpy.abs(calibrate_list[0][0].vis -
self.blockvis_list[0].vis)) < 2e-6
