/usr/include/dolfin/fem/BasisFunction.h is in libdolfin-dev 2017.2.0.post0-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 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 | // Copyright (C) 2013 Anders Logg
//
// This file is part of DOLFIN.
//
// DOLFIN is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// DOLFIN is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with DOLFIN. If not, see <http://www.gnu.org/licenses/>.
//
// First added: 2009-01-01
// Last changed: 2013-03-04
#ifndef __BASIS_FUNCTION_H
#define __BASIS_FUNCTION_H
#include <memory>
#include <vector>
#include <ufc.h>
#include <dolfin/fem/FiniteElement.h>
namespace dolfin
{
/// Represention of a finite element basis function.
/// It can be used for computation of basis function values
/// and derivatives.
///
/// Evaluation of basis functions is also possible through the use
/// of the functions ``evaluate_basis`` and
/// ``evaluate_basis_derivatives`` available in the _FiniteElement_
/// class. The BasisFunction class relies on these functions for
/// evaluation but also implements the ufc::function interface which
/// allows evaluate_dof to be evaluated for a basis function (on a
/// possibly different element).
class BasisFunction : public ufc::function
{
public:
/// Create basis function with given index on element on given cell
///
/// @param index (std::size_t)
/// The index of the basis function.
/// @param element (_FiniteElement_)
/// The element to create basis function on.
/// @param coordinate_dofs (std::vector<double>&)
/// The coordinate dofs of the cell
BasisFunction(std::size_t index,
std::shared_ptr<const FiniteElement> element,
const std::vector<double>& coordinate_dofs)
: _index(index), _element(element), _coordinate_dofs(coordinate_dofs) {}
/// Destructor
~BasisFunction() {}
/// Update the basis function index
///
/// @param index (std::size_t)
/// The index of the basis function.
void update_index(std::size_t index)
{ _index = index; }
/// Evaluate basis function at given point
///
/// @param values (double)
/// The values of the function at the point.
/// @param x (double)
/// The coordinates of the point.
void eval(double* values, const double* x) const
{
// Note: assuming cell_orientation = 0
dolfin_assert(_element);
_element->evaluate_basis(_index, values, x, _coordinate_dofs.data(), 0);
}
/// Evaluate all order n derivatives at given point
///
/// @param values (double)
/// The values of derivatives at the point.
/// @param x (double)
/// The coordinates of the point.
/// @param n (std::size_t)
/// The order of derivation.
void eval_derivatives(double* values, const double* x, std::size_t n) const
{
// Note: assuming cell_orientation = 0
dolfin_assert(_element);
_element->evaluate_basis_derivatives(_index, n, values, x,
_coordinate_dofs.data(), 0);
}
//--- Implementation of ufc::function interface ---
/// Evaluate function at given point in cell
///
/// @param values (double)
/// The values of the function at the point..
/// @param coordinates (double)
/// The coordinates of the point.
/// @param cell (ufc::cell)
/// The cell.
void evaluate(double* values, const double* coordinates,
const ufc::cell& cell) const
{ eval(values, coordinates); }
private:
// The index
std::size_t _index;
// The finite element
std::shared_ptr<const FiniteElement> _element;
// Cell coordinates
const std::vector<double> _coordinate_dofs;
};
}
#endif
|