3 #ifndef DUNE_PDELAB_GRIDOPERATOR_FASTDG_LOCALASSEMBLER_HH
4 #define DUNE_PDELAB_GRIDOPERATOR_FASTDG_LOCALASSEMBLER_HH
8 #include <dune/common/shared_ptr.hh>
10 #include <dune/typetree/typetree.hh>
36 template<
typename GO,
typename LOP,
bool nonoverlapping_mode = false>
39 typename GO::Traits::TrialGridFunctionSpaceConstraints,
40 typename GO::Traits::TestGridFunctionSpaceConstraints>
48 typedef typename Traits::Residual::ElementType
RangeField;
88 : lop(stackobject_to_shared_ptr(lop_)),
90 doPreProcessing_(true),
91 doPostProcessing_(true),
92 pattern_engine(*this,border_dof_exchanger), residual_engine(*this), jacobian_engine(*this)
93 , jacobian_apply_engine(*this)
99 std::shared_ptr<typename GO::BorderDOFExchanger> border_dof_exchanger)
101 lop(stackobject_to_shared_ptr(lop_)),
103 doPreProcessing_(true),
104 doPostProcessing_(true),
105 pattern_engine(*this,border_dof_exchanger), residual_engine(*this), jacobian_engine(*this)
106 , jacobian_apply_engine(*this)
132 models<Impl::HasDoSkipEntity, LOP>(),
133 "Your local operator does not provide the 'doSkipEntity' flag. "
134 "If you are porting a previous implementation, set the flag to "
137 LOP, LOP::doSkipEntity>::skip_entity(*lop, eg);
149 static_assert(models<Impl::HasDoSkipIntersection, LOP>(),
150 "Your local operator does not provide the "
151 "'doSkipIntersection' flag. "
152 "If you are porting a previous implementation, set the "
155 LOP, LOP::doSkipIntersection>::skip_intersection(*lop,
ig);
180 void preStep (
Real time_,
Real dt_, std::size_t stages_){ lop->preStep(time_,dt_,stages_); }
188 return _reconstruct_border_entries;
200 return pattern_engine;
210 return residual_engine;
220 return jacobian_engine;
230 return jacobian_apply_engine;
241 return jacobian_apply_engine;
263 static constexpr
bool isLinear() {
return LOP::isLinear;}
278 doPreProcessing_ = v;
293 doPostProcessing_ = v;
296 template<
typename GFSV,
typename GC,
typename C>
299 typedef typename C::const_iterator global_row_iterator;
300 for (global_row_iterator cit = c.begin(); cit != c.end(); ++cit)
301 globalcontainer.clear_row_block(cit->first,1);
304 template<
typename GFSV,
typename GC>
309 template<
typename GFSV,
typename GC>
312 globalcontainer.flush();
314 globalcontainer.finalize();
320 std::shared_ptr<LOP> lop;
327 bool doPreProcessing_;
331 bool doPostProcessing_;
341 bool _reconstruct_border_entries;
const P & p
Definition: constraints.hh:148
const IG & ig
Definition: constraints.hh:149
For backward compatibility – Do not use this!
Definition: adaptivity.hh:28
Impl::LocalAssemblerCallSwitchHelper< LOP, doIt > LocalAssemblerCallSwitch
Definition: callswitch.hh:407
Definition: constraintstransformation.hh:112
Definition: lfsindexcache.hh:979
Create a local function space from a global function space.
Definition: localfunctionspace.hh:725
Definition: assemblerutilities.hh:51
GO::Traits::Range Residual
The type of the range (residual).
Definition: assemblerutilities.hh:88
MatrixBackend::template Pattern< Jacobian, TestGridFunctionSpace, TrialGridFunctionSpace > MatrixPattern
The matrix pattern.
Definition: assemblerutilities.hh:102
GO::Traits::TrialGridFunctionSpace TrialGridFunctionSpace
The trial grid function space.
Definition: assemblerutilities.hh:54
GO::Traits::Jacobian Jacobian
The type of the jacobian.
Definition: assemblerutilities.hh:95
GO::Traits::TestGridFunctionSpace TestGridFunctionSpace
The test grid function space.
Definition: assemblerutilities.hh:57
GO::Traits::Domain Solution
The type of the domain (solution).
Definition: assemblerutilities.hh:78
GO::Traits::Range Range
The type of the range (residual).
Definition: assemblerutilities.hh:85
GO::Traits::Domain Domain
The type of the domain (solution).
Definition: assemblerutilities.hh:75
GO::Traits::TrialGridFunctionSpaceConstraints TrialGridFunctionSpaceConstraints
The type of the trial grid function space constraints.
Definition: assemblerutilities.hh:61
GO::Traits::TestGridFunctionSpaceConstraints TestGridFunctionSpaceConstraints
The type of the test grid function space constraints.
Definition: assemblerutilities.hh:64
Base class for local assembler.
Definition: assemblerutilities.hh:217
const GO::Traits::TestGridFunctionSpaceConstraints * pconstraintsv
Definition: assemblerutilities.hh:688
void setResult(Range &result_)
Definition: fastdg/jacobianapplyengine.hh:125
void setSolution(const Domain &solution_)
Definition: fastdg/jacobianapplyengine.hh:107
void setUpdate(const Domain &update_)
Definition: fastdg/jacobianapplyengine.hh:117
void setSolution(const Solution &solution_)
Definition: fastdg/jacobianengine.hh:142
void setJacobian(Jacobian &jacobian_)
Definition: fastdg/jacobianengine.hh:132
The local assembler for DUNE grids.
Definition: fastdg/localassembler.hh:41
static constexpr bool doAlphaVolumePostSkeleton()
Definition: fastdg/localassembler.hh:256
static constexpr bool doPatternVolume()
Definition: fastdg/localassembler.hh:259
static constexpr bool doLambdaBoundary()
Definition: fastdg/localassembler.hh:255
Traits::TrialGridFunctionSpaceConstraints CU
Definition: fastdg/localassembler.hh:54
FastDGLocalAssembler(LOP &lop_, const CU &cu_, const CV &cv_, std::shared_ptr< typename GO::BorderDOFExchanger > border_dof_exchanger)
Constructor for non trivial constraints.
Definition: fastdg/localassembler.hh:98
FastDGLocalResidualAssemblerEngine< FastDGLocalAssembler > LocalResidualAssemblerEngine
Definition: fastdg/localassembler.hh:78
const LOP & localOperator() const
get a reference to the local operator
Definition: fastdg/localassembler.hh:117
static const bool isNonOverlapping
Definition: fastdg/localassembler.hh:63
void postStage()
Definition: fastdg/localassembler.hh:182
LFSIndexCache< LFSV, CV, true > LFSVCache
Definition: fastdg/localassembler.hh:71
bool reconstructBorderEntries() const
Definition: fastdg/localassembler.hh:186
LOP LocalOperator
The local operator.
Definition: fastdg/localassembler.hh:61
Traits::Residual::ElementType RangeField
The local operators type for real numbers e.g. time.
Definition: fastdg/localassembler.hh:48
bool skipEntity(const EG &eg) const
Definition: fastdg/localassembler.hh:129
Dune::PDELab::LocalFunctionSpace< GFSU, Dune::PDELab::TrialSpaceTag > LFSU
Definition: fastdg/localassembler.hh:68
static constexpr bool doLambdaVolume()
Definition: fastdg/localassembler.hh:251
void postStep()
Definition: fastdg/localassembler.hh:181
LOP & localOperator()
get a reference to the local operator
Definition: fastdg/localassembler.hh:111
static constexpr bool doPatternBoundary()
Definition: fastdg/localassembler.hh:261
Traits::TestGridFunctionSpaceConstraints CV
Definition: fastdg/localassembler.hh:55
void postProcessing(bool v)
Definition: fastdg/localassembler.hh:291
static constexpr bool doAlphaSkeleton()
Definition: fastdg/localassembler.hh:252
Dune::PDELab::LocalFunctionSpace< GFSV, Dune::PDELab::TestSpaceTag > LFSV
Definition: fastdg/localassembler.hh:69
static constexpr bool doAlphaVolume()
Query methods for the assembler engines. Theses methods do not belong to the assembler interface,...
Definition: fastdg/localassembler.hh:250
static constexpr bool doLambdaVolumePostSkeleton()
Definition: fastdg/localassembler.hh:257
RangeField Real
Definition: fastdg/localassembler.hh:49
void preStage(Real time_, int r_)
Definition: fastdg/localassembler.hh:179
bool skipIntersection(const IG &ig) const
Definition: fastdg/localassembler.hh:147
static constexpr bool doAlphaBoundary()
Definition: fastdg/localassembler.hh:254
void preProcessing(bool v)
Definition: fastdg/localassembler.hh:276
FastDGLocalAssembler(LOP &lop_, std::shared_ptr< typename GO::BorderDOFExchanger > border_dof_exchanger)
Constructor with empty constraints.
Definition: fastdg/localassembler.hh:87
LocalJacobianApplyAssemblerEngine & localJacobianApplyAssemblerEngine(const typename Traits::Domain &update, typename Traits::Range &result)
Definition: fastdg/localassembler.hh:226
static constexpr bool doLambdaSkeleton()
Definition: fastdg/localassembler.hh:253
void setTime(Real time_)
Definition: fastdg/localassembler.hh:161
bool doPreProcessing() const
Query whether to do preprocessing in the engines.
Definition: fastdg/localassembler.hh:270
LocalResidualAssemblerEngine & localResidualAssemblerEngine(typename Traits::Residual &r, const typename Traits::Solution &x)
Definition: fastdg/localassembler.hh:206
LFSIndexCache< LFSU, CU, true > LFSUCache
Definition: fastdg/localassembler.hh:70
bool doPostProcessing() const
Query whether to do postprocessing in the engines.
Definition: fastdg/localassembler.hh:285
Traits::TestGridFunctionSpace GFSV
Definition: fastdg/localassembler.hh:52
Traits::TrialGridFunctionSpace GFSU
Definition: fastdg/localassembler.hh:51
LocalJacobianAssemblerEngine & localJacobianAssemblerEngine(typename Traits::Jacobian &a, const typename Traits::Solution &x)
Definition: fastdg/localassembler.hh:216
LocalPatternAssemblerEngine & localPatternAssemblerEngine(typename Traits::MatrixPattern &p)
Definition: fastdg/localassembler.hh:197
static constexpr bool doPatternSkeleton()
Definition: fastdg/localassembler.hh:260
RangeField weight() const
Obtain the weight that was set last.
Definition: fastdg/localassembler.hh:167
void preStep(Real time_, Real dt_, std::size_t stages_)
Definition: fastdg/localassembler.hh:180
Dune::PDELab::LocalAssemblerBase< typename Traits::MatrixBackend, CU, CV > Base
The base class of this local assembler.
Definition: fastdg/localassembler.hh:58
void set_trivial_rows(const GFSV &gfsv, GC &globalcontainer, const C &c) const
Definition: fastdg/localassembler.hh:297
FastDGLocalPatternAssemblerEngine< FastDGLocalAssembler > LocalPatternAssemblerEngine
Definition: fastdg/localassembler.hh:77
FastDGLocalJacobianAssemblerEngine< FastDGLocalAssembler > LocalJacobianAssemblerEngine
Definition: fastdg/localassembler.hh:79
Real suggestTimestep(Real dt) const
Definition: fastdg/localassembler.hh:183
void setWeight(RangeField weight)
Notifies the assembler about the current weight of assembling.
Definition: fastdg/localassembler.hh:173
static constexpr bool doPatternVolumePostSkeleton()
Definition: fastdg/localassembler.hh:262
Dune::PDELab::LocalAssemblerTraits< GO > Traits
The traits class.
Definition: fastdg/localassembler.hh:45
static constexpr bool doSkeletonTwoSided()
Definition: fastdg/localassembler.hh:258
void set_trivial_rows(const GFSV &gfsv, GC &globalcontainer, const EmptyTransformation &c) const
Definition: fastdg/localassembler.hh:305
FastDGLocalJacobianApplyAssemblerEngine< FastDGLocalAssembler > LocalJacobianApplyAssemblerEngine
Definition: fastdg/localassembler.hh:80
LocalJacobianApplyAssemblerEngine & localJacobianApplyAssemblerEngine(const typename Traits::Domain &solution, const typename Traits::Domain &update, typename Traits::Range &result)
Definition: fastdg/localassembler.hh:236
static constexpr bool isLinear()
Definition: fastdg/localassembler.hh:263
void handle_dirichlet_constraints(const GFSV &gfsv, GC &globalcontainer) const
Definition: fastdg/localassembler.hh:310
void setPattern(Pattern &pattern_)
Definition: fastdg/patternengine.hh:92
void setSolution(const Solution &solution_)
Definition: fastdg/residualengine.hh:148
void setResidual(Residual &residual_)
Definition: fastdg/residualengine.hh:140