/usr/include/dune/localfunctions/whitney/edges0.5/interpolation.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 | // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_LOCALFUNCTIONS_WHITNEY_EDGES0_5_INTERPOLATION_HH
#define DUNE_LOCALFUNCTIONS_WHITNEY_EDGES0_5_INTERPOLATION_HH
#include <cstddef>
#include <vector>
#include <dune/localfunctions/whitney/edges0.5/common.hh>
namespace Dune {
//////////////////////////////////////////////////////////////////////
//
// Interpolation
//
//! Interpolation for lowest order edge elements on simplices
/**
* \tparam Geometry Type of the local-to-global map.
* \tparam RF Type to represent the field in the range.
*
* \nosubgrouping
*/
template<class Geometry, class Traits_>
class EdgeS0_5Interpolation :
private EdgeS0_5Common<Traits_::dimDomainLocal,
typename Traits_::DomainField>
{
public:
typedef Traits_ Traits;
private:
static const std::size_t dim = Traits::dimDomainLocal;
typedef EdgeS0_5Common<dim, typename Traits::DomainField> Base;
using Base::refelem;
using Base::s;
std::vector<typename Traits::DomainGlobal> edgev;
public:
//! constructor
/**
* \param geo Geometry of the element to contruct a local basis
* for.
* \param vertexOrder Vertex ordering information. Only the vertex order
* on the dim=1 sub-entities (edges) is required.
*/
template<typename VertexOrder>
EdgeS0_5Interpolation(const Geometry& geo,
const VertexOrder& vertexOrder) :
edgev(s)
{
for(std::size_t i = 0; i < s; ++i) {
const std::size_t i0 = refelem.subEntity(i,dim-1,0,dim);
const std::size_t i1 = refelem.subEntity(i,dim-1,1,dim);
edgev[i] = geo.corner(i1);
edgev[i] -= geo.corner(i0);
edgev[i] /= edgev[i].two_norm();
const typename VertexOrder::iterator& edgeVertexOrder =
vertexOrder.begin(dim-1, i);
if(edgeVertexOrder[0] > edgeVertexOrder[1])
edgev[i] *= -1;
}
}
//! Interpolation of a function
template<typename F, typename C>
void interpolate(const F& f, std::vector<C>& out) const {
typename Traits::Range y;
out.resize(s);
for(std::size_t i = 0; i < s; ++i) {
f.evaluate(refelem.position(i,dim-1), y);
out[i] = y * edgev[i];
}
}
};
} // namespace Dune
#endif // DUNE_LOCALFUNCTIONS_WHITNEY_EDGES0_5_INTERPOLATION_HH
|