Add dtype to matrix

* linsolver/solvers.py (gaussian_elimination): Specify dtype of xx

1 files changed, 1 insertions, 1 deletions

diff --git a/solvers.py b/solvers.py
index da38e2b..933bc2f 100644
--- a/solvers.py
+++ b/solvers.py
@@ -52,7 +52,7 @@ def gaussian_eliminate(
         return None
 
     # Back substitution
-    xx = np.zeros((nn, 1))
+    xx = np.zeros((nn, 1), dtype=np.float_)
     for i in range(nn - 1, -1, -1):
         xx[i] = (ee[i, nn] - np.dot(ee[i, i:nn], xx[i:nn])) / ee[i, i]