diff --git a/src/MVPSolver.cc b/src/MVPSolver.cc
index 45237e36..7feb564d 100644
--- a/src/MVPSolver.cc
+++ b/src/MVPSolver.cc
@@ -347,6 +347,7 @@ int MVPSolver<OrbitalsType, MatrixType>::solve(OrbitalsType& orbitals)
             = dynamic_cast<ProjectedMatrices<MatrixType>*>(
                 orbitals.getProjMatrices());
         projmatrices->setDM(*work_, orbitals.getIterativeIndex());
+        projmatrices->setEigenvalues(proj_mat_work_->getEigenvalues());
     }
 
     // Generate new density
diff --git a/src/ProjectedMatrices.h b/src/ProjectedMatrices.h
index 702637f4..6725b827 100644
--- a/src/ProjectedMatrices.h
+++ b/src/ProjectedMatrices.h
@@ -409,6 +409,12 @@ class ProjectedMatrices : public ProjectedMatricesInterface
     {
         dm_->setMatrix(mat, orbitals_index);
     }
+    void setEigenvalues(const std::vector<double>& eigenvalues)
+    {
+        memcpy(eigenvalues_.data(), eigenvalues.data(),
+            eigenvalues.size() * sizeof(double));
+    }
+    const std::vector<double>& getEigenvalues() const { return eigenvalues_; }
     void computeGenEigenInterval(std::vector<double>& interval,
         const int maxits, const double padding = 0.01);
     DensityMatrix<MatrixType>& getDM() { return *dm_; }
diff --git a/tests/MVP/test.py b/tests/MVP/test.py
index 708fd45c..bd5d708e 100644
--- a/tests/MVP/test.py
+++ b/tests/MVP/test.py
@@ -70,6 +70,30 @@
     print("Found HDF5 error")
     sys.exit(1)
 
+flag = 0
+eigenvalues=[]
+for line in lines:
+  if line.count(b'FERMI'):
+    flag = 0
+  if flag==1:
+    words=line.split()
+    for w in words:
+      eigenvalues.append(eval(w))
+  if line.count(b'Eigenvalues'):
+    flag = 1
+    eigenvalues=[]
+
+print(eigenvalues)
+tol = 1.e-4
+eigenvalue0 = -0.208
+if abs(eigenvalues[0]-eigenvalue0)>tol:
+  print("Expected eigenvalue 0 to be {}".format(eigenvalue0))
+  sys.exit(1)
+eigenvalue50 = 0.208
+if abs(eigenvalues[50]-eigenvalue50)>tol:
+  print("Expected eigenvalue 50 to be {}".format(eigenvalue50))
+  sys.exit(1)
+
 niterations = len(energies)
 print("MVP solver ran for {} iterations".format(niterations))
 if niterations>180: