kernel_matrix() and fill_lisf() funcions instead of the many for loops
they previously had.  In particular, using fill_lisf() makes the reducers
generally a lot faster and a little more accurate.

--HG--
extra : convert_revision : svn%3Afdd8eb12-d10e-0410-9acb-85c331704f74/trunk%403680