/usr/include/dune/localfunctions/meta/power/basis.hh is in libdune-localfunctions-dev 2.5.0-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 | // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
#define DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
#include <numeric>
#include <cstddef>
#include <vector>
#include <dune/common/fmatrix.hh>
#include <dune/common/fvector.hh>
namespace Dune {
//! Meta-basis turning a scalar basis into vector-valued basis
/**
* @ingroup BasisImplementation
*
* \tparam Backend Type of basis to take the power of.
* \tparam dimR Power to raise the basis to.
*/
template<class Backend, std::size_t dimR>
class PowerBasis {
static_assert(Backend::Traits::dimRange == 1,
"PowerBasis works only with scalar backends");
// don't use a reference here so this class stays copyable
const Backend *backend;
public:
//! types of domain and range
struct Traits : public Backend::Traits
{
//! Dimension of the range values
static const std::size_t dimRange = dimR;
//! Type used for range values
typedef FieldVector<typename Traits::RangeField, dimR> Range;
//! Jacobian properties
/**
* \note The Jacobian should be some matrix type with \c dimRange x
* \c dimDomainGlobal components of type \c RangeField.
*/
typedef FieldMatrix<typename Traits::RangeField, dimR,
Traits::dimDomainGlobal> Jacobian;
//! maximum number of partial derivatives supported
static const std::size_t diffOrder =
Backend::Traits::diffOrder > 0 ? 1 : 0;
};
//! Construct a PowerBasis
/**
* \param backend_ Backend basis object to construct this object from.
* This object holds a reference to the backend object.
* This reference is also copied when this object is
* copied.
*/
PowerBasis(const Backend &backend_) : backend(&backend_) { }
//! Number of shape functions
std::size_t size () const { return backend->size()*dimR; }
//! Polynomial order of the shape functions for quadrature
std::size_t order () const { return backend->order(); }
//! Evaluate all shape functions at given position
void evaluateFunction(const typename Traits::DomainLocal& in,
std::vector<typename Traits::Range>& out) const
{
std::vector<typename Backend::Traits::Range> backendValues;
backend->evaluateFunction(in, backendValues);
out.assign(size(), typename Traits::Range(0));
for(std::size_t d = 0; d < dimR; ++d)
for(std::size_t i = 0; i < backend->size(); ++i)
out[d*backend->size()+i][d] = backendValues[i][0];
}
//! Evaluate Jacobian of all shape functions at given position
void evaluateJacobian(const typename Traits::DomainLocal& in,
std::vector<typename Traits::Jacobian>& out) const
{
std::vector<typename Backend::Traits::Jacobian> backendValues;
backend->evaluateJacobian(in, backendValues);
out.assign(size(), typename Traits::Jacobian(0));
for(std::size_t d = 0; d < dimR; ++d)
for(std::size_t i = 0; i < backend->size(); ++i)
out[d*backend->size()+i][d] = backendValues[i][0];
}
//! \brief Evaluate partial derivatives of all shape functions
void partial (const std::array<unsigned int, Backend::Traits::dimDomainGlobal>& order,
const typename Traits::DomainLocal& in, // position
std::vector<typename Traits::Range>& out) const // return value
{
auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
if (totalOrder == 0) {
evaluateFunction(in, out);
} else {
DUNE_THROW(NotImplemented, "Desired derivative order is not implemented");
}
}
};
} // namespace Dune
#endif // DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
|