Removed overwriting of RHS which may affect solution

jeffrey-hokanson · Nov 12, 2020 · c149ced · c149ced
1 parent 6e5581b
commit c149ced
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/polyrat/vecfit.py b/polyrat/vecfit.py
@@ -14,7 +14,7 @@
 def _solve_linearized_vecfit(num_basis, denom_basis, y):
 	A = np.hstack([ num_basis, -(denom_basis[:,1:].T * y).T ])
 	b = y
-	x, res, rank, s = scipy.linalg.lstsq(A, b, overwrite_a = True, overwrite_b = True)
+	x, res, rank, s = scipy.linalg.lstsq(A, b, overwrite_a = True, overwrite_b = False)
 	a = x[0:num_basis.shape[1]]
 	b = np.hstack([1, x[num_basis.shape[1]:]])
 	return a, b, s[0]/s[-1]