[Python] What is wrong with this NumPy linear algebra matrix rank calculator?

Consider:

import numpy as np
from numpy.linalg import matrix_rank

for i in range(1, 70):
y = np.linspace(0,i, i+1)
y = 2 ** y
x = np.ones(i+1)
matrix = np.stack([x,y])
print(matrix_rank(matrix))


After 48, I see the rank shifting from 2 to 1. How do I fix this?

Continue reading...