1 #ifndef DUNE_PDELAB_GRIDOPERATOR_FASTDG_NONLINEARJACOBIANAPPLYENGINE_HH 2 #define DUNE_PDELAB_GRIDOPERATOR_FASTDG_NONLINEARJACOBIANAPPLYENGINE_HH 26 template<
typename TrialConstra
intsContainer,
typename TestConstra
intsContainer>
36 typedef typename LA::LocalOperator
LOP;
39 typedef typename LA::Traits::Residual
Residual;
43 typedef typename LA::Traits::Solution
Solution;
47 typedef typename LA::LFSU
LFSU;
49 typedef typename LFSU::Traits::GridFunctionSpace
GFSU;
50 typedef typename LA::LFSV
LFSV;
52 typedef typename LFSV::Traits::GridFunctionSpace
GFSV;
54 typedef typename Solution::template ConstAliasedLocalView<LFSUCache>
SolutionView;
55 typedef typename Residual::template AliasedLocalView<LFSVCache>
ResidualView;
64 : local_assembler(local_assembler_),
65 lop(local_assembler_.localOperator())
71 {
return local_assembler.doAlphaSkeleton(); }
73 {
return local_assembler.doSkeletonTwoSided(); }
75 {
return local_assembler.doAlphaVolume(); }
77 {
return local_assembler.doAlphaSkeleton(); }
79 {
return local_assembler.doAlphaBoundary(); }
81 {
return local_assembler.doAlphaVolumePostSkeleton(); }
87 return local_assembler;
91 const typename LocalAssembler::Traits::TrialGridFunctionSpaceConstraints&
trialConstraints()
const 97 const typename LocalAssembler::Traits::TestGridFunctionSpaceConstraints&
testConstraints()
const 106 global_rl_view.attach(residual_);
107 global_rn_view.attach(residual_);
114 global_sl_view.attach(solution_);
115 global_sn_view.attach(solution_);
122 global_zl_view.attach(update_);
123 global_zn_view.attach(update_);
129 template<
typename EG,
typename LFSUC,
typename LFSVC>
130 void onBindLFSUV(
const EG & eg,
const LFSUC & lfsu_cache,
const LFSVC & lfsv_cache)
132 global_sl_view.bind(lfsu_cache);
133 global_zl_view.bind(lfsu_cache);
136 template<
typename EG,
typename LFSVC>
139 global_rl_view.bind(lfsv_cache);
142 template<
typename IG,
typename LFSUC,
typename LFSVC>
145 global_sl_view.bind(lfsu_cache);
146 global_zl_view.bind(lfsu_cache);
149 template<
typename IG,
typename LFSUC,
typename LFSVC>
151 const LFSUC & lfsu_s_cache,
const LFSVC & lfsv_s_cache,
152 const LFSUC & lfsu_n_cache,
const LFSVC & lfsv_n_cache)
154 global_sn_view.bind(lfsu_n_cache);
155 global_zn_view.bind(lfsu_n_cache);
158 template<
typename IG,
typename LFSVC>
161 global_rl_view.bind(lfsv_cache);
164 template<
typename IG,
typename LFSVC>
166 const LFSVC & lfsv_s_cache,
167 const LFSVC & lfsv_n_cache)
169 global_rn_view.bind(lfsv_n_cache);
177 template<
typename EG,
typename LFSVC>
180 global_rl_view.commit();
181 global_rl_view.unbind();
184 template<
typename IG,
typename LFSVC>
187 global_rl_view.commit();
188 global_rl_view.unbind();
191 template<
typename IG,
typename LFSVC>
193 const LFSVC & lfsv_s_cache,
194 const LFSVC & lfsv_n_cache)
196 global_rn_view.commit();
197 global_rn_view.unbind();
203 template<
typename LFSUC>
207 template<
typename LFSUC>
211 template<
typename LFSUC>
214 DUNE_THROW(Dune::NotImplemented,
"No coupling lfsu available for ");
223 if(local_assembler.doPostProcessing())
225 global_rl_view.container());
240 template<
typename EG>
243 return LocalAssembler::isNonOverlapping && eg.entity().partitionType() != Dune::InteriorEntity;
246 template<
typename EG,
typename LFSUC,
typename LFSVC>
249 global_rl_view.setWeight(local_assembler.weight());
254 template<
typename IG,
typename LFSUC,
typename LFSVC>
256 const LFSUC & lfsu_n_cache,
const LFSVC & lfsv_n_cache)
258 global_rl_view.setWeight(local_assembler.weight());
259 global_rn_view.setWeight(local_assembler.weight());
262 lfsu_s_cache.localFunctionSpace(),global_sl_view,global_zl_view,lfsv_s_cache.localFunctionSpace(),
263 lfsu_n_cache.localFunctionSpace(),global_sn_view,global_zn_view,lfsv_n_cache.localFunctionSpace(),
264 global_rl_view,global_rn_view);
267 template<
typename IG,
typename LFSUC,
typename LFSVC>
270 global_rl_view.setWeight(local_assembler.weight());
275 template<
typename IG,
typename LFSUC,
typename LFSVC>
277 const LFSUC & lfsu_s_cache,
const LFSVC & lfsv_s_cache,
278 const LFSUC & lfsu_n_cache,
const LFSVC & lfsv_n_cache,
279 const LFSUC & lfsu_coupling_cache,
const LFSVC & lfsv_coupling_cache)
280 {DUNE_THROW(Dune::NotImplemented,
"Assembling of coupling spaces is not implemented for ");}
282 template<
typename EG,
typename LFSUC,
typename LFSVC>
285 global_rl_view.setWeight(local_assembler.weight());
295 const LocalAssembler & local_assembler;
301 ResidualView global_rl_view;
302 ResidualView global_rn_view;
305 SolutionView global_sl_view;
306 SolutionView global_sn_view;
309 SolutionView global_zl_view;
310 SolutionView global_zn_view;
const LocalAssembler & localAssembler() const
Public access to the wrapping local assembler.
Definition: fastdg/nonlinearjacobianapplyengine.hh:85
void assembleUVVolumePostSkeleton(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:283
static void nonlinear_jacobian_apply_volume_post_skeleton(const LA &la, const EG &eg, const LFSU &lfsu, const X &x, const X &z, const LFSV &lfsv, Y &y)
Definition: callswitch.hh:159
void onBindLFSV(const EG &eg, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:137
const LocalAssembler::Traits::TestGridFunctionSpaceConstraints & testConstraints() const
Test space constraints.
Definition: fastdg/nonlinearjacobianapplyengine.hh:97
Residual::template AliasedLocalView< LFSVCache > ResidualView
Definition: fastdg/nonlinearjacobianapplyengine.hh:55
void loadCoefficientsLFSUInside(const LFSUC &lfsu_s_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:204
bool requireSkeletonTwoSided() const
Definition: fastdg/nonlinearjacobianapplyengine.hh:72
void setUpdate(const Solution &update_)
Definition: fastdg/nonlinearjacobianapplyengine.hh:120
void constrain_residual(const CG &cg, XG &xg)
transform residual into transformed basis: r -> r~
Definition: constraints.hh:906
LA::LocalOperator LOP
The type of the local operator.
Definition: fastdg/nonlinearjacobianapplyengine.hh:36
Solution::ElementType SolutionElement
Definition: fastdg/nonlinearjacobianapplyengine.hh:44
void onBindLFSUV(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:130
LA::LFSVCache LFSVCache
Definition: fastdg/nonlinearjacobianapplyengine.hh:51
Definition: localfunctionspacetags.hh:54
void onBindLFSVInside(const IG &ig, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:159
bool requireUVSkeleton() const
Definition: fastdg/nonlinearjacobianapplyengine.hh:76
void setResidual(Residual &residual_)
Definition: fastdg/nonlinearjacobianapplyengine.hh:104
void postAssembly(const GFSU &gfsu, const GFSV &gfsv)
Definition: fastdg/nonlinearjacobianapplyengine.hh:221
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/nonlinearjacobianapplyengine.hh:276
bool requireUVVolume() const
Definition: fastdg/nonlinearjacobianapplyengine.hh:74
bool requireUVBoundary() const
Definition: fastdg/nonlinearjacobianapplyengine.hh:78
Base class for LocalAssemblerEngine implementations to avoid boilerplate code.
Definition: localassemblerenginebase.hh:21
static void nonlinear_jacobian_apply_boundary(const LA &la, const IG &ig, const LFSU &lfsu_s, const X &x_s, const X &z_s, const LFSV &lfsv_s, Y &y_s)
Definition: callswitch.hh:170
void onUnbindLFSVInside(const IG &ig, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:185
LA::Traits::Solution Solution
The type of the solution vector.
Definition: fastdg/nonlinearjacobianapplyengine.hh:43
void onBindLFSUVInside(const IG &ig, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:143
void setSolution(const Solution &solution_)
Definition: fastdg/nonlinearjacobianapplyengine.hh:112
static void nonlinear_jacobian_apply_volume(const LA &la, const EG &eg, const LFSU &lfsu, const X &x, const X &z, const LFSV &lfsv, Y &y)
Definition: callswitch.hh:155
void loadCoefficientsLFSUOutside(const LFSUC &lfsu_n_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:208
LFSU::Traits::GridFunctionSpace GFSU
Definition: fastdg/nonlinearjacobianapplyengine.hh:49
For backward compatibility – Do not use this!
Definition: adaptivity.hh:28
Solution::template ConstAliasedLocalView< LFSUCache > SolutionView
Definition: fastdg/nonlinearjacobianapplyengine.hh:54
void assembleUVBoundary(const IG &ig, const LFSUC &lfsu_s_cache, const LFSVC &lfsv_s_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:268
LFSV::Traits::GridFunctionSpace GFSV
Definition: fastdg/nonlinearjacobianapplyengine.hh:52
const LocalAssembler::Traits::TrialGridFunctionSpaceConstraints & trialConstraints() const
Trial space constraints.
Definition: fastdg/nonlinearjacobianapplyengine.hh:91
void onUnbindLFSV(const EG &eg, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:178
const IG & ig
Definition: constraints.hh:148
bool assembleCell(const EG &eg)
Definition: fastdg/nonlinearjacobianapplyengine.hh:241
LA::Traits::Residual Residual
The type of the residual vector.
Definition: fastdg/nonlinearjacobianapplyengine.hh:39
bool requireUVVolumePostSkeleton() const
Definition: fastdg/nonlinearjacobianapplyengine.hh:80
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/nonlinearjacobianapplyengine.hh:255
bool needsConstraintsCaching(const TrialConstraintsContainer &cu, const TestConstraintsContainer &cv) const
Definition: fastdg/nonlinearjacobianapplyengine.hh:27
LA::LFSUCache LFSUCache
Definition: fastdg/nonlinearjacobianapplyengine.hh:48
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/nonlinearjacobianapplyengine.hh:150
LA LocalAssembler
The type of the wrapping local assembler.
Definition: fastdg/nonlinearjacobianapplyengine.hh:33
Residual::ElementType ResidualElement
Definition: fastdg/nonlinearjacobianapplyengine.hh:40
void loadCoefficientsLFSUCoupling(const LFSUC &lfsu_c_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:212
bool requireSkeleton() const
Definition: fastdg/nonlinearjacobianapplyengine.hh:70
LA::LFSU LFSU
The local function spaces.
Definition: fastdg/nonlinearjacobianapplyengine.hh:47
void onUnbindLFSVOutside(const IG &ig, const LFSVC &lfsv_s_cache, const LFSVC &lfsv_n_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:192
void assembleUVVolume(const EG &eg, const LFSUC &lfsu_cache, const LFSVC &lfsv_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:247
void onBindLFSVOutside(const IG &ig, const LFSVC &lfsv_s_cache, const LFSVC &lfsv_n_cache)
Definition: fastdg/nonlinearjacobianapplyengine.hh:165
FastDGLocalNonlinearJacobianApplyAssemblerEngine(const LocalAssembler &local_assembler_)
Constructor.
Definition: fastdg/nonlinearjacobianapplyengine.hh:63
Definition: localfunctionspacetags.hh:48
LA::LFSV LFSV
Definition: fastdg/nonlinearjacobianapplyengine.hh:50
The fast DG local assembler engine for DUNE grids which assembles the local application of the Jacobi...
Definition: fastdg/nonlinearjacobianapplyengine.hh:21
static void nonlinear_jacobian_apply_skeleton(const LA &la, const IG &ig, const LFSU &lfsu_s, const X &x_s, const X &z_s, const LFSV &lfsv_s, const LFSU &lfsu_n, const X &x_n, const X &z_n, const LFSV &lfsv_n, Y &y_s, Y &y_n)
Definition: callswitch.hh:163