From 87efacb3690ae72b32572dc1a82aaaf8b2e8b954 Mon Sep 17 00:00:00 2001
From: Paul Dechamps <paulzer@Pauls-MacBook-Pro.local>
Date: Tue, 19 Apr 2022 12:12:52 +0200
Subject: [PATCH] Correction Eigen BL fluxes
---
dart/src/wTimeSolver.cpp | 4 ++--
dart/src/wViscFluxes.cpp | 6 +++---
2 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/dart/src/wTimeSolver.cpp b/dart/src/wTimeSolver.cpp
index fcd3078..ae69b24 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 8626978..5d039e5 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")
--
GitLab