use BLAS 
more complete helpers.h
Half dynamic slice?
fix irls?
iterators
