} // namespace internal
+
template <int dim, typename PolynomialType>
void
TensorProductPolynomials<dim, PolynomialType>::evaluate(
std::vector<Tensor<3, dim>> &third_derivatives,
std::vector<Tensor<4, dim>> &fourth_derivatives) const
{
- if constexpr (dim == 0)
- return;
-
Assert(dim <= 3, ExcNotImplemented());
Assert(values.size() == this->n() || values.empty(),
ExcDimensionMismatch2(values.size(), this->n(), 0));
+template <>
+void
+TensorProductPolynomials<0, Polynomials::Polynomial<double>>::evaluate(
+ const Point<0> &,
+ std::vector<double> &,
+ std::vector<Tensor<1, 0>> &,
+ std::vector<Tensor<2, 0>> &,
+ std::vector<Tensor<3, 0>> &,
+ std::vector<Tensor<4, 0>> &) const
+{
+ constexpr int dim = 0;
+ AssertThrow(dim > 0, ExcNotImplemented());
+}
+
+
+
template <int dim, typename PolynomialType>
std::unique_ptr<ScalarPolynomialsBase<dim>>
TensorProductPolynomials<dim, PolynomialType>::clone() const