Rozwiązywanie układu liniowego przez liczby całkowite z numpy
Próbuję rozwiązać naddeterminowany układ liniowy równań z numpy. Obecnie robię coś takiego (jako prosty przykład):
a = np.array([[1,0], [0,1], [-1,1]])
b = np.array([1,1,0])
print np.linalg.lstsq(a,b)[0]
[ 1. 1.]
To działa, ale używa pływaków. Czy jest jakiś sposób na rozwiązanie systemu tylko przez liczby całkowite? Próbowałem czegoś na wzór
print map(int, np.linalg.lstsq(a,b)[0])
[0, 1]
w celu przekonwertowania rozwiązania na tablicę intów[1, 1]
, ale wyraźnie coś mi brakuje. Czy ktoś mógłby wskazać mi właściwy kierunek?