diff --git a/.gitignore b/.gitignore
index 64cb8dbc..fdcb032f 100644
--- a/.gitignore
+++ b/.gitignore
@@ -50,7 +50,7 @@ config.status
 .deps
 Make.inc
 eigen.inc
-Eigen.inc
+# Eigen.inc
 
 # http://www.gnu.org/software/autoconf #
 ########################################
diff --git a/lib/Eigen.inc b/lib/Eigen.inc
new file mode 100644
index 00000000..3b3328fe
--- /dev/null
+++ b/lib/Eigen.inc
@@ -0,0 +1,301 @@
+eigen_files =\
+  Eigen/Eigenvalues \
+  Eigen/Geometry \
+  Eigen/Eigen \
+  Eigen/StdList \
+  Eigen/StdVector \
+  Eigen/Householder \
+  Eigen/StdDeque \
+  Eigen/OrderingMethods \
+  Eigen/PardisoSupport \
+  Eigen/SparseQR \
+  Eigen/Dense \
+  Eigen/Cholesky \
+  Eigen/Jacobi \
+  Eigen/QtAlignedMalloc \
+  Eigen/SPQRSupport \
+  Eigen/UmfPackSupport \
+  Eigen/SparseCholesky \
+  Eigen/Core \
+  Eigen/MetisSupport \
+  Eigen/src/Eigenvalues/RealSchur.h \
+  Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h \
+  Eigen/src/Eigenvalues/GeneralizedSelfAdjointEigenSolver.h \
+  Eigen/src/Eigenvalues/EigenSolver.h \
+  Eigen/src/Eigenvalues/SelfAdjointEigenSolver_LAPACKE.h \
+  Eigen/src/Eigenvalues/Tridiagonalization.h \
+  Eigen/src/Eigenvalues/GeneralizedEigenSolver.h \
+  Eigen/src/Eigenvalues/RealQZ.h \
+  Eigen/src/Eigenvalues/MatrixBaseEigenvalues.h \
+  Eigen/src/Eigenvalues/ComplexEigenSolver.h \
+  Eigen/src/Eigenvalues/ComplexSchur_LAPACKE.h \
+  Eigen/src/Eigenvalues/RealSchur_LAPACKE.h \
+  Eigen/src/Eigenvalues/HessenbergDecomposition.h \
+  Eigen/src/Eigenvalues/ComplexSchur.h \
+  Eigen/src/Geometry/Rotation2D.h \
+  Eigen/src/Geometry/OrthoMethods.h \
+  Eigen/src/Geometry/Transform.h \
+  Eigen/src/Geometry/ParametrizedLine.h \
+  Eigen/src/Geometry/Scaling.h \
+  Eigen/src/Geometry/Homogeneous.h \
+  Eigen/src/Geometry/Umeyama.h \
+  Eigen/src/Geometry/Quaternion.h \
+  Eigen/src/Geometry/AngleAxis.h \
+  Eigen/src/Geometry/Hyperplane.h \
+  Eigen/src/Geometry/AlignedBox.h \
+  Eigen/src/Geometry/EulerAngles.h \
+  Eigen/src/Geometry/RotationBase.h \
+  Eigen/src/Geometry/arch/Geometry_SSE.h \
+  Eigen/src/Geometry/Translation.h \
+  Eigen/src/Householder/HouseholderSequence.h \
+  Eigen/src/Householder/Householder.h \
+  Eigen/src/Householder/BlockHouseholder.h \
+  Eigen/src/OrderingMethods/Amd.h \
+  Eigen/src/OrderingMethods/Ordering.h \
+  Eigen/src/OrderingMethods/Eigen_Colamd.h \
+  Eigen/src/PardisoSupport/PardisoSupport.h \
+  Eigen/src/SparseQR/SparseQR.h \
+  Eigen/src/Cholesky/LLT.h \
+  Eigen/src/Cholesky/LLT_LAPACKE.h \
+  Eigen/src/Cholesky/LDLT.h \
+  Eigen/src/Jacobi/Jacobi.h \
+  Eigen/src/SPQRSupport/SuiteSparseQRSupport.h \
+  Eigen/src/UmfPackSupport/UmfPackSupport.h \
+  Eigen/src/SparseCholesky/SimplicialCholesky_impl.h \
+  Eigen/src/SparseCholesky/SimplicialCholesky.h \
+  Eigen/src/Core/GlobalFunctions.h \
+  Eigen/src/Core/AssignEvaluator.h \
+  Eigen/src/Core/PermutationMatrix.h \
+  Eigen/src/Core/Fuzzy.h \
+  Eigen/src/Core/SelfCwiseBinaryOp.h \
+  Eigen/src/Core/Reverse.h \
+  Eigen/src/Core/Select.h \
+  Eigen/src/Core/VectorBlock.h \
+  Eigen/src/Core/Assign_MKL.h \
+  Eigen/src/Core/CoreIterators.h \
+  Eigen/src/Core/CwiseBinaryOp.h \
+  Eigen/src/Core/functors/TernaryFunctors.h \
+  Eigen/src/Core/functors/NullaryFunctors.h \
+  Eigen/src/Core/functors/UnaryFunctors.h \
+  Eigen/src/Core/functors/StlFunctors.h \
+  Eigen/src/Core/functors/BinaryFunctors.h \
+  Eigen/src/Core/functors/AssignmentFunctors.h \
+  Eigen/src/Core/SelfAdjointView.h \
+  Eigen/src/Core/CoreEvaluators.h \
+  Eigen/src/Core/Swap.h \
+  Eigen/src/Core/TriangularMatrix.h \
+  Eigen/src/Core/Ref.h \
+  Eigen/src/Core/CwiseNullaryOp.h \
+  Eigen/src/Core/PlainObjectBase.h \
+  Eigen/src/Core/Array.h \
+  Eigen/src/Core/Random.h \
+  Eigen/src/Core/Transpositions.h \
+  Eigen/src/Core/Block.h \
+  Eigen/src/Core/NestByValue.h \
+  Eigen/src/Core/ArrayWrapper.h \
+  Eigen/src/Core/CwiseUnaryOp.h \
+  Eigen/src/Core/MathFunctionsImpl.h \
+  Eigen/src/Core/CommaInitializer.h \
+  Eigen/src/Core/DenseStorage.h \
+  Eigen/src/Core/GeneralProduct.h \
+  Eigen/src/Core/MapBase.h \
+  Eigen/src/Core/Inverse.h \
+  Eigen/src/Core/Stride.h \
+  Eigen/src/Core/Visitor.h \
+  Eigen/src/Core/Solve.h \
+  Eigen/src/Core/BandMatrix.h \
+  Eigen/src/Core/util/Memory.h \
+  Eigen/src/Core/util/XprHelper.h \
+  Eigen/src/Core/util/ReenableStupidWarnings.h \
+  Eigen/src/Core/util/DisableStupidWarnings.h \
+  Eigen/src/Core/util/Meta.h \
+  Eigen/src/Core/util/MKL_support.h \
+  Eigen/src/Core/util/Constants.h \
+  Eigen/src/Core/util/StaticAssert.h \
+  Eigen/src/Core/util/NonMPL2.h \
+  Eigen/src/Core/util/ForwardDeclarations.h \
+  Eigen/src/Core/util/Macros.h \
+  Eigen/src/Core/util/BlasUtil.h \
+  Eigen/src/Core/Dot.h \
+  Eigen/src/Core/NoAlias.h \
+  Eigen/src/Core/StableNorm.h \
+  Eigen/src/Core/Transpose.h \
+  Eigen/src/Core/ReturnByValue.h \
+  Eigen/src/Core/SolverBase.h \
+  Eigen/src/Core/ConditionEstimator.h \
+  Eigen/src/Core/Product.h \
+  Eigen/src/Core/DenseCoeffsBase.h \
+  Eigen/src/Core/SolveTriangular.h \
+  Eigen/src/Core/Matrix.h \
+  Eigen/src/Core/GenericPacketMath.h \
+  Eigen/src/Core/MathFunctions.h \
+  Eigen/src/Core/Diagonal.h \
+  Eigen/src/Core/CwiseUnaryView.h \
+  Eigen/src/Core/ProductEvaluators.h \
+  Eigen/src/Core/DiagonalMatrix.h \
+  Eigen/src/Core/NumTraits.h \
+  Eigen/src/Core/Redux.h \
+  Eigen/src/Core/ArrayBase.h \
+  Eigen/src/Core/products/SelfadjointMatrixMatrix.h \
+  Eigen/src/Core/products/GeneralMatrixMatrix.h \
+  Eigen/src/Core/products/SelfadjointProduct.h \
+  Eigen/src/Core/products/SelfadjointMatrixMatrix_BLAS.h \
+  Eigen/src/Core/products/SelfadjointRank2Update.h \
+  Eigen/src/Core/products/SelfadjointMatrixVector_BLAS.h \
+  Eigen/src/Core/products/Parallelizer.h \
+  Eigen/src/Core/products/TriangularMatrixMatrix.h \
+  Eigen/src/Core/products/SelfadjointMatrixVector.h \
+  Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h \
+  Eigen/src/Core/products/TriangularSolverVector.h \
+  Eigen/src/Core/products/TriangularSolverMatrix.h \
+  Eigen/src/Core/products/TriangularMatrixMatrix_BLAS.h \
+  Eigen/src/Core/products/TriangularMatrixVector.h \
+  Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h \
+  Eigen/src/Core/products/TriangularSolverMatrix_BLAS.h \
+  Eigen/src/Core/products/TriangularMatrixVector_BLAS.h \
+  Eigen/src/Core/products/GeneralMatrixMatrix_BLAS.h \
+  Eigen/src/Core/products/GeneralMatrixVector_BLAS.h \
+  Eigen/src/Core/products/GeneralBlockPanelKernel.h \
+  Eigen/src/Core/products/GeneralMatrixVector.h \
+  Eigen/src/Core/ForceAlignedAccess.h \
+  Eigen/src/Core/EigenBase.h \
+  Eigen/src/Core/Map.h \
+  Eigen/src/Core/IO.h \
+  Eigen/src/Core/BooleanRedux.h \
+  Eigen/src/Core/Assign.h \
+  Eigen/src/Core/arch/CUDA/PacketMath.h \
+  Eigen/src/Core/arch/CUDA/Complex.h \
+  Eigen/src/Core/arch/CUDA/MathFunctions.h \
+  Eigen/src/Core/arch/CUDA/TypeCasting.h \
+  Eigen/src/Core/arch/CUDA/Half.h \
+  Eigen/src/Core/arch/CUDA/PacketMathHalf.h \
+  Eigen/src/Core/arch/AVX/PacketMath.h \
+  Eigen/src/Core/arch/AVX/Complex.h \
+  Eigen/src/Core/arch/AVX/MathFunctions.h \
+  Eigen/src/Core/arch/AVX/TypeCasting.h \
+  Eigen/src/Core/arch/NEON/PacketMath.h \
+  Eigen/src/Core/arch/NEON/Complex.h \
+  Eigen/src/Core/arch/NEON/MathFunctions.h \
+  Eigen/src/Core/arch/AltiVec/PacketMath.h \
+  Eigen/src/Core/arch/AltiVec/Complex.h \
+  Eigen/src/Core/arch/AltiVec/MathFunctions.h \
+  Eigen/src/Core/arch/AVX512/PacketMath.h \
+  Eigen/src/Core/arch/AVX512/MathFunctions.h \
+  Eigen/src/Core/arch/SSE/PacketMath.h \
+  Eigen/src/Core/arch/SSE/Complex.h \
+  Eigen/src/Core/arch/SSE/MathFunctions.h \
+  Eigen/src/Core/arch/SSE/TypeCasting.h \
+  Eigen/src/Core/arch/Default/Settings.h \
+  Eigen/src/Core/arch/ZVector/PacketMath.h \
+  Eigen/src/Core/arch/ZVector/Complex.h \
+  Eigen/src/Core/arch/ZVector/MathFunctions.h \
+  Eigen/src/Core/CwiseTernaryOp.h \
+  Eigen/src/Core/DenseBase.h \
+  Eigen/src/Core/MatrixBase.h \
+  Eigen/src/Core/Replicate.h \
+  Eigen/src/Core/DiagonalProduct.h \
+  Eigen/src/Core/VectorwiseOp.h \
+  Eigen/src/MetisSupport/MetisSupport.h \
+  Eigen/src/plugins/ArrayCwiseBinaryOps.h \
+  Eigen/src/plugins/MatrixCwiseUnaryOps.h \
+  Eigen/src/plugins/BlockMethods.h \
+  Eigen/src/plugins/CommonCwiseBinaryOps.h \
+  Eigen/src/plugins/MatrixCwiseBinaryOps.h \
+  Eigen/src/plugins/CommonCwiseUnaryOps.h \
+  Eigen/src/plugins/ArrayCwiseUnaryOps.h \
+  Eigen/src/CholmodSupport/CholmodSupport.h \
+  Eigen/src/SVD/UpperBidiagonalization.h \
+  Eigen/src/SVD/JacobiSVD_LAPACKE.h \
+  Eigen/src/SVD/BDCSVD.h \
+  Eigen/src/SVD/SVDBase.h \
+  Eigen/src/SVD/JacobiSVD.h \
+  Eigen/src/SparseLU/SparseLUImpl.h \
+  Eigen/src/SparseLU/SparseLU_column_bmod.h \
+  Eigen/src/SparseLU/SparseLU_relax_snode.h \
+  Eigen/src/SparseLU/SparseLU_pivotL.h \
+  Eigen/src/SparseLU/SparseLU_Memory.h \
+  Eigen/src/SparseLU/SparseLU_column_dfs.h \
+  Eigen/src/SparseLU/SparseLU_panel_dfs.h \
+  Eigen/src/SparseLU/SparseLU_kernel_bmod.h \
+  Eigen/src/SparseLU/SparseLU_copy_to_ucol.h \
+  Eigen/src/SparseLU/SparseLU_panel_bmod.h \
+  Eigen/src/SparseLU/SparseLU_pruneL.h \
+  Eigen/src/SparseLU/SparseLU_Utils.h \
+  Eigen/src/SparseLU/SparseLU.h \
+  Eigen/src/SparseLU/SparseLU_Structs.h \
+  Eigen/src/SparseLU/SparseLU_gemm_kernel.h \
+  Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h \
+  Eigen/src/SparseLU/SparseLU_heap_relax_snode.h \
+  Eigen/src/SparseCore/SparseColEtree.h \
+  Eigen/src/SparseCore/SparseDenseProduct.h \
+  Eigen/src/SparseCore/SparseAssign.h \
+  Eigen/src/SparseCore/SparseTriangularView.h \
+  Eigen/src/SparseCore/AmbiVector.h \
+  Eigen/src/SparseCore/SparseSelfAdjointView.h \
+  Eigen/src/SparseCore/SparseDiagonalProduct.h \
+  Eigen/src/SparseCore/SparseCompressedBase.h \
+  Eigen/src/SparseCore/SparseMatrixBase.h \
+  Eigen/src/SparseCore/MappedSparseMatrix.h \
+  Eigen/src/SparseCore/TriangularSolver.h \
+  Eigen/src/SparseCore/SparseCwiseUnaryOp.h \
+  Eigen/src/SparseCore/SparseProduct.h \
+  Eigen/src/SparseCore/SparseSparseProductWithPruning.h \
+  Eigen/src/SparseCore/SparseCwiseBinaryOp.h \
+  Eigen/src/SparseCore/SparseRef.h \
+  Eigen/src/SparseCore/SparsePermutation.h \
+  Eigen/src/SparseCore/SparseMap.h \
+  Eigen/src/SparseCore/SparseRedux.h \
+  Eigen/src/SparseCore/SparseSolverBase.h \
+  Eigen/src/SparseCore/SparseView.h \
+  Eigen/src/SparseCore/CompressedStorage.h \
+  Eigen/src/SparseCore/SparseFuzzy.h \
+  Eigen/src/SparseCore/SparseMatrix.h \
+  Eigen/src/SparseCore/SparseBlock.h \
+  Eigen/src/SparseCore/SparseUtil.h \
+  Eigen/src/SparseCore/SparseVector.h \
+  Eigen/src/SparseCore/SparseTranspose.h \
+  Eigen/src/SparseCore/ConservativeSparseSparseProduct.h \
+  Eigen/src/SparseCore/SparseDot.h \
+  Eigen/src/PaStiXSupport/PaStiXSupport.h \
+  Eigen/src/StlSupport/StdList.h \
+  Eigen/src/StlSupport/details.h \
+  Eigen/src/StlSupport/StdVector.h \
+  Eigen/src/StlSupport/StdDeque.h \
+  Eigen/src/LU/InverseImpl.h \
+  Eigen/src/LU/FullPivLU.h \
+  Eigen/src/LU/PartialPivLU_LAPACKE.h \
+  Eigen/src/LU/PartialPivLU.h \
+  Eigen/src/LU/arch/Inverse_SSE.h \
+  Eigen/src/LU/Determinant.h \
+  Eigen/src/SuperLUSupport/SuperLUSupport.h \
+  Eigen/src/QR/FullPivHouseholderQR.h \
+  Eigen/src/QR/HouseholderQR.h \
+  Eigen/src/QR/ColPivHouseholderQR.h \
+  Eigen/src/QR/HouseholderQR_LAPACKE.h \
+  Eigen/src/QR/ColPivHouseholderQR_LAPACKE.h \
+  Eigen/src/QR/CompleteOrthogonalDecomposition.h \
+  Eigen/src/misc/lapacke_mangling.h \
+  Eigen/src/misc/Kernel.h \
+  Eigen/src/misc/Image.h \
+  Eigen/src/misc/lapacke.h \
+  Eigen/src/misc/lapack.h \
+  Eigen/src/misc/RealSvd2x2.h \
+  Eigen/src/misc/blas.h \
+  Eigen/src/IterativeLinearSolvers/BiCGSTAB.h \
+  Eigen/src/IterativeLinearSolvers/SolveWithGuess.h \
+  Eigen/src/IterativeLinearSolvers/BasicPreconditioners.h \
+  Eigen/src/IterativeLinearSolvers/IncompleteCholesky.h \
+  Eigen/src/IterativeLinearSolvers/LeastSquareConjugateGradient.h \
+  Eigen/src/IterativeLinearSolvers/IterativeSolverBase.h \
+  Eigen/src/IterativeLinearSolvers/ConjugateGradient.h \
+  Eigen/src/IterativeLinearSolvers/IncompleteLUT.h \
+  Eigen/CholmodSupport \
+  Eigen/SVD \
+  Eigen/SparseLU \
+  Eigen/SparseCore \
+  Eigen/Sparse \
+  Eigen/PaStiXSupport \
+  Eigen/LU \
+  Eigen/SuperLUSupport \
+  Eigen/QR \
+  Eigen/IterativeLinearSolvers