'Test the bmat atom.'

def test_bmat(self):

    v_np = np.ones((3, 1))
    expr = bmat([[v_np, v_np], [[0, 0], [1, 2]]])
    self.assertEqual(expr.size, (5, 2))
    const = np.bmat([[v_np, v_np], [np.zeros((2, 1)), np.mat([1, 2]).T]])
    self.assertItemsAlmostEqual(expr.value, const)
