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

Fix boundary condition memory allocation

Browse Source
This commit is contained in:
Max Luebke 2022-01-27 09:54:13 +01:00
parent 6f0fd990da
commit 0828da7a09
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/BTCSDiffusion.cpp
View File

@ -46,10 +46,10 @@ void BTCSDiffusion::updateInternals() {
switch (grid_dim) {
case 1:
bc.resize(2, {BTCSDiffusion::BC_CLOSED, 0});
bc.resize(grid_cells[0] + 2, {BTCSDiffusion::BC_CLOSED, 0});
break;
case 2:
bc.resize(2 * grid_cells[0] + 2 * grid_cells[1],
bc.resize((grid_cells[0] + 2) * (grid_cells[1] + 2),
{BTCSDiffusion::BC_CLOSED, 0});
break;
case 3:
@ -163,7 +163,8 @@ void BTCSDiffusion::simulate(std::vector<double> &c,
// double bc_right =
// getBCFromTuple(1, c[c.size() - 1], alpha[alpha.size() - 1]);
simulate1D(c, bc[0], bc[1], alpha, this->deltas[0], this->grid_cells[0]);
simulate1D(c, bc[0], bc[grid_cells[0] + 1], alpha, this->deltas[0],
this->grid_cells[0]);
}
}