def test(): map_nl = mapper_nl(tod, model, hypers=2*[1.], norms=[norm2_ellipsoid(invntt)] + 2*[norm2], tol=1.e-4, maxiter=1000, callback=None if tamasis.var.verbose else Callback(), ) print 'Elapsed time: ' + str(map_nl.header['TIME']) + ' after ' + \ str(map_nl.header['NITER']) + ' iterations.' if map_nl.header['NITER'] > 150: raise TestFailure()
masking_tod = Masking(tod.mask) model = masking_tod * projection # iterative map, taking all map pixels class Callback(): def __init__(self): self.niterations = 0 def __call__(self, x): self.niterations += 1 invntt = Diagonal(1 / obs.get_detector_stddev(100)**2) invntt = Identity() map_nl = mapper_nl( tod, model, hypers=2 * [1.], norms=[norm2_ellipsoid(invntt)] + 2 * [norm2], tol=1.e-4, maxiter=1000, callback=None if tamasis.var.verbose else Callback(), ) print 'Elapsed time: ' + str(map_nl.header['TIME']) + ' after ' + \ str(map_nl.header['NITER']) + ' iterations.' #if map_nl.header['NITER'] > 48: # raise TestFailure()
obs = PacsObservation(filename=data_dir+'frames_blue.fits', fine_sampling_factor=1) obs.pointing.chop[:] = 0 tod = obs.get_tod(flatfielding=False) projection = Projection(obs, resolution=3.2, oversampling=False, npixels_per_sample=6) masking_tod = Masking(tod.mask) model = masking_tod * projection # iterative map, taking all map pixels class Callback(): def __init__(self): self.niterations = 0 def __call__(self, x): self.niterations += 1 invntt = Diagonal(1/obs.get_detector_stddev(100)**2) invntt = Identity() map_nl = mapper_nl(tod, model, hypers=2*[1.], norms=[norm2_ellipsoid(invntt)] + 2*[norm2], tol=1.e-4, maxiter=1000, callback=None if tamasis.var.verbose else Callback(), ) print 'Elapsed time: ' + str(map_nl.header['TIME']) + ' after ' + \ str(map_nl.header['NITER']) + ' iterations.' #if map_nl.header['NITER'] > 48: # raise TestFailure()