dune-pdelab  2.7-git
fastdg/jacobianengine.hh
Go to the documentation of this file.
1 #ifndef DUNE_PDELAB_GRIDOPERATOR_FASTDG_JACOBIANENGINE_HH
2 #define DUNE_PDELAB_GRIDOPERATOR_FASTDG_JACOBIANENGINE_HH
3 
13 
14 namespace Dune{
15  namespace PDELab{
16 
24  template<typename LA>
27  {
28  public:
29 
30  template<typename TrialConstraintsContainer, typename TestConstraintsContainer>
31  bool needsConstraintsCaching(const TrialConstraintsContainer& cu, const TestConstraintsContainer& cv)
32  {
33  return cu.containsNonDirichletConstraints() || cv.containsNonDirichletConstraints();
34  }
35 
37  typedef LA LocalAssembler;
38 
40  typedef typename LA::LocalOperator LOP;
41 
43  typedef typename LA::LFSU LFSU;
44  typedef typename LA::LFSUCache LFSUCache;
45  typedef typename LFSU::Traits::GridFunctionSpace GFSU;
46  typedef typename LA::LFSV LFSV;
47  typedef typename LA::LFSVCache LFSVCache;
48  typedef typename LFSV::Traits::GridFunctionSpace GFSV;
49 
51  typedef typename LA::Traits::Jacobian Jacobian;
52  typedef typename Jacobian::ElementType JacobianElement;
53  typedef typename Jacobian::template AliasedLocalView<LFSVCache,LFSUCache> JacobianView;
54 
56  typedef typename LA::Traits::Solution Solution;
57  typedef typename Solution::ElementType SolutionElement;
58  typedef typename Solution::template ConstAliasedLocalView<LFSUCache> SolutionView;
59 
67  : local_assembler(local_assembler_),
68  lop(local_assembler_.localOperator())
69  {}
70 
72 
87  local_assembler(other.local_assembler), lop(other.lop),
88  global_s_s_view(other.global_s_s_view),
89  global_s_n_view(other.global_s_n_view),
90  global_a_ss_view(other.global_a_ss_view),
91  global_a_sn_view(other.global_a_sn_view),
92  global_a_ns_view(other.global_a_ns_view),
93  global_a_nn_view(other.global_a_nn_view)
94  { }
95 
98  bool requireSkeleton() const
99  { return local_assembler.doAlphaSkeleton(); }
101  { return local_assembler.doSkeletonTwoSided(); }
102  bool requireUVVolume() const
103  { return local_assembler.doAlphaVolume(); }
104  bool requireUVSkeleton() const
105  { return local_assembler.doAlphaSkeleton(); }
106  bool requireUVBoundary() const
107  { return local_assembler.doAlphaBoundary(); }
109  { return local_assembler.doAlphaVolumePostSkeleton(); }
111 
114  {
115  return local_assembler;
116  }
117 
119  const typename LocalAssembler::Traits::TrialGridFunctionSpaceConstraints& trialConstraints() const
120  {
121  return localAssembler().trialConstraints();
122  }
123 
125  const typename LocalAssembler::Traits::TestGridFunctionSpaceConstraints& testConstraints() const
126  {
127  return localAssembler().testConstraints();
128  }
129 
132  void setJacobian(Jacobian & jacobian_)
133  {
134  global_a_ss_view.attach(jacobian_);
135  global_a_sn_view.attach(jacobian_);
136  global_a_ns_view.attach(jacobian_);
137  global_a_nn_view.attach(jacobian_);
138  }
139 
142  void setSolution(const Solution & solution_)
143  {
144  global_s_s_view.attach(solution_);
145  global_s_n_view.attach(solution_);
146  }
147 
151  template<typename EG, typename LFSUC, typename LFSVC>
152  void onBindLFSUV(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache)
153  {
154  global_s_s_view.bind(lfsu_cache);
155  global_a_ss_view.bind(lfsv_cache,lfsu_cache);
156  }
157 
158  template<typename IG, typename LFSUC, typename LFSVC>
159  void onBindLFSUVOutside(const IG & ig,
160  const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
161  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache)
162  {
163  global_s_n_view.bind(lfsu_n_cache);
164  global_a_sn_view.bind(lfsv_s_cache,lfsu_n_cache);
165  global_a_ns_view.bind(lfsv_n_cache,lfsu_s_cache);
166  global_a_nn_view.bind(lfsv_n_cache,lfsu_n_cache);
167  }
168 
170 
174  template<typename EG, typename LFSUC, typename LFSVC>
175  void onUnbindLFSUV(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache)
176  {
177  global_a_ss_view.unbind();
178  }
179 
180  template<typename IG, typename LFSUC, typename LFSVC>
181  void onUnbindLFSUVOutside(const IG & ig,
182  const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
183  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache)
184  {
185  global_a_sn_view.unbind();
186  global_a_ns_view.unbind();
187  global_a_nn_view.unbind();
188  }
189 
191 
194  template<typename LFSUC>
195  void loadCoefficientsLFSUInside(const LFSUC & lfsu_cache)
196  {}
197  template<typename LFSUC>
198  void loadCoefficientsLFSUOutside(const LFSUC & lfsu_n_cache)
199  {}
200  template<typename LFSUC>
201  void loadCoefficientsLFSUCoupling(const LFSUC & lfsu_c_cache)
202  {
203  DUNE_THROW(Dune::NotImplemented,"No coupling lfsu_cache available for ");
204  }
206 
209  void postAssembly(const GFSU& gfsu, const GFSV& gfsv)
210  {
211  Jacobian& jacobian = global_a_ss_view.container();
212  global_s_s_view.detach();
213  global_s_n_view.detach();
214  global_a_ss_view.detach();
215  global_a_sn_view.detach();
216  global_a_ns_view.detach();
217  global_a_nn_view.detach();
218 
219  if(local_assembler.doPostProcessing())
220  local_assembler.handle_dirichlet_constraints(gfsv,jacobian);
221  }
223 
226 
233  template<typename EG>
234  bool assembleCell(const EG & eg)
235  {
236  return LocalAssembler::isNonOverlapping && eg.entity().partitionType() != Dune::InteriorEntity;
237  }
238 
239  template<typename EG, typename LFSUC, typename LFSVC>
240  void assembleUVVolume(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache)
241  {
242  global_a_ss_view.setWeight(local_assembler.weight());
244  jacobian_volume(lop,eg,lfsu_cache.localFunctionSpace(),global_s_s_view,lfsv_cache.localFunctionSpace(),global_a_ss_view);
245  }
246 
247  template<typename IG, typename LFSUC, typename LFSVC>
248  void assembleUVSkeleton(const IG & ig, const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
249  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache)
250  {
251  global_a_ss_view.setWeight(local_assembler.weight());
252  global_a_sn_view.setWeight(local_assembler.weight());
253  global_a_ns_view.setWeight(local_assembler.weight());
254  global_a_nn_view.setWeight(local_assembler.weight());
257  lfsu_s_cache.localFunctionSpace(),global_s_s_view,lfsv_s_cache.localFunctionSpace(),
258  lfsu_n_cache.localFunctionSpace(),global_s_n_view,lfsv_n_cache.localFunctionSpace(),
259  global_a_ss_view, global_a_sn_view,
260  global_a_ns_view, global_a_nn_view);
261  }
262 
263  template<typename IG, typename LFSUC, typename LFSVC>
264  void assembleUVBoundary(const IG & ig, const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache)
265  {
266  global_a_ss_view.setWeight(local_assembler.weight());
268  jacobian_boundary(lop,ig,lfsu_s_cache.localFunctionSpace(),global_s_s_view,lfsv_s_cache.localFunctionSpace(),global_a_ss_view);
269  }
270 
271  template<typename IG, typename LFSUC, typename LFSVC>
272  static void assembleUVEnrichedCoupling(const IG & ig,
273  const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
274  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache,
275  const LFSUC & lfsu_coupling_cache, const LFSVC & lfsv_coupling_cache)
276  {
277  DUNE_THROW(Dune::NotImplemented,"Assembling of coupling spaces is not implemented for ");
278  }
279 
280  template<typename IG, typename LFSVC>
281  static void assembleVEnrichedCoupling(const IG & ig,
282  const LFSVC & lfsv_s_cache,
283  const LFSVC & lfsv_n_cache,
284  const LFSVC & lfsv_coupling_cache)
285  {
286  DUNE_THROW(Dune::NotImplemented,"Assembling of coupling spaces is not implemented for ");
287  }
288 
289  template<typename EG, typename LFSUC, typename LFSVC>
290  void assembleUVVolumePostSkeleton(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache)
291  {
292  global_a_ss_view.setWeight(local_assembler.weight());
294  jacobian_volume_post_skeleton(lop,eg,lfsu_cache.localFunctionSpace(),global_s_s_view,lfsv_cache.localFunctionSpace(),global_a_ss_view);
295  }
296 
298 
299  private:
302  const LocalAssembler & local_assembler;
303 
305  const LOP & lop;
306 
308  SolutionView global_s_s_view;
309  SolutionView global_s_n_view;
310 
312  JacobianView global_a_ss_view;
313  JacobianView global_a_sn_view;
314  JacobianView global_a_ns_view;
315  JacobianView global_a_nn_view;
316 
319  typedef Dune::PDELab::TrialSpaceTag LocalTrialSpaceTag;
320  typedef Dune::PDELab::TestSpaceTag LocalTestSpaceTag;
321 
323  typedef typename std::conditional<
324  std::is_base_of<
326  LOP
327  >::value,
330  >::type JacobianMatrix;
331 
333 
334  }; // End of class FastDGLocalJacobianAssemblerEngine
335 
336  }
337 }
338 #endif
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::FastDGLocalJacobianAssemblerEngine
FastDGLocalJacobianAssemblerEngine(const LocalAssembler &local_assembler_)
Constructor.
Definition: fastdg/jacobianengine.hh:66
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::LocalAssembler
LA LocalAssembler
The type of the wrapping local assembler.
Definition: fastdg/jacobianengine.hh:37
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::assembleUVSkeleton
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/jacobianengine.hh:248
gridoperatorutilities.hh
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::SolutionElement
Solution::ElementType SolutionElement
Definition: fastdg/jacobianengine.hh:57
Dune::PDELab::TrialSpaceTag
Definition: localfunctionspacetags.hh:48
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::assembleUVBoundary
void assembleUVBoundary(const IG &ig, const LFSUC &lfsu_s_cache, const LFSVC &lfsv_s_cache)
Definition: fastdg/jacobianengine.hh:264
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::requireSkeleton
bool requireSkeleton() const
Definition: fastdg/jacobianengine.hh:98
Dune::PDELab::DiagonalLocalMatrix
A dense matrix for storing data associated with the degrees of freedom of a pair of LocalFunctionSpac...
Definition: diagonallocalmatrix.hh:28
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::assembleVEnrichedCoupling
static void assembleVEnrichedCoupling(const IG &ig, const LFSVC &lfsv_s_cache, const LFSVC &lfsv_n_cache, const LFSVC &lfsv_coupling_cache)
Definition: fastdg/jacobianengine.hh:281
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::postAssembly
void postAssembly(const GFSU &gfsu, const GFSV &gfsv)
Definition: fastdg/jacobianengine.hh:209
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::onBindLFSUV
void onBindLFSUV(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/jacobianengine.hh:152
Dune::PDELab::FastDGLocalJacobianAssemblerEngine
The fast DG local assembler engine for DUNE grids which assembles the jacobian matrix.
Definition: fastdg/jacobianengine.hh:25
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::GFSU
LFSU::Traits::GridFunctionSpace GFSU
Definition: fastdg/jacobianengine.hh:45
Dune
For backward compatibility – Do not use this!
Definition: adaptivity.hh:28
localassemblerenginebase.hh
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::assembleUVVolume
void assembleUVVolume(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/jacobianengine.hh:240
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::requireUVVolume
bool requireUVVolume() const
Definition: fastdg/jacobianengine.hh:102
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::loadCoefficientsLFSUCoupling
void loadCoefficientsLFSUCoupling(const LFSUC &lfsu_c_cache)
Definition: fastdg/jacobianengine.hh:201
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::LFSU
LA::LFSU LFSU
The local function spaces.
Definition: fastdg/jacobianengine.hh:43
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::needsConstraintsCaching
bool needsConstraintsCaching(const TrialConstraintsContainer &cu, const TestConstraintsContainer &cv)
Definition: fastdg/jacobianengine.hh:31
value
static const unsigned int value
Definition: gridfunctionspace/tags.hh:139
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::localAssembler
const LocalAssembler & localAssembler() const
Public access to the wrapping local assembler.
Definition: fastdg/jacobianengine.hh:113
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::SolutionView
Solution::template ConstAliasedLocalView< LFSUCache > SolutionView
Definition: fastdg/jacobianengine.hh:58
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::JacobianView
Jacobian::template AliasedLocalView< LFSVCache, LFSUCache > JacobianView
Definition: fastdg/jacobianengine.hh:53
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::LOP
LA::LocalOperator LOP
The type of the local operator.
Definition: fastdg/jacobianengine.hh:40
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::LFSVCache
LA::LFSVCache LFSVCache
Definition: fastdg/jacobianengine.hh:47
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::requireUVSkeleton
bool requireUVSkeleton() const
Definition: fastdg/jacobianengine.hh:104
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::assembleCell
bool assembleCell(const EG &eg)
Definition: fastdg/jacobianengine.hh:234
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::LFSV
LA::LFSV LFSV
Definition: fastdg/jacobianengine.hh:46
localvector.hh
Dune::PDELab::TestSpaceTag
Definition: localfunctionspacetags.hh:54
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::requireUVVolumePostSkeleton
bool requireUVVolumePostSkeleton() const
Definition: fastdg/jacobianengine.hh:108
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::loadCoefficientsLFSUInside
void loadCoefficientsLFSUInside(const LFSUC &lfsu_cache)
Definition: fastdg/jacobianengine.hh:195
callswitch.hh
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::setSolution
void setSolution(const Solution &solution_)
Definition: fastdg/jacobianengine.hh:142
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::requireSkeletonTwoSided
bool requireSkeletonTwoSided() const
Definition: fastdg/jacobianengine.hh:100
Dune::PDELab::LocalMatrix
A dense matrix for storing data associated with the degrees of freedom of a pair of LocalFunctionSpac...
Definition: localmatrix.hh:183
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::trialConstraints
const LocalAssembler::Traits::TrialGridFunctionSpaceConstraints & trialConstraints() const
Trial space constraints.
Definition: fastdg/jacobianengine.hh:119
Dune::PDELab::LocalAssemblerCallSwitch
Impl::LocalAssemblerCallSwitchHelper< LOP, doIt > LocalAssemblerCallSwitch
Definition: callswitch.hh:344
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::GFSV
LFSV::Traits::GridFunctionSpace GFSV
Definition: fastdg/jacobianengine.hh:48
Dune::PDELab::LocalAssemblerEngineBase
Base class for LocalAssemblerEngine implementations to avoid boilerplate code.
Definition: localassemblerenginebase.hh:21
flags.hh
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::loadCoefficientsLFSUOutside
void loadCoefficientsLFSUOutside(const LFSUC &lfsu_n_cache)
Definition: fastdg/jacobianengine.hh:198
Dune::PDELab::FastDGLocalAssembler
The local assembler for DUNE grids.
Definition: fastdg/localassembler.hh:37
Dune::PDELab::LocalVector< SolutionElement, LocalTrialSpaceTag >
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::assembleUVVolumePostSkeleton
void assembleUVVolumePostSkeleton(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/jacobianengine.hh:290
constraints.hh
diagonallocalmatrix.hh
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::requireUVBoundary
bool requireUVBoundary() const
Definition: fastdg/jacobianengine.hh:106
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::testConstraints
const LocalAssembler::Traits::TestGridFunctionSpaceConstraints & testConstraints() const
Test space constraints.
Definition: fastdg/jacobianengine.hh:125
Dune::PDELab::lop::DiagonalJacobian
Decorator base class for local operators that have a diagonal jacobian matrix.
Definition: flags.hh:101
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::JacobianElement
Jacobian::ElementType JacobianElement
Definition: fastdg/jacobianengine.hh:52
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::Solution
LA::Traits::Solution Solution
The type of the solution vector.
Definition: fastdg/jacobianengine.hh:56
ig
const IG & ig
Definition: constraints.hh:149
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::assembleUVEnrichedCoupling
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/jacobianengine.hh:272
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::setJacobian
void setJacobian(Jacobian &jacobian_)
Definition: fastdg/jacobianengine.hh:132
localmatrix.hh
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::LFSUCache
LA::LFSUCache LFSUCache
Definition: fastdg/jacobianengine.hh:44
assemblerutilities.hh
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::onUnbindLFSUV
void onUnbindLFSUV(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/jacobianengine.hh:175
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::Jacobian
LA::Traits::Jacobian Jacobian
The type of the jacobian matrix.
Definition: fastdg/jacobianengine.hh:51
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::onUnbindLFSUVOutside
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/jacobianengine.hh:181
Dune::PDELab::FastDGLocalJacobianAssemblerEngine::onBindLFSUVOutside
void onBindLFSUVOutside(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/jacobianengine.hh:159