Test tclean_VLASS_1v2

#	Test Name	Description	Run Time
1	test_j1302_awproject	J1302 (southern source), awproject gridder, Stokes I imaging	20870.99s
2	test_j1302_mosaic_cube	J1302 (southern source), coarse cube (continum image per spw), mosaic gridder, Stokes IQUV imaging	9084.27s
3	test_j1302_mtmfs	J1302 (southern source), mosaic gridder, Stokes I imaging	908.88s
4	test_j1302_ql	J1302 (southern source), QL, mosaic gridder, Stokes I imaging	835.86s
5	test_j1927_mosaic_cube	J1927 (northern source), coarse cube (continuum image per spw), mosaic gridder, Stokes IQUV imaging	15295.42s
6	test_j1927_mtmfs	J1927 (northern source), mosaic gridder, Stokes I imaging	2291.55s
7	test_j1927_ql	J1927 (northern source) QL, mosaic gridder, Stokes I imaging	2523.87s
8	Total	Total test runtime	51810.84s

test_j1302_awproject

_{J1302 (southern source), awproject gridder, Stokes I imaging}

Elapsed Time: 20870.986520051956 Seconds

Status: PASS

Task Executions:

immath(imagename=['J1302_awproject_secondmask.mask', 'J1302_awproject_QLcatmask.mask'], expr='IM0+IM1', outfile='J1302_awproject_sum_of_masks.mask')
tclean('J1302-12fields.ms', cfcache='cache0d.cf', niter=0, calcres=False, calcpsf=True, wbawp=False, datacolumn='corrected', mask='', usemask='user', pbmask=0.0, nsigma=2.0, cycleniter=5000, scales=[0], savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 5250, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'awproject', 'wprojplanes': 32, 'conjbeams': True, 'mosweight': False, 'usepointing': True, 'rotatepastep': 5.0, 'pointingoffsetsigdev': [300, 30], 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', cfcache='cache2.cf', niter=0, calcres=True, calcpsf=True, datacolumn='corrected', mask='', usemask='user', pbmask=0.0, wbawp=True, nsigma=2.0, cycleniter=5000, scales=[0], savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 5250, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'awproject', 'wprojplanes': 32, 'conjbeams': True, 'mosweight': False, 'usepointing': True, 'rotatepastep': 5.0, 'pointingoffsetsigdev': [300, 30], 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', cfcache='cache2.cf', niter=20000, scales=[0, 5, 12], nsigma=3.0, cycleniter=3000, mask='J1302_awproject_QLcatmask.mask', datacolumn='corrected', usemask='user', pbmask=0.0, wbawp=True, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 5250, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'awproject', 'wprojplanes': 32, 'conjbeams': True, 'mosweight': False, 'usepointing': True, 'rotatepastep': 5.0, 'pointingoffsetsigdev': [300, 30], 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', cfcache='cache2.cf', niter=0, datacolumn='data', savemodel='modelcolumn', mask='', usemask='user', pbmask=0.0, wbawp=True, nsigma=2.0, cycleniter=5000, scales=[0], calcres=False, calcpsf=False, common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 5250, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'awproject', 'wprojplanes': 32, 'conjbeams': True, 'mosweight': False, 'usepointing': True, 'rotatepastep': 5.0, 'pointingoffsetsigdev': [300, 30], 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', cfcache='cache2.cf', niter=20000, scales=[0, 5, 12], nsigma=3.0, cycleniter=3000, mask='J1302_awproject_combined.mask', datacolumn='corrected', usemask='user', pbmask=0.0, wbawp=True, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 5250, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'awproject', 'wprojplanes': 32, 'conjbeams': True, 'mosweight': False, 'usepointing': True, 'rotatepastep': 5.0, 'pointingoffsetsigdev': [300, 30], 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', cfcache='cache2.cf', niter=20000, scales=[0, 5, 12], nsigma=4.5, cycleniter=500, mask='', usemask='pb', pbmask=0.4, datacolumn='corrected', wbawp=True, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 5250, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'awproject', 'wprojplanes': 32, 'conjbeams': True, 'mosweight': False, 'usepointing': True, 'rotatepastep': 5.0, 'pointingoffsetsigdev': [300, 30], 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tt0stats=imstat(imagename='J1302_awproject_iter2.image.tt0', box='2625,2625,2625,2625')
tt1stats=imstat(imagename='J1302_awproject_iter2.image.tt1', box='2625,2625,2625,2625')
alphastats=imstat(imagename='J1302_awproject_iter2.alpha', box='2625,2625,2625,2625')
restbeam = imhead('J1302_awproject_iter2.image.tt0')['restoringbeam']

images:

j1302_awproject_iter2_tt0.png

j1302_awproject_iter2_tt1.png

Serial Mode: MPI Environment Not Enabled

Number of processors available in the system: 24

report: [ test_j1302_awproject ] Frac Diff tt0 vs. 6.2.1 image is 0.005605126507764093 ( Pass : should be < 0.1) (raw actual/expected values: 0.3156070113182068/0.317386) [ test_j1302_awproject ] Frac Diff tt1 vs. 6.2.1 image is 0.006482147306811382 ( Pass : should be < 0.2751310082438085) (raw actual/expected values: -0.01527363806962967/-0.01537329) [ test_j1302_awproject ] Abs Diff alpha vs. 6.2.1 image is 4.272741719246048e-05 ( Pass : should be < 0.2) (raw actual/expected values: -0.04839448258280754/-0.04843721) [ test_j1302_awproject ] Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(7.572259882794729e-06), np.float64(2.0896352852411895e-05), np.float64(6.714091586307467e-05)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(3.0541276931762695), np.float64(2.4986119270324707), np.float64(10.82421588897705)]/[np.float64(3.05415082), np.float64(2.49866414), np.float64(10.82348919)]) [ test_j1302_awproject ] count('MODEL_DATA') is 1 ( Pass : should be 1, Exact: True )

Mem Use: 2505.3046875 MiB

Re-Run: /export/cvpost134-nvme/casa/bamboohome/xml-data/build-dir/CASASPECS-MASTERTESTSTKPY312-V1STE/work/pkg/casa-6.7.0-11-py3.12.el8/bin/python3 /export/cvpost134-nvme/casa/bamboohome/xml-data/build-dir/CASASPECS-MASTERTESTSTKPY312-V1STE/work/casasources//casa6/casatests/stakeholder/test_stk_vlass_pipeline_imaging_1v2.py test_j1302.test_j1302_awproject

test_j1302_mosaic_cube

_{J1302 (southern source), coarse cube (continum image per spw), mosaic gridder, Stokes IQUV imaging}

Elapsed Time: 9084.272276878357 Seconds

Status: PASS

Task Executions:

r[0] = tclean('J1302-12fields.ms', datacolumn='corrected', niter=0, spw='2', stokes='IQUV', reffreq='3.564GHz', mask='', usemask='user', pbmask=0.0, nsigma=2.0, cycleniter=500, scales=[0], calcres=True, calcpsf=True, psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'usepointing': False, 'uvtaper': [''], 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
r[1] = tclean('J1302-12fields.ms', datacolumn='corrected', scales=[0, 5, 12], nsigma=3.0, niter=20000, cycleniter=500, mask='J1302_mosaic_cube_QLcatmask.mask', calcres=False, calcpsf=False, spw='2', stokes='IQUV', reffreq='3.564GHz', usemask='user', pbmask=0.0, psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'usepointing': False, 'uvtaper': [''], 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
r[2] = tclean('J1302-12fields.ms', datacolumn='corrected', scales=[0, 5, 12], nsigma=3.0, niter=20000, cycleniter=500, mask='J1302_mosaic_cube_combined.mask', calcres=False, calcpsf=False, spw='2', stokes='IQUV', reffreq='3.564GHz', usemask='user', pbmask=0.0, psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'usepointing': False, 'uvtaper': [''], 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
r[3] = tclean('J1302-12fields.ms', datacolumn='corrected', scales=[0, 5, 12], nsigma=4.5, niter=20000, cycleniter=100, mask='', calcres=False, calcpsf=False, usemask='pb', pbmask=0.4, spw='2', stokes='IQUV', reffreq='3.564GHz', psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'usepointing': False, 'uvtaper': [''], 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tt0statsI=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='I')
tt0statsQ=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='Q')
tt0statsU=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='U')
tt0statsV=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='V')
rmstt0statsI=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.residual.tt0', stokes='I')
rmstt0statsQ=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.residual.tt0', stokes='Q')
rmstt0statsU=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.residual.tt0', stokes='U')
rmstt0statsV=imstat(imagename='J1302_mosaic_cube_iter2_2_IQUV.residual.tt0', stokes='V')
header = imhead('J1302_mosaic_cube_iter2_2_IQUV.psf.tt0')

images:

j1302_mosaic_cube_spw2.png

j1302_mosaic_cube_spw8.png

j1302_mosaic_cube_spw14.png

Serial Mode: MPI Environment Not Enabled

Number of processors available in the system: 24

report: [ test_j1302_mosaic_cube ] Frac Diff F_nu vs. 6.2.1 image is 0.0003249264090835648 ( Pass : should be < 0.1) (raw actual/expected values: 0.31241716961565213/0.31231569) [ test_j1302_mosaic_cube ] Abs Diff alpha vs. on-axis is 0.0852119669706424 ( Pass : should be < 0.24880476650159936) (raw actual/expected values: 0.037611966970642396/-0.0476) [ test_j1302_mosaic_cube ] Abs Diff alpha vs. 6.2.1 image is 0.0009761269706423936 ( Pass : should be < 0.24880476650159936) (raw actual/expected values: 0.037611966970642396/0.03663584) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff IQUV vs 6.2.1 is [np.float64(8.194903638711692e-06), np.float64(1.7627668441545275e-06), np.float64(0.000591511809575387), np.float64(0.008206980317289033)] ( Pass : should be < [0.1 0.75293858 1. 0.73512655]) (raw actual/expected values: [np.float64(0.3023480176925659), np.float64(-0.0016958370106294751), np.float64(-0.0004073947202414274), np.float64(-0.0017226750496774912)]/[np.float64(0.30234554), np.float64(-0.00169584), np.float64(-0.00040815), np.float64(-0.00173693)]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(0.005293841013337406), np.float64(0.0022214445370190164), np.float64(0.03097696138265046)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(4.653586387634277), np.float64(3.7662594318389893), np.float64(8.387100219726562)]/[np.float64(4.62908077), np.float64(3.75791144), np.float64(8.6552124)]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff IQUV vs 6.2.1 is [np.float64(0.00010676764633294085), np.float64(1.9945999188831424e-11), np.float64(2.0500930016558436e-07), np.float64(7.202567402098499e-06)] ( Pass : should be < [0.1 1. 1. 1. ]) (raw actual/expected values: [np.float64(0.32466351985931396), np.float64(-0.00010136422497453168), np.float64(-0.0002730565902311355), np.float64(-0.000981407705694437)]/[np.float64(0.32462886), np.float64(-0.000101364225), np.float64(-0.000273056852), np.float64(-0.000981398509)]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(0.0051022469247153455), np.float64(8.956149494319215e-05), np.float64(0.08865688123329074)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(3.334955930709839), np.float64(2.8108577728271484), np.float64(10.876094818115234)]/[np.float64(3.31802654), np.float64(2.81110954), np.float64(11.9341383)]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff IQUV vs 6.2.1 is [np.float64(0.00058517705936094), np.float64(6.953201046220804e-07), np.float64(4.103472707746564e-07), np.float64(1.5983402780352728e-06)] ( Pass : should be < [0.1 1. 0.54531621 1. ]) (raw actual/expected values: [np.float64(0.3073750436306), np.float64(0.0011263691121712327), np.float64(-0.002341510960832238), np.float64(-0.000724207959137857)]/[np.float64(0.30719528), np.float64(0.00112637), np.float64(-0.00234151), np.float64(-0.00072421)]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(0.014410739880193054), np.float64(0.0030767082376187992), np.float64(0.06023197020698886)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(2.026454210281372), np.float64(1.619074821472168), np.float64(7.9184465408325195)]/[np.float64(1.99766636), np.float64(1.61410868), np.float64(8.42595863)])

Mem Use: -424.421875 MiB

test_j1302_mtmfs

_{J1302 (southern source), mosaic gridder, Stokes I imaging}

Elapsed Time: 908.880274772644 Seconds

Status: PASS

Task Executions:

immath(imagename=['J1302_mtmfs_secondmask.mask', 'J1302_mtmfs_QLcatmask.mask'], expr='IM0+IM1', outfile='J1302_mtmfs_sum_of_masks.mask')
#records.append('tstobj.run_tclean(**params') )
tclean('J1302-12fields.ms', calcres=True, calcpsf=True, datacolumn='corrected', mask='', usemask='user', pbmask=0.0, nsigma=2.0, cycleniter=500, scales=[0], savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imagename': 'J1302_mtmfs_iter2', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', mask='J1302_mtmfs_QLcatmask.mask', nsigma=3.0, scales=[0, 5, 12], datacolumn='corrected', usemask='user', pbmask=0.0, cycleniter=500, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imagename': 'J1302_mtmfs_iter2', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', datacolumn='data', savemodel='modelcolumn', mask='', usemask='user', pbmask=0.0, nsigma=2.0, cycleniter=500, scales=[0], calcres=False, calcpsf=False, common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imagename': 'J1302_mtmfs_iter2', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', mask='J1302_mtmfs_combined.mask', nsigma=3.0, scales=[0, 5, 12], datacolumn='corrected', usemask='user', pbmask=0.0, cycleniter=500, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imagename': 'J1302_mtmfs_iter2', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', usemask='pb', mask='', pbmask=0.4, nsigma=4.5, scales=[0, 5, 12], cycleniter=100, datacolumn='corrected', calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imagename': 'J1302_mtmfs_iter2', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tt0stats = imstat(imagename='J1302_mtmfs_iter2.image.tt0', box='2000,2000,2000,2000')
tt1stats = imstat(imagename='J1302_mtmfs_iter2.image.tt1', box='2000,2000,2000,2000')
alphastats = imstat(imagename='J1302_mtmfs_iter2.alpha', box='2000,2000,2000,2000')
restbeam = imhead('J1302_mtmfs_iter2.image.tt0')['restoringbeam']

images:

j1302_mtmfs_tt0.png

j1302_mtmfs_tt1.png

Serial Mode: MPI Environment Not Enabled

Number of processors available in the system: 24

report: [ test_j1302_mtmfs ] Frac Diff tt0 vs. 6.2.1 image is 1.275879967226669e-05 ( Pass : should be < 0.1) (raw actual/expected values: 0.3196171820163727/0.31962126) [ test_j1302_mtmfs ] Frac Diff tt1 vs. 6.2.1 image is 0.003656827251761227 ( Pass : should be < 0.13582807312418335) (raw actual/expected values: 0.019359225407242775/0.01928869) [ test_j1302_mtmfs ] Abs Diff alpha vs. on-axis is 0.10817003885507584 ( Pass : should be < 0.2) (raw actual/expected values: 0.060570038855075836/-0.0476) [ test_j1302_mtmfs ] Abs Diff alpha vs. 6.2.1 image is 0.00022146885507583897 ( Pass : should be < 0.2) (raw actual/expected values: 0.060570038855075836/0.06034857) [ test_j1302_mtmfs ] Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(1.2747511586951134e-09), np.float64(1.2464764565995781e-09), np.float64(1.274943184034465e-10)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(3.150028944015503), np.float64(2.591778516769409), np.float64(11.412883758544922)]/[np.float64(3.15002894), np.float64(2.59177852), np.float64(11.41288376)]) [ test_j1302_mtmfs ] count('MODEL_DATA') is 1 ( Pass : should be 1, Exact: True )

Mem Use: 2.56640625 MiB

test_j1302_ql

_{J1302 (southern source), QL, mosaic gridder, Stokes I imaging}

Elapsed Time: 835.8621916770935 Seconds

Status: PASS

Task Executions:

tclean('J1302-12fields.ms', niter=0, restoration=False, mask='', usemask='user', pbmask=0.0, nsigma=0.0, cycleniter=-1, cyclefactor=1.0, scales=[0], calcres=True, calcpsf=True, savemodel='none', restoringbeam=None, common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 7290, 'cell': '1.0arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'datacolumn': 'data', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1302-12fields.ms', niter=20000, restoration=True, restoringbeam='common', nsigma=4.5, cycleniter=500, cyclefactor=2.0, calcres=False, calcpsf=False, mask='', usemask='user', pbmask=0.0, scales=[0], savemodel='none', common_args={'vis': 'J1302-12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 7290, 'cell': '1.0arcsec', 'phasecenter': '13:03:13.874 -10.51.16.73', 'datacolumn': 'data', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'pblimit': 0.02, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'restart': True, 'interactive': False, 'parallel': False})
impbcor(imagename='J1302_ql_iter1.residual.tt0', pbimage='J1302_ql_iter1.pb.tt0', outfile='J1302_ql_iter1.image.residual.pbcor.tt0', mode='divide', cutoff=-1.0, stretch=False)
imdev(imagename='J1302_ql_iter1.image.pbcor.tt0', outfile='J1302_ql_iter1.image.pbcor.tt0.rms', overwrite=True, stretch=False, grid=[10, 10], anchor='ref', xlength='60arcsec', ylength='60arcsec', interp='cubic', stattype='xmadm', statalg='chauvenet', zscore=-1, maxiter=-1)
imhead(imagename='J1302_ql_iter1.pb.tt0')
imsubimage(imagename='J1302_ql_iter1.pb.tt0', outfile='J1302_ql_iter1.pb.tt0.subim', box='1785,1785,5506,5506')
imstat_vals = imstat(imagename='J1302_ql_iter1.image.pbcor.tt0.subim', box='1860,1860,1860,1860')
restbeam = imhead('J1302_ql_iter1.image.pbcor.tt0.subim')['restoringbeam']

images:

j1302_ql.png

Serial Mode: MPI Environment Not Enabled

Number of processors available in the system: 24

report: [ test_j1302_ql ] image_exists('J1302_ql_iter1.image.pbcor.tt0') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.image.residual.pbcor.tt0') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.image.pbcor.tt0.rms') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.image.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.residual.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.image.pbcor.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.image.pbcor.tt0.rms.subim') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.psf.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.image.residual.pbcor.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] image_exists('J1302_ql_iter1.pb.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1302_ql ] Frac Diff F_nu vs. on-axis is [np.float64(0.037053740806911095)] ( Pass : should be < [0.1]) (raw actual/expected values: 0.32133516669273376/[np.float64(0.3337)]) [ test_j1302_ql ] Frac Diff F_nu vs. 6.2.1 image is [np.float64(0.0024528077173352476)] ( Pass : should be < [0.1]) (raw actual/expected values: 0.32133516669273376/[np.float64(0.3205489218235016)]) [ test_j1302_ql ] Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(4.567121531468285e-07), np.float64(2.000210684915838e-06), np.float64(9.680161446864748e-06)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(3.123652696609497), np.float64(2.620107650756836), np.float64(13.793045997619629)]/[np.float64(3.12365127), np.float64(2.62010241), np.float64(13.79291248)])

Mem Use: 25.73046875 MiB

test_j1927_mosaic_cube

_{J1927 (northern source), coarse cube (continuum image per spw), mosaic gridder, Stokes IQUV imaging}

Elapsed Time: 15295.41998577118 Seconds

Status: PASS

Task Executions:

r[0] = tclean('J1927-12fields.ms', datacolumn='corrected', niter=0, spw='2', stokes='IQUV', reffreq='3.564GHz', mask='', usemask='user', pbmask=0.0, cfcache='', wbawp=True, nsigma=2.0, cycleniter=5000, scales=[0], calcres=True, calcpsf=True, psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
r[1] = tclean('J1927-12fields.ms', datacolumn='corrected', scales=[0, 5, 12], nsigma=3.0, niter=20000, cycleniter=500, mask='J1927_mosaic_cube_QLcatmask.mask', calcres=False, calcpsf=False, spw='2', stokes='IQUV', reffreq='3.564GHz', usemask='user', pbmask=0.0, cfcache='', wbawp=True, psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
r[2] = tclean('J1927-12fields.ms', datacolumn='corrected', scales=[0, 5, 12], nsigma=3.0, niter=20000, cycleniter=500, mask='J1927_mosaic_cube_combined.mask', calcres=False, calcpsf=False, spw='2', stokes='IQUV', reffreq='3.564GHz', usemask='user', pbmask=0.0, cfcache='', wbawp=True, psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
r[3] = tclean('J1927-12fields.ms', datacolumn='corrected', scales=[0, 5, 12], nsigma=4.5, niter=20000, cycleniter=100, mask='', calcres=False, calcpsf=False, usemask='pb', pbmask=0.4, spw='2', stokes='IQUV', reffreq='3.564GHz', cfcache='', wbawp=True, psfcutoff=0.5, savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 1, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tt0statsI=imstat(imagename='J1927_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='I')
tt0statsQ=imstat(imagename='J1927_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='Q')
tt0statsU=imstat(imagename='J1927_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='U')
tt0statsV=imstat(imagename='J1927_mosaic_cube_iter2_2_IQUV.image.tt0', box='2000,2000,2000,2000', stokes='V')
header=imhead('J1927_mosaic_cube_iter2_2_IQUV.psf.tt0')['perplanebeams']['beams']['*0']['*0']

images:

j1927_mosaic_cube_spw2.png

j1927_mosaic_cube_spw8.png

j1927_mosaic_cube_spw14.png

Serial Mode: MPI Environment Not Enabled

Number of processors available in the system: 24

report: [ test_j1927_mosaic_cube ] Frac Diff F_nu vs. on-axis is 0.06586175976714789 ( Pass : should be < 0.1) (raw actual/expected values: 0.888272052637419/0.9509) [ test_j1927_mosaic_cube ] Frac Diff F_nu vs. 6.2.1 image is 9.513900773083878e-05 ( Pass : should be < 0.1) (raw actual/expected values: 0.888272052637419/0.88835657) [ test_j1927_mosaic_cube ] Abs Diff alpha vs. on-axis is 0.051217442554740256 ( Pass : should be < 0.2) (raw actual/expected values: 0.41131744255474023/0.3601) [ test_j1927_mosaic_cube ] Abs Diff alpha vs. 6.2.1 image is 0.0007050874452597666 ( Pass : should be < 0.2) (raw actual/expected values: 0.41131744255474023/0.41202253) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff IQUV vs 6.2.1 is [np.float64(0.00021162304355661548), np.float64(0.000270046717227037), np.float64(0.0009542526851552822), np.float64(0.009819722844876116)] ( Pass : should be < [0.1 0.16326569 1. 0.14952827]) (raw actual/expected values: [np.float64(0.7555659413337708), np.float64(0.0057603451423347), np.float64(0.000794297200627625), np.float64(-0.006349605042487383)]/[np.float64(0.75540608), np.float64(0.00575879), np.float64(0.0007934), np.float64(-0.00628786)]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(0.0030857301278356913), np.float64(0.0007222850070653711), np.float64(0.04683829777110915)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(3.665053606033325), np.float64(3.029024839401245), np.float64(-23.843441009521484)]/[np.float64(3.65377903), np.float64(3.03121424), np.float64(-22.77662277)]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff IQUV vs 6.2.1 is [np.float64(0.0001153497613042239), np.float64(4.613063307735884e-05), np.float64(8.178594781789252e-06), np.float64(0.04347684735705346)] ( Pass : should be < [0.1 0.17989271 0.15549112 0.12771143]) (raw actual/expected values: [np.float64(0.8643968105316162), np.float64(0.005226761102676392), np.float64(-0.006046779453754425), np.float64(-0.007041933014988899)]/[np.float64(0.86449653), np.float64(0.00522652), np.float64(-0.00604673), np.float64(-0.00736201)]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(0.000252757728101065), np.float64(0.0033029517209369062), np.float64(0.04842347672051004)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(2.655961275100708), np.float64(2.23101544380188), np.float64(-28.002300262451172)]/[np.float64(2.65529013), np.float64(2.2384088), np.float64(-26.70895958)]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff IQUV vs 6.2.1 is [np.float64(0.00018329475508388157), np.float64(0.0011830922710417853), np.float64(0.007697839054717622), np.float64(0.02464209567842369)] ( Pass : should be < [0.1 0.17300434 0.1 0.26293629]) (raw actual/expected values: [np.float64(0.9558669924736023), np.float64(0.005428190343081951), np.float64(-0.009582056663930416), np.float64(-0.003663935698568821)]/[np.float64(0.95604223), np.float64(0.00543462), np.float64(-0.00965639), np.float64(-0.00357582)]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(0.0013304395597801193), np.float64(0.005160215228592125), np.float64(0.01120997146896644)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(2.0644712448120117), np.float64(1.6848171949386597), np.float64(-25.778289794921875)]/[np.float64(2.06172824), np.float64(1.69355631), np.float64(-25.49251938)])

Mem Use: 0.6953125 MiB

test_j1927_mtmfs

_{J1927 (northern source), mosaic gridder, Stokes I imaging}

Elapsed Time: 2291.5477945804596 Seconds

Status: PASS

Task Executions:

immath(imagename=['J1927_mtmfs_secondmask.mask', 'J1927_mtmfs_QLcatmask.mask'], expr='IM0+IM1', outfile='J1927_mtmfs_sum_of_masks.mask')
tclean('J1927-12fields.ms', niter=0, datacolumn='corrected', calcres=True, calcpsf=True, mask='', usemask='user', pbmask=0.0, nsigma=2.0, cycleniter=500, scales=[0], savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1927-12fields.ms', niter=20000, datacolumn='corrected', mask='J1927_mtmfs_QLcatmask.mask', nsigma=3.0, scales=[0, 5, 12], usemask='user', pbmask=0.0, cycleniter=500, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1927-12fields.ms', niter=0, datacolumn='data', savemodel='modelcolumn', mask='', usemask='user', pbmask=0.0, nsigma=2.0, cycleniter=500, scales=[0], calcres=False, calcpsf=False, common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1927-12fields.ms', niter=20000, datacolumn='corrected', mask='J1927_mtmfs_combined.mask', nsigma=3.0, scales=[0, 5, 12], usemask='user', pbmask=0.0, cycleniter=500, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1927-12fields.ms', niter=20000, datacolumn='corrected', usemask='pb', mask='', pbmask=0.4, nsigma=4.5, scales=[0, 5, 12], cycleniter=100, calcres=False, calcpsf=False, savemodel='none', common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 4000, 'cell': '0.6arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'spw': '', 'reffreq': '3.0GHz', 'gridder': 'mosaic', 'mosweight': False, 'pblimit': 0.1, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.4, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'cyclefactor': 3.0, 'restart': True, 'interactive': False, 'parallel': False})
tt0stats = imstat(imagename='J1927_mtmfs_iter2.image.tt0', box='2000,2000,2000,2000')
tt1stats = imstat(imagename='J1927_mtmfs_iter2.image.tt1', box='2000,2000,2000,2000')
alphastats = imstat(imagename='J1927_mtmfs_iter2.alpha', box='2000,2000,2000,2000')
restbeam = imhead('J1927_mtmfs_iter2.image.tt0')['restoringbeam']

images:

j1927_mtmfs_tt0.png

j1927_mtmfs_tt1.png

Serial Mode: MPI Environment Not Enabled

Number of processors available in the system: 24

report: [ test_j1927_mtmfs ] Frac Diff tt0 vs. on-axis is 0.06627108138346471 ( Pass : should be < 0.1) (raw actual/expected values: 0.8878828287124634/0.9509) [ test_j1927_mtmfs ] Frac Diff tt0 vs. 6.2.1 image is 9.867519453105647e-06 ( Pass : should be < 0.1) (raw actual/expected values: 0.8878828287124634/0.88789159) [ test_j1927_mtmfs ] Frac Diff tt1 vs. 6.2.1 image is 0.00025937317315127465 ( Pass : should be < 0.1) (raw actual/expected values: 0.4129445254802704/0.41305166) [ test_j1927_mtmfs ] Abs Diff alpha vs. on-axis is 0.08548899331092835 ( Pass : should be < 0.2) (raw actual/expected values: 0.46508899331092834/0.3796) [ test_j1927_mtmfs ] Abs Diff alpha vs. 6.2.1 image is 0.00011607668907165358 ( Pass : should be < 0.2) (raw actual/expected values: 0.46508899331092834/0.46520507) [ test_j1927_mtmfs ] Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(1.8401506739554768e-09), np.float64(2.285955441074116e-09), np.float64(7.435039114081441e-12)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(2.4488351345062256), np.float64(2.0478949546813965), np.float64(-23.64226531982422)]/[np.float64(2.44883513), np.float64(2.04789495), np.float64(-23.64226532)]) [ test_j1927_mtmfs ] count('MODEL_DATA') is 1 ( Pass : should be 1, Exact: True )

Mem Use: 0.45703125 MiB

test_j1927_ql

_{J1927 (northern source) QL, mosaic gridder, Stokes I imaging}

Elapsed Time: 2523.8681139945984 Seconds

Status: PASS

Task Executions:

tclean('J1927-12fields.ms', niter=0, restoration=False, mask='', usemask='user', pbmask=0.0, nsigma=0.0, cycleniter=-1, cyclefactor=1.0, scales=[0], calcres=True, calcpsf=True, common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 7290, 'cell': '1.0arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'reffreq': '3.0GHz', 'datacolumn': 'data', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'restoringbeam': 'common', 'pblimit': 0.2, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.0, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'restart': True, 'interactive': False, 'parallel': False})
tclean('J1927-12fields.ms', niter=20000, restoration=True, nsigma=4.5, cycleniter=500, cyclefactor=2.0, calcres=False, calcpsf=False, mask='', usemask='user', pbmask=0.0, scales=[0], common_args={'vis': 'J1927_12fields.ms', 'uvrange': '<12km', 'intent': 'OBSERVE_TARGET#UNSPECIFIED', 'imsize': 7290, 'cell': '1.0arcsec', 'phasecenter': '19:27:30.443 +61.17.32.898', 'reffreq': '3.0GHz', 'datacolumn': 'data', 'gridder': 'mosaic', 'conjbeams': False, 'mosweight': False, 'restoringbeam': 'common', 'pblimit': 0.2, 'deconvolver': 'mtmfs', 'nterms': 2, 'smallscalebias': 0.0, 'weighting': 'briggs', 'robust': 1.0, 'threshold': 0.0, 'restart': True, 'interactive': False, 'parallel': False})
impbcor(imagename='J1927_ql_iter1.residual.tt0', pbimage='J1927_ql_iter1.pb.tt0', outfile='J1927_ql_iter1.image.residual.pbcor.tt0', mode='divide', cutoff=-1.0, stretch=False)
imdev(imagename='J1927_ql_iter1.image.pbcor.tt0', outfile='J1927_ql_iter1.image.pbcor.tt0.rms', overwrite=True, stretch=False, grid=[10, 10], anchor='ref', xlength='60arcsec', ylength='60arcsec', interp='cubic', stattype='xmadm', statalg='chauvenet', zscore=-1, maxiter=-1)
imhead(imagename='J1927_ql_iter1.pb.tt0')
imsubimage(imagename='J1927_ql_iter1.pb.tt0', outfile='J1927_ql_iter1.pb.tt0.subim', box='1785,1785,5506,5506')
imstat_vals = imstat(imagename='J1927_ql_iter1.image.pbcor.tt0.subim', box='1860,1860,1860,1860')
restbeam = imhead('J1927_ql_iter1.image.pbcor.tt0.subim')['restoringbeam']

images:

j1927_ql.png

Serial Mode: MPI Environment Not Enabled

Number of processors available in the system: 24

report: [ test_j1927_ql ] image_exists('J1927_ql_iter1.image.pbcor.tt0') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.image.residual.pbcor.tt0') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.image.pbcor.tt0.rms') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.image.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.residual.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.image.pbcor.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.image.pbcor.tt0.rms.subim') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.psf.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.image.residual.pbcor.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] image_exists('J1927_ql_iter1.pb.tt0.subim') is True ( Pass : should be True, Exact: True ) [ test_j1927_ql ] Frac Diff F_nu vs. on-axis is [np.float64(0.0481328156884755)] ( Pass : should be < [0.1]) (raw actual/expected values: 0.9051305055618286/[np.float64(0.9509)]) [ test_j1927_ql ] Frac Diff F_nu vs. 6.2.1 image is [[np.float64(3.9575479559659906e-05)]] ( Pass : should be < [[0.1]]) (raw actual/expected values: 0.9051305055618286/[[np.float64(0.9051663279533386)]]) [ test_j1927_ql ] Frac Diff Maj, Min, PA vs 6.2.1 is [np.float64(4.1805116728073546e-06), np.float64(6.423501762312727e-06), np.float64(1.6162444857662526e-05)] ( Pass : should be < 0.1) (raw actual/expected values: [np.float64(2.4521117210388184), np.float64(2.0791704654693604), np.float64(-18.881975173950195)]/[np.float64(2.45210147), np.float64(2.07915711), np.float64(-18.88167)])

Mem Use: 5.98828125 MiB

Test tclean_VLASS_1v2_pipeline

test_j1302_awproject

test_j1302_mosaic_cube

test_j1302_mtmfs

test_j1302_ql

test_j1927_mosaic_cube

test_j1927_mtmfs

test_j1927_ql