Lacosmicx

I have finished the optimized LA Cosmic implementation. It uses Cython for the main part of the algorithm and uses C for the individual routines. Everything is multithreaded using OpenMP (the Cython standard). For a 4K x 2K image, this implementation is ~35 times faster than cosmics.py written by Malte Tewes (which is already significantly faster than the original IRAF version). 

Given the complexity of this code, should this be included in imageutils or should it be its own package? 

Any other comments are welcome.
Ah ... you have to add functions / classes you want to export to the end-user (in your case just one function that was called `run` I think) to the `__all__` list in `imageutils/__init__.py`:
https://github.com/astropy/imageutils/blob/master/imageutils/__init__.py#L18

Then `python setup.py build_sphinx` will pick it up and generate API docs because of the `.. automodapi:: imageutils` statement here:
https://github.com/astropy/imageutils/blame/master/docs/imageutils/index.rst#L55