def test_min_at(): k = np.array([3, 2, 2, 2, 1, 1, 1]) values = np.array([10., 1, 2, 3, 31, 21, 11]) i = np.zeros(k.shape, int) j = np.zeros(k.shape, int) a = np.zeros([1, 1, 4]) + 100. _min_at(a, (i, j, k), values) npt.assert_array_equal(a, [[[100, 11, 1, 10]]])