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>
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;
76 using LinAlgTraits =
typename Traits::LinAlgTraits;
77 using Mat =
typename LinAlgTraits::template MatrixImpl<typename Traits::CoefficientType>;
78 using Vec =
typename LinAlgTraits::template VectorImpl<typename Traits::CoefficientType>;
79 using PartitionSet =
typename LinAlgTraits::PartitionSet;
98 template <
class Gr
id_>
107 template <
class Gr
id_,
class Basis_,
class B_ = Underlying_t<Basis_>,
108 REQUIRES(Concepts::GlobalBasis<B_>)>
117 template <class Grid_, class PBF_, class GV_ = typename Underlying_t<Grid_>::LeafGridView,
118 REQUIRES(Concepts::PreBasisFactory<PBF_, GV_, MultiIndex_t<PBF_>>)>
122 adoptGlobalBasis(makeSharedPtr(GlobalBasis{
grid_->leafGridView(), preBasisFactory}));
160 template <
class Operator,
class RowTreePath = RootTreePath,
class ColTreePath = RootTreePath>
163 static constexpr
bool isValidTreePath =
164 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, RowTreePath> &&
165 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, ColTreePath>;
166 static_assert(isValidTreePath,
"Invalid row and/or col treepath passed to addMatrixOperator!");
168 if constexpr (isValidTreePath)
191 template <
class Operator,
class RowTreePath = RootTreePath,
class ColTreePath = RootTreePath>
194 using I =
typename GridView::Intersection;
195 static constexpr
bool isValidTreePath =
196 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, RowTreePath> &&
197 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, ColTreePath>;
198 static_assert(isValidTreePath,
"Invalid row and/or col treepath passed to addMatrixOperator!");
200 if constexpr (isValidTreePath)
223 template <
class Operator,
class TreePath = RootTreePath>
226 static constexpr
bool isValidTreePath =
227 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, TreePath>;
228 static_assert(isValidTreePath,
"Invalid treepath passed to addVectorOperator!");
230 if constexpr (isValidTreePath)
251 template <
class Operator,
class TreePath = RootTreePath>
254 using I =
typename GridView::Intersection;
255 static constexpr
bool isValidTreePath =
256 Concepts::ValidTreePath<typename GlobalBasis::LocalView::Tree, TreePath>;
257 static_assert(isValidTreePath,
"Invalid treepath passed to addVectorOperator!");
259 if constexpr (isValidTreePath)
287 template <
class Predicate,
class RowTreePath,
class ColTreePath,
class Values>
289 RowTreePath row, ColTreePath col,
290 Values
const& values);
292 template <
class RowTreePath,
class ColTreePath,
class Values>
294 RowTreePath row, ColTreePath col,
295 Values
const& values);
297 template <
class Identifier,
class Values>
300 addDirichletBC(FWD(
id), RootTreePath{}, RootTreePath{}, FWD(values));
324 bool asmMatrix =
true,
325 bool asmVector =
true)
override;
336 bool createMatrixData =
true,
337 bool storeMatrixData =
false)
override;
361 std::string
const&
name()
const override {
return name_; }
365 std::shared_ptr<Grid const>
grid()
const {
return grid_; }
392 std::shared_ptr<SystemVector const>
rhsVector()
const {
return rhs_; }
400 template <
class Range = void,
class...
Indices>
403 assert(
bool(
solution_) &&
"You have to call initialize() before.");
404 return valueOf<Range>(*
solution_, ii...);
412 template <
class Range = void,
class...
Indices>
415 assert(
bool(
solution_) &&
"You have to call initialize() before.");
416 return valueOf<Range>(*
solution_, ii...);
423 template <
class Solver_>
433 template <
class Gr
id_>
438 createMatricesAndVectors();
448 template <
class Marker_>
452 auto it =
marker_.emplace(marker->name(), marker);
454 it.first->second->setMaximumMarking(
true);
460 std::size_t num =
marker_.erase(name);
461 test_warning(num == 1,
"A marker with the given name '{}' does not exist.", name);
473 void createGlobalBasis();
475 void createMatricesAndVectors();
478 void createFileWriter();
480 void adoptGlobalBasis(std::shared_ptr<GlobalBasis>
globalBasis)
486 void adoptGrid(std::shared_ptr<Grid>
const&
grid,
494 void adoptGrid(std::shared_ptr<Grid>
const& grid)
499 void adoptGrid(std::shared_ptr<typename Grid::HostGrid>
const& hostGrid)
501 auto grid = std::make_shared<Grid>(hostGrid);
507 void createGlobalBasisImpl(std::true_type);
508 void createGlobalBasisImpl(std::false_type);
510 void initGlobalBasis();
533 std::map<std::string, std::shared_ptr<Marker<Grid>>>
marker_;
549 std::shared_ptr<SystemVector>
rhs_;
556 std::list<BoundaryCondition<SystemMatrix, SolutionVector, SystemVector>>
constraints_;
562 template <
class Gr
id,
class B,
class =
void>
563 struct DeducedProblemTraits;
565 template <
class Gr
id,
class PB>
566 struct DeducedProblemTraits<Grid,GlobalBasis<PB>,void>
571 template <
class G,
class PBF>
572 struct DeducedProblemTraits<G,PBF,
573 std::enable_if_t<Concepts::PreBasisFactory<PBF, typename G::LeafGridView, MultiIndex_t<PBF>>>>
576 using GridView =
typename Grid::LeafGridView;
577 using Basis = decltype(GlobalBasis{std::declval<GridView>(),std::declval<PBF>()});
582 template <
class Gr
id,
class Basis>
583 using DeducedProblemTraits_t =
typename DeducedProblemTraits<Grid,Basis>::type;
588 template <
class Gr
id,
class Basis>
599 #include "ProblemStat.inc.hpp" Flag oneIteration(AdaptInfo &adaptInfo, Flag toDo=FULL_ITERATION) override
Implementation of StandardProblemIteration::oneIteration.
Definition: ProblemStat.hpp:316
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:413
std::string const & name() const override
Implementation of ProblemStatBase::name.
Definition: ProblemStat.hpp:361
void addPeriodicBC(BoundaryType id, WorldMatrix const &A, WorldVector const &b)
Definition: ProblemStat.inc.hpp:338
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:99
std::map< std::string, std::vector< double > > estimates_
Definition: ProblemStat.hpp:553
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:192
void addDirichletBC(Predicate const &predicate, RowTreePath row, ColTreePath col, Values const &values)
Add boundary conditions to the system.
Definition: ProblemStat.inc.hpp:294
std::index_sequence< I... > Indices
class that represents a sequence of indices
Definition: Index.hpp:40
void addMarker(Marker_ &&m)
Store the shared_ptr and the name of the marker in the problem.
Definition: ProblemStat.hpp:449
typename Impl::AdaptiveGridImpl< HostGrid >::type AdaptiveGrid_t
Definition: AdaptiveGrid.hpp:367
void addMatrixOperator(Operator const &op, RowTreePath row={}, ColTreePath col={})
Add an operator to A.
Definition: ProblemStat.hpp:161
std::shared_ptr< SystemMatrix > systemMatrix_
Matrix that is filled during assembling.
Definition: ProblemStat.hpp:542
void estimate(AdaptInfo &) override
Implementation of ProblemStatBase::estimate.
Definition: ProblemStat.hpp:340
constexpr bool Predicate
A predicate is a function that returns a boolean.
Definition: Concepts.hpp:142
void setGrid(Grid_ &&grid)
Definition: ProblemStat.hpp:434
std::shared_ptr< SystemVector > rhs_
Definition: ProblemStat.hpp:549
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:92
std::shared_ptr< SolutionVector > solution_
Vector with the solution components.
Definition: ProblemStat.hpp:545
The basic container that stores a base vector and a corresponding basis.
Definition: DOFVector.hpp:27
Contains all classes needed for solving linear and non linear equation systems.
Definition: AdaptBase.hpp:6
StandardProblemIteration when derived from ProblemStat.
Definition: StandardProblemIteration.hpp:71
Definition: OperatorList.hpp:16
Base class for all markers.
Definition: Marker.hpp:28
auto solution(Indices... ii)
Return a mutable view to a solution component.
Definition: ProblemStat.hpp:401
std::shared_ptr< LinearSolver > solver()
Return a reference to the linear solver, linearSolver.
Definition: ProblemStat.hpp:379
std::shared_ptr< SolutionVector > solutionVector()
Returns a reference to the solution vector, solution_.
Definition: ProblemStat.hpp:387
Wrapper around a global basis providing default traits.
Definition: ProblemStatTraits.hpp:77
std::shared_ptr< GlobalBasis > globalBasis()
Return the globalBasis_.
Definition: ProblemStat.hpp:375
std::list< std::shared_ptr< FileWriterInterface > > filewriter_
A FileWriter object.
Definition: ProblemStat.hpp:530
Definition: ProblemStat.hpp:52
void initialize(Flag initFlag, Self *adoptProblem=nullptr, Flag adoptFlag=INIT_NOTHING)
Initialisation of the problem.
Definition: ProblemStat.inc.hpp:23
Standard implementation of ProblemTimeInterface for a time dependent problems.
Definition: ProblemInstat.hpp:22
std::shared_ptr< GlobalBasis > globalBasis_
FE space of this problem.
Definition: ProblemStat.hpp:527
void restore(Flag initFlag)
Read the grid and solution from backup files and initialize the problem.
Definition: ProblemStat.inc.hpp:108
std::map< std::string, std::shared_ptr< Marker< Grid > > > marker_
Pointer to the adaptation markers.
Definition: ProblemStat.hpp:533
ProblemStat(std::string const &name, Grid_ &&grid, Basis_ &&globalBasis)
Constructor taking a grid and basis. Wraps both in shared pointers.
Definition: ProblemStat.hpp:109
std::shared_ptr< Grid > grid()
Return the grid_.
Definition: ProblemStat.hpp:364
std::shared_ptr< SystemMatrix > systemMatrix()
Returns a reference to system-matrix, systemMatrix_.
Definition: ProblemStat.hpp:383
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:224
std::shared_ptr< SystemVector > rhsVector()
Return a reference to the rhs system-vector, rhs.
Definition: ProblemStat.hpp:391
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:458
Abstract base class for linear solvers.
Definition: LinearSolverInterface.hpp:27
Flag adaptGrid(AdaptInfo &adaptInfo) override
Implementation of ProblemStatBase::refineMesh.
Definition: ProblemStat.inc.hpp:441
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:119
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:512
std::shared_ptr< Grid > grid_
Grid of this problem.
Definition: ProblemStat.hpp:518
void solve(AdaptInfo &adaptInfo, bool createMatrixData=true, bool storeMatrixData=false) override
Implementation of ProblemStatBase::solve.
Definition: ProblemStat.inc.hpp:350
typename PreBasis::GridView GridView
The grid view that the FE space is defined on.
Definition: GlobalBasis.hpp:66
GridView gridView() const
Return the gridView of the basis.
Definition: ProblemStat.hpp:368
Flag globalRefine(int n) override
Uniform global refinement by n level.
Definition: ProblemStat.inc.hpp:426
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:252
std::shared_ptr< T > wrap_or_share(T &t)
Definition: SharedPtr.hpp:22
std::shared_ptr< BoundaryManager< Grid > > boundaryManager()
Return the boundary manager to identify boundary segments.
Definition: ProblemStat.hpp:371
std::shared_ptr< BoundaryManager< Grid > > boundaryManager_
Management of boundary conditions.
Definition: ProblemStat.hpp:524
void buildAfterAdapt(AdaptInfo &adaptInfo, Flag flag, bool asmMatrix=true, bool asmVector=true) override
Implementation of ProblemStatBase::buildAfterCoarse.
Definition: ProblemStat.inc.hpp:456
Definition: OperatorList.hpp:14
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:556
void assemble(AdaptInfo &adaptInfo)
Assemble the linear system by calling buildAfterAdapt with asmMatrix and asmVector set to true...
Definition: ProblemStat.hpp:329
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:395
void setSolver(Solver_ &&solver)
Set a new linear solver for the problem.
Definition: ProblemStat.hpp:424
void removeMarker(Marker< Grid > const &marker)
Remove a marker from the problem.
Definition: ProblemStat.hpp:465
Interface class for boundary conditions.
Definition: BoundaryCondition.hpp:61
std::string gridName_
Name of the grid.
Definition: ProblemStat.hpp:521
void test_warning(bool condition, std::string const &str, Args &&... args)
test for condition and in case of failure print message
Definition: Output.hpp:183
std::string name_
Name of this problem.
Definition: ProblemStat.hpp:515
std::shared_ptr< LinearSolver > linearSolver_
Pointer to the estimators for this problem.
Definition: ProblemStat.hpp:539
Flag oneIteration(AdaptInfo &adaptInfo, Flag toDo) override
Implementation of ProblemIterationInterface::oneIteration()
Definition: StandardProblemIteration.cpp:35