def test_primary_beam():
def primary_beam(theta, phi):
import numpy as np
with settingerr(invalid='ignore'):
return theta <= np.radians(9)
a = QubicAcquisition(150, s, primary_beam=primary_beam)
p = a.get_projection_operator()
assert_equal(p.matrix.ncolmax, 5)
def test_primary_beam():
def primary_beam(theta, phi):
import numpy as np
with settingerr(invalid='ignore'):
return theta <= np.radians(9)
a = QubicAcquisition(150, s, primary_beam=primary_beam)
p = a.get_projection_operator()
assert_equal(p.matrix.ncolmax, 5)
def func(instrument, expected):
op = instrument.get_polarizer_operator(sampling, scene)
sky = np.ones((len(instrument), 1))
assert_equal(op(sky).ravel(), expected)
def func(instrument, expected):
op = instrument.get_polarizer_operator(sampling, scene)
sky = np.ones((len(instrument), 1))
assert_equal(op(sky).ravel(), expected)
