max/tug
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/tug.git synced 2025-12-14 09:58:22 +01:00
Code Issues Packages Projects Releases Wiki Activity

fix: conversion of inner index when using inner boundary conditions

Browse Source
This commit is contained in:
Max Luebke 2023-01-24 12:01:56 +01:00
parent d4e3ab8544
commit 79d7a32fc2
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/BoundaryCondition.cpp
View File

@ -175,7 +175,7 @@ void tug::bc::BoundaryCondition::setInnerBC(boundary_condition bc, int x,
if (x >= this->sizes[X_DIM] || y >= this->sizes[Y_DIM]) { if (x >= this->sizes[X_DIM] || y >= this->sizes[Y_DIM]) {
throw_out_of_range("One input parameter is out of range"); throw_out_of_range("One input parameter is out of range");
} }
uint32_t index = x * this->sizes[Y_DIM] + y; uint32_t index = y * this->sizes[X_DIM] + x;
auto it = this->inner_cells.find(index); auto it = this->inner_cells.find(index);
if (it != this->inner_cells.end()) { if (it != this->inner_cells.end()) {
@ -187,7 +187,7 @@ void tug::bc::BoundaryCondition::setInnerBC(boundary_condition bc, int x,
} }
void tug::bc::BoundaryCondition::unsetInnerBC(int x, int y) { void tug::bc::BoundaryCondition::unsetInnerBC(int x, int y) {
uint32_t index = x * this->sizes[Y_DIM] + y; uint32_t index = y * this->sizes[X_DIM] + x;
this->inner_cells.erase(index); this->inner_cells.erase(index);
} }
@ -197,7 +197,7 @@ auto tug::bc::BoundaryCondition::getInnerBC(int x, int y = 0)
throw_out_of_range("One input parameter is out of range"); throw_out_of_range("One input parameter is out of range");
} }
uint32_t index = x * this->sizes[Y_DIM] + y; uint32_t index = y * this->sizes[X_DIM] + x;
auto it = this->inner_cells.find(index); auto it = this->inner_cells.find(index);