Build: #124 was successful
Job: Build and Package RHEL6 was successful
Build log
The build generated 3,221 lines of output.The output is too long and has been truncated to the last 1,000 lines. Download or view full build log.
| 30-Jul-2019 15:10:09 | + if numpy.isnan(val) or numpy.isinf(val): |
| 30-Jul-2019 15:10:09 | + out=False |
| 30-Jul-2019 15:10:09 | + if bound == None and val != None: |
| 30-Jul-2019 15:10:09 | + out = False |
| 30-Jul-2019 15:10:09 | + if bound != None and val == None: |
| 30-Jul-2019 15:10:09 | + out = False |
| 30-Jul-2019 15:10:09 | + if out==True and val != None: |
| 30-Jul-2019 15:10:09 | + if val > bound: |
| 30-Jul-2019 15:10:09 | + out=False |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + pstr = "[ {} ] {} is {} ( {} : should be less than {} )\n".format(testname, valname, str(val), verdict(out), str(bound)) |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + logging.info(pstr) |
| 30-Jul-2019 15:10:09 | + return out, pstr |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | +def check_val_greater_than(val, bound, valname='Value'): |
| 30-Jul-2019 15:10:09 | + pstr = '' |
| 30-Jul-2019 15:10:09 | + if casa5: |
| 30-Jul-2019 15:10:09 | + testname = inspect.stack()[2][3] # Make Sure this is correct |
| 30-Jul-2019 15:10:09 | + else: |
| 30-Jul-2019 15:10:09 | + testname = "TODO" |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + out = True |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + if numpy.isnan(val) or numpy.isinf(val): |
| 30-Jul-2019 15:10:09 | + out=False |
| 30-Jul-2019 15:10:09 | + if bound == None and val != None: |
| 30-Jul-2019 15:10:09 | + out = False |
| 30-Jul-2019 15:10:09 | + if bound != None and val == None: |
| 30-Jul-2019 15:10:09 | + out = False |
| 30-Jul-2019 15:10:09 | + if out==True and val != None: |
| 30-Jul-2019 15:10:09 | + if val < bound: |
| 30-Jul-2019 15:10:09 | + out=False |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + pstr = "[ {} ] {} is {} ( {} : should be greater than {} )\n".format(testname, valname, str(val), verdict(out), str(bound)) |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + logging.info(pstr) |
| 30-Jul-2019 15:10:09 | + return out, pstr |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | def check_ims(imlist,truth): |
| 30-Jul-2019 15:10:09 | if casa5: |
| 30-Jul-2019 15:10:09 | testname = inspect.stack()[2][3] |
| 30-Jul-2019 15:10:09 | @@ -1466,17 +1514,17 @@ def check_im_keywords(imname, check_misc=True, check_extended=True): |
| 30-Jul-2019 15:10:09 | # Looks like a refconcat image, nothing to check |
| 30-Jul-2019 15:10:09 | #return '' |
| 30-Jul-2019 15:10:09 | # make a bit more informative |
| 30-Jul-2019 15:10:09 | - pstr = 'Looks like it is a refconcat image. Skipping the imageinfo keywords check.' |
| 30-Jul-2019 15:10:09 | + pstr = 'Looks like it is a refconcat image. Skipping the imageinfo keywords check.\n' |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | else: |
| 30-Jul-2019 15:10:09 | - pstr = 'Cannot open image table to check keywords: {0}'.format(imname) |
| 30-Jul-2019 15:10:09 | + pstr = 'Cannot open image table to check keywords: {0}\n'.format(imname) |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | finally: |
| 30-Jul-2019 15:10:09 | tbt.close() |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | pstr = '' |
| 30-Jul-2019 15:10:09 | if len(keys) <= 0: |
| 30-Jul-2019 15:10:09 | - pstr += ('No keywords found ({0})'.format(verdict(False))) |
| 30-Jul-2019 15:10:09 | + pstr += ('No keywords found ({0})\n'.format(verdict(False))) |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | # Records that need to be present |
| 30-Jul-2019 15:10:09 | @@ -1605,46 +1653,46 @@ def check_ref_freq(imname,theval=0, epsilon=0.05): |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | ################################### |
| 30-Jul-2019 15:10:09 | -def check_imexist(imexist): |
| 30-Jul-2019 15:10:09 | +def check_imexist(imgexist): |
| 30-Jul-2019 15:10:09 | pstr = '' |
| 30-Jul-2019 15:10:09 | - if imexist != None: |
| 30-Jul-2019 15:10:09 | - if type(imexist)==list: |
| 30-Jul-2019 15:10:09 | - pstr += check_ims(imexist, True) |
| 30-Jul-2019 15:10:09 | - print "pstr after checkims=",pstr |
| 30-Jul-2019 15:10:09 | - pstr += check_keywords(imexist) |
| 30-Jul-2019 15:10:09 | - print "pstr after check_keywords=",pstr |
| 30-Jul-2019 15:10:09 | + if imgexist != None: |
| 30-Jul-2019 15:10:09 | + if type(imgexist)==list: |
| 30-Jul-2019 15:10:09 | + pstr += check_ims(imgexist, True) |
| 30-Jul-2019 15:10:09 | + print("pstr after checkims = {}".format(pstr)) |
| 30-Jul-2019 15:10:09 | + pstr += check_keywords(imgexist) |
| 30-Jul-2019 15:10:09 | + print("pstr after check_keywords = {}".format(pstr)) |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | -def check_imexistnot(imexistnot): |
| 30-Jul-2019 15:10:09 | +def check_imexistnot(imgexistnot): |
| 30-Jul-2019 15:10:09 | pstr = '' |
| 30-Jul-2019 15:10:09 | - if imexistnot != None: |
| 30-Jul-2019 15:10:09 | - if type(imexistnot)==list: |
| 30-Jul-2019 15:10:09 | - pstr += check_ims(imexistnot, False) |
| 30-Jul-2019 15:10:09 | + if imgexistnot != None: |
| 30-Jul-2019 15:10:09 | + if type(imgexistnot)==list: |
| 30-Jul-2019 15:10:09 | + pstr += check_ims(imgexistnot, False) |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | -def check_imval(imval): |
| 30-Jul-2019 15:10:09 | +def check_imval(imgval, epsilon=0.05): |
| 30-Jul-2019 15:10:09 | pstr = '' |
| 30-Jul-2019 15:10:09 | - if imval != None: |
| 30-Jul-2019 15:10:09 | - if type(imval)==list: |
| 30-Jul-2019 15:10:09 | - for ii in imval: |
| 30-Jul-2019 15:10:09 | + if imgval != None: |
| 30-Jul-2019 15:10:09 | + if type(imgval)==list: |
| 30-Jul-2019 15:10:09 | + for ii in imgval: |
| 30-Jul-2019 15:10:09 | if type(ii)==tuple and len(ii)==3: |
| 30-Jul-2019 15:10:09 | - pstr += check_pix_val(ii[0],ii[1],ii[2]) |
| 30-Jul-2019 15:10:09 | + pstr += check_pix_val(ii[0],ii[1],ii[2],epsilon=epsilon) |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | -def check_imvalexact(imvalexact): |
| 30-Jul-2019 15:10:09 | +def check_imvalexact(imgvalexact, epsilon=0.05): |
| 30-Jul-2019 15:10:09 | pstr = '' |
| 30-Jul-2019 15:10:09 | - if imvalexact != None: |
| 30-Jul-2019 15:10:09 | - if type(imvalexact)==list: |
| 30-Jul-2019 15:10:09 | - for ii in imvalexact: |
| 30-Jul-2019 15:10:09 | + if imgvalexact != None: |
| 30-Jul-2019 15:10:09 | + if type(imgvalexact)==list: |
| 30-Jul-2019 15:10:09 | + for ii in imgvalexact: |
| 30-Jul-2019 15:10:09 | if type(ii)==tuple and len(ii)==3: |
| 30-Jul-2019 15:10:09 | - pstr += check_pix_val(ii[0],ii[1],ii[2], exact=True) |
| 30-Jul-2019 15:10:09 | + pstr += check_pix_val(ii[0],ii[1],ii[2], exact=True,epsilon=epsilon) |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | -def check_immask(immask): |
| 30-Jul-2019 15:10:09 | +def check_immask(imgmask): |
| 30-Jul-2019 15:10:09 | pstr = '' |
| 30-Jul-2019 15:10:09 | - if immask != None: |
| 30-Jul-2019 15:10:09 | - if type(immask)==list: |
| 30-Jul-2019 15:10:09 | - for ii in immask: |
| 30-Jul-2019 15:10:09 | + if imgmask != None: |
| 30-Jul-2019 15:10:09 | + if type(imgmask)==list: |
| 30-Jul-2019 15:10:09 | + for ii in imgmask: |
| 30-Jul-2019 15:10:09 | if type(ii)==tuple and len(ii)==3: |
| 30-Jul-2019 15:10:09 | pstr += check_pixmask(ii[0],ii[1],ii[2]) |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | @@ -1654,7 +1702,7 @@ def check_tabcache(tabcache): |
| 30-Jul-2019 15:10:09 | if tabcache==True: |
| 30-Jul-2019 15:10:09 | opentabs = tb.showcache() |
| 30-Jul-2019 15:10:09 | if len(opentabs)>0 : |
| 30-Jul-2019 15:10:09 | - pstr += "["+inspect.stack()[1][3]+"] " + verdict(False) + ": Found open tables after run " |
| 30-Jul-2019 15:10:09 | + pstr += "["+inspect.stack()[1][3]+"] " + verdict(False) + ": Found open tables after run \n" |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def check_stopcode(stopcode): |
| 30-Jul-2019 15:10:09 | @@ -1662,7 +1710,7 @@ def check_stopcode(stopcode): |
| 30-Jul-2019 15:10:09 | if stopcode != None: |
| 30-Jul-2019 15:10:09 | if type(stopcode)==int: |
| 30-Jul-2019 15:10:09 | stopstr = "["+inspect.stack()[1][3]+"] Stopcode is " + str(ret['stopcode']) + " (" + verdict(ret['stopcode']==stopcode) + " : should be " + str(stopcode) + ")\n" |
| 30-Jul-2019 15:10:09 | - print stopstr |
| 30-Jul-2019 15:10:09 | + print(stopstr) |
| 30-Jul-2019 15:10:09 | pstr += stopstr |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | @@ -1676,18 +1724,18 @@ def check_reffreq(reffreq): |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | -def checkall( ret=None, peakres=None, modflux=None, iterdone=None, nmajordone=None, imexist=None, imexistnot=None, imval=None, imvalexact=None, immask=None, tabcache=True, stopcode=None, reffreq=None ): |
| 30-Jul-2019 15:10:09 | +def checkall( ret=None, peakres=None, modflux=None, iterdone=None, nmajordone=None, imgexist=None, imgexistnot=None, imgval=None, imgvalexact=None, imgmask=None, tabcache=True, stopcode=None, reffreq=None, epsilon=0.05 ): |
| 30-Jul-2019 15:10:09 | """ |
| 30-Jul-2019 15:10:09 | ret=None, |
| 30-Jul-2019 15:10:09 | peakres=None, # a float |
| 30-Jul-2019 15:10:09 | modflux=None, # a float |
| 30-Jul-2019 15:10:09 | iterdone=None, # an int |
| 30-Jul-2019 15:10:09 | nmajordone=None, # an int |
| 30-Jul-2019 15:10:09 | - imexist=None, # list of image names |
| 30-Jul-2019 15:10:09 | - imexistnot=None, # list of image names |
| 30-Jul-2019 15:10:09 | - imval=None, # list of tuples of (imagename,val,pos) |
| 30-Jul-2019 15:10:09 | - imvalexact=None, # list of tuples of (imagename,val,pos) |
| 30-Jul-2019 15:10:09 | - immask=None, #list of tuples to check mask value |
| 30-Jul-2019 15:10:09 | + imgexist=None, # list of image names |
| 30-Jul-2019 15:10:09 | + imgexistnot=None, # list of image names |
| 30-Jul-2019 15:10:09 | + imgval=None, # list of tuples of (imagename,val,pos) |
| 30-Jul-2019 15:10:09 | + imgvalexact=None, # list of tuples of (imagename,val,pos) |
| 30-Jul-2019 15:10:09 | + imgmask=None, #list of tuples to check mask value |
| 30-Jul-2019 15:10:09 | tabcache=True, |
| 30-Jul-2019 15:10:09 | stopcode=None, |
| 30-Jul-2019 15:10:09 | reffreq=None # list of tuples of (imagename, reffreq) |
| 30-Jul-2019 15:10:09 | @@ -1712,18 +1760,25 @@ def checkall( ret=None, peakres=None, modflux=None, iterdone=None, nmajordone=No |
| 30-Jul-2019 15:10:09 | except Exception as e: |
| 30-Jul-2019 15:10:09 | logging.info(ret) |
| 30-Jul-2019 15:10:09 | raise |
| 30-Jul-2019 15:10:09 | - |
| 30-Jul-2019 15:10:09 | - pstr += check_imexist(imexist) |
| 30-Jul-2019 15:10:09 | - pstr += check_imexistnot(imexistnot) |
| 30-Jul-2019 15:10:09 | - pstr += check_imval(imval) |
| 30-Jul-2019 15:10:09 | - pstr += check_imvalexact(imvalexact) |
| 30-Jul-2019 15:10:09 | - pstr += check_immask(immask) |
| 30-Jul-2019 15:10:09 | + logging.info("Epsilon: {}".format(epsilon)) |
| 30-Jul-2019 15:10:09 | + pstr += check_imexist(imgexist) |
| 30-Jul-2019 15:10:09 | + pstr += check_imexistnot(imgexistnot) |
| 30-Jul-2019 15:10:09 | + pstr += check_imval(imgval,epsilon=epsilon) |
| 30-Jul-2019 15:10:09 | + pstr += check_imvalexact(imgvalexact,epsilon=epsilon) |
| 30-Jul-2019 15:10:09 | + pstr += check_immask(imgmask) |
| 30-Jul-2019 15:10:09 | pstr += check_tabcache(tabcache) |
| 30-Jul-2019 15:10:09 | pstr += check_stopcode(stopcode) |
| 30-Jul-2019 15:10:09 | pstr += check_reffreq(reffreq) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | return pstr |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | +def check_final(pstr=""): |
| 30-Jul-2019 15:10:09 | + if not isinstance(pstr, six.string_types): |
| 30-Jul-2019 15:10:09 | + return False |
| 30-Jul-2019 15:10:09 | + casalog.post(pstr,'INFO') |
| 30-Jul-2019 15:10:09 | + if( pstr.count("Fail") > 0 ): |
| 30-Jul-2019 15:10:09 | + return False |
| 30-Jul-2019 15:10:09 | + return True |
| 30-Jul-2019 15:10:09 | ############################################################################################ |
| 30-Jul-2019 15:10:09 | ################################## Decorators ################################## |
| 30-Jul-2019 15:10:09 | ############################################################################################ |
| 30-Jul-2019 15:10:09 | @@ -1765,29 +1820,9 @@ def skipIfMissingModule(required_module,strict=False): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | #import casaTestHelper |
| 30-Jul-2019 15:10:09 | #@casaTestHelper.time_execution |
| 30-Jul-2019 15:10:09 | -def time_execution(out_dict): |
| 30-Jul-2019 15:10:09 | - def time_decorator(function): |
| 30-Jul-2019 15:10:09 | - ''' |
| 30-Jul-2019 15:10:09 | - Decorator: time execution of test |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | - Example: |
| 30-Jul-2019 15:10:09 | - @casaTestHelper.time_execution |
| 30-Jul-2019 15:10:09 | - def test_test(self): |
| 30-Jul-2019 15:10:09 | - ''' |
| 30-Jul-2019 15:10:09 | - @wraps(function) |
| 30-Jul-2019 15:10:09 | - def function_timer(*args, **kwargs): |
| 30-Jul-2019 15:10:09 | - t0 = time.time() |
| 30-Jul-2019 15:10:09 | - result = function(*args, **kwargs) |
| 30-Jul-2019 15:10:09 | - t1 = time.time() |
| 30-Jul-2019 15:10:09 | - #print ("Total time running %s: %s seconds" % (function.__name__, str(t1-t0))) |
| 30-Jul-2019 15:10:09 | - casalog.post("Total time running {}: {} seconds".format(function.__name__, str(t1-t0))) |
| 30-Jul-2019 15:10:09 | - out_dict[function.__name__]['runtime'] = t1-t0 |
| 30-Jul-2019 15:10:09 | - return result |
| 30-Jul-2019 15:10:09 | - |
| 30-Jul-2019 15:10:09 | - return function_timer |
| 30-Jul-2019 15:10:09 | - return time_decorator |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | -def time_execution_alternative(out_dict): |
| 30-Jul-2019 15:10:09 | +def time_execution(out_dict): |
| 30-Jul-2019 15:10:09 | # TODO Ver if this is the better option |
| 30-Jul-2019 15:10:09 | def time_decorator(function): |
| 30-Jul-2019 15:10:09 | ''' |
| 30-Jul-2019 15:10:09 | @@ -1802,6 +1837,7 @@ def time_execution_alternative(out_dict): |
| 30-Jul-2019 15:10:09 | failed = False |
| 30-Jul-2019 15:10:09 | result = None |
| 30-Jul-2019 15:10:09 | t0 = time.time() |
| 30-Jul-2019 15:10:09 | + print(out_dict) |
| 30-Jul-2019 15:10:09 | try: |
| 30-Jul-2019 15:10:09 | result = function(*args, **kwargs) |
| 30-Jul-2019 15:10:09 | except: |
| 30-Jul-2019 15:10:09 | @@ -1815,8 +1851,8 @@ def time_execution_alternative(out_dict): |
| 30-Jul-2019 15:10:09 | t1 = time.time() |
| 30-Jul-2019 15:10:09 | #print ("Total time running %s: %s seconds" % (function.__name__, str(t1-t0))) |
| 30-Jul-2019 15:10:09 | casalog.post("Total time running {}: {} seconds".format(function.__name__, str(t1-t0))) |
| 30-Jul-2019 15:10:09 | - print('======================================================') |
| 30-Jul-2019 15:10:09 | - print(function.__name__) |
| 30-Jul-2019 15:10:09 | + #print('======================================================') |
| 30-Jul-2019 15:10:09 | + #print(function.__name__) |
| 30-Jul-2019 15:10:09 | out_dict[function.__name__]['runtime'] = t1-t0 |
| 30-Jul-2019 15:10:09 | out_dict[function.__name__]['status'] = True |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | @@ -1928,7 +1964,7 @@ def mem_use_deco(out_dict): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def stats_dict(out_dict): |
| 30-Jul-2019 15:10:09 | def stats_decorator(function): |
| 30-Jul-2019 15:10:09 | - @time_execution_alternative(out_dict) |
| 30-Jul-2019 15:10:09 | + @time_execution(out_dict) |
| 30-Jul-2019 15:10:09 | #@cpu_usage(out_dict) |
| 30-Jul-2019 15:10:09 | #@peakmem(out_dict) |
| 30-Jul-2019 15:10:09 | @mem_use_deco(out_dict) |
| 30-Jul-2019 15:10:09 | diff --git a/tests/tasks/test_simobserve.py b/tests/tasks/test_simobserve.py |
| 30-Jul-2019 15:10:09 | index 8588757..8ab3e7e 100644 |
| 30-Jul-2019 15:10:09 | --- a/tests/tasks/test_simobserve.py |
| 30-Jul-2019 15:10:09 | +++ b/tests/tasks/test_simobserve.py |
| 30-Jul-2019 15:10:09 | @@ -55,7 +55,7 @@ class simobserve_unittest_base(unittest.TestCase): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | thistask = "simobserve" |
| 30-Jul-2019 15:10:09 | imkeys=['max','mean','min','npts','rms','blc','blcf','trc','trcf','sigma','sum','sumsq'] |
| 30-Jul-2019 15:10:09 | - # relative and ablsolute tolerance |
| 30-Jul-2019 15:10:09 | + # relative and absolute tolerance |
| 30-Jul-2019 15:10:09 | # (atol=0. means to ignore absolute tolerance) |
| 30-Jul-2019 15:10:09 | rtol = 5.0e-3 |
| 30-Jul-2019 15:10:09 | atol = 0. |
| 30-Jul-2019 15:10:09 | @@ -83,7 +83,6 @@ class simobserve_unittest_base(unittest.TestCase): |
| 30-Jul-2019 15:10:09 | _ms.close() |
| 30-Jul-2019 15:10:09 | return stats[stats.keys()[0]] |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | - # TODO: need to check image axes |
| 30-Jul-2019 15:10:09 | def _check_imstats(self, name, ref, rtol=None, atol=None): |
| 30-Jul-2019 15:10:09 | # ref: a dictionary of reference statistics or reference image name |
| 30-Jul-2019 15:10:09 | # name: the name of image to compare statistics |
| 30-Jul-2019 15:10:09 | @@ -110,10 +109,27 @@ class simobserve_unittest_base(unittest.TestCase): |
| 30-Jul-2019 15:10:09 | (key, type(stats[key]), str(stats[key]), str(ref[key])))) |
| 30-Jul-2019 15:10:09 | message="image statistic '%s' does not match: %s (expected: %s)" % \ |
| 30-Jul-2019 15:10:09 | (key, str(stats[key]), str(ref[key])) |
| 30-Jul-2019 15:10:09 | - if type(stats[key])==str: |
| 30-Jul-2019 15:10:09 | - self.assertEqual(stats[key],ref[key], |
| 30-Jul-2019 15:10:09 | - msg=message) |
| 30-Jul-2019 15:10:09 | + if type(stats[key])==str: |
| 30-Jul-2019 15:10:09 | + # only maxposf, minposf, blcf, trcf return <type 'str'> |
| 30-Jul-2019 15:10:09 | + # these are actually all lists |
| 30-Jul-2019 15:10:09 | + ax_stats = [x.strip() for x in stats[key].split(',')] |
| 30-Jul-2019 15:10:09 | + ax_ref = [x.strip() for x in ref[key].split(',')] |
| 30-Jul-2019 15:10:09 | + # compare dimension of image axes |
| 30-Jul-2019 15:10:09 | + self.assertEqual(len(ax_stats),len(ax_ref),msg=message) |
| 30-Jul-2019 15:10:09 | + # extract, compare numerical data from axis world coordinates |
| 30-Jul-2019 15:10:09 | + for kk in zip(ax_stats, ax_ref): |
| 30-Jul-2019 15:10:09 | + # only check the first element in tuple |
| 30-Jul-2019 15:10:09 | + if qa.isquantity(kk[0]): |
| 30-Jul-2019 15:10:09 | + # test and reference numbers |
| 30-Jul-2019 15:10:09 | + s_val = qa.quantity(kk[0])['value'] |
| 30-Jul-2019 15:10:09 | + r_val = qa.quantity(kk[1])['value'] |
| 30-Jul-2019 15:10:09 | + ret=numpy.allclose(s_val,r_val, |
| 30-Jul-2019 15:10:09 | + rtol=rtol,atol=atol) |
| 30-Jul-2019 15:10:09 | + self.assertEqual(ret,True,msg=message) |
| 30-Jul-2019 15:10:09 | + else: # should only be Stokes axis |
| 30-Jul-2019 15:10:09 | + self.assertEqual(kk[0],kk[1],msg=message) |
| 30-Jul-2019 15:10:09 | else: |
| 30-Jul-2019 15:10:09 | + # not a string so expect numpy arrays |
| 30-Jul-2019 15:10:09 | ret=numpy.allclose(stats[key],ref[key], |
| 30-Jul-2019 15:10:09 | rtol=rtol,atol=atol) |
| 30-Jul-2019 15:10:09 | self.assertEqual(ret,True,msg=message) |
| 30-Jul-2019 15:10:09 | @@ -261,6 +277,7 @@ class simobserve_sky(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | #pass |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | # Tests of skymodel simulations |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"'int' object is not subscriptable. In task_simobserve. - needs to be fixed") |
| 30-Jul-2019 15:10:09 | def testSky_skymodel(self): |
| 30-Jul-2019 15:10:09 | """Test skymodel simulation: only modify model""" |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | @@ -287,6 +304,7 @@ class simobserve_sky(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | self._get_data_prefix(antennalist,self.project)+".skymodel" |
| 30-Jul-2019 15:10:09 | self._check_imstats(currmodel, self.refmodel) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"'int' object is not subscriptable. In task_simobserve. - needs to be fixed") |
| 30-Jul-2019 15:10:09 | def testSky_almaptg(self): |
| 30-Jul-2019 15:10:09 | """Test skymodel simulation: only setpointing (maptype='ALMA')""" |
| 30-Jul-2019 15:10:09 | skymodel = self.refmodel |
| 30-Jul-2019 15:10:09 | @@ -307,6 +325,7 @@ class simobserve_sky(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | refptg = self.refpref + "alma.alma.out01.ptg.txt" |
| 30-Jul-2019 15:10:09 | self._check_ptgfile(currptg, refptg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"'int' object is not subscriptable. In task_simobserve. - needs to be fixed") |
| 30-Jul-2019 15:10:09 | def testSky_hexptg(self): |
| 30-Jul-2019 15:10:09 | """Test skymodel simulation: only setpointing (maptype='hexagonal')""" |
| 30-Jul-2019 15:10:09 | skymodel = self.refmodel |
| 30-Jul-2019 15:10:09 | @@ -327,6 +346,7 @@ class simobserve_sky(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | refptg = self.refpref + "hex.aca.i.ptg.txt" |
| 30-Jul-2019 15:10:09 | self._check_ptgfile(currptg, refptg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Both antennnalist and sdantlist are defined. - needs to be fixed") |
| 30-Jul-2019 15:10:09 | def testSky_sqptg(self): |
| 30-Jul-2019 15:10:09 | """Test skymodel simulation: only setpointing (maptype='square')""" |
| 30-Jul-2019 15:10:09 | skymodel = self.refmodel |
| 30-Jul-2019 15:10:09 | @@ -481,6 +501,7 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | project = simobserve_unittest_base.thistask+"_comp" |
| 30-Jul-2019 15:10:09 | incomp = "core5ps.clist" |
| 30-Jul-2019 15:10:09 | compwidth = "10MHz" |
| 30-Jul-2019 15:10:09 | + comp_nchan = 1 |
| 30-Jul-2019 15:10:09 | direction = "J2000 19h00m00 -23d00m00" |
| 30-Jul-2019 15:10:09 | sdantlist = "aca.tp.cfg" |
| 30-Jul-2019 15:10:09 | antlist = "alma.out01.cfg" |
| 30-Jul-2019 15:10:09 | @@ -505,6 +526,11 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | # reference simulated MS |
| 30-Jul-2019 15:10:09 | self.refms_sd = self.refpref_sd+".sd.ms" |
| 30-Jul-2019 15:10:09 | self.refms_int = self.refpref_int+".ms" |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + # new data for comp_nchan > 1 |
| 30-Jul-2019 15:10:09 | + self.refmodel_int_8ch = self.refpref_int+".8ch.compskymodel" |
| 30-Jul-2019 15:10:09 | + self.refms_int_8ch = self.refpref_int+".8ch.ms" |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | # copy input components list |
| 30-Jul-2019 15:10:09 | self._copy_input(self.incomp) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | @@ -518,11 +544,12 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | """Test complist simulation: only generating input model""" |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | antennalist="alma.out01.cfg" # necessary even if only modifymodel |
| 30-Jul-2019 15:10:09 | sdantlist = "" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=True,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | antennalist=antennalist,sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | thermalnoise="",graphics=self.graphics) |
| 30-Jul-2019 15:10:09 | @@ -536,13 +563,14 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | """Test complist simulation: only setpointing (maptype='ALMA')""" |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = True |
| 30-Jul-2019 15:10:09 | maptype = "ALMA" |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | antennalist = "alma.out01.cfg" |
| 30-Jul-2019 15:10:09 | sdantlist = "" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,maptype=maptype, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -558,13 +586,14 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | """Test complist simulation: only setpointing (maptype='hexagonal')""" |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = True |
| 30-Jul-2019 15:10:09 | maptype = "hexagonal" |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | antennalist = "aca.i.cfg" |
| 30-Jul-2019 15:10:09 | sdantlist = "" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,maptype=maptype, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -576,17 +605,19 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | refptg = self.refpref + "hex.aca.i.ptg.txt" |
| 30-Jul-2019 15:10:09 | self._check_ptgfile(currptg, refptg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | def testComp_sqptg(self): |
| 30-Jul-2019 15:10:09 | """Test complist simulation: only setpointing (maptype='square')""" |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = True |
| 30-Jul-2019 15:10:09 | maptype = "square" |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | antennalist = "" |
| 30-Jul-2019 15:10:09 | sdantlist = self.sdantlist |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,maptype=maptype, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -610,7 +641,7 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | sdantlist = self.sdantlist |
| 30-Jul-2019 15:10:09 | totaltime = "144s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,ptgfile=ptgfile, |
| 30-Jul-2019 15:10:09 | integration=integration,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist,totaltime=totaltime, |
| 30-Jul-2019 15:10:09 | @@ -626,6 +657,7 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | """Test complist simulation: only observation (INT)""" |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = False |
| 30-Jul-2019 15:10:09 | ptgfile = self.refpref_int + ".ptg.txt" |
| 30-Jul-2019 15:10:09 | integration = "4s" |
| 30-Jul-2019 15:10:09 | @@ -633,7 +665,7 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | antennalist = 'alma.out01.cfg' |
| 30-Jul-2019 15:10:09 | totaltime = "28s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,ptgfile=ptgfile, |
| 30-Jul-2019 15:10:09 | integration=integration,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | antennalist=antennalist,totaltime=totaltime, |
| 30-Jul-2019 15:10:09 | @@ -654,7 +686,7 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | leakage = 0.5 |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,ptgfile=ptgfile, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,thermalnoise="", |
| 30-Jul-2019 15:10:09 | leakage=leakage,graphics=self.graphics) |
| 30-Jul-2019 15:10:09 | @@ -664,6 +696,7 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | """Test complist simulation: single dish""" |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | integration = "4s" |
| 30-Jul-2019 15:10:09 | direction = self.direction |
| 30-Jul-2019 15:10:09 | mapsize = ["60arcsec", "60arcsec"] |
| 30-Jul-2019 15:10:09 | @@ -672,7 +705,8 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | sdantlist = "aca.tp.cfg" |
| 30-Jul-2019 15:10:09 | totaltime = "144s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth,setpointings=True, |
| 30-Jul-2019 15:10:09 | + compwidth = |
| 30-Jul-2019 15:10:09 | + compwidth,comp_nchan=comp_nchan,setpointings=True, |
| 30-Jul-2019 15:10:09 | integration=integration,direction=direction, |
| 30-Jul-2019 15:10:09 | mapsize=mapsize,maptype=maptype,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | totaltime=totaltime,antennalist="",sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -690,6 +724,7 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | """Test complist simulation: interferometer""" |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | integration = "4s" |
| 30-Jul-2019 15:10:09 | direction = self.direction |
| 30-Jul-2019 15:10:09 | mapsize = ['20arcsec', '20arcsec'] |
| 30-Jul-2019 15:10:09 | @@ -698,11 +733,12 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | antennalist = 'alma.out01.cfg' |
| 30-Jul-2019 15:10:09 | totaltime = "28s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | - compwidth = compwidth,setpointings=True, |
| 30-Jul-2019 15:10:09 | - integration=integration,direction=direction, |
| 30-Jul-2019 15:10:09 | - mapsize=mapsize,maptype=maptype,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | - totaltime=totaltime,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | - thermalnoise="",graphics=self.graphics) |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | + setpointings=True,integration=integration, |
| 30-Jul-2019 15:10:09 | + direction=direction,mapsize=mapsize,maptype=maptype, |
| 30-Jul-2019 15:10:09 | + obsmode=obsmode,totaltime=totaltime, |
| 30-Jul-2019 15:10:09 | + antennalist=antennalist,thermalnoise="", |
| 30-Jul-2019 15:10:09 | + graphics=self.graphics) |
| 30-Jul-2019 15:10:09 | self.assertTrue(res) |
| 30-Jul-2019 15:10:09 | # compare outputs |
| 30-Jul-2019 15:10:09 | currpref = self.project + "/" + \ |
| 30-Jul-2019 15:10:09 | @@ -711,6 +747,36 @@ class simobserve_comp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | self._check_ptgfile(currpref+".ptg.txt", self.refpref_int+".ptg.txt") |
| 30-Jul-2019 15:10:09 | self._check_msstats(currpref+".ms",self.refms_int) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Failed AlwaysAssert sdirections.nelements() == nptg - needs to be fixed.") |
| 30-Jul-2019 15:10:09 | + def testComp_intNchan(self): |
| 30-Jul-2019 15:10:09 | + """Test complist simulation: interferometer, but with comp_nchan > 1""" |
| 30-Jul-2019 15:10:09 | + complist = self.incomp |
| 30-Jul-2019 15:10:09 | + compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = 8 |
| 30-Jul-2019 15:10:09 | + integration = "4s" |
| 30-Jul-2019 15:10:09 | + direction = self.direction |
| 30-Jul-2019 15:10:09 | + mapsize = ['20arcsec', '20arcsec'] |
| 30-Jul-2019 15:10:09 | + maptype = "ALMA" |
| 30-Jul-2019 15:10:09 | + obsmode = 'int' |
| 30-Jul-2019 15:10:09 | + antennalist = 'alma.out01.cfg' |
| 30-Jul-2019 15:10:09 | + totaltime = "28s" |
| 30-Jul-2019 15:10:09 | + res = simobserve(project=self.project,complist=complist, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | + setpointings=True, |
| 30-Jul-2019 15:10:09 | + integration=integration,direction=direction, |
| 30-Jul-2019 15:10:09 | + mapsize=mapsize,maptype=maptype,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | + totaltime=totaltime,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | + thermalnoise="",graphics=self.graphics) |
| 30-Jul-2019 15:10:09 | + self.assertTrue(res) |
| 30-Jul-2019 15:10:09 | + # compare outputs |
| 30-Jul-2019 15:10:09 | + currpref = self.project + "/" + \ |
| 30-Jul-2019 15:10:09 | + self._get_data_prefix(antennalist,self.project) |
| 30-Jul-2019 15:10:09 | + self._check_imstats(currpref+".compskymodel", |
| 30-Jul-2019 15:10:09 | + self.refmodel_int_8ch) |
| 30-Jul-2019 15:10:09 | + self._check_ptgfile(currpref+".ptg.txt", |
| 30-Jul-2019 15:10:09 | + self.refpref_int+".8ch.ptg.txt") |
| 30-Jul-2019 15:10:09 | + self._check_msstats(currpref+".ms", |
| 30-Jul-2019 15:10:09 | + self.refms_int_8ch) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | ######################################################################## |
| 30-Jul-2019 15:10:09 | @@ -728,6 +794,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | indirection = "J2000 19h00m00 -23d00m00" |
| 30-Jul-2019 15:10:09 | incomp = "ps5.clist" |
| 30-Jul-2019 15:10:09 | compwidth = "10MHz" |
| 30-Jul-2019 15:10:09 | + comp_nchan = 1 |
| 30-Jul-2019 15:10:09 | sdantlist = "aca.tp.cfg" |
| 30-Jul-2019 15:10:09 | antlist = "alma.out01.cfg" |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | @@ -770,6 +837,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | #setpointings = False |
| 30-Jul-2019 15:10:09 | #ptgfile = # necessary even if only modifymodel |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | @@ -777,6 +845,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | sdantlist = "" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=True,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | antennalist=antennalist,sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | thermalnoise="",graphics=self.graphics) |
| 30-Jul-2019 15:10:09 | @@ -792,6 +861,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = True |
| 30-Jul-2019 15:10:09 | maptype = "ALMA" |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | @@ -799,6 +869,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | sdantlist = "" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,maptype=maptype, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -815,6 +886,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = True |
| 30-Jul-2019 15:10:09 | maptype = "hexagonal" |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | @@ -822,6 +894,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | sdantlist = "" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,maptype=maptype, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -838,6 +911,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = True |
| 30-Jul-2019 15:10:09 | maptype = "square" |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | @@ -845,6 +919,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | sdantlist = self.sdantlist |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,maptype=maptype, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -862,6 +937,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = False |
| 30-Jul-2019 15:10:09 | ptgfile = self.refpref_sd + ".ptg.txt" |
| 30-Jul-2019 15:10:09 | integration = "4s" |
| 30-Jul-2019 15:10:09 | @@ -870,6 +946,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | totaltime = "144s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,ptgfile=ptgfile, |
| 30-Jul-2019 15:10:09 | integration=integration,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | sdantlist=sdantlist,totaltime=totaltime, |
| 30-Jul-2019 15:10:09 | @@ -886,6 +963,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = False |
| 30-Jul-2019 15:10:09 | ptgfile = self.refpref_int + ".ptg.txt" |
| 30-Jul-2019 15:10:09 | integration = "4s" |
| 30-Jul-2019 15:10:09 | @@ -894,6 +972,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | totaltime = "28s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,ptgfile=ptgfile, |
| 30-Jul-2019 15:10:09 | integration=integration,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | antennalist=antennalist,totaltime=totaltime, |
| 30-Jul-2019 15:10:09 | @@ -910,11 +989,13 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | setpointings = False |
| 30-Jul-2019 15:10:09 | obsmode = "" |
| 30-Jul-2019 15:10:09 | leakage = 0.5 |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=setpointings,ptgfile=ptgfile, |
| 30-Jul-2019 15:10:09 | obsmode=obsmode,thermalnoise="", |
| 30-Jul-2019 15:10:09 | leakage=leakage,graphics=self.graphics) |
| 30-Jul-2019 15:10:09 | @@ -925,6 +1006,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | integration = "4s" |
| 30-Jul-2019 15:10:09 | mapsize = ["60arcsec", "60arcsec"] |
| 30-Jul-2019 15:10:09 | maptype = "square" |
| 30-Jul-2019 15:10:09 | @@ -933,6 +1015,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | totaltime = "144s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=True,integration=integration, |
| 30-Jul-2019 15:10:09 | mapsize=mapsize,maptype=maptype,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | totaltime=totaltime,antennalist="",sdantlist=sdantlist, |
| 30-Jul-2019 15:10:09 | @@ -951,6 +1034,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | skymodel = self.inmodel |
| 30-Jul-2019 15:10:09 | complist = self.incomp |
| 30-Jul-2019 15:10:09 | compwidth = self.compwidth |
| 30-Jul-2019 15:10:09 | + comp_nchan = self.comp_nchan |
| 30-Jul-2019 15:10:09 | integration = "4s" |
| 30-Jul-2019 15:10:09 | mapsize = ['20arcsec', '20arcsec'] |
| 30-Jul-2019 15:10:09 | maptype = "ALMA" |
| 30-Jul-2019 15:10:09 | @@ -959,6 +1043,7 @@ class simobserve_skycomp(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | totaltime = "28s" |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,skymodel=skymodel, |
| 30-Jul-2019 15:10:09 | complist=complist,compwidth=compwidth, |
| 30-Jul-2019 15:10:09 | + comp_nchan=comp_nchan, |
| 30-Jul-2019 15:10:09 | setpointings=True,integration=integration, |
| 30-Jul-2019 15:10:09 | mapsize=mapsize,maptype=maptype,obsmode=obsmode, |
| 30-Jul-2019 15:10:09 | totaltime=totaltime,antennalist=antennalist, |
| 30-Jul-2019 15:10:09 | @@ -1024,6 +1109,7 @@ class simobserve_noise(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | #-----------------------------------------------------------------# |
| 30-Jul-2019 15:10:09 | # thermalnoise = "tsys-manual" |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"int object is not subscriptable' in task_simobserve - needs to be fixed") |
| 30-Jul-2019 15:10:09 | def testNZ_intMan(self): |
| 30-Jul-2019 15:10:09 | """Test INT thermal noise (tsys-manual)""" |
| 30-Jul-2019 15:10:09 | project = self.project_int |
| 30-Jul-2019 15:10:09 | @@ -1054,6 +1140,7 @@ class simobserve_noise(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | self.assertTrue(abs((msnoise-ananoise)/ananoise) < 1.e-1, \ |
| 30-Jul-2019 15:10:09 | msg=self.anamsg % (msnoise, ananoise)) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Both antennalist and sdantlist are defined. Define one of them. - needs to be fixed") |
| 30-Jul-2019 15:10:09 | def testNZ_sdMan(self): |
| 30-Jul-2019 15:10:09 | """Test SD thermal noise (tsys-manual): standard parameter set""" |
| 30-Jul-2019 15:10:09 | thermalnoise="tsys-manual" |
| 30-Jul-2019 15:10:09 | @@ -1079,6 +1166,7 @@ class simobserve_noise(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | self.assertTrue(abs((msnoise-ananoise)/ananoise) < 1.e-1, \ |
| 30-Jul-2019 15:10:09 | msg=self.anamsg % (msnoise, ananoise)) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Both antennalist and sdantlist are defined. Define one of them. - needs to be fixed") |
| 30-Jul-2019 15:10:09 | def testNZ_sdMan_tau(self): |
| 30-Jul-2019 15:10:09 | """Test SD thermal noise (tsys-manual): tau0=1.5""" |
| 30-Jul-2019 15:10:09 | thermalnoise="tsys-manual" |
| 30-Jul-2019 15:10:09 | @@ -1572,6 +1660,7 @@ class simobserve_badinputs(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Exception message thrown is different - antennnalist/sdantennnalist - fix this") |
| 30-Jul-2019 15:10:09 | def test_noProject(self): |
| 30-Jul-2019 15:10:09 | """Test no project name""" |
| 30-Jul-2019 15:10:09 | project = '' |
| 30-Jul-2019 15:10:09 | @@ -1704,15 +1793,31 @@ class simobserve_badinputs(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | """Test bad compwidth""" |
| 30-Jul-2019 15:10:09 | # not frequency |
| 30-Jul-2019 15:10:09 | compwidth="2arcsec" |
| 30-Jul-2019 15:10:09 | + comp_nchan=1 |
| 30-Jul-2019 15:10:09 | try: |
| 30-Jul-2019 15:10:09 | res = simobserve(project=self.project,complist=self.incomp, |
| 30-Jul-2019 15:10:09 | totaltime=self.tottime,mapsize=self.mapsize, |
| 30-Jul-2019 15:10:09 | - compwidth=compwidth) |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan) |
| 30-Jul-2019 15:10:09 | self.fail(self.failmsg) |
| 30-Jul-2019 15:10:09 | except Exception as e: |
| 30-Jul-2019 15:10:09 | pos=str(e).find("Quantum::operator- unequal units 'GHz, 'arcsec'") |
| 30-Jul-2019 15:10:09 | msg = self.errmsg % str(e) |
| 30-Jul-2019 15:10:09 | self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Exception message thrown is different - unexpected comp_nchan keyword - fix this") |
| 30-Jul-2019 15:10:09 | + def testBad_comp_nchan(self): |
| 30-Jul-2019 15:10:09 | + """Test bad comp_nchan""" |
| 30-Jul-2019 15:10:09 | + compwidth="2arcsec" |
| 30-Jul-2019 15:10:09 | + comp_nchan=self.badnum |
| 30-Jul-2019 15:10:09 | + try: |
| 30-Jul-2019 15:10:09 | + res = simobserve(project=self.project,complist=self.incomp, |
| 30-Jul-2019 15:10:09 | + totaltime=self.tottime,mapsize=self.mapsize, |
| 30-Jul-2019 15:10:09 | + compwidth=compwidth,comp_nchan=comp_nchan) |
| 30-Jul-2019 15:10:09 | + self.fail(self.failmsg) |
| 30-Jul-2019 15:10:09 | + except Exception as e: |
| 30-Jul-2019 15:10:09 | + pos=str(e).find("Parameter verification failed") |
| 30-Jul-2019 15:10:09 | + msg = self.errmsg % str(e) |
| 30-Jul-2019 15:10:09 | + self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | @unittest.skipIf(is_CASA6,"Exception message thrown is different - antennnalist/sdantennnalist - fix this") |
| 30-Jul-2019 15:10:09 | def testBad_ptgfile(self): |
| 30-Jul-2019 15:10:09 | @@ -1835,6 +1940,7 @@ class simobserve_badinputs(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | msg = self.errmsg % str(e) |
| 30-Jul-2019 15:10:09 | self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Exception message thrown is different - antennnalist/sdantennnalist - fix this") |
| 30-Jul-2019 15:10:09 | def testBad_antennalist(self): |
| 30-Jul-2019 15:10:09 | """Test bad antennalist name""" |
| 30-Jul-2019 15:10:09 | antennalist = self.badname |
| 30-Jul-2019 15:10:09 | @@ -1876,6 +1982,7 @@ class simobserve_badinputs(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | msg = self.errmsg % str(e) |
| 30-Jul-2019 15:10:09 | self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Exception message thrown is different - antennnalist/sdantennnalist - fix this") |
| 30-Jul-2019 15:10:09 | def testBad_sdantlist(self): |
| 30-Jul-2019 15:10:09 | """Test bad sdantlist name""" |
| 30-Jul-2019 15:10:09 | obsmode = "sd" |
| 30-Jul-2019 15:10:09 | @@ -1911,6 +2018,7 @@ class simobserve_badinputs(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | msg = self.errmsg % str(e) |
| 30-Jul-2019 15:10:09 | self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Exception message thrown is different - antennnalist/sdantennnalist - fix this") |
| 30-Jul-2019 15:10:09 | def testBad_refdate(self): |
| 30-Jul-2019 15:10:09 | """Test bad refdate""" |
| 30-Jul-2019 15:10:09 | obsmode = "sd" |
| 30-Jul-2019 15:10:09 | @@ -1928,6 +2036,7 @@ class simobserve_badinputs(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | msg = self.errmsg % str(e) |
| 30-Jul-2019 15:10:09 | self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Exception message thrown is different - antennnalist/sdantennnalist - fix this") |
| 30-Jul-2019 15:10:09 | def testBad_hourangle(self): |
| 30-Jul-2019 15:10:09 | """Test bad hourangle""" |
| 30-Jul-2019 15:10:09 | obsmode = "sd" |
| 30-Jul-2019 15:10:09 | @@ -1946,6 +2055,7 @@ class simobserve_badinputs(simobserve_unittest_base): |
| 30-Jul-2019 15:10:09 | self.assertNotEqual(pos,-1,msg=msg) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | # casapy crashes for totaltime < 0 |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(is_CASA6,"Exception message thrown is different - antennnalist/sdantennnalist - fix this") |
| 30-Jul-2019 15:10:09 | def testBad_totaltime(self): |
| 30-Jul-2019 15:10:09 | """Test bad totaltime""" |
| 30-Jul-2019 15:10:09 | obsmode = "sd" |
| 30-Jul-2019 15:10:09 | diff --git a/tests/tasks/test_tclean.py b/tests/tasks/test_tclean.py |
| 30-Jul-2019 15:10:09 | index 64ec8b4..97ca6f4 100644 |
| 30-Jul-2019 15:10:09 | --- a/tests/tasks/test_tclean.py |
| 30-Jul-2019 15:10:09 | +++ b/tests/tasks/test_tclean.py |
| 30-Jul-2019 15:10:09 | @@ -231,7 +231,7 @@ class test_onefield(testref_base): |
| 30-Jul-2019 15:10:09 | """ [onefield] Test_Onefield_multiscale : mfs with multiscale minor cycle """ |
| 30-Jul-2019 15:10:09 | self.prepData('refim_eptwochan.ms') |
| 30-Jul-2019 15:10:09 | ret = tclean(vis=self.msfile,imagename=self.img,imsize=200,cell='8.0arcsec',niter=10,deconvolver='multiscale',scales=[0,20,40,100],interactive=0,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | - report=self.th.checkall(ret=ret, peakres=1.28, modflux=2.87, iterdone=10, imexist=[self.img+'.psf', self.img+'.residual', self.img+'.image',self.img+'.model'], imval=[(self.img+'.psf',1.0,[100,100,0,0])]) |
| 30-Jul-2019 15:10:09 | + report=self.th.checkall(ret=ret, peakres=0.823, modflux=3.816, iterdone=10, imexist=[self.img+'.psf', self.img+'.residual', self.img+'.image',self.img+'.model'], imval=[(self.img+'.psf',1.0,[100,100,0,0])]) |
| 30-Jul-2019 15:10:09 | self.checkfinal(pstr=report) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def test_onefield_mtmfs(self): |
| 30-Jul-2019 15:10:09 | @@ -525,8 +525,14 @@ class test_onefield(testref_base): |
| 30-Jul-2019 15:10:09 | report=self.th.checkall(ret=ret, peakres=0.369, modflux=0.689, iterdone=10, imexist=[self.img+'.psf.tt0', self.img+'.residual.tt0', self.img+'.image.tt0', self.img+'.model.tt0'], imval=[(self.img+'.psf.tt0',1.0,[50,50,0,0]),(self.img+'.image.tt0',1.05,[50,50,0,0])]) |
| 30-Jul-2019 15:10:09 | ## iterdone=11 only because of the return (iterdone_p+1) in MultiTermMatrixCleaner::mtclean() ! |
| 30-Jul-2019 15:10:09 | self.checkfinal(pstr=report) |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + def test_onefield_mtmfs_smallscalebias(self): |
| 30-Jul-2019 15:10:09 | + """ [onefield] Test_Onefield_mtmfs : mt-mfs with minor cycle iterations and smallscalebias = 0.9 """ |
| 30-Jul-2019 15:10:09 | + self.prepData('refim_eptwochan.ms') |
| 30-Jul-2019 15:10:09 | + ret = tclean(vis=self.msfile,imagename=self.img,imsize=200,cell='8.0arcsec',niter=10,deconvolver='mtmfs',nterms=1,interactive=0,parallel=self.parallel,smallscalebias=0.9,scales=[0,20,40,100]) |
| 30-Jul-2019 15:10:09 | + report=self.th.checkall(ret=ret, peakres=0.73153, modflux=2.9194, iterdone=10, imexist=[self.img+'.psf.tt0', self.img+'.residual.tt0', self.img+'.image.tt0', self.img+'.model.tt0'],imval=[(self.img+'.image.tt0',0.526,[100,100,0,0])]) |
| 30-Jul-2019 15:10:09 | + self.checkfinal(pstr=report) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | - |
| 30-Jul-2019 15:10:09 | def test_onefield_gridders(self): |
| 30-Jul-2019 15:10:09 | """ [onefield] Test_Onefield_gridders : Check all single field gridder equivalent names are accepted """ |
| 30-Jul-2019 15:10:09 | self.prepData('refim_twochan.ms') |
| 30-Jul-2019 15:10:09 | @@ -786,7 +792,7 @@ class test_iterbot(testref_base): |
| 30-Jul-2019 15:10:09 | ret2 = tclean(vis=self.msfile,imagename=self.img+'2',imsize=100,cell='8.0arcsec',niter=10,threshold='0.1Jy', interactive=0,deconvolver='hogbom',parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | report2=self.th.checkall(ret=ret2, peakres=0.3530, modflux=0.7719, iterdone=10, nmajordone=2,imexist=[self.img+'2.psf', self.img+'2.residual', self.img+'2.image']) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | - ret3 = tclean(vis=self.msfile,imagename=self.img+'3',imsize=100,cell='8.0arcsec',niter=10,threshold='0.1Jy', interactive=0,deconvolver='multiscale',parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | + ret3 = tclean(vis=self.msfile,imagename=self.img+'3',imsize=100,cell='8.0arcsec',niter=10,threshold='0.1Jy', interactive=0,deconvolver='multiscale',parallel=self.parallel,smallscalebias=0.6) |
| 30-Jul-2019 15:10:09 | report3=self.th.checkall(ret=ret3, peakres=0.3922, modflux=0.7327, iterdone=10, nmajordone=2,imexist=[self.img+'3.psf', self.img+'3.residual', self.img+'3.image']) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | @@ -2057,6 +2063,22 @@ class test_cube(testref_base): |
| 30-Jul-2019 15:10:09 | report=self.th.checkall(imexist=[self.img+'cc.image'],imval=[(self.img+'cc.image',1.5002,[50,50,0,0]) , (self.img+'cc.image',0.769,[50,50,0,19]) ]) |
| 30-Jul-2019 15:10:09 | self.assertTrue( self.th.checkmodelchan(self.msfile,5) > 0.0 and self.th.checkmodelchan(self.msfile,18) > 0.0 ) |
| 30-Jul-2019 15:10:09 | self.checkfinal(report) |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(ParallelTaskHelper.isMPIEnabled(), "Skip the test temporarily") |
| 30-Jul-2019 15:10:09 | + def test_cube_mtmfs_nterms1(self): |
| 30-Jul-2019 15:10:09 | + """ [cube] Test mtmfs with cube and nterms = 1 """ |
| 30-Jul-2019 15:10:09 | + self.prepData('refim_eptwochan.ms') |
| 30-Jul-2019 15:10:09 | + ret = tclean(vis=self.msfile,imagename=self.img+'cc', specmode='cube', imsize=200,cell='8.0arcsec',niter=10,deconvolver='mtmfs',nterms=1,interactive=0,parallel=self.parallel,scales=[0,20,40,100]) |
| 30-Jul-2019 15:10:09 | + report=self.th.checkall(ret=ret, imexist=[self.img+'cc.psf.tt0', self.img+'cc.residual.tt0', self.img+'cc.image.tt0', self.img+'cc.model.tt0'],imval=[(self.img+'cc.image.tt0',1.0,[100,100,0,0]),(self.img+'cc.image.tt0',0.492,[100,100,0,1]),(self.img+'cc.image.tt0',0.281,[100,100,0,2])]) |
| 30-Jul-2019 15:10:09 | + self.checkfinal(report) |
| 30-Jul-2019 15:10:09 | + |
| 30-Jul-2019 15:10:09 | + @unittest.skipIf(ParallelTaskHelper.isMPIEnabled(), "Skip the test temporarily") |
| 30-Jul-2019 15:10:09 | + def test_cubedata_mtmfs_nterms1(self): |
| 30-Jul-2019 15:10:09 | + """ [cube] Test mtmfs with cube data and nterms = 1 """ |
| 30-Jul-2019 15:10:09 | + self.prepData('refim_eptwochan.ms') |
| 30-Jul-2019 15:10:09 | + ret = tclean(vis=self.msfile,imagename=self.img+'cc', specmode='cubedata', imsize=200,cell='8.0arcsec',niter=10,deconvolver='mtmfs',nterms=1,interactive=0,parallel=self.parallel,scales=[0,20,40,100]) |
| 30-Jul-2019 15:10:09 | + report=self.th.checkall(ret=ret, imexist=[self.img+'cc.psf.tt0', self.img+'cc.residual.tt0', self.img+'cc.image.tt0', self.img+'cc.model.tt0'],imval=[(self.img+'cc.image.tt0',1.0,[100,100,0,0]),(self.img+'cc.image.tt0',0.492,[100,100,0,1]),(self.img+'cc.image.tt0',0.281,[100,100,0,2])]) |
| 30-Jul-2019 15:10:09 | + self.checkfinal(report) |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | ############################################## |
| 30-Jul-2019 15:10:09 | ############################################## |
| 30-Jul-2019 15:10:09 | @@ -2747,35 +2769,6 @@ class test_widefield(testref_base): |
| 30-Jul-2019 15:10:09 | #do stokes V too.. |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | -class test_widefield_failing(testref_base): |
| 30-Jul-2019 15:10:09 | - |
| 30-Jul-2019 15:10:09 | - @unittest.skip("Skip failing widefield test") |
| 30-Jul-2019 15:10:09 | - def test_widefield_imagemosaic(self): |
| 30-Jul-2019 15:10:09 | - """ [widefield] Test_Widefield_imagemosaic : Image domain mosaic for single-term mfs (or narrowband) """ |
| 30-Jul-2019 15:10:09 | - self.prepData("refim_mawproject.ms") |
| 30-Jul-2019 15:10:09 | - ret = tclean(vis=self.msfile,spw='1',field='*',imagename=self.img,imsize=512,cell='10.0arcsec',phasecenter="J2000 19:59:28.500 +40.44.01.50", |
| 30-Jul-2019 15:10:09 | - niter=30,gridder='imagemosaic',deconvolver='hogbom',parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | - report=self.th.checkall(imexist=[self.img+'.image', self.img+'.psf', self.img+'.weight'],imval=[(self.img+'.image',1.0,[256,256,0,0]),(self.img+'.weight',0.493,[256,256,0,0]) ] ) |
| 30-Jul-2019 15:10:09 | - self.checkfinal(report) |
| 30-Jul-2019 15:10:09 | - |
| 30-Jul-2019 15:10:09 | - @unittest.skip("Skip failing widefield test") |
| 30-Jul-2019 15:10:09 | - def test_widefield_mosaic_outlier(self): |
| 30-Jul-2019 15:10:09 | - """ [multifield] Test_widefield_mosaic_outlier : Mosaic with an outlier field """ |
| 30-Jul-2019 15:10:09 | - #### Need another dataset for this. |
| 30-Jul-2019 15:10:09 | - self.prepData("refim_mawproject.ms") |
| 30-Jul-2019 15:10:09 | - ## Outlier uses gridft |
| 30-Jul-2019 15:10:09 | - self.th.write_file(self.img+'.out.txt', 'imagename='+self.img+'1\nimsize=[80,80]\ncell=[8.0arcsec,8.0arcsec]\nphasecenter=J2000 19:58:42.0 +40.55.58.543\nmask=circle[[40pix,40pix],10pix]\ngridder=gridft') |
| 30-Jul-2019 15:10:09 | - ret = tclean(vis=self.msfile,imagename=self.img,imsize=512,cell='10.0arcsec',phasecenter="J2000 19:59:28.500 +40.44.01.50", |
| 30-Jul-2019 15:10:09 | - outlierfile=self.img+'.out.txt',niter=10,deconvolver='hogbom',gridder='mosaicft',interactive=0,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | - report=self.th.checkall(ret=ret, |
| 30-Jul-2019 15:10:09 | - iterdone=10, # outlier field has nothing in it :). |
| 30-Jul-2019 15:10:09 | - nmajordone=2, |
| 30-Jul-2019 15:10:09 | - imexist=[self.img+'.image', self.img+'1.image'], |
| 30-Jul-2019 15:10:09 | - imval=[(self.img+'.image',0.933,[256,256,0,0]), |
| 30-Jul-2019 15:10:09 | - (self.img+'1.image',0.0,[40,40,0,0])]) |
| 30-Jul-2019 15:10:09 | - self.checkfinal(report) |
| 30-Jul-2019 15:10:09 | - #### NOT WORKING as the model is being picked from tst1, so residual after 0.6 model and one major cycle is same as dirty and the output image is 1.6 instead of 0.93. |
| 30-Jul-2019 15:10:09 | - |
| 30-Jul-2019 15:10:09 | ############################################## |
| 30-Jul-2019 15:10:09 | ############################################## |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | @@ -3575,7 +3568,7 @@ class test_mosaic_cube(testref_base): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def test_cube_mosaic_cbFalse_mwFalse_twofield(self): |
| 30-Jul-2019 15:10:09 | self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 30-Jul-2019 15:10:09 | - phasecenter = '' |
| 30-Jul-2019 15:10:09 | + phasecenter = 'J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 30-Jul-2019 15:10:09 | field='0,1' |
| 30-Jul-2019 15:10:09 | tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.10354316235,[512,596,0,0]),(self.img+'.image.pbcor',0.981979846954,[512,596,0,1]),(self.img+'.image.pbcor', 0.895015060902,[512,596,0,2])]) |
| 30-Jul-2019 15:10:09 | @@ -3626,7 +3619,7 @@ class test_mosaic_cube(testref_base): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def test_cube_mosaic_cbFalse_mwTrue_twofield(self): |
| 30-Jul-2019 15:10:09 | self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 30-Jul-2019 15:10:09 | - phasecenter = '' |
| 30-Jul-2019 15:10:09 | + phasecenter = 'J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 30-Jul-2019 15:10:09 | field='0,1' |
| 30-Jul-2019 15:10:09 | tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.10354316235,[512,596,0,0]),(self.img+'.image.pbcor', 0.981979727745,[512,596,0,1]),(self.img+'.image.pbcor', 0.895014822483,[512,596,0,2])]) |
| 30-Jul-2019 15:10:09 | @@ -3678,7 +3671,7 @@ class test_mosaic_cube(testref_base): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def test_cube_awproject_cbFalse_mwFalse_twofield(self): |
| 30-Jul-2019 15:10:09 | self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 30-Jul-2019 15:10:09 | - phasecenter = '' |
| 30-Jul-2019 15:10:09 | + phasecenter = 'J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 30-Jul-2019 15:10:09 | field='0,1' |
| 30-Jul-2019 15:10:09 | cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 30-Jul-2019 15:10:09 | tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | @@ -3731,7 +3724,7 @@ class test_mosaic_cube(testref_base): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def test_cube_awproject_cbFalse_mwTrue_twofield(self): |
| 30-Jul-2019 15:10:09 | self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 30-Jul-2019 15:10:09 | - phasecenter = '' |
| 30-Jul-2019 15:10:09 | + phasecenter = 'J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 30-Jul-2019 15:10:09 | field='0,1' |
| 30-Jul-2019 15:10:09 | cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 30-Jul-2019 15:10:09 | tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | @@ -3785,7 +3778,7 @@ class test_mosaic_cube(testref_base): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def test_cube_mosaic_cbFalse_mwFalse_twofield_upTrue(self): |
| 30-Jul-2019 15:10:09 | self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 30-Jul-2019 15:10:09 | - phasecenter = '' |
| 30-Jul-2019 15:10:09 | + phasecenter = 'J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 30-Jul-2019 15:10:09 | field='0,1' |
| 30-Jul-2019 15:10:09 | tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.10354316235,[512,596,0,0]),(self.img+'.image.pbcor',0.981979846954,[512,596,0,1]),(self.img+'.image.pbcor', 0.895015060902,[512,596,0,2])]) |
| 30-Jul-2019 15:10:09 | @@ -3837,7 +3830,7 @@ class test_mosaic_cube(testref_base): |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | def test_cube_awproject_cbFalse_mwFalse_twofield_upTrue(self): |
| 30-Jul-2019 15:10:09 | self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 30-Jul-2019 15:10:09 | - phasecenter = '' |
| 30-Jul-2019 15:10:09 | + phasecenter = 'J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 30-Jul-2019 15:10:09 | field='0,1' |
| 30-Jul-2019 15:10:09 | cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 30-Jul-2019 15:10:09 | tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 30-Jul-2019 15:10:09 | Argument values: |
| 30-Jul-2019 15:10:09 | casatasks |
| 30-Jul-2019 15:10:09 | /export/home/cbt-el7-6/bamboohome/xml-data/build-dir/CASA-CTB-JOB1/casatasks |
| 30-Jul-2019 15:10:09 | master |
| 30-Jul-2019 15:10:09 | None |
| 30-Jul-2019 15:10:09 | Last Tag: 2019.126 |
| 30-Jul-2019 15:10:09 | New Tag: 2019.127 |
| 30-Jul-2019 15:10:09 | Creating a new tag |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | None |
| 30-Jul-2019 15:10:09 | Push enabled. Pushing new tag upstream. |
| 30-Jul-2019 15:10:09 | To ssh://git@open-bitbucket.nrao.edu:7999/casa/casatasks.git |
| 30-Jul-2019 15:10:09 | * [new tag] 2019.127 -> 2019.127 |
| 30-Jul-2019 15:10:09 | |
| 30-Jul-2019 15:10:09 | None |
| 30-Jul-2019 15:10:09 | Finished task 'Create wheel' with result: Success |
| 30-Jul-2019 15:10:09 | Running post build plugin 'Docker Container Cleanup' |
| 30-Jul-2019 15:10:09 | Beginning to execute external process for build 'CASA - Casa6 Tasks Wheel Build - Build and Package RHEL6 #124 (CASA-CTB-JOB1-124)' ... running command line: /bin/docker rm -f task-wheel-container ... in: /export/home/cbt-el7-6/bamboohome/xml-data/build-dir/CASA-CTB-JOB1 |
| 30-Jul-2019 15:10:11 | task-wheel-container |
| 30-Jul-2019 15:10:11 | Running post build plugin 'NCover Results Collector' |
| 30-Jul-2019 15:10:11 | Running post build plugin 'Clover Results Collector' |
| 30-Jul-2019 15:10:11 | Running post build plugin 'npm Cache Cleanup' |
| 30-Jul-2019 15:10:11 | Running post build plugin 'Artifact Copier' |
| 30-Jul-2019 15:10:11 | Publishing an artifact: EL6 wheel |
| 30-Jul-2019 15:10:11 | Finished publishing of artifact Non required shared artifact: [EL6 wheel], pattern: [*-any.whl] anchored at: [wheeldirectory/] in 264.7 ms |
| 30-Jul-2019 15:10:11 | Finalising the build... |
| 30-Jul-2019 15:10:11 | Stopping timer. |
| 30-Jul-2019 15:10:11 | Build CASA-CTB-JOB1-124 completed. |
| 30-Jul-2019 15:10:11 | Running on server: post build plugin 'NCover Results Collector' |
| 30-Jul-2019 15:10:11 | Running on server: post build plugin 'Build Hanging Detection Configuration' |
| 30-Jul-2019 15:10:11 | Running on server: post build plugin 'Clover Delta Calculator' |
| 30-Jul-2019 15:10:11 | Running on server: post build plugin 'Maven Dependencies Postprocessor' |
| 30-Jul-2019 15:10:11 | All post build plugins have finished |
| 30-Jul-2019 15:10:11 | Generating build results summary... |
| 30-Jul-2019 15:10:11 | Saving build results to disk... |
| 30-Jul-2019 15:10:11 | Logging substituted variables... |
| 30-Jul-2019 15:10:11 | Indexing build results... |
| 30-Jul-2019 15:10:11 | Finished building CASA-CTB-JOB1-124. |