Correction Eigen BL fluxes

dart/src/wTimeSolver.cpp

@@ -92,8 +92,8 @@ int TimeSolver::Integration(size_t iPoint, BLRegion *bl)
     double normSysRes0 = SysRes.norm();
     double normSysRes = normSysRes0;
 
-    MatrixXd JacMatrix(5,5);
-    VectorXd slnIncr(5);
+    MatrixXd JacMatrix = MatrixXd::Zero(5, 5);
+    VectorXd slnIncr = VectorXd::Zero(5);
 
     nErrors = 0;                 // Number of errors encountered
     unsigned int innerIters = 0; // Inner (non-linear) iterations

dart/src/wViscFluxes.cpp

@@ -38,7 +38,7 @@ MatrixXd ViscFluxes::ComputeJacobian(size_t iPoint, BLRegion *bl, VectorXd &sysR
 {
 
     unsigned int nVar = bl->GetnVar();
-    Matrix<double, 5, 5> JacMatrix;
+    MatrixXd JacMatrix = MatrixXd::Zero(5,5);
     std::vector<double> uUp(nVar, 0.);
     for (size_t k = 0; k < nVar; ++k)
     {
@@ -62,8 +62,8 @@ VectorXd ViscFluxes::BLlaws(size_t iPoint, BLRegion *bl, std::vector<double> u)
 {
     unsigned int nVar = bl->GetnVar();
 
-    Matrix<double, 5, 5> timeMatrix;
-    Vector<double, 5> spaceVector;
+    MatrixXd timeMatrix = MatrixXd::Zero(5, 5);
+    VectorXd spaceVector = VectorXd::Zero(5);
 
     double dissipRatio;
     if (bl->name == "wake")