Go to the documentation of this file. 1 #ifndef DUNE_PDELAB_GRIDOPERATOR_FASTDG_PATTERNENGINE_HH
2 #define DUNE_PDELAB_GRIDOPERATOR_FASTDG_PATTERNENGINE_HH
27 template<
typename TrialConstra
intsContainer,
typename TestConstra
intsContainer>
30 return cu.containsNonDirichletConstraints() || cv.containsNonDirichletConstraints();
37 typedef typename LA::LocalOperator
LOP;
40 typedef typename LA::LFSU
LFSU;
42 typedef typename LFSU::Traits::GridFunctionSpace
GFSU;
43 typedef typename LA::LFSV
LFSV;
45 typedef typename LFSV::Traits::GridFunctionSpace
GFSV;
52 typedef typename LA::Traits::MatrixPattern
Pattern;
65 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();
107 return local_assembler.doPatternVolume();
112 return local_assembler.doPatternSkeleton();
117 return local_assembler.doPatternBoundary();
122 return local_assembler.doPatternVolumePostSkeleton();
127 return Partitions::all;
134 template<
typename LFSVC,
typename LFSUC>
138 if (local_assembler.reconstructBorderEntries() &&
139 !communicationCache().initialized())
141 communicationCache().addEntries(lfsv_cache,lfsu_cache,
p);
145 template<
typename LFSVC,
typename LFSUC>
150 template<
typename LFSVC,
typename LFSUC>
154 local_assembler.add_entry(*pattern,
169 template<
typename EG,
typename LFSUC,
typename LFSVC>
170 void onUnbindLFSUV(
const EG & eg,
const LFSUC & lfsu_cache,
const LFSVC & lfsv_cache)
173 localpattern.clear();
176 template<
typename IG,
typename LFSUC,
typename LFSVC>
178 const LFSUC& lfsu_s_cache,
const LFSVC& lfsv_s_cache,
179 const LFSUC& lfsu_n_cache,
const LFSVC& lfsv_n_cache)
181 add_pattern(lfsv_s_cache,lfsu_n_cache,localpattern_sn);
182 localpattern_sn.clear();
183 add_pattern(lfsv_n_cache,lfsu_s_cache,localpattern_ns);
184 localpattern_ns.clear();
192 template<
typename EG,
typename LFSUC,
typename LFSVC>
196 pattern_volume(lop,lfsu_cache.localFunctionSpace(),lfsv_cache.localFunctionSpace(),localpattern);
199 template<
typename IG,
typename LFSUC,
typename LFSVC>
201 const LFSUC & lfsu_n_cache,
const LFSVC & lfsv_n_cache)
205 lfsu_s_cache.localFunctionSpace(),lfsv_s_cache.localFunctionSpace(),
206 lfsu_n_cache.localFunctionSpace(),lfsv_n_cache.localFunctionSpace(),
207 localpattern_sn, localpattern_ns);
210 template<
typename IG,
typename LFSUC,
typename LFSVC>
214 pattern_boundary(lop,lfsu_s_cache.localFunctionSpace(),lfsv_s_cache.localFunctionSpace(),localpattern);
217 template<
typename IG,
typename LFSUC,
typename LFSVC>
219 const LFSUC & lfsu_s_cache,
const LFSVC & lfsv_s_cache,
220 const LFSUC & lfsu_n_cache,
const LFSVC & lfsv_n_cache,
221 const LFSUC & lfsu_coupling_cache,
const LFSVC & lfsv_coupling_cache)
223 DUNE_THROW(Dune::NotImplemented,
"Assembling of coupling spaces is not implemented for ");
226 template<
typename IG,
typename LFSVC>
228 const LFSVC & lfsv_s_cache,
229 const LFSVC & lfsv_n_cache,
230 const LFSVC & lfsv_coupling_cache)
232 DUNE_THROW(Dune::NotImplemented,
"Assembling of coupling spaces is not implemented for ");
235 template<
typename EG,
typename LFSUC,
typename LFSVC>
251 if(local_assembler.doPostProcessing() &&
252 local_assembler.reconstructBorderEntries())
254 communicationCache().finishInitialization();
256 typename LA::Traits::BorderDOFExchanger::template PatternExtender<Pattern>
257 data_handle(*_border_dof_exchanger,gfsu,gfsv,*pattern);
258 gfsv.gridView().communicate(data_handle,
259 InteriorBorder_InteriorBorder_Interface,
260 ForwardCommunication);
272 typename LA::Traits::BorderDOFExchanger::CommunicationCache&
275 return _border_dof_exchanger->communicationCache();
278 const typename LA::Traits::BorderDOFExchanger::CommunicationCache&
279 communicationCache()
const
281 return _border_dof_exchanger->communicationCache();
299 shared_ptr<BorderDOFExchanger> _border_dof_exchanger;
const LocalAssembler & localAssembler() const
Public access to the wrapping local assembler.
Definition: fastdg/patternengine.hh:73
FastDGLocalPatternAssemblerEngine(const LocalAssembler &local_assembler_, shared_ptr< typename LA::Traits::BorderDOFExchanger > border_dof_exchanger)
Constructor.
Definition: fastdg/patternengine.hh:64
const LocalAssembler::Traits::TestGridFunctionSpaceConstraints & testConstraints() const
Test space constraints.
Definition: fastdg/patternengine.hh:85
void postAssembly(const GFSU &gfsu, const GFSV &gfsv)
Definition: fastdg/patternengine.hh:243
LA::Traits::MatrixPattern Pattern
The type of the solution vector.
Definition: fastdg/patternengine.hh:52
The local assembler for DUNE grids.
Definition: fastdg/localassembler.hh:38
static void pattern_skeleton(const LA &la, const LFSU &lfsu_s, const LFSV &lfsv_s, const LFSU &lfsu_n, const LFSV &lfsv_n, LocalPattern &pattern_sn, LocalPattern &pattern_ns)
Definition: callswitch.hh:78
void assembleUVVolume(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/patternengine.hh:193
static void pattern_volume(const LA &la, const LFSU &lfsu, const LFSV &lfsv, LocalPattern &pattern)
Definition: callswitch.hh:67
auto partition() const
Definition: fastdg/patternengine.hh:125
bool requireUVVolume() const
Definition: fastdg/patternengine.hh:105
For backward compatibility – Do not use this!
Definition: adaptivity.hh:28
void add_border_pattern(std::false_type, const LFSVC &lfsv_cache, const LFSUC &lfsu_cache, const LocalPattern &p)
Definition: fastdg/patternengine.hh:146
The fast DG local assembler engine for DUNE grids which creates the matrix pattern.
Definition: fastdg/patternengine.hh:22
static void assembleUVEnrichedCoupling(const IG &ig, const LFSUC &lfsu_s_cache, const LFSVC &lfsv_s_cache, const LFSUC &lfsu_n_cache, const LFSVC &lfsv_n_cache, const LFSUC &lfsu_coupling_cache, const LFSVC &lfsv_coupling_cache)
Definition: fastdg/patternengine.hh:218
std::size_t size_type
Definition: fastdg/patternengine.hh:56
const LocalAssembler::Traits::TrialGridFunctionSpaceConstraints & trialConstraints() const
Trial space constraints.
Definition: fastdg/patternengine.hh:79
LA::LFSUCache LFSUCache
Definition: fastdg/patternengine.hh:41
static void pattern_volume_post_skeleton(const LA &la, const LFSU &lfsu, const LFSV &lfsv, LocalPattern &pattern)
Definition: callswitch.hh:72
static void assembleVEnrichedCoupling(const IG &ig, const LFSVC &lfsv_s_cache, const LFSVC &lfsv_n_cache, const LFSVC &lfsv_coupling_cache)
Definition: fastdg/patternengine.hh:227
void assembleUVBoundary(const IG &ig, const LFSUC &lfsu_s_cache, const LFSVC &lfsv_s_cache)
Definition: fastdg/patternengine.hh:211
void onUnbindLFSUVOutside(const IG &ig, const LFSUC &lfsu_s_cache, const LFSVC &lfsv_s_cache, const LFSUC &lfsu_n_cache, const LFSVC &lfsv_n_cache)
Definition: fastdg/patternengine.hh:177
LA::LFSU LFSU
The local function spaces.
Definition: fastdg/patternengine.hh:40
LA::LocalOperator LOP
The type of the local operator.
Definition: fastdg/patternengine.hh:37
void post_border_pattern_assembly(std::false_type, const GFSU &gfsu, const GFSV &gfsv)
Definition: fastdg/patternengine.hh:264
LFSU::Traits::GridFunctionSpace GFSU
Definition: fastdg/patternengine.hh:42
LFSV::Traits::GridFunctionSpace GFSV
Definition: fastdg/patternengine.hh:45
void add_border_pattern(std::true_type, const LFSVC &lfsv_cache, const LFSUC &lfsu_cache, const LocalPattern &p)
Definition: fastdg/patternengine.hh:135
const IG & ig
Definition: constraints.hh:148
BorderDOFExchanger::BorderPattern BorderPattern
Definition: fastdg/patternengine.hh:49
LA LocalAssembler
The type of the wrapping local assembler.
Definition: fastdg/patternengine.hh:34
void setPattern(Pattern &pattern_)
Definition: fastdg/patternengine.hh:92
Base class for LocalAssemblerEngine implementations to avoid boilerplate code.
Definition: localassemblerenginebase.hh:21
Layout description for a sparse linear operator.
Definition: assemblerutilities.hh:129
bool requireUVSkeleton() const
Definition: fastdg/patternengine.hh:110
LA::LFSVCache LFSVCache
Definition: fastdg/patternengine.hh:44
void post_border_pattern_assembly(std::true_type, const GFSU &gfsu, const GFSV &gfsv)
Definition: fastdg/patternengine.hh:249
Dune::PDELab::LocalSparsityPattern LocalPattern
Definition: fastdg/patternengine.hh:54
LA::Traits::BorderDOFExchanger BorderDOFExchanger
helper classes
Definition: fastdg/patternengine.hh:48
bool requireUVBoundary() const
Definition: fastdg/patternengine.hh:115
bool requireSkeleton() const
Definition: fastdg/patternengine.hh:100
bool needsConstraintsCaching(const TrialConstraintsContainer &cu, const TestConstraintsContainer &cv) const
Definition: fastdg/patternengine.hh:28
LA::LFSV LFSV
Definition: fastdg/patternengine.hh:43
static void pattern_boundary(const LA &la, const LFSU &lfsu_s, const LFSV &lfsv_s, LocalPattern &pattern_ss)
Definition: callswitch.hh:85
void assembleUVVolumePostSkeleton(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/patternengine.hh:236
void onUnbindLFSUV(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/patternengine.hh:170
bool requireUVVolumePostSkeleton() const
Definition: fastdg/patternengine.hh:120
void add_pattern(const LFSVC &lfsv_cache, const LFSUC &lfsu_cache, const LocalPattern &p)
Definition: fastdg/patternengine.hh:151
void assembleUVSkeleton(const IG &ig, const LFSUC &lfsu_s_cache, const LFSVC &lfsv_s_cache, const LFSUC &lfsu_n_cache, const LFSVC &lfsv_n_cache)
Definition: fastdg/patternengine.hh:200
const P & p
Definition: constraints.hh:147