Commit 69c02886 authored by Santiago Ospina's avatar Santiago Ospina
Browse files

[Flux Reconstruction] Fix erroneus symmetry term switch in Richards operator

parent af8f93b7
......@@ -792,7 +792,7 @@ public:
normal_flux *= std::abs( ig.centerUnitOuterNormal() * param.gravity() );
if (not diff_term)
return;
continue;
// update residual (flux term)
if constexpr (lopcase != LOPCase::RTVolume)
for (unsigned int i = 0; i<lfsv_s.size(); i++)
......@@ -864,7 +864,7 @@ public:
if constexpr (lopcase != LOPCase::RTSkeleton)
{
if (sym_term)
if (not sym_term)
continue;
// update residual (symmetry term)
// (non-)symmetric IP term
......@@ -958,21 +958,25 @@ public:
time = t;
}
//! Switch symmetry term on/off in alpha methods.
void symmetry_term(bool t)
{
sym_term = t;
}
//! Switch diffusion term on/off in alpha methods.
void diffusion_term(bool t)
{
diff_term = t;
}
//! Switch consistency term on/off in alpha methods.
void consistency_term(bool t)
{
cs_term = t;
}
//! Switch penalty term on/off in alpha methods.
void penalty_term(bool t)
{
pty_term = t;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment