From 94d3eb9804f6dea3f6a4c3dc9cb652e0dda2f5fa Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?St=C3=A9phane=20Del=20Pino?= <stephane.delpino44@gmail.com>
Date: Fri, 22 Dec 2023 16:39:52 +0100
Subject: [PATCH] Compute geometric error in implicit cells only
---
src/scheme/ImplicitAcousticSolver.cpp | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/scheme/ImplicitAcousticSolver.cpp b/src/scheme/ImplicitAcousticSolver.cpp
index f823c1b8f..350f34d9b 100644
--- a/src/scheme/ImplicitAcousticSolver.cpp
+++ b/src/scheme/ImplicitAcousticSolver.cpp
@@ -1950,7 +1950,9 @@ class ImplicitAcousticSolverHandler::ImplicitAcousticSolver final
// const auto& cell_nodes = cell_to_node_matrix[j];
// V_j^n+1/tau_j^n+1 - M_j
// sum_tau_rho += std::abs(new_tau[j] - 1. / new_rho[j]);
- max_tau_error = std::max(max_tau_error, std::abs(1 - new_tau[j] * new_rho[j]));
+ if (m_is_implicit_cell[j]) {
+ max_tau_error = std::max(max_tau_error, std::abs(1 - new_tau[j] * new_rho[j]));
+ }
}
// std::cout << "sum_tau_rho =" << sum_tau_rho << '\n';
// std::cout << "max_tau_error=" << max_tau_error << '\n';
--
GitLab