diff --git a/Utils/Davidson.ml b/Utils/Davidson.ml index 1450fcf..d0e078f 100644 --- a/Utils/Davidson.ml +++ b/Utils/Davidson.ml @@ -4,7 +4,7 @@ type t let make ?guess - ?(n_states=8) + ?(n_states=1) ?(n_iter=10) ?(threshold=1.e-6) diagonal @@ -37,7 +37,7 @@ let make ) |> Util.normalize in - List.init m (fun i -> random_vector i) + List.init m (fun i -> random_vector (i+1)) in let pick_new u = @@ -73,7 +73,7 @@ let make matrix_prod ( u_new_ortho |> Mat.of_col_vecs_list - |> Matrix.sparse_of_mat ) + |> Matrix.dense_of_mat ) |> Matrix.to_mat |> Mat.to_col_vecs_list in