dune-localfunctions 2.10
Loading...
Searching...
No Matches
meta/power/basis.hh
Go to the documentation of this file.
1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2// vi: set et ts=4 sw=2 sts=2:
3// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file LICENSE.md in module root
4// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
5
6#ifndef DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
7#define DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
8
9#include <numeric>
10#include <cstddef>
11#include <vector>
12
13#include <dune/common/fmatrix.hh>
14#include <dune/common/fvector.hh>
15
16namespace Dune {
17
19
25 template<class Backend, std::size_t dimR>
26 class PowerBasis {
27 static_assert(Backend::Traits::dimRange == 1,
28 "PowerBasis works only with scalar backends");
29
30 // don't use a reference here so this class stays copyable
31 const Backend *backend;
32
33 public:
35 struct Traits : public Backend::Traits
36 {
38 static const std::size_t dimRange = dimR;
40 typedef FieldVector<typename Traits::RangeField, dimR> Range;
41
43
47 typedef FieldMatrix<typename Traits::RangeField, dimR,
48 Traits::dimDomainGlobal> Jacobian;
49 };
50
52
58 PowerBasis(const Backend &backend_) : backend(&backend_) { }
59
61 std::size_t size () const { return backend->size()*dimR; }
63 std::size_t order () const { return backend->order(); }
64
66 void evaluateFunction(const typename Traits::DomainLocal& in,
67 std::vector<typename Traits::Range>& out) const
68 {
69 std::vector<typename Backend::Traits::Range> backendValues;
70 backend->evaluateFunction(in, backendValues);
71 out.assign(size(), typename Traits::Range(0));
72 for(std::size_t d = 0; d < dimR; ++d)
73 for(std::size_t i = 0; i < backend->size(); ++i)
74 out[d*backend->size()+i][d] = backendValues[i][0];
75 }
76
78 void evaluateJacobian(const typename Traits::DomainLocal& in,
79 std::vector<typename Traits::Jacobian>& out) const
80 {
81 std::vector<typename Backend::Traits::Jacobian> backendValues;
82 backend->evaluateJacobian(in, backendValues);
83 out.assign(size(), typename Traits::Jacobian(0));
84 for(std::size_t d = 0; d < dimR; ++d)
85 for(std::size_t i = 0; i < backend->size(); ++i)
86 out[d*backend->size()+i][d] = backendValues[i][0];
87 }
88
90 void partial (const std::array<unsigned int, Backend::Traits::dimDomainGlobal>& order,
91 const typename Traits::DomainLocal& in, // position
92 std::vector<typename Traits::Range>& out) const // return value
93 {
94 auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
95 if (totalOrder == 0) {
96 evaluateFunction(in, out);
97 } else {
98 DUNE_THROW(NotImplemented, "Desired derivative order is not implemented");
99 }
100 }
101 };
102
103} // namespace Dune
104
105#endif // DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
Definition bdfmcube.hh:18
Meta-basis turning a scalar basis into vector-valued basis.
Definition meta/power/basis.hh:26
PowerBasis(const Backend &backend_)
Construct a PowerBasis.
Definition meta/power/basis.hh:58
std::size_t size() const
Number of shape functions.
Definition meta/power/basis.hh:61
void partial(const std::array< unsigned int, Backend::Traits::dimDomainGlobal > &order, const typename Traits::DomainLocal &in, std::vector< typename Traits::Range > &out) const
Evaluate partial derivatives of all shape functions.
Definition meta/power/basis.hh:90
void evaluateJacobian(const typename Traits::DomainLocal &in, std::vector< typename Traits::Jacobian > &out) const
Evaluate Jacobian of all shape functions at given position.
Definition meta/power/basis.hh:78
void evaluateFunction(const typename Traits::DomainLocal &in, std::vector< typename Traits::Range > &out) const
Evaluate all shape functions at given position.
Definition meta/power/basis.hh:66
std::size_t order() const
Polynomial order of the shape functions for quadrature.
Definition meta/power/basis.hh:63
types of domain and range
Definition meta/power/basis.hh:36
FieldMatrix< typename Traits::RangeField, dimR, Traits::dimDomainGlobal > Jacobian
Jacobian properties.
Definition meta/power/basis.hh:48
FieldVector< typename Traits::RangeField, dimR > Range
Type used for range values.
Definition meta/power/basis.hh:40
static const std::size_t dimRange
Dimension of the range values.
Definition meta/power/basis.hh:38