max/tug
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/tug.git synced 2025-12-13 17:38:23 +01:00
Code Issues Packages Projects Releases Wiki Activity

Refactor loop breakout to same syntax as in fillMatrix

Browse Source
This commit is contained in:
Max Luebke 2022-02-11 17:38:39 +01:00
parent 6ea3bf0182
commit 08dba0975b
1 changed files with 13 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
src/BTCSDiffusion.cpp
View File

@ -269,22 +269,22 @@ void BTCSDiffusion::fillVectorFromRow2D(Eigen::Map<DMatrixRowMajor> &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;
}
}