#encoding:latin-1
import nnls
print "module __doc__ says:"
print
print nnls.nnls.__doc__
print
print 70*"="
print 
print "_nnls.py __doc__ says:"
print 
print nnls._nnls.__doc__

from numpy import *
import time


A=random.rand(2000,800)-.6
x=arange(800)
b=dot(A,x)

print
print 70*"="
print
print "solve 2000 x 800 system, this will need some seconds"
print

start = time.time()
x, resid = nnls.nnls(A,b)

print 70*"="
print
print "residual is", resid
print
print "needed %.2f seconds" % (time.time()-start)
print