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

Change loop index variable

Browse Source
This commit is contained in:
Max Luebke 2022-02-17 16:43:04 +01:00
parent b3e0f7987b
commit 0dad4f93fd
1 changed files with 9 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/BTCSDiffusion.cpp
View File

@ -257,27 +257,26 @@ void BTCSDiffusion::fillVectorFromRowADI(Eigen::Map<DMatrixRowMajor> &c,
getBCFromFlux(right, c(row, ncol - 1), alpha[ncol - 1]);
}
for (int j = 1; j < offset - 1; j++) {
boundary_condition tmp_bc = this->bc(row, j-1);
for (int j = 0; j < ncol; j++) {
boundary_condition tmp_bc = this->bc(row, j);
if (tmp_bc.type == BTCSDiffusion::BC_CONSTANT) {
b_vector[offset * row + j] = tmp_bc.value;
b_vector[offset * row + (j+1)] = tmp_bc.value;
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);
(row != 0 ? c(row - 1, j) : getBCFromFlux(tmp_bc, c(row, j), alpha[j]));
y_values[1] = c(row, j);
y_values[2] =
(row != nrow - 1 ? c(row + 1, j - 1)
: getBCFromFlux(tmp_bc, c(row, j - 1), alpha[j - 1]));
(row != nrow - 1 ? c(row + 1, j)
: getBCFromFlux(tmp_bc, c(row, j), alpha[j]));
double t0_c =
alpha[j - 1] *
alpha[j] *
((y_values[0] - 2 * y_values[1] + y_values[2]) / (delta * delta));
b_vector[offset * row + j] = -c(row, j - 1) - t0_c;
b_vector[offset * row + (j + 1)] = -c(row, j) - t0_c;
}
}