11 #include <dune/common/fvector.hh> 12 #include <dune/common/fmatrix.hh> 13 #include <dune/common/shared_ptr.hh> 15 #include <dune/grid/common/grid.hh> 17 #include <amdis/AdaptInfo.hpp> 18 #include <amdis/AdaptiveGrid.hpp> 19 #include <amdis/BiLinearForm.hpp> 20 #include <amdis/CreatorInterface.hpp> 21 #include <amdis/CreatorMap.hpp> 22 #include <amdis/DirichletBC.hpp> 23 #include <amdis/DOFVector.hpp> 25 #include <amdis/Flag.hpp> 26 #include <amdis/Initfile.hpp> 27 #include <amdis/LinearAlgebra.hpp> 28 #include <amdis/LinearForm.hpp> 29 #include <amdis/Marker.hpp> 30 #include <amdis/MeshCreator.hpp> 31 #include <amdis/OperatorList.hpp> 32 #include <amdis/PeriodicBC.hpp> 33 #include <amdis/ProblemStatBase.hpp> 34 #include <amdis/ProblemStatTraits.hpp> 35 #include <amdis/StandardProblemIteration.hpp> 36 #include <amdis/common/SharedPtr.hpp> 37 #include <amdis/common/TupleUtility.hpp> 38 #include <amdis/common/TypeTraits.hpp> 39 #include <amdis/GridFunctions.hpp> 40 #include <amdis/gridfunctions/DiscreteFunction.hpp> 41 #include <amdis/io/FileWriterBase.hpp> 42 #include <amdis/typetree/Concepts.hpp> 43 #include <amdis/typetree/TreePath.hpp> 48 template <
class Traits>
51 template <
class Traits>
64 using Grid = AdaptiveGrid_t<typename GridView::Grid>;
65 using Element =
typename GridView::template Codim<0>::Entity;
66 using WorldVector =
typename Element::Geometry::GlobalCoordinate;
67 using WorldMatrix = FieldMatrix<typename WorldVector::field_type, WorldVector::dimension, WorldVector::dimension>;
70 static constexpr
int dim = Grid::dimension;
73 static constexpr
int dow = Grid::dimensionworld;
75 using value_type =
typename Traits::CoefficientType;
77 using LinAlgTraits =
typename Traits::LinAlgTraits;
78 using Mat =
typename LinAlgTraits::template MatrixImpl<value_type>;
79 using Vec =
typename LinAlgTraits::template VectorImpl<value_type>;
80 using PartitionSet =
typename LinAlgTraits::PartitionSet;
100 template <
class Gr
id_>
104 adoptGrid(wrap_or_share(FWD(
grid)));
109 template <
class Gr
id_,
class Basis_,
class B_ = Underlying_t<Basis_>,
110 REQUIRES(Concepts::GlobalBasis<B_>)>
119 template <class Grid_, class PBF_, class GV_ = typename Underlying_t<Grid_>::LeafGridView,
120 REQUIRES(Concepts::PreBasisFactory<PBF_, GV_, MultiIndex_t<PBF_>>)>
124 adoptGlobalBasis(makeSharedPtr(GlobalBasis{
grid_->leafGridView(), preBasisFactory}));
162 template <
class Operator,
class RowTreePath = RootTreePath,
class ColTreePath = RootTreePath>
165 static constexpr
bool isValidTreePath =
166 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, RowTreePath> &&
167 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, ColTreePath>;
168 static_assert(isValidTreePath,
"Invalid row and/or col treepath passed to addMatrixOperator!");
170 if constexpr (isValidTreePath)
193 template <
class Operator,
class RowTreePath = RootTreePath,
class ColTreePath = RootTreePath>
196 using I =
typename GridView::Intersection;
197 static constexpr
bool isValidTreePath =
198 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, RowTreePath> &&
199 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, ColTreePath>;
200 static_assert(isValidTreePath,
"Invalid row and/or col treepath passed to addMatrixOperator!");
202 if constexpr (isValidTreePath)
225 template <
class Operator,
class TreePath = RootTreePath>
228 static constexpr
bool isValidTreePath =
229 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, TreePath>;
230 static_assert(isValidTreePath,
"Invalid treepath passed to addVectorOperator!");
232 if constexpr (isValidTreePath)
253 template <
class Operator,
class TreePath = RootTreePath>
256 using I =
typename GridView::Intersection;
257 static constexpr
bool isValidTreePath =
258 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, TreePath>;
259 static_assert(isValidTreePath,
"Invalid treepath passed to addVectorOperator!");
261 if constexpr (isValidTreePath)
289 template <
class Predicate,
class RowTreePath,
class ColTreePath,
class Values>
291 RowTreePath row, ColTreePath col,
292 Values
const& values);
294 template <
class RowTreePath,
class ColTreePath,
class Values>
296 RowTreePath row, ColTreePath col,
297 Values
const& values);
299 template <
class Identifier,
class Values>
302 addDirichletBC(FWD(
id), RootTreePath{}, RootTreePath{}, FWD(values));
326 bool asmMatrix =
true,
327 bool asmVector =
true)
override;
338 bool createMatrixData =
true,
339 bool storeMatrixData =
false)
override;
363 std::string
const&
name()
const override {
return name_; }
367 std::shared_ptr<Grid const>
grid()
const {
return grid_; }
394 std::shared_ptr<SystemVector const>
rhsVector()
const {
return rhs_; }
402 template <
class Range = void,
class... Indices>
405 assert(
bool(
solution_) &&
"You have to call initialize() before.");
406 return valueOf<Range>(*
solution_, ii...);
414 template <
class Range = void,
class... Indices>
417 assert(
bool(
solution_) &&
"You have to call initialize() before.");
418 return valueOf<Range>(*
solution_, ii...);
425 template <
class Solver_>
435 template <
class Gr
id_>
438 adoptGrid(wrap_or_share(FWD(
grid)));
440 createMatricesAndVectors();
450 template <
class Marker_>
453 auto marker = wrap_or_share(FWD(m));
454 auto it =
marker_.emplace(marker->name(), marker);
456 it.first->second->setMaximumMarking(
true);
462 std::size_t num =
marker_.erase(name);
463 test_warning(num == 1,
"A marker with the given name '{}' does not exist.", name);
475 void createGlobalBasis();
477 void createMatricesAndVectors();
480 void createFileWriter();
482 void adoptGlobalBasis(std::shared_ptr<GlobalBasis>
globalBasis)
488 void adoptGrid(std::shared_ptr<Grid>
const&
grid,
496 void adoptGrid(std::shared_ptr<Grid>
const& grid)
501 void adoptGrid(std::shared_ptr<typename Grid::HostGrid>
const& hostGrid)
503 auto grid = std::make_shared<Grid>(hostGrid);
509 void createGlobalBasisImpl(std::true_type);
510 void createGlobalBasisImpl(std::false_type);
512 void initGlobalBasis();
535 std::map<std::string, std::shared_ptr<Marker<Grid>>>
marker_;
551 std::shared_ptr<SystemVector>
rhs_;
558 std::list<BoundaryCondition<SystemMatrix, SolutionVector, SystemVector>>
constraints_;
564 template <
class Gr
id,
class B,
class =
void>
565 struct DeducedProblemTraits;
567 template <
class Gr
id,
class PB>
568 struct DeducedProblemTraits<Grid,GlobalBasis<PB>,void>
573 template <
class G,
class PBF>
574 struct DeducedProblemTraits<G,PBF,
575 std::enable_if_t<Concepts::PreBasisFactory<PBF, typename G::LeafGridView, MultiIndex_t<PBF>>>>
577 using Grid = AdaptiveGrid_t<G>;
578 using GridView =
typename Grid::LeafGridView;
579 using Basis = decltype(GlobalBasis{std::declval<GridView>(),std::declval<PBF>()});
584 template <
class Gr
id,
class Basis>
585 using DeducedProblemTraits_t =
typename DeducedProblemTraits<Grid,Basis>::type;
590 template <
class Gr
id,
class Basis>
601 #include "ProblemStat.inc.hpp" Flag oneIteration(AdaptInfo &adaptInfo, Flag toDo=FULL_ITERATION) override
Implementation of StandardProblemIteration::oneIteration.
Definition: ProblemStat.hpp:318
static constexpr int dow
Dimension of the world.
Definition: ProblemStat.hpp:73
auto solution(Indices... ii) const
Return a const view to a solution component.
Definition: ProblemStat.hpp:415
std::string const & name() const override
Implementation of ProblemStatBase::name.
Definition: ProblemStat.hpp:363
void addPeriodicBC(BoundaryType id, WorldMatrix const &A, WorldVector const &b)
Definition: ProblemStat.inc.hpp:334
The Flag class encapsulates flags which represents simple information. Used e.g. while mesh traversal...
Definition: Flag.hpp:13
ProblemStat(std::string const &name, Grid_ &&grid)
Definition: ProblemStat.hpp:101
std::map< std::string, std::vector< double > > estimates_
Definition: ProblemStat.hpp:555
void addMatrixOperator(BoundaryType b, Operator const &op, RowTreePath row={}, ColTreePath col={})
Operator evaluated on the boundary of the domain with boundary index b
Definition: ProblemStat.hpp:194
void addDirichletBC(Predicate const &predicate, RowTreePath row, ColTreePath col, Values const &values)
Add boundary conditions to the system.
Definition: ProblemStat.inc.hpp:290
void addMarker(Marker_ &&m)
Store the shared_ptr and the name of the marker in the problem.
Definition: ProblemStat.hpp:451
The base class for an operator to be used in an Assembler.
Definition: Operator.hpp:79
void addMatrixOperator(Operator const &op, RowTreePath row={}, ColTreePath col={})
Add an operator to A.
Definition: ProblemStat.hpp:163
std::shared_ptr< SystemMatrix > systemMatrix_
Matrix that is filled during assembling.
Definition: ProblemStat.hpp:544
void estimate(AdaptInfo &) override
Implementation of ProblemStatBase::estimate.
Definition: ProblemStat.hpp:342
constexpr bool Predicate
A predicate is a function that returns a boolean.
Definition: Concepts.hpp:142
void setGrid(Grid_ &&grid)
Definition: ProblemStat.hpp:436
std::shared_ptr< SystemVector > rhs_
Definition: ProblemStat.hpp:551
ProblemStat(std::string const &name)
Constructor. Takes the name of the problem that is used to access values corresponding to this proble...
Definition: ProblemStat.hpp:94
std::shared_ptr< SolutionVector > solution_
Vector with the solution components.
Definition: ProblemStat.hpp:547
The basic container that stores a base vector and a corresponding basis.
Definition: DOFVector.hpp:27
Definition: AdaptBase.hpp:6
StandardProblemIteration when derived from ProblemStat.
Definition: StandardProblemIteration.hpp:71
Definition: OperatorList.hpp:19
Base class for all markers.
Definition: Marker.hpp:28
auto solution(Indices... ii)
Return a mutable view to a solution component.
Definition: ProblemStat.hpp:403
std::shared_ptr< SolutionVector > solutionVector()
Returns a reference to the solution vector, solution_.
Definition: ProblemStat.hpp:389
Wrapper around a global basis providing default traits.
Definition: ProblemStatTraits.hpp:77
std::shared_ptr< GlobalBasis > globalBasis()
Return the globalBasis_.
Definition: ProblemStat.hpp:377
std::list< std::shared_ptr< FileWriterInterface > > filewriter_
A FileWriter object.
Definition: ProblemStat.hpp:532
Definition: ProblemStat.hpp:52
void initialize(Flag initFlag, Self *adoptProblem=nullptr, Flag adoptFlag=INIT_NOTHING)
Initialisation of the problem.
Definition: ProblemStat.inc.hpp:22
Standard implementation of ProblemTimeInterface for a time dependent problems.
Definition: ProblemInstat.hpp:24
std::shared_ptr< GlobalBasis > globalBasis_
FE space of this problem.
Definition: ProblemStat.hpp:529
void restore(Flag initFlag)
Read the grid and solution from backup files and initialize the problem.
Definition: ProblemStat.inc.hpp:107
std::map< std::string, std::shared_ptr< Marker< Grid > > > marker_
Pointer to the adaptation markers.
Definition: ProblemStat.hpp:535
ProblemStat(std::string const &name, Grid_ &&grid, Basis_ &&globalBasis)
Constructor taking a grid and basis. Wraps both in shared pointers.
Definition: ProblemStat.hpp:111
std::shared_ptr< Grid > grid()
Return the grid_.
Definition: ProblemStat.hpp:366
std::shared_ptr< SystemMatrix > systemMatrix()
Returns a reference to system-matrix, systemMatrix_.
Definition: ProblemStat.hpp:385
static std::optional< T > get(std::string const &key)
Get parameter-values from parameter-tree.
Definition: Initfile.hpp:25
static constexpr int dim
Dimension of the grid.
Definition: ProblemStat.hpp:70
void addVectorOperator(Operator const &op, TreePath path={})
Add an operator to rhs.
Definition: ProblemStat.hpp:226
std::shared_ptr< SystemVector > rhsVector()
Return a reference to the rhs system-vector, rhs.
Definition: ProblemStat.hpp:393
Flag markElements(AdaptInfo &adaptInfo) override
Implementation of ProblemStatBase::markElements.
Definition: ProblemStat.inc.hpp:379
void removeMarker(std::string name)
Remove a marker with the given name from the problem.
Definition: ProblemStat.hpp:460
Definition: LinearSolverInterface.hpp:8
Flag adaptGrid(AdaptInfo &adaptInfo) override
Implementation of ProblemStatBase::refineMesh.
Definition: ProblemStat.inc.hpp:446
ProblemStat(std::string const &name, Grid_ &&grid, PBF_ const &preBasisFactory)
Constructor taking a grid and pre-basis factory to create a global basis on the fly.
Definition: ProblemStat.hpp:121
Interface for time independent problems. Concrete problems must override all pure virtual methods...
Definition: ProblemStatBase.hpp:58
void writeFiles(AdaptInfo &adaptInfo, bool force=false)
Writes output files. If force=true write even if timestep out of write rhythm.
Definition: ProblemStat.inc.hpp:520
Implementation of LinearSolverInterface for EIGEN solvers.
Definition: LinearSolver.hpp:14
std::shared_ptr< Grid > grid_
Grid of this problem.
Definition: ProblemStat.hpp:520
void solve(AdaptInfo &adaptInfo, bool createMatrixData=true, bool storeMatrixData=false) override
Implementation of ProblemStatBase::solve.
Definition: ProblemStat.inc.hpp:346
typename PreBasis::GridView GridView
The grid view that the FE space is defined on.
Definition: GlobalBasis.hpp:63
GridView gridView() const
Return the gridView of the basis.
Definition: ProblemStat.hpp:370
Flag globalRefine(int n) override
Uniform global refinement by n level.
Definition: ProblemStat.inc.hpp:431
Holds adapt parameters and infos about the problem.
Definition: AdaptInfo.hpp:25
void addVectorOperator(BoundaryType b, Operator const &op, TreePath path={})
Operator evaluated on the boundary of the domain with boundary index b
Definition: ProblemStat.hpp:254
std::shared_ptr< LinearSolverInterface > solver()
Return a reference to the linear solver, linearSolver.
Definition: ProblemStat.hpp:381
std::shared_ptr< BoundaryManager< Grid > > boundaryManager()
Return the boundary manager to identify boundary segments.
Definition: ProblemStat.hpp:373
std::shared_ptr< BoundaryManager< Grid > > boundaryManager_
Management of boundary conditions.
Definition: ProblemStat.hpp:526
void buildAfterAdapt(AdaptInfo &adaptInfo, Flag flag, bool asmMatrix=true, bool asmVector=true) override
Implementation of ProblemStatBase::buildAfterCoarse.
Definition: ProblemStat.inc.hpp:461
Definition: OperatorList.hpp:15
std::shared_ptr< LinearSolverInterface > linearSolver_
Pointer to the estimators for this problem.
Definition: ProblemStat.hpp:541
constexpr bool GlobalBasis
A Dune::Functions::GlobalBasis type.
Definition: Concepts.hpp:190
std::list< BoundaryCondition< SystemMatrix, SolutionVector, SystemVector > > constraints_
List of constraints to apply to matrix, solution and rhs.
Definition: ProblemStat.hpp:558
void assemble(AdaptInfo &adaptInfo)
Assemble the linear system by calling buildAfterAdapt with asmMatrix and asmVector set to true...
Definition: ProblemStat.hpp:331
Manage boundary ids of boundary segments in a grid.
Definition: BoundaryManager.hpp:52
std::string const & name() const
Returns name_ of the Marker.
Definition: Marker.hpp:81
Flag globalCoarsen(int n) override
Uniform global grid coarsening by up to n level.
Definition: ProblemStat.inc.hpp:400
void setSolver(Solver_ &&solver)
Set a new linear solver for the problem.
Definition: ProblemStat.hpp:426
void removeMarker(Marker< Grid > const &marker)
Remove a marker from the problem.
Definition: ProblemStat.hpp:467
Interface class for boundary conditions.
Definition: BoundaryCondition.hpp:61
std::string gridName_
Name of the grid.
Definition: ProblemStat.hpp:523
std::string name_
Name of this problem.
Definition: ProblemStat.hpp:517
Flag oneIteration(AdaptInfo &adaptInfo, Flag toDo) override
Implementation of ProblemIterationInterface::oneIteration()
Definition: StandardProblemIteration.cpp:35