try closed boundary

This commit is contained in:
Max Lübke 2023-09-05 14:26:16 +02:00 committed by Max Lübke
parent 115dd87ba3
commit 32b2c5c741

View File

@ -40,7 +40,8 @@ inline double getFluxApplyConc(bool inbound, std::uint32_t curr_index,
// On inbound flux and non-boundary condition
if (inbound) {
if (neighbor_index == -1) {
return bc;
// HACK: assume closed boundaries for benchmark
return conc[curr_index];
}
return conc[neighbor_index];
}