Optimizing Orthogonal Matching Pursuit code in Numpy, part 2

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 …

Optimizing Orthogonal Matching Pursuit code in Numpy, part 1

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 …