Intermediate possibly buggy state
Updated generate_sparse_precision_matrix

Added an option to get an corresponding covariance matrix with
only ones on the diagonal.
Optimized Newton-Raphson step.

Dropped dependence on scipy for Newton-Raphson.
Code cleanup