Merge pull request #62 from kamesy/kmsy/perf

Fix initial dual edge case
jondeuce · Apr 15, 2024 · e9d8eb2 · e9d8eb2
2 parents 266312d + 30c6cd6
commit e9d8eb2
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/lsqnonneg.jl b/src/lsqnonneg.jl
@@ -67,6 +67,7 @@ function solve!(
  w[j] = wj
  end
  @inbounds w[end] = 0
+ @inbounds w[end] = all(<=(0), w)
 
  # initialize nnls workspace
  @inbounds for i in 1:m
@@ -130,6 +131,7 @@ function solve!(
  w[j] = wj
  end
  @inbounds w[end] = 0
+ @inbounds w[end] = all(<=(0), w)
 
  # initialize nnls workspace
  @inbounds for i in 1:m