Hello! I have supported constrained solvers for linear matrix problems for about 10 years in C++, but have now switched to Python. I am going to submit a couple of new routines for linalg called autoreg(A,b) and autoregnn(A,b). They work just like lstsq(A,b) normally, but when they detect that the problem is dominated by noise they revert to an automatic regularization scheme that returns a better behaved result than one gets from lstsq. In addition, autoregnn enforces a nonnegativity constraint on the solution. I have put on my web site a slightly fuller featured version of these same two algorithms, using a Class implementation to facilitate retuning several diagnostic or other artifacts. The web site contains tutorials on these methods and a number of examples of their use. See http://www.rejones7.net/autorej/ . I hope this community can take a look at these routines and see whether they are appropriate for linalg or should be in another location. Ron Jones
Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10
_______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@python.org https://mail.python.org/mailman/listinfo/numpy-discussion
