precision is being lost in computations of P and P^(-1); crs_matrix values member storage plays a role...