diff --git a/dart/src/wTimeSolver.cpp b/dart/src/wTimeSolver.cpp
index fcd307833bb2c5580030b6e7fed1e706e58fe462..ae69b24551bd5ad07b5d0a2af09b21372bf420a4 100644
--- a/dart/src/wTimeSolver.cpp
+++ b/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
diff --git a/dart/src/wViscFluxes.cpp b/dart/src/wViscFluxes.cpp
index 86269781ba7ab35260ee7d3970705db84bb8f188..5d039e5559060e3564959445dae84c0f212fab16 100644
--- a/dart/src/wViscFluxes.cpp
+++ b/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")