>>> 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]])