Added a C function for calculating a difference matrix between genotypes.

This function is mostly complete, but needs optimization, documentation, and an R wrapper function.
This commit also contains a cleaned up version of mlg_clustering.c, but no functional changes.