def test_getkernel(self):
(ns, M2, M1) = self.fmap.shape
N2 = M2
fmap = np.swapaxes(self.fmap, -1, -2)
chi = np.swapaxes(self.chi, -1, -2)
b = np.arange(M2) - M2 / 2
n2 = np.arange(M2) - M2 / 2
Tl = 0.000742
K = GU.get_kernel(M2, Tl, b, n2, fmap, chi)
assert_array_almost_equal(K, self.ref_kernel)
def test_getkernel(self):
(ns, M2, M1) = self.fmap.shape
N2 = M2
fmap = np.swapaxes(self.fmap, -1, -2)
chi = np.swapaxes(self.chi, -1, -2)
b = np.arange(M2)-M2/2
n2 = np.arange(M2)-M2/2
Tl = 0.000742
K = GU.get_kernel(M2, Tl, b, n2, fmap, chi)
assert_array_almost_equal(K, self.ref_kernel)
