Added 'entire_loss' style loss functions, which take an entire clustering. These are of course much less efficient than Stats based losses.

Added support for using function handles as loss functions, so custom losses can be passed to the optimizer. These get called with the entire clustering. Octave only for now.