scipy.optimize.nnls

scipy.optimize.nnls(A, b)

Solve || Ax - b ||_2 -> min with x>=0

Inputs:

A – matrix as above b – vector as above

Outputs:

x – solution vector rnorm – residual || Ax-b ||_2

wrapper around NNLS.F code below nnls/ directory