Commit 1db0866a authored by Dominic Kempf's avatar Dominic Kempf

[bugfix][fastdg] Fix pattern assembly in FastDG

This fixes a bug that originate from a partial exadune port.
In Exadune, the pattern base classes have a fast version, which
was not ported to PDELab. However, in the absence of these, the
PatternEngine of the FastDGGridOperator needs to be more clever
and only operate on one entry of each block.
parent bb8579aa
Pipeline #15308 passed with stage
in 14 minutes and 48 seconds
......@@ -150,12 +150,15 @@ namespace Dune{
template<typename LFSVC, typename LFSUC>
void add_pattern(const LFSVC& lfsv_cache, const LFSUC& lfsu_cache, const LocalPattern& p)
{
for (size_type k=0; k<p.size(); ++k)
if (p.size() > 0)
{
local_assembler.add_entry(*pattern,
lfsv_cache,p[k].i(),
lfsu_cache,p[k].j()
lfsv_cache,p[0].i(),
lfsu_cache,p[0].j()
);
}
add_border_pattern(std::integral_constant<bool,LocalAssembler::isNonOverlapping>(),
lfsv_cache,
lfsu_cache,
......
Markdown is supported
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