>>> sample = np.matrix(np.ones((50,50)) * np.inf)
>>> sample
matrix([[ inf, inf, inf, ..., inf, inf, inf],
[ inf, inf, inf, ..., inf, inf, inf],
[ inf, inf, inf, ..., inf, inf, inf],
...,
[ inf, inf, inf, ..., inf, inf, inf],
[ inf, inf, inf, ..., inf, inf, inf],
[ inf, inf, inf, ..., inf, inf, inf]])