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)