Category: scikit-learn
#bug #omp #orthogonal matching pursuit #dictionary learning #scikit-learn
A while back, Bob L. Sturm blogged about a similar implementation of OMP to the one in scikit-learn. Instead of using the Cholesky decomposition like we did, his Matlab code uses the QR decomposition, to a similar (or maybe even identical) outcome, in theory. So lucky that Alejandro pointed out …
Category: scikit-learn
#blas #efficient #lapack #numpy #omp #orthogonal matching pursuit #potrs #scipy #dictionary learning #python #scikit-learn
EDIT: There was a bug in the final version of the code presented here. It is fixed now, for its backstory, check out my blog post on it.
When we last saw our hero, he was fighting with the dreaded implementation of least-angle regression, knowing full well that it was …
Category: scikit-learn
#efficient #numpy #omp #orthogonal matching pursuit #scipy #dictionary learning #scikit-learn #python
After intense code optimization work, my implementation of OMP finally beat least-angle regression! This was the primary issue discussed during the pull request, so once performance was taken care of, the code was ready for merge. Orthogonal matching pursuit is now available in scikits.learn as a sparse linear regression …