dune-pdelab  2.4-dev
default/jacobianengine.hh
Go to the documentation of this file.
1 #ifndef DUNE_PDELAB_DEFAULT_JACOBIANENGINE_HH
2 #define DUNE_PDELAB_DEFAULT_JACOBIANENGINE_HH
3 
12 
13 namespace Dune{
14  namespace PDELab{
15 
23  template<typename LA>
26  {
27  public:
28 
29  template<typename TrialConstraintsContainer, typename TestConstraintsContainer>
30  bool needsConstraintsCaching(const TrialConstraintsContainer& cu, const TestConstraintsContainer& cv)
31  {
32  return cu.containsNonDirichletConstraints() || cv.containsNonDirichletConstraints();
33  }
34 
36  typedef LA LocalAssembler;
37 
39  typedef typename LA::LocalOperator LOP;
40 
42  typedef typename LA::LFSU LFSU;
43  typedef typename LA::LFSUCache LFSUCache;
44  typedef typename LFSU::Traits::GridFunctionSpace GFSU;
45  typedef typename LA::LFSV LFSV;
46  typedef typename LA::LFSVCache LFSVCache;
47  typedef typename LFSV::Traits::GridFunctionSpace GFSV;
48 
50  typedef typename LA::Traits::Jacobian Jacobian;
51  typedef typename Jacobian::ElementType JacobianElement;
52  typedef typename Jacobian::template LocalView<LFSVCache,LFSUCache> JacobianView;
53 
55  typedef typename LA::Traits::Solution Solution;
56  typedef typename Solution::ElementType SolutionElement;
57  typedef typename Solution::template ConstLocalView<LFSUCache> SolutionView;
58 
65  DefaultLocalJacobianAssemblerEngine(const LocalAssembler & local_assembler_)
66  : local_assembler(local_assembler_), lop(local_assembler_.lop),
67  al_view(al,1.0),
68  al_sn_view(al_sn,1.0),
69  al_ns_view(al_ns,1.0),
70  al_nn_view(al_nn,1.0)
71  {}
72 
75  bool requireSkeleton() const
76  { return local_assembler.doAlphaSkeleton(); }
78  { return local_assembler.doSkeletonTwoSided(); }
79  bool requireUVVolume() const
80  { return local_assembler.doAlphaVolume(); }
81  bool requireUVSkeleton() const
82  { return local_assembler.doAlphaSkeleton(); }
83  bool requireUVBoundary() const
84  { return local_assembler.doAlphaBoundary(); }
86  { return local_assembler.doAlphaVolumePostSkeleton(); }
88 
90  const LocalAssembler & localAssembler() const { return local_assembler; }
91 
93  const typename LocalAssembler::Traits::TrialGridFunctionSpaceConstraints& trialConstraints() const
94  {
95  return localAssembler().trialConstraints();
96  }
97 
99  const typename LocalAssembler::Traits::TestGridFunctionSpaceConstraints& testConstraints() const
100  {
101  return localAssembler().testConstraints();
102  }
103 
106  void setJacobian(Jacobian & jacobian_){
107  global_a_ss_view.attach(jacobian_);
108  global_a_sn_view.attach(jacobian_);
109  global_a_ns_view.attach(jacobian_);
110  global_a_nn_view.attach(jacobian_);
111  }
112 
115  void setSolution(const Solution & solution_){
116  global_s_s_view.attach(solution_);
117  global_s_n_view.attach(solution_);
118  }
119 
123  template<typename EG, typename LFSUC, typename LFSVC>
124  void onBindLFSUV(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache){
125  global_s_s_view.bind(lfsu_cache);
126  xl.resize(lfsu_cache.size());
127  global_a_ss_view.bind(lfsv_cache,lfsu_cache);
128  al.assign(lfsv_cache.size(),lfsu_cache.size(),0.0);
129  }
130 
131  template<typename IG, typename LFSUC, typename LFSVC>
132  void onBindLFSUVOutside(const IG & ig,
133  const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
134  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache)
135  {
136  global_s_n_view.bind(lfsu_n_cache);
137  xn.resize(lfsu_n_cache.size());
138  global_a_sn_view.bind(lfsv_s_cache,lfsu_n_cache);
139  al_sn.assign(lfsv_s_cache.size(),lfsu_n_cache.size(),0.0);
140  global_a_ns_view.bind(lfsv_n_cache,lfsu_s_cache);
141  al_ns.assign(lfsv_n_cache.size(),lfsu_s_cache.size(),0.0);
142  global_a_nn_view.bind(lfsv_n_cache,lfsu_n_cache);
143  al_nn.assign(lfsv_n_cache.size(),lfsu_n_cache.size(),0.0);
144  }
145 
147 
151  template<typename EG, typename LFSUC, typename LFSVC>
152  void onUnbindLFSUV(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache){
153  local_assembler.scatter_jacobian(al,global_a_ss_view,false);
154  }
155 
156  template<typename IG, typename LFSUC, typename LFSVC>
157  void onUnbindLFSUVOutside(const IG & ig,
158  const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
159  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache)
160  {
161  local_assembler.scatter_jacobian(al_sn,global_a_sn_view,false);
162  local_assembler.scatter_jacobian(al_ns,global_a_ns_view,false);
163  local_assembler.scatter_jacobian(al_nn,global_a_nn_view,false);
164  }
165 
167 
170  template<typename LFSUC>
171  void loadCoefficientsLFSUInside(const LFSUC & lfsu_cache){
172  global_s_s_view.read(xl);
173  }
174  template<typename LFSUC>
175  void loadCoefficientsLFSUOutside(const LFSUC & lfsu_n_cache){
176  global_s_n_view.read(xn);
177  }
178  template<typename LFSUC>
179  void loadCoefficientsLFSUCoupling(const LFSUC & lfsu_c_cache)
180  {DUNE_THROW(Dune::NotImplemented,"No coupling lfsu_cache available for ");}
182 
185  void postAssembly(const GFSU& gfsu, const GFSV& gfsv){
186  Jacobian& jacobian = global_a_ss_view.container();
187  global_s_s_view.detach();
188  global_s_n_view.detach();
189  global_a_ss_view.detach();
190  global_a_sn_view.detach();
191  global_a_ns_view.detach();
192  global_a_nn_view.detach();
193 
194  if(local_assembler.doPostProcessing){
195  local_assembler.handle_dirichlet_constraints(gfsv,jacobian);
196  }
197  }
199 
202 
209  template<typename EG>
210  bool assembleCell(const EG & eg)
211  {
212  return LocalAssembler::isNonOverlapping && eg.entity().partitionType() != Dune::InteriorEntity;
213  }
214 
215  template<typename EG, typename LFSUC, typename LFSVC>
216  void assembleUVVolume(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache)
217  {
218  al_view.setWeight(local_assembler.weight);
220  jacobian_volume(lop,eg,lfsu_cache.localFunctionSpace(),xl,lfsv_cache.localFunctionSpace(),al_view);
221  }
222 
223  template<typename IG, typename LFSUC, typename LFSVC>
224  void assembleUVSkeleton(const IG & ig, const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
225  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache)
226  {
227  al_view.setWeight(local_assembler.weight);
228  al_sn_view.setWeight(local_assembler.weight);
229  al_ns_view.setWeight(local_assembler.weight);
230  al_nn_view.setWeight(local_assembler.weight);
231 
233  jacobian_skeleton(lop,ig,lfsu_s_cache.localFunctionSpace(),xl,lfsv_s_cache.localFunctionSpace(),lfsu_n_cache.localFunctionSpace(),xn,lfsv_n_cache.localFunctionSpace(),al_view,al_sn_view,al_ns_view,al_nn_view);
234  }
235 
236  template<typename IG, typename LFSUC, typename LFSVC>
237  void assembleUVBoundary(const IG & ig, const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache)
238  {
239  al_view.setWeight(local_assembler.weight);
241  jacobian_boundary(lop,ig,lfsu_s_cache.localFunctionSpace(),xl,lfsv_s_cache.localFunctionSpace(),al_view);
242  }
243 
244  template<typename IG, typename LFSUC, typename LFSVC>
245  static void assembleUVEnrichedCoupling(const IG & ig,
246  const LFSUC & lfsu_s_cache, const LFSVC & lfsv_s_cache,
247  const LFSUC & lfsu_n_cache, const LFSVC & lfsv_n_cache,
248  const LFSUC & lfsu_coupling_cache, const LFSVC & lfsv_coupling_cache)
249  {DUNE_THROW(Dune::NotImplemented,"Assembling of coupling spaces is not implemented for ");}
250 
251  template<typename IG, typename LFSVC>
252  static void assembleVEnrichedCoupling(const IG & ig,
253  const LFSVC & lfsv_s_cache,
254  const LFSVC & lfsv_n_cache,
255  const LFSVC & lfsv_coupling_cache)
256  {DUNE_THROW(Dune::NotImplemented,"Assembling of coupling spaces is not implemented for ");}
257 
258  template<typename EG, typename LFSUC, typename LFSVC>
259  void assembleUVVolumePostSkeleton(const EG & eg, const LFSUC & lfsu_cache, const LFSVC & lfsv_cache)
260  {
261  al_view.setWeight(local_assembler.weight);
263  jacobian_volume_post_skeleton(lop,eg,lfsu_cache.localFunctionSpace(),xl,lfsv_cache.localFunctionSpace(),al_view);
264  }
265 
267 
268  private:
271  const LocalAssembler & local_assembler;
272 
274  const LOP & lop;
275 
277  SolutionView global_s_s_view;
278  SolutionView global_s_n_view;
279 
281  JacobianView global_a_ss_view;
282  JacobianView global_a_sn_view;
283  JacobianView global_a_ns_view;
284  JacobianView global_a_nn_view;
285 
288  typedef Dune::PDELab::TrialSpaceTag LocalTrialSpaceTag;
289  typedef Dune::PDELab::TestSpaceTag LocalTestSpaceTag;
290 
292  typedef typename std::conditional<
293  std::is_base_of<
295  LOP
296  >::value,
299  >::type JacobianMatrix;
300 
301  SolutionVector xl;
302  SolutionVector xn;
303 
304  JacobianMatrix al;
305  JacobianMatrix al_sn;
306  JacobianMatrix al_ns;
307  JacobianMatrix al_nn;
308 
309  typename JacobianMatrix::WeightedAccumulationView al_view;
310  typename JacobianMatrix::WeightedAccumulationView al_sn_view;
311  typename JacobianMatrix::WeightedAccumulationView al_ns_view;
312  typename JacobianMatrix::WeightedAccumulationView al_nn_view;
313 
315 
316  }; // End of class DefaultLocalJacobianAssemblerEngine
317 
318  }
319 }
320 #endif
bool assembleCell(const EG &eg)
Definition: default/jacobianengine.hh:210
LFSV::Traits::GridFunctionSpace GFSV
Definition: default/jacobianengine.hh:47
void onBindLFSUV(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: default/jacobianengine.hh:124
void resize(size_type size)
Resize the container.
Definition: localvector.hh:251
bool requireUVBoundary() const
Definition: default/jacobianengine.hh:83
static void assembleVEnrichedCoupling(const IG &ig, const LFSVC &lfsv_s_cache, const LFSVC &lfsv_n_cache, const LFSVC &lfsv_coupling_cache)
Definition: default/jacobianengine.hh:252
A dense matrix for storing data associated with the degrees of freedom of a pair of LocalFunctionSpac...
Definition: diagonallocalmatrix.hh:28
The local assembler engine for DUNE grids which assembles the jacobian matrix.
Definition: default/jacobianengine.hh:24
A dense matrix for storing data associated with the degrees of freedom of a pair of LocalFunctionSpac...
Definition: localmatrix.hh:184
static void jacobian_volume_post_skeleton(const LA &la, const EG &eg, const LFSU &lfsu, const X &x, const LFSV &lfsv, M &mat)
Definition: callswitch.hh:105
Base class for LocalAssemblerEngine implementations to avoid boilerplate code.
Definition: localassemblerenginebase.hh:21
void setJacobian(Jacobian &jacobian_)
Definition: default/jacobianengine.hh:106
void loadCoefficientsLFSUOutside(const LFSUC &lfsu_n_cache)
Definition: default/jacobianengine.hh:175
void loadCoefficientsLFSUCoupling(const LFSUC &lfsu_c_cache)
Definition: default/jacobianengine.hh:179
LA::LFSV LFSV
Definition: default/jacobianengine.hh:45
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: default/jacobianengine.hh:224
Jacobian::template LocalView< LFSVCache, LFSUCache > JacobianView
Definition: default/jacobianengine.hh:52
const IG & ig
Definition: constraints.hh:147
bool needsConstraintsCaching(const TrialConstraintsContainer &cu, const TestConstraintsContainer &cv)
Definition: default/jacobianengine.hh:30
static const unsigned int value
Definition: gridfunctionspace/tags.hh:175
LA LocalAssembler
The type of the wrapping local assembler.
Definition: default/jacobianengine.hh:36
void setSolution(const Solution &solution_)
Definition: default/jacobianengine.hh:115
static void jacobian_boundary(const LA &la, const IG &ig, const LFSU &lfsu_s, const X &x_s, const LFSV &lfsv_s, M &mat_ss)
Definition: callswitch.hh:117
Decorator base class for local operators that have a diagonal jacobian matrix.
Definition: flags.hh:98
Definition: localfunctionspacetags.hh:48
void postAssembly(const GFSU &gfsu, const GFSV &gfsv)
Definition: default/jacobianengine.hh:185
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: default/jacobianengine.hh:157
Definition: localfunctionspacetags.hh:54
LA::LFSU LFSU
The local function spaces.
Definition: default/jacobianengine.hh:42
LA::Traits::Jacobian Jacobian
The type of the jacobian matrix.
Definition: default/jacobianengine.hh:50
Jacobian::ElementType JacobianElement
Definition: default/jacobianengine.hh:51
LFSU::Traits::GridFunctionSpace GFSU
Definition: default/jacobianengine.hh:44
static void jacobian_skeleton(const LA &la, const IG &ig, const LFSU &lfsu_s, const X &x_s, const LFSV &lfsv_s, const LFSU &lfsu_n, const X &x_n, const LFSV &lfsv_n, M &mat_ss, M &mat_sn, M &mat_ns, M &mat_nn)
Definition: callswitch.hh:109
Solution::template ConstLocalView< LFSUCache > SolutionView
Definition: default/jacobianengine.hh:57
bool requireSkeleton() const
Definition: default/jacobianengine.hh:75
LA::LFSVCache LFSVCache
Definition: default/jacobianengine.hh:46
const LocalAssembler::Traits::TestGridFunctionSpaceConstraints & testConstraints() const
Test space constraints.
Definition: default/jacobianengine.hh:99
Definition: adaptivity.hh:27
void assembleUVVolumePostSkeleton(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: default/jacobianengine.hh:259
static void jacobian_volume(const LA &la, const EG &eg, const LFSU &lfsu, const X &x, const LFSV &lfsv, M &mat)
Definition: callswitch.hh:101
bool requireUVSkeleton() const
Definition: default/jacobianengine.hh:81
void onUnbindLFSUV(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: default/jacobianengine.hh:152
const LocalAssembler & localAssembler() const
Public access to the wrapping local assembler.
Definition: default/jacobianengine.hh:90
bool requireUVVolume() const
Definition: default/jacobianengine.hh:79
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: default/jacobianengine.hh:132
void assembleUVBoundary(const IG &ig, const LFSUC &lfsu_s_cache, const LFSVC &lfsv_s_cache)
Definition: default/jacobianengine.hh:237
LA::Traits::Solution Solution
The type of the solution vector.
Definition: default/jacobianengine.hh:55
Solution::ElementType SolutionElement
Definition: default/jacobianengine.hh:56
LA::LocalOperator LOP
The type of the local operator.
Definition: default/jacobianengine.hh:39
bool requireUVVolumePostSkeleton() const
Definition: default/jacobianengine.hh:85
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: default/jacobianengine.hh:245
void loadCoefficientsLFSUInside(const LFSUC &lfsu_cache)
Definition: default/jacobianengine.hh:171
bool requireSkeletonTwoSided() const
Definition: default/jacobianengine.hh:77
LA::LFSUCache LFSUCache
Definition: default/jacobianengine.hh:43
const LocalAssembler::Traits::TrialGridFunctionSpaceConstraints & trialConstraints() const
Trial space constraints.
Definition: default/jacobianengine.hh:93
const EG & eg
Definition: constraints.hh:280
void assembleUVVolume(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: default/jacobianengine.hh:216
DefaultLocalJacobianAssemblerEngine(const LocalAssembler &local_assembler_)
Constructor.
Definition: default/jacobianengine.hh:65