/usr/include/dune/localfunctions/meta/power/basis.hh is in libdune-localfunctions-dev 2.3.1-1.
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 | // -*- 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 <cstddef>
#include <vector>
#include <dune/common/fmatrix.hh>
#include <dune/common/fvector.hh>
#include <dune/common/static_assert.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 {
dune_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];
}
};
} // namespace Dune
#endif // DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
|