From 08dba0975b7eefcad79c0145ddb5bdeaa77776b0 Mon Sep 17 00:00:00 2001 From: Max Luebke Date: Fri, 11 Feb 2022 17:38:39 +0100 Subject: [PATCH] Refactor loop breakout to same syntax as in fillMatrix --- src/BTCSDiffusion.cpp | 26 +++++++++++++------------- 1 file changed, 13 insertions(+), 13 deletions(-) diff --git a/src/BTCSDiffusion.cpp b/src/BTCSDiffusion.cpp index 507ce21..dc68093 100644 --- a/src/BTCSDiffusion.cpp +++ b/src/BTCSDiffusion.cpp @@ -269,22 +269,22 @@ void BTCSDiffusion::fillVectorFromRow2D(Eigen::Map &c, if (tmp_bc.type == BTCSDiffusion::BC_CONSTANT) { b_vector[offset * row + j] = tmp_bc.value; - } else { + continue; + } - double y_values[3]; - y_values[0] = - (row != 0 ? c(row - 1, j - 1) - : getBCFromFlux(tmp_bc, c(row, j - 1), alpha[j - 1])); - y_values[1] = c(row, j - 1); - y_values[2] = (row != nrow - 1 - ? c(row + 1, j - 1) + double y_values[3]; + y_values[0] = + (row != 0 ? c(row - 1, j - 1) + : getBCFromFlux(tmp_bc, c(row, j - 1), alpha[j - 1])); + y_values[1] = c(row, j - 1); + y_values[2] = + (row != nrow - 1 ? c(row + 1, j - 1) : getBCFromFlux(tmp_bc, c(row, j - 1), alpha[j - 1])); - double t0_c = - alpha[j - 1] * - ((y_values[0] - 2 * y_values[1] + y_values[2]) / (delta * delta)); - b_vector[offset * row + j] = -c(row, j - 1) - t0_c; - } + double t0_c = + alpha[j - 1] * + ((y_values[0] - 2 * y_values[1] + y_values[2]) / (delta * delta)); + b_vector[offset * row + j] = -c(row, j - 1) - t0_c; } }