avoid another float coeff

shashi · shashi · commit 614dfcd1c0da · 2020-08-31T15:37:50.000-04:00
diff --git a/src/solve.jl b/src/solve.jl
@@ -93,7 +93,7 @@ function ldiv(A::LU, b)
         sub = simplifying_dot(b[i+1:end], U[i,i+1:end])
         den = U[i,i]
         x[i] = _iszero(sub) ? x[i] : x[i] - sub
-        x[i] = _isone(U[i,i]) ? x[i] : x[i] / U[i,i]
+        x[i] = _isone(U[i,i]) ? x[i] : _isone(-U[i,i]) ? -x[i] : x[i] / U[i,i]
     end
     x
 end

Original file line number	Diff line number	Diff line change
`@@ -93,7 +93,7 @@ function ldiv(A::LU, b)`
`93`	`93`	`sub = simplifying_dot(b[i+1:end], U[i,i+1:end])`
`94`	`94`	`den = U[i,i]`
`95`	`95`	`x[i] = _iszero(sub) ? x[i] : x[i] - sub`
`96`		`- x[i] = _isone(U[i,i]) ? x[i] : x[i] / U[i,i]`
	`96`	`+ x[i] = _isone(U[i,i]) ? x[i] : _isone(-U[i,i]) ? -x[i] : x[i] / U[i,i]`
`97`	`97`	`end`
`98`	`98`	`x`
`99`	`99`	`end`