irls_gaussian_cpp_fast function

Fast Iterative Reweighed Least Square algorithm for Gaussians