You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Ah, I see! what clang-tidy says makes a lot of sense!
So I think the origin of that first if (std::abs(cfz-csz) <= 1e-14) comes from Stone et.al 2018 paper (https://iopscience.iop.org/article/10.1086/588755). On page 171, it says "In the degenerate case in which Ca=Cax=a, so cfz =csz , then equation (A16) becomes alf= 1 and als= 0."
So the question is, if that first if is removed, then for that case when cfz=csz , would these terms alf and als approach the "theoretical" 1 and 0 or numerically have values far from it? or should the second and third ifs only be done if cfz is "significantly" different than csz and then avoid overwriting?
I'm not sure if trying to not overwrite the first if matters, since the tests in #2880 are passing
In computing the eigensystem of the MHD equations in
mhd_eigen.H
, we have logic like:So the values set in the first test are never used, since they are overwritten by the
else
in the next 2 tests.This was identified by
clang-tidy
The text was updated successfully, but these errors were encountered: