fix standardize_col
clip calc_mu to avoid impossible mu
Better mean calculation
Streamlined gradient descent
identify sufficient NLL progress
Alternate gradient descent method