max/tug
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/tug.git synced 2025-12-15 18:38:23 +01:00
Code Issues Packages Projects Releases Wiki Activity

Fix alpha intercell calculation

Browse Source
This commit is contained in:
Max Lübke 2024-03-05 10:26:15 +01:00
parent 48000710c7
commit 2f737ce09e
1 changed files with 3 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
include/tug/Core/BTCS.hpp
View File

@ -138,10 +138,9 @@ template <typename T>
constexpr T calcExplicitConcentrationsBoundaryConstant(T conc_center, T conc_bc, constexpr T calcExplicitConcentrationsBoundaryConstant(T conc_center, T conc_bc,
T alpha_center, T alpha_center,
T alpha_neighbor, T sy) { T alpha_neighbor, T sy) {
return sy * calcAlphaIntercell(alpha_center, alpha_neighbor) * conc_center + const T inter_cell = calcAlphaIntercell(alpha_center, alpha_neighbor);
(1 - sy * (calcAlphaIntercell(alpha_center, alpha_center) + return sy * inter_cell * conc_center +
alpha_center)) * (1 - sy * (inter_cell + alpha_center)) * conc_center +
conc_center +
sy * alpha_center * conc_bc; sy * alpha_center * conc_bc;
} }