diff options
Diffstat (limited to 'src/box_obstacle.cc')
-rw-r--r-- | src/box_obstacle.cc | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/src/box_obstacle.cc b/src/box_obstacle.cc new file mode 100644 index 0000000..9cb6ad0 --- /dev/null +++ b/src/box_obstacle.cc @@ -0,0 +1,28 @@ +#include "box_obstacle.h" + +#include "boundary_conditions.h" + +BoxObstacle::BoxObstacle(std::size_t lX, std::size_t lY, std::size_t uX, std::size_t uY): + lower_x_(lX), lower_y_(lY), upper_x_(uX), upper_y_(uY) { } + +bool BoxObstacle::isInside(std::size_t x, std::size_t y) const { + return x > lower_x_ + && x < upper_x_ + && y > lower_y_ + && y < upper_y_; +} + +void BoxObstacle::applyBoundary(DataCellBuffer& pop) const { + for ( std::size_t x = lower_x_+1; x < upper_x_; ++x ) { + computeWallCell(pop, {x, lower_y_}, { 0,-1}); + computeWallCell(pop, {x, upper_y_}, { 0, 1}); + } + for ( std::size_t y = lower_y_+1; y < upper_y_; ++y ) { + computeWallCell(pop, {lower_x_, y}, {-1, 0}); + computeWallCell(pop, {upper_x_, y}, { 1, 0}); + } + computeWallCell(pop, {lower_x_, lower_y_}, {-1,-1}); + computeWallCell(pop, {upper_x_, lower_y_}, { 1,-1}); + computeWallCell(pop, {upper_x_, upper_y_}, { 1, 1}); + computeWallCell(pop, {lower_x_, upper_y_}, {-1, 1}); +} |