ENH pairwise L1 distances for sparse matrices

Simple, O(n_features) temp space algorithm: densify row by row, then subtract
and compute L1 norm.

Added BLAS support code (cblas_dasum) to speed this up by a factor of two on
x86-64 w/ GCC and ATLAS for pair of 93% sparse matrices of shape 1000*3000.
That's an order of magnitude faster than the dense version.

Also cleaned up chi2 kernel code while I was at it and added a nogil decl.