IGD for elastic net is working now

Is much better now,

but still needs a threshold in __gaussian_igd_result() to remove all
tiny coefficients, and calculate the corresponding loss and log-
likelihood. Standard deviations for all features will be needed to
compare the relative magnitudes of all coefficients.