DUNE PDELab (2.7)

enginebase.hh
1#ifndef DUNE_PDELAB_GRIDOPERATOR_ONESTEP_ENGINEBASE_HH
2#define DUNE_PDELAB_GRIDOPERATOR_ONESTEP_ENGINEBASE_HH
3
4namespace Dune{
5 namespace PDELab{
6
14 template<typename OSLA, typename LAE0, typename LAE1>
15 class OneStepLocalAssemblerEngineBase
16 {
17 public:
19 typedef OSLA OneStepLocalAssembler;
20
21 typedef typename LAE0::Traits Traits;
22
23 template<typename TrialConstraintsContainer, typename TestConstraintsContainer>
24 bool needsConstraintsCaching(const TrialConstraintsContainer& cu, const TestConstraintsContainer& cv) const
25 {
26 return (lae0->needsConstraintsCaching(cu,cv) or lae1->needsConstraintsCaching(cu,cv));
27 }
28
29
32 typedef typename OSLA::LocalAssemblerDT0 LocalAssemblerDT0;
33 typedef typename OSLA::LocalAssemblerDT1 LocalAssemblerDT1;
34
35 typedef LAE0 LocalAssemblerEngineDT0;
36 typedef LAE1 LocalAssemblerEngineDT1;
38
40 typedef typename OSLA::Real Real;
41
42 typedef OSLA LocalAssembler;
43
50 OneStepLocalAssemblerEngineBase(const LocalAssembler & local_assembler_)
51 : invalid_lae0(nullptr),
52 invalid_lae1(nullptr),
53 la(local_assembler_),
54 lae0(invalid_lae0), lae1(invalid_lae1),
55 implicit(true)
56 {}
57
60 bool requireSkeleton() const
61 { return implicit && (lae0->requireSkeleton() || lae1->requireSkeleton()); }
62 bool requireSkeletonTwoSided() const
63 { return lae0->requireSkeletonTwoSided() || lae1->requireSkeletonTwoSided(); }
64 bool requireUVVolume() const
65 { return lae0->requireUVVolume() || lae1->requireUVVolume(); }
66 bool requireVVolume() const
67 { return lae0->requireVVolume() || lae1->requireVVolume(); }
68 bool requireUVSkeleton() const
69 { return lae0->requireUVSkeleton() || lae1->requireUVSkeleton(); }
70 bool requireVSkeleton() const
71 { return lae0->requireVSkeleton() || lae1->requireVSkeleton(); }
72 bool requireUVBoundary() const
73 { return lae0->requireUVBoundary() || lae1->requireUVBoundary(); }
74 bool requireVBoundary() const
75 { return lae0->requireVBoundary() || lae1->requireVBoundary(); }
76 bool requireUVProcessor() const
77 { return lae0->requireUVProcessor() || lae1->requireUVProcessor(); }
78 bool requireVProcessor() const
79 { return lae0->requireVProcessor() || lae1->requireVProcessor(); }
80 bool requireUVEnrichedCoupling() const
81 { return lae0->requireUVEnrichedCoupling() || lae1->requireUVEnrichedCoupling(); }
82 bool requireVEnrichedCoupling() const
83 { return lae0->requireVEnrichedCoupling() || lae1->requireVEnrichedCoupling(); }
84 bool requireUVVolumePostSkeleton() const
85 { return lae0->requireUVVolumePostSkeleton() || lae1->requireUVVolumePostSkeleton();}
86 bool requireVVolumePostSkeleton() const
87 { return lae0->requireVVolumePostSkeleton() || lae1->requireVVolumePostSkeleton(); }
89
90
92 const LocalAssembler & localAssembler(){ return la; }
93
94 LocalAssemblerEngineDT0& localAssemblerEngineDT0()
95 {
96 return *lae0;
97 }
98
99 const LocalAssemblerEngineDT0& localAssemblerEngineDT0() const
100 {
101 return *lae0;
102 }
103
104 LocalAssemblerEngineDT1& localAssemblerEngineDT1()
105 {
106 return *lae1;
107 }
108
109 const LocalAssemblerEngineDT1& localAssemblerEngineDT1() const
110 {
111 return *lae1;
112 }
113
114 auto partition() const
115 {
116 return localAssemblerEngineDT0().partition();
117 }
118
119 void setLocalAssemblerEngineDT0(LocalAssemblerEngineDT0& lae0_)
120 {
121 lae0 = &lae0_;
122 }
123
124 void setLocalAssemblerEngineDT1(LocalAssemblerEngineDT1& lae1_)
125 {
126 lae1 = &lae1_;
127 }
128
129 const typename OneStepLocalAssembler::Traits::TrialGridFunctionSpaceConstraints& trialConstraints() const
130 {
131 return localAssemblerEngineDT0().trialConstraints();
132 }
133
134 const typename OneStepLocalAssembler::Traits::TestGridFunctionSpaceConstraints& testConstraints() const
135 {
136 return localAssemblerEngineDT0().testConstraints();
137 }
138
142 template<typename EG, typename LFSU, typename LFSV>
143 void onBindLFSUV(const EG & eg, const LFSU & lfsu, const LFSV & lfsv)
144 {
145 lae0->onBindLFSUV(eg,lfsu,lfsv);
146 lae1->onBindLFSUV(eg,lfsu,lfsv);
147 }
148
149 template<typename EG, typename LFSV>
150 void onBindLFSV(const EG & eg, const LFSV & lfsv)
151 {
152 lae0->onBindLFSV(eg,lfsv);
153 lae1->onBindLFSV(eg,lfsv);
154 }
155
156 template<typename IG, typename LFSU, typename LFSV>
157 void onBindLFSUVInside(const IG & ig, const LFSU & lfsu, const LFSV & lfsv)
158 {
159 lae0->onBindLFSUVInside(ig,lfsu,lfsv);
160 lae1->onBindLFSUVInside(ig,lfsu,lfsv);
161 }
162
163 template<typename IG,
164 typename LFSU_S, typename LFSV_S,
165 typename LFSU_N, typename LFSV_N>
166 void onBindLFSUVOutside(const IG & ig,
167 const LFSU_S & lfsu_s, const LFSV_S & lfsv_s,
168 const LFSU_N & lfsu_n, const LFSV_N & lfsv_n)
169 {
170 lae0->onBindLFSUVOutside(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n);
171 lae1->onBindLFSUVOutside(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n);
172 }
173
174 template<typename IG, typename LFSV>
175 void onBindLFSVInside(const IG & ig, const LFSV & lfsv)
176 {
177 lae0->onBindLFSVInside(ig,lfsv);
178 lae1->onBindLFSVInside(ig,lfsv);
179 }
180
181 template<typename IG,
182 typename LFSV_S,
183 typename LFSV_N>
184 void onBindLFSVOutside(const IG & ig,
185 const LFSV_S & lfsv_s,
186 const LFSV_N & lfsv_n)
187 {
188 lae0->onBindLFSVOutside(ig,lfsv_s,lfsv_n);
189 lae1->onBindLFSVOutside(ig,lfsv_s,lfsv_n);
190 }
191
192 template<typename IG,
193 typename LFSU_S, typename LFSV_S,
194 typename LFSU_N, typename LFSV_N,
195 typename LFSU_C, typename LFSV_C>
196 void onBindLFSUVCoupling(const IG & ig,
197 const LFSU_S & lfsu_s, const LFSV_S & lfsv_s,
198 const LFSU_N & lfsu_n, const LFSV_N & lfsv_n,
199 const LFSU_C & lfsu_c, const LFSV_C & lfsv_c)
200 {
201 lae0->onBindLFSUVCoupling(ig,
202 lfsu_s,lfsv_s,
203 lfsu_n,lfsv_n,
204 lfsu_c,lfsv_c);
205 lae1->onBindLFSUVCoupling(ig,
206 lfsu_s,lfsv_s,
207 lfsu_n,lfsv_n,
208 lfsu_c,lfsv_c);
209 }
210
211 template<typename IG,
212 typename LFSV_S,
213 typename LFSV_N,
214 typename LFSV_C>
215 void onBindLFSVCoupling(const IG & ig,
216 const LFSV_S & lfsv_s,
217 const LFSV_N & lfsv_n,
218 const LFSV_C & lfsv_c)
219 {
220 lae0->onBindLFSVCoupling(ig,lfsv_s,lfsv_n,lfsv_c);
221 lae1->onBindLFSVCoupling(ig,lfsv_s,lfsv_n,lfsv_c);
222 }
223
225
229 template<typename EG, typename LFSU, typename LFSV>
230 void onUnbindLFSUV(const EG & eg, const LFSU & lfsu, const LFSV & lfsv)
231 {
232 lae0->onUnbindLFSUV(eg,lfsu, lfsv);
233 lae1->onUnbindLFSUV(eg,lfsu, lfsv);
234 }
235
236 template<typename EG, typename LFSV>
237 void onUnbindLFSV(const EG & eg, const LFSV & lfsv)
238 {
239 lae0->onUnbindLFSV(eg,lfsv);
240 lae1->onUnbindLFSV(eg,lfsv);
241 }
242
243 template<typename IG, typename LFSU, typename LFSV>
244 void onUnbindLFSUVInside(const IG & ig, const LFSU & lfsu, const LFSV & lfsv)
245 {
246 lae0->onUnbindLFSUVInside(ig,lfsu, lfsv);
247 lae1->onUnbindLFSUVInside(ig,lfsu, lfsv);
248 }
249
250 template<typename IG,
251 typename LFSU_S, typename LFSV_S,
252 typename LFSU_N, typename LFSV_N>
253 void onUnbindLFSUVOutside(const IG & ig,
254 const LFSU_S & lfsu_s, const LFSV_S & lfsv_s,
255 const LFSU_N & lfsu_n, const LFSV_N & lfsv_n)
256 {
257 lae0->onUnbindLFSUVOutside(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n);
258 lae1->onUnbindLFSUVOutside(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n);
259 }
260
261 template<typename IG, typename LFSV>
262 void onUnbindLFSVInside(const IG & ig, const LFSV & lfsv)
263 {
264 lae0->onUnbindLFSVInside(ig,lfsv);
265 lae1->onUnbindLFSVInside(ig,lfsv);
266 }
267
268 template<typename IG,
269 typename LFSV_S,
270 typename LFSV_N>
271 void onUnbindLFSVOutside(const IG & ig,
272 const LFSV_S & lfsv_s,
273 const LFSV_N & lfsv_n)
274 {
275 lae0->onUnbindLFSVOutside(ig,lfsv_s,lfsv_n);
276 lae1->onUnbindLFSVOutside(ig,lfsv_s,lfsv_n);
277 }
278
279 template<typename IG,
280 typename LFSU_S, typename LFSV_S,
281 typename LFSU_N, typename LFSV_N,
282 typename LFSU_C, typename LFSV_C>
283 void onUnbindLFSUVCoupling(const IG & ig,
284 const LFSU_S & lfsu_s, const LFSV_S & lfsv_s,
285 const LFSU_N & lfsu_n, const LFSV_N & lfsv_n,
286 const LFSU_C & lfsu_c, const LFSV_C & lfsv_c)
287 {
288 lae0->onUnbindLFSUVCoupling(ig,
289 lfsu_s,lfsv_s,
290 lfsu_n,lfsv_n,
291 lfsu_c,lfsv_c);
292 lae1->onUnbindLFSUVCoupling(ig,
293 lfsu_s,lfsv_s,
294 lfsu_n,lfsv_n,
295 lfsu_c,lfsv_c);
296 }
297
298 template<typename IG,
299 typename LFSV_S,
300 typename LFSV_N,
301 typename LFSV_C>
302 void onUnbindLFSVCoupling(const IG & ig,
303 const LFSV_S & lfsv_s,
304 const LFSV_N & lfsv_n,
305 const LFSV_C & lfsv_c)
306 {
307 lae0->onUnbindLFSVCoupling(ig,lfsv_s,lfsv_n,lfsv_c);
308 lae1->onUnbindLFSVCoupling(ig,lfsv_s,lfsv_n,lfsv_c);
309 }
310
312
316 template<typename LFSU>
317 void loadCoefficientsLFSUInside(const LFSU & lfsu_s)
318 {
319 lae0->loadCoefficientsLFSUInside(lfsu_s);
320 lae1->loadCoefficientsLFSUInside(lfsu_s);
321 }
322 template<typename LFSU>
323 void loadCoefficientsLFSUOutside(const LFSU & lfsu_n)
324 {
325 lae0->loadCoefficientsLFSUOutside(lfsu_n);
326 lae1->loadCoefficientsLFSUOutside(lfsu_n);
327 }
328 template<typename LFSU>
329 void loadCoefficientsLFSUCoupling(const LFSU & lfsu_c)
330 {
331 lae0->loadCoefficientsLFSUCoupling(lfsu_c);
332 lae1->loadCoefficientsLFSUCoupling(lfsu_c);
333 }
335
338
339 template<typename EG>
340 bool assembleCell(const EG & eg)
341 {
342 bool rv = true;
343 rv &= lae0->assembleCell(eg);
344 rv &= lae1->assembleCell(eg);
345 return rv;
346 }
347
348 template<typename EG, typename LFSU, typename LFSV>
349 void assembleUVVolume(const EG & eg, const LFSU & lfsu, const LFSV & lfsv)
350 {
351 if(implicit)
352 lae0->assembleUVVolume(eg,lfsu,lfsv);
353 lae1->assembleUVVolume(eg,lfsu,lfsv);
354 }
355
356 template<typename EG, typename LFSV>
357 void assembleVVolume(const EG & eg, const LFSV & lfsv)
358 {
359 if(implicit)
360 lae0->assembleVVolume(eg,lfsv);
361 lae1->assembleVVolume(eg,lfsv);
362 }
363
364 template<typename IG, typename LFSU_S, typename LFSV_S, typename LFSU_N, typename LFSV_N>
365 void assembleUVSkeleton(const IG & ig, const LFSU_S & lfsu_s, const LFSV_S & lfsv_s,
366 const LFSU_N & lfsu_n, const LFSV_N & lfsv_n)
367 {
368 if(implicit)
369 lae0->assembleUVSkeleton(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n);
370 lae1->assembleUVSkeleton(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n);
371 }
372
373 template<typename IG, typename LFSV_S, typename LFSV_N>
374 void assembleVSkeleton(const IG & ig, const LFSV_S & lfsv_s, const LFSV_N & lfsv_n)
375 {
376 if(implicit)
377 lae0->assembleVSkeleton(ig,lfsv_s,lfsv_n);
378 lae1->assembleVSkeleton(ig,lfsv_s,lfsv_n);
379 }
380
381 template<typename IG, typename LFSU_S, typename LFSV_S>
382 void assembleUVBoundary(const IG & ig, const LFSU_S & lfsu_s, const LFSV_S & lfsv_s)
383 {
384 if(implicit)
385 lae0->assembleUVBoundary(ig,lfsu_s,lfsv_s);
386 lae1->assembleUVBoundary(ig,lfsu_s,lfsv_s);
387 }
388
389 template<typename IG, typename LFSV_S>
390 void assembleVBoundary(const IG & ig, const LFSV_S & lfsv_s)
391 {
392 if(implicit)
393 lae0->assembleVBoundary(ig,lfsv_s);
394 lae1->assembleVBoundary(ig,lfsv_s);
395 }
396
397 template<typename IG, typename LFSU_S, typename LFSV_S>
398 void assembleUVProcessor(const IG & ig, const LFSU_S & lfsu_s, const LFSV_S & lfsv_s)
399 {
400 if(implicit)
401 lae0->assembleUVProcessor(ig,lfsu_s,lfsv_s);
402 lae1->assembleUVProcessor(ig,lfsu_s,lfsv_s);
403 }
404
405 template<typename IG, typename LFSV_S>
406 void assembleVProcessor(const IG & ig, const LFSV_S & lfsv_s)
407 {
408 if(implicit)
409 lae0->assembleVProcessor(ig,lfsv_s);
410 lae1->assembleVProcessor(ig,lfsv_s);
411 }
412
413 template<typename IG, typename LFSU_S, typename LFSV_S, typename LFSU_N, typename LFSV_N,
414 typename LFSU_C, typename LFSV_C>
415 void assembleUVEnrichedCoupling(const IG & ig,
416 const LFSU_S & lfsu_s, const LFSV_S & lfsv_s,
417 const LFSU_N & lfsu_n, const LFSV_N & lfsv_n,
418 const LFSU_C & lfsu_c, const LFSV_C & lfsv_c)
419 {
420 if(implicit)
421 lae0->assembleUVEnrichedCoupling(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n,lfsu_c,lfsv_c);
422 lae1->assembleUVEnrichedCoupling(ig,lfsu_s,lfsv_s,lfsu_n,lfsv_n,lfsu_c,lfsv_c);
423 }
424
425 template<typename IG, typename LFSV_S, typename LFSV_N, typename LFSV_C>
426 void assembleVEnrichedCoupling(const IG & ig,
427 const LFSV_S & lfsv_s,
428 const LFSV_N & lfsv_n,
429 const LFSV_C & lfsv_c)
430 {
431 if(implicit)
432 lae0->assembleVEnrichedCoupling(ig,lfsv_s,lfsv_n,lfsv_c);
433 }
434
435 template<typename EG, typename LFSU, typename LFSV>
436 void assembleUVVolumePostSkeleton(const EG & eg, const LFSU & lfsu, const LFSV & lfsv)
437 {
438 if(implicit)
439 lae0->assembleUVVolumePostSkeleton(eg,lfsu,lfsv);
440 }
441
442 template<typename EG, typename LFSV>
443 void assembleVVolumePostSkeleton(const EG & eg, const LFSV & lfsv)
444 {
445 if(implicit)
446 lae0->assembleVVolumePostSkeleton(eg,lfsv);
447 }
449
450 private:
451
452 LocalAssemblerEngineDT0 * const invalid_lae0;
453 LocalAssemblerEngineDT1 * const invalid_lae1;
454
455 protected:
456
459 const LocalAssembler & la;
460
461 LocalAssemblerEngineDT0 * lae0;
462 LocalAssemblerEngineDT1 * lae1;
463
464 bool implicit;
465
466 }; // End of class OneStepLocalAssemblerEngineBase
467
468 }
469}
470
471#endif // DUNE_PDELAB_GRIDOPERATOR_ONESTEP_ENGINEBASE_HH
Dune::PDELab::OneStepLocalAssemblerEngineBase
The local assembler engine for UDG sub triangulations which assembles the residual vector.
Definition: enginebase.hh:16
Dune::PDELab::OneStepLocalAssemblerEngineBase::onUnbindLFSUV
void onUnbindLFSUV(const EG &eg, const LFSU &lfsu, const LFSV &lfsv)
Definition: enginebase.hh:230
Dune::PDELab::OneStepLocalAssemblerEngineBase::la
const LocalAssembler & la
Definition: enginebase.hh:459
Dune::PDELab::OneStepLocalAssemblerEngineBase::LocalAssemblerDT0
OSLA::LocalAssemblerDT0 LocalAssemblerDT0
Definition: enginebase.hh:32
Dune::PDELab::OneStepLocalAssemblerEngineBase::localAssembler
const LocalAssembler & localAssembler()
Public access to the wrapping local assembler.
Definition: enginebase.hh:92
Dune::PDELab::OneStepLocalAssemblerEngineBase::OneStepLocalAssemblerEngineBase
OneStepLocalAssemblerEngineBase(const LocalAssembler &local_assembler_)
Constructor.
Definition: enginebase.hh:50
Dune::PDELab::OneStepLocalAssemblerEngineBase::loadCoefficientsLFSUInside
void loadCoefficientsLFSUInside(const LFSU &lfsu_s)
Definition: enginebase.hh:317
Dune::PDELab::OneStepLocalAssemblerEngineBase::onBindLFSUV
void onBindLFSUV(const EG &eg, const LFSU &lfsu, const LFSV &lfsv)
Definition: enginebase.hh:143
Dune::PDELab::OneStepLocalAssemblerEngineBase::Real
OSLA::Real Real
The type for real numbers.
Definition: enginebase.hh:40
Dune::PDELab::OneStepLocalAssemblerEngineBase::requireSkeleton
bool requireSkeleton() const
Definition: enginebase.hh:60
Dune::PDELab::OneStepLocalAssemblerEngineBase::OneStepLocalAssembler
OSLA OneStepLocalAssembler
The type of the wrapping local assembler.
Definition: enginebase.hh:19
Dune
Dune namespace.
Definition: alignedallocator.hh:14
Dune::PDELab::LocalAssemblerTraits::TrialGridFunctionSpaceConstraints
GO::Traits::TrialGridFunctionSpaceConstraints TrialGridFunctionSpaceConstraints
The type of the trial grid function space constraints.
Definition: assemblerutilities.hh:33
Dune::PDELab::LocalAssemblerTraits::TestGridFunctionSpaceConstraints
GO::Traits::TestGridFunctionSpaceConstraints TestGridFunctionSpaceConstraints
The type of the test grid function space constraints.
Definition: assemblerutilities.hh:36
Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.111.3 (Jul 15, 22:36, 2024)