Commit cc255a6e authored by Oliver Sander's avatar Oliver Sander Committed by Christoph Grüninger

Remove the 'evaluate<order>' method the the QkLocalBasis class

parent 5e458f7b
......@@ -6,7 +6,6 @@
#include <numeric>
#include <dune/common/deprecated.hh>
#include <dune/common/fvector.hh>
#include <dune/common/fmatrix.hh>
#include <dune/common/power.hh>
......@@ -173,41 +172,6 @@ namespace Dune
}
}
/** \brief Evaluate derivative in a given direction
* \param [in] direction The direction to derive in
* \param [in] in Position where to evaluate
* \param [out] out The return value
*/
template<int diffOrder>
inline void DUNE_DEPRECATED_MSG("Use method 'partial' instead!")
evaluate(
const std::array<int,1>& direction,
const typename Traits::DomainType& in,
std::vector<typename Traits::RangeType>& out) const
{
static_assert(diffOrder == 1, "We only can compute first derivatives");
out.resize(size());
// Loop over all shape functions
for (size_t i=0; i<size(); i++)
{
// convert index i to multiindex
Dune::FieldVector<int,d> alpha(multiindex(i));
// Loop over all coordinate directions
std::size_t j = direction[0];
// Initialize: the overall expression is a product
// if j-th bit of i is set to -1, else 1
out[i][0] = dp(alpha[j],in[j]);
// rest of the product
for (std::size_t l=0; l<d; l++)
if (l!=j)
out[i][0] *= p(alpha[l],in[l]);
}
}
//! \brief Polynomial order of the shape functions
unsigned int order () const
{
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment