Test tclean_VLASS_1v2

#	Test Name	Description	Run Time
1	test_j1302_awproject	J1302 (southern source), awproject gridder, Stokes I imaging	20473.59s
2	test_j1302_mosaic_cube	J1302 (southern source), coarse cube (continum image per spw), mosaic gridder, Stokes IQUV imaging	9294.68s
3	test_j1302_mtmfs	J1302 (southern source), mosaic gridder, Stokes I imaging	877.86s
4	test_j1302_ql	J1302 (southern source), QL, mosaic gridder, Stokes I imaging	804.15s
5	test_j1927_mosaic_cube	J1927 (northern source), coarse cube (continuum image per spw), mosaic gridder, Stokes IQUV imaging	15556.42s
6	test_j1927_mtmfs	J1927 (northern source), mosaic gridder, Stokes I imaging	2198.69s
7	Total	Total test runtime	49205.4s

test_j1302_awproject

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

Elapsed Time: 20473.587675094604 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.006482268467938845 ( Pass : should be < 0.2751310082438085) (raw actual/expected values: -0.01527363620698452/-0.01537329) [ test_j1302_awproject ] Abs Diff alpha vs. 6.2.1 image is 4.2734867773057406e-05 ( Pass : should be < 0.2) (raw actual/expected values: -0.048394475132226944/-0.04843721) [ test_j1302_awproject ] Frac Diff Maj, Min, PA vs 6.2.1 is [7.572259882794729e-06, 2.0800934434512177e-05, 6.72290274045254e-05] ( Pass : should be < 0.1) (raw actual/expected values: [3.0541276931762695, 2.49861216545105, 10.824216842651367]/[3.05415082, 2.49866414, 10.82348919]) [ test_j1302_awproject ] count('MODEL_DATA') is 1 ( Pass : should be 1, Exact: True )

Mem Use: 1439.3359375 MiB

Re-Run: /export/cvpost135-nvme/casa/bamboohome/xml-data/build-dir/CASA-CMTS-V1STEM/work/pkg/casa-6.6.4-19-py3.8.el7/bin/python3 /export/cvpost135-nvme/casa/bamboohome/xml-data/build-dir/CASA-CMTS-V1STEM/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: 9294.681391954422 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.000324889384938603 ( Pass : should be < 0.1) (raw actual/expected values: 0.31241715805243075/0.31231569) [ test_j1302_mosaic_cube ] Abs Diff alpha vs. on-axis is 0.0852119572524668 ( Pass : should be < 0.2488044293728314) (raw actual/expected values: 0.037611957252466806/-0.0476) [ test_j1302_mosaic_cube ] Abs Diff alpha vs. 6.2.1 image is 0.000976117252466803 ( Pass : should be < 0.2488044293728314) (raw actual/expected values: 0.037611957252466806/0.03663584) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff IQUV vs 6.2.1 is [8.194903638711692e-06, 1.7627668441545275e-06, 0.0005915573960199425, 0.008207114364512384] ( Pass : should be < [0.1 0.75293858 1. 0.73512655]) (raw actual/expected values: [0.3023480176925659, -0.0016958370106294751, -0.0004073946620337665, -0.0017226748168468475]/[0.30234554, -0.00169584, -0.00040815, -0.00173693]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff Maj, Min, PA vs 6.2.1 is [0.005293841013337406, 0.0022214445370190164, 0.03097696138265046] ( Pass : should be < 0.1) (raw actual/expected values: [4.653586387634277, 3.7662594318389893, 8.387100219726562]/[4.62908077, 3.75791144, 8.6552124]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff IQUV vs 6.2.1 is [0.00010667584204176449, 1.9945999188831424e-11, 2.0500930016558436e-07, 7.1113945129874556e-06] ( Pass : should be < [0.1 1. 1. 1. ]) (raw actual/expected values: [0.3246634900569916, -0.00010136422497453168, -0.0002730565902311355, -0.0009814075892791152]/[0.32462886, -0.000101364225, -0.000273056852, -0.000981398509]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff Maj, Min, PA vs 6.2.1 is [0.0051022469247153455, 8.956149494319215e-05, 0.08865688123329074] ( Pass : should be < 0.1) (raw actual/expected values: [3.334955930709839, 2.8108577728271484, 10.876094818115234]/[3.31802654, 2.81110954, 11.9341383]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff IQUV vs 6.2.1 is [0.00058517705936094, 6.953201046220804e-07, 4.103472707746564e-07, 1.5071673889242288e-06] ( Pass : should be < [0.1 1. 0.54531621 1. ]) (raw actual/expected values: [0.3073750436306, 0.0011263691121712327, -0.002341510960832238, -0.0007242080755531788]/[0.30719528, 0.00112637, -0.00234151, -0.00072421]) [ test_j1302_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff Maj, Min, PA vs 6.2.1 is [0.014410739880193054, 0.0030767082376187992, 0.06023197020698886] ( Pass : should be < 0.1) (raw actual/expected values: [2.026454210281372, 1.619074821472168, 7.9184465408325195]/[1.99766636, 1.61410868, 8.42595863])

Mem Use: 315.59375 MiB

test_j1302_mtmfs

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

Elapsed Time: 877.8628330230713 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.2665557056341532e-05 ( Pass : should be < 0.1) (raw actual/expected values: 0.31961721181869507/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.10817003512978554 ( Pass : should be < 0.2) (raw actual/expected values: 0.06057003512978554/-0.0476) [ test_j1302_mtmfs ] Abs Diff alpha vs. 6.2.1 image is 0.0002214651297855405 ( Pass : should be < 0.2) (raw actual/expected values: 0.06057003512978554/0.06034857) [ test_j1302_mtmfs ] Frac Diff Maj, Min, PA vs 6.2.1 is [1.2747511586951134e-09, 1.2464764565995781e-09, 1.274943184034465e-10] ( Pass : should be < 0.1) (raw actual/expected values: [3.150028944015503, 2.591778516769409, 11.412883758544922]/[3.15002894, 2.59177852, 11.41288376]) [ test_j1302_mtmfs ] count('MODEL_DATA') is 1 ( Pass : should be 1, Exact: True )

Mem Use: 0.62109375 MiB

test_j1302_ql

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

Elapsed Time: 804.152197599411 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 [0.03704900744435027] ( Pass : should be < [0.1]) (raw actual/expected values: 0.3213367462158203/[0.3337]) [ test_j1302_ql ] Frac Diff F_nu vs. 6.2.1 image is [0.002457735274344526] ( Pass : should be < [0.1]) (raw actual/expected values: 0.3213367462158203/[0.3205489218235016]) [ test_j1302_ql ] Frac Diff Maj, Min, PA vs 6.2.1 is [6.118642719743648e-07, 2.367592555597574e-06, 1.2168819645515459e-05] ( Pass : should be < 0.1) (raw actual/expected values: [3.1236493587493896, 2.620096206665039, 13.792744636535645]/[3.12365127, 2.62010241, 13.79291248])

Mem Use: 1.4453125 MiB

test_j1927_mosaic_cube

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

Elapsed Time: 15556.42256808281 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.0658617027220473 ( Pass : should be < 0.1) (raw actual/expected values: 0.8882721068816052/0.9509) [ test_j1927_mosaic_cube ] Frac Diff F_nu vs. 6.2.1 image is 9.507794645434835e-05 ( Pass : should be < 0.1) (raw actual/expected values: 0.8882721068816052/0.88835657) [ test_j1927_mosaic_cube ] Abs Diff alpha vs. on-axis is 0.05121765053167787 ( Pass : should be < 0.2) (raw actual/expected values: 0.41131765053167785/0.3601) [ test_j1927_mosaic_cube ] Abs Diff alpha vs. 6.2.1 image is 0.0007048794683221504 ( Pass : should be < 0.2) (raw actual/expected values: 0.41131765053167785/0.41202253) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff IQUV vs 6.2.1 is [0.00021162304355661548, 0.000270046717227037, 0.0009542526851552822, 0.009819722844876116] ( Pass : should be < [0.1 0.16326569 1. 0.14952827]) (raw actual/expected values: [0.7555659413337708, 0.0057603451423347, 0.000794297200627625, -0.006349605042487383]/[0.75540608, 0.00575879, 0.0007934, -0.00628786]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 0), Frac Diff Maj, Min, PA vs 6.2.1 is [0.0030857301278356913, 0.0007222850070653711, 0.04683829777110915] ( Pass : should be < 0.1) (raw actual/expected values: [3.665053606033325, 3.029024839401245, -23.843441009521484]/[3.65377903, 3.03121424, -22.77662277]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff IQUV vs 6.2.1 is [0.0001153497613042239, 4.613063307735884e-05, 8.178594781789252e-06, 0.04347684735705346] ( Pass : should be < [0.1 0.17989271 0.15549112 0.12771143]) (raw actual/expected values: [0.8643968105316162, 0.005226761102676392, -0.006046779453754425, -0.007041933014988899]/[0.86449653, 0.00522652, -0.00604673, -0.00736201]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 1), Frac Diff Maj, Min, PA vs 6.2.1 is [0.000252757728101065, 0.0033029517209369062, 0.04842347672051004] ( Pass : should be < 0.1) (raw actual/expected values: [2.655961275100708, 2.23101544380188, -28.002300262451172]/[2.65529013, 2.2384088, -26.70895958]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff IQUV vs 6.2.1 is [0.00018317006468233856, 0.0011831779552823225, 0.007697839054717622, 0.02464209567842369] ( Pass : should be < [0.1 0.17300434 0.1 0.26293629]) (raw actual/expected values: [0.9558671116828918, 0.005428189877420664, -0.009582056663930416, -0.003663935698568821]/[0.95604223, 0.00543462, -0.00965639, -0.00357582]) [ test_j1927_mosaic_cube ] Stokes Comparison (spw 2), Frac Diff Maj, Min, PA vs 6.2.1 is [0.0013304395597801193, 0.005160215228592125, 0.01120997146896644] ( Pass : should be < 0.1) (raw actual/expected values: [2.0644712448120117, 1.6848171949386597, -25.778289794921875]/[2.06172824, 1.69355631, -25.49251938])

Mem Use: 8.8828125 MiB

test_j1927_mtmfs

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

Elapsed Time: 2198.6927392482758 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.00025930102158943925 ( Pass : should be < 0.1) (raw actual/expected values: 0.4129445552825928/0.41305166) [ test_j1927_mtmfs ] Abs Diff alpha vs. on-axis is 0.08548902311325074 ( Pass : should be < 0.2) (raw actual/expected values: 0.46508902311325073/0.3796) [ test_j1927_mtmfs ] Abs Diff alpha vs. 6.2.1 image is 0.00011604688674926589 ( Pass : should be < 0.2) (raw actual/expected values: 0.46508902311325073/0.46520507) [ test_j1927_mtmfs ] Frac Diff Maj, Min, PA vs 6.2.1 is [1.8401506739554768e-09, 2.285955441074116e-09, 7.435039114081441e-12] ( Pass : should be < 0.1) (raw actual/expected values: [2.4488351345062256, 2.0478949546813965, -23.64226531982422]/[2.44883513, 2.04789495, -23.64226532]) [ test_j1927_mtmfs ] count('MODEL_DATA') is 1 ( Pass : should be 1, Exact: True )

Mem Use: 0.421875 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