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)