def test_raises(self) :
Data = self.blocks[0]
rotate_pol.rotate(Data, (1,))
self.assertRaises(ce.DataError, pol_cal.calibrate_pol, Data,self.mueler)
Data = self.blocks[1]
rotate_pol.rotate(Data, (1,2,3,4))
Data.field['CRVAL4'][0] = 3
Data.field['CRVAL4'][2] = 1
self.assertRaises(ce.DataError, pol_cal.calibrate_pol, Data,
self.mueler)
def setUp(self) :
Reader = fitsGBT.Reader(test_file, feedback=0)
blocks = Reader.read((),())
for Data in blocks :
rotate_pol.rotate(Data, (1,2,3,4))
Writer = fitsGBT.Writer(blocks, feedback=0)
Writer.write('test_rot.testout.fits')
self.params = {'pc_file_middles' : ("test",),
'pc_input_end' : "_rot.testout.fits",
'pc_output_end' : "_polcal.testout.fits",
'pc_mueler_file' : test_mueler_file
}
def setUp(self):
Reader = fitsGBT.Reader(test_file, feedback=0)
self.blocks = Reader.read((), ())
b = self.blocks[0]
b.calc_freq()
map = sp.zeros((150, 10, 11))
map = algebra.make_vect(map, axis_names=('freq', 'ra', 'dec'))
map.set_axis_info('freq', 700.0e6, -2.0e6)
map.set_axis_info('ra', 325.6, -0.2)
map.set_axis_info('dec', 0, 0.2)
self.map = map
for Data in self.blocks:
rotate_pol.rotate(Data, (1, 2, 3, 4))
def setUp(self) :
Reader = fitsGBT.Reader(test_file, feedback=0)
self.blocks = Reader.read((),())
b = self.blocks[0]
b.calc_freq()
map = sp.zeros((150, 10, 11))
map = algebra.make_vect(map, axis_names=('freq', 'ra', 'dec'))
map.set_axis_info('freq', 700.0e6, -2.0e6)
map.set_axis_info('ra', 325.6, -0.2)
map.set_axis_info('dec', 0, 0.2)
self.map = map
for Data in self.blocks :
rotate_pol.rotate(Data, (1, 2, 3, 4))
def setUp(self) :
Reader = fitsGBT.Reader(test_file, feedback=0)
blocks = Reader.read((),())
for Data in blocks :
rotate_pol.rotate(Data, (1,2,3,4))
Writer = fitsGBT.Writer(blocks, feedback=0)
Writer.write('test_rot.testout.fits')
map = sp.zeros((150, 10, 11))
map = algebra.make_vect(map, axis_names=('freq', 'ra', 'dec'))
map.set_axis_info('freq', 700.0e6, -2.0e6)
map.set_axis_info('ra', 325.6, -0.2)
map.set_axis_info('dec', 0, 0.2)
algebra.save('test_inmap.npy', map)
self.params = {'sm_file_middles' : ("test",),
'sm_input_end' : "_rot.testout.fits",
'sm_output_end' : "_sub.testout.fits",
'sm_solve_for_gain' : True,
'sm_gain_output_end' : '_gain.pickle',
'sm_map_file' : 'test_inmap.npy'
}
def test_runs(self) :
for Data in self.blocks :
rotate_pol.rotate(Data, (1,2,3,4))
pol_cal.calibrate_pol(Data, self.mueler)
def setUp(self) :
Reader = fitsGBT.Reader(test_file, feedback=0)
self.Data = Reader.read(1,0)
rotate_pol.rotate(self.Data, (1,))
self.DataSub = Reader.read(1,0)
rotate_pol.rotate(self.DataSub, (1,))
def setUp(self):
Reader = fitsGBT.Reader(test_file, feedback=0)
self.Data = Reader.read(1, 0)
rotate_pol.rotate(self.Data, (1, ))
self.DataSub = Reader.read(1, 0)
rotate_pol.rotate(self.DataSub, (1, ))