The flux variables should be simpler by using ifs of constexprs instead of template specializations.
if
constexpr
