1#ifndef DUNE_PDELAB_GRIDOPERATOR_DEFAULT_PATTERNENGINE_HH
2#define DUNE_PDELAB_GRIDOPERATOR_DEFAULT_PATTERNENGINE_HH
5#include <dune/pdelab/gridoperator/common/localmatrix.hh>
6#include <dune/pdelab/gridoperator/common/assemblerutilities.hh>
7#include <dune/pdelab/gridoperator/common/localassemblerenginebase.hh>
8#include <dune/pdelab/constraints/common/constraints.hh>
9#include <dune/pdelab/localoperator/callswitch.hh>
27 template<
typename TrialConstra
intsContainer,
typename TestConstra
intsContainer>
28 bool needsConstraintsCaching(
const TrialConstraintsContainer& cu,
const TestConstraintsContainer& cv)
const
30 return cu.containsNonDirichletConstraints() || cv.containsNonDirichletConstraints();
37 typedef typename LA::LocalOperator
LOP;
40 typedef typename LA::LFSU
LFSU;
41 typedef typename LA::LFSUCache LFSUCache;
42 typedef typename LFSU::Traits::GridFunctionSpace GFSU;
43 typedef typename LA::LFSV LFSV;
44 typedef typename LA::LFSVCache LFSVCache;
45 typedef typename LFSV::Traits::GridFunctionSpace GFSV;
49 typedef typename BorderDOFExchanger::BorderPattern BorderPattern;
52 typedef typename LA::Traits::MatrixPattern
Pattern;
56 typedef std::size_t size_type;
65 std::shared_ptr<typename LA::Traits::BorderDOFExchanger> border_dof_exchanger)
66 : local_assembler(local_assembler_)
67 , lop(local_assembler.localOperator())
69 , _border_dof_exchanger(border_dof_exchanger)
75 return local_assembler;
79 const typename LocalAssembler::Traits::TrialGridFunctionSpaceConstraints&
trialConstraints()
const
85 const typename LocalAssembler::Traits::TestGridFunctionSpaceConstraints&
testConstraints()
const
102 return local_assembler.doPatternSkeleton();
105 bool requireUVVolume()
const
107 return local_assembler.doPatternVolume();
110 bool requireUVSkeleton()
const
112 return local_assembler.doPatternSkeleton();
115 bool requireUVBoundary()
const
117 return local_assembler.doPatternBoundary();
120 bool requireUVVolumePostSkeleton()
const
122 return local_assembler.doPatternVolumePostSkeleton();
129 template<
typename LFSVC,
typename LFSUC>
130 void add_border_pattern(std::true_type,
const LFSVC& lfsv_cache,
const LFSUC& lfsu_cache,
131 const LocalPattern& p)
133 if (local_assembler.reconstructBorderEntries() &&
134 !communicationCache().initialized())
136 communicationCache().addEntries(lfsv_cache,lfsu_cache,p);
140 template<
typename LFSVC,
typename LFSUC>
141 void add_border_pattern(std::false_type,
const LFSVC& lfsv_cache,
const LFSUC& lfsu_cache,
142 const LocalPattern& p)
145 template<
typename LFSVC,
typename LFSUC>
146 void add_pattern(
const LFSVC& lfsv_cache,
const LFSUC& lfsu_cache,
const LocalPattern& p)
148 for (size_type k=0; k<p.size(); ++k)
149 local_assembler.add_entry(*pattern,
154 add_border_pattern(std::integral_constant<bool,LocalAssembler::isNonOverlapping>(),
164 template<
typename EG,
typename LFSUC,
typename LFSVC>
165 void onUnbindLFSUV(
const EG & eg,
const LFSUC & lfsu_cache,
const LFSVC & lfsv_cache)
167 add_pattern(lfsv_cache,lfsu_cache,localpattern);
168 localpattern.clear();
171 template<
typename IG,
typename LFSUC,
typename LFSVC>
172 void onUnbindLFSUVOutside(
const IG& ig,
173 const LFSUC& lfsu_s_cache,
const LFSVC& lfsv_s_cache,
174 const LFSUC& lfsu_n_cache,
const LFSVC& lfsv_n_cache)
176 add_pattern(lfsv_s_cache,lfsu_n_cache,localpattern_sn);
177 localpattern_sn.clear();
178 add_pattern(lfsv_n_cache,lfsu_s_cache,localpattern_ns);
179 localpattern_ns.clear();
187 template<
typename EG,
typename LFSUC,
typename LFSVC>
190 Dune::PDELab::LocalAssemblerCallSwitch<LOP,LOP::doPatternVolume>::
191 pattern_volume(lop,lfsu_cache.localFunctionSpace(),lfsv_cache.localFunctionSpace(),localpattern);
194 template<
typename IG,
typename LFSUC,
typename LFSVC>
195 void assembleUVSkeleton(
const IG & ig,
const LFSUC & lfsu_s_cache,
const LFSVC & lfsv_s_cache,
196 const LFSUC & lfsu_n_cache,
const LFSVC & lfsv_n_cache)
198 Dune::PDELab::LocalAssemblerCallSwitch<LOP,LOP::doPatternSkeleton>::
199 pattern_skeleton(lop,
200 lfsu_s_cache.localFunctionSpace(),lfsv_s_cache.localFunctionSpace(),
201 lfsu_n_cache.localFunctionSpace(),lfsv_n_cache.localFunctionSpace(),
202 localpattern_sn, localpattern_ns);
205 template<
typename IG,
typename LFSUC,
typename LFSVC>
206 void assembleUVBoundary(
const IG & ig,
const LFSUC & lfsu_s_cache,
const LFSVC & lfsv_s_cache)
208 Dune::PDELab::LocalAssemblerCallSwitch<LOP,LOP::doPatternBoundary>::
209 pattern_boundary(lop,lfsu_s_cache.localFunctionSpace(),lfsv_s_cache.localFunctionSpace(),localpattern);
212 template<
typename IG,
typename LFSUC,
typename LFSVC>
213 static void assembleUVEnrichedCoupling(
const IG & ig,
214 const LFSUC & lfsu_s_cache,
const LFSVC & lfsv_s_cache,
215 const LFSUC & lfsu_n_cache,
const LFSVC & lfsv_n_cache,
216 const LFSUC & lfsu_coupling_cache,
const LFSVC & lfsv_coupling_cache)
221 template<
typename IG,
typename LFSVC>
222 static void assembleVEnrichedCoupling(
const IG & ig,
223 const LFSVC & lfsv_s_cache,
224 const LFSVC & lfsv_n_cache,
225 const LFSVC & lfsv_coupling_cache)
230 template<
typename EG,
typename LFSUC,
typename LFSVC>
231 void assembleUVVolumePostSkeleton(
const EG & eg,
const LFSUC & lfsu_cache,
const LFSVC & lfsv_cache)
233 Dune::PDELab::LocalAssemblerCallSwitch<LOP,LOP::doPatternVolumePostSkeleton>::
234 pattern_volume_post_skeleton(lop,lfsu_cache.localFunctionSpace(),lfsv_cache.localFunctionSpace(),localpattern);
238 void postAssembly(
const GFSU& gfsu,
const GFSV& gfsv)
240 post_border_pattern_assembly(std::integral_constant<bool,LocalAssembler::isNonOverlapping>(),
245 void post_border_pattern_assembly(std::true_type,
const GFSU& gfsu,
const GFSV& gfsv)
247 if(local_assembler.doPostProcessing() and
248 local_assembler.reconstructBorderEntries())
250 communicationCache().finishInitialization();
252 typename LA::Traits::BorderDOFExchanger::template PatternExtender<Pattern>
253 data_handle(*_border_dof_exchanger,gfsu,gfsv,*pattern);
254 gfsv.gridView().communicate(data_handle,
260 void post_border_pattern_assembly(std::false_type,
const GFSU& gfsu,
const GFSV& gfsv)
268 typename LA::Traits::BorderDOFExchanger::CommunicationCache&
271 return _border_dof_exchanger->communicationCache();
274 const typename LA::Traits::BorderDOFExchanger::CommunicationCache&
275 communicationCache()
const
277 return _border_dof_exchanger->communicationCache();
290 LocalPattern localpattern;
291 LocalPattern localpattern_sn, localpattern_ns;
293 BorderPattern _border_pattern;
295 std::shared_ptr<BorderDOFExchanger> _border_dof_exchanger;
Default exception for dummy implementations.
Definition: exceptions.hh:261
The local assembler engine for DUNE grids which creates the matrix pattern.
Definition: patternengine.hh:24
LA::Traits::MatrixPattern Pattern
The type of the solution vector.
Definition: patternengine.hh:52
DefaultLocalPatternAssemblerEngine(const LocalAssembler &local_assembler_, std::shared_ptr< typename LA::Traits::BorderDOFExchanger > border_dof_exchanger)
Constructor.
Definition: patternengine.hh:64
LA::LocalOperator LOP
The type of the local operator.
Definition: patternengine.hh:37
LA LocalAssembler
The type of the wrapping local assembler.
Definition: patternengine.hh:34
bool requireSkeleton() const
Definition: patternengine.hh:100
const LocalAssembler & localAssembler() const
Public access to the wrapping local assembler.
Definition: patternengine.hh:73
void setPattern(Pattern &pattern_)
Definition: patternengine.hh:92
void onUnbindLFSUV(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: patternengine.hh:165
LA::LFSU LFSU
The local function spaces.
Definition: patternengine.hh:40
const LocalAssembler::Traits::TestGridFunctionSpaceConstraints & testConstraints() const
Test space constraints.
Definition: patternengine.hh:85
const LocalAssembler::Traits::TrialGridFunctionSpaceConstraints & trialConstraints() const
Trial space constraints.
Definition: patternengine.hh:79
LA::Traits::BorderDOFExchanger BorderDOFExchanger
helper classes
Definition: patternengine.hh:48
void assembleUVVolume(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: patternengine.hh:188
Base class for LocalAssemblerEngine implementations to avoid boilerplate code.
Definition: localassemblerenginebase.hh:22
Layout description for a sparse linear operator.
Definition: assemblerutilities.hh:137
#define DUNE_THROW(E, m)
Definition: exceptions.hh:216
Dune namespace.
Definition: alignedallocator.hh:14