133 for (
unsigned int i = 0; i < dim; ++i)
134 for (
unsigned int j = i;
j < dim; ++
j)
148 for (
unsigned int i = 0; i < dim; ++i)
149 for (
unsigned int j = i;
j < dim; ++
j)
150 out[i][
j] += val * in[i] * in[
j];
158 const unsigned int n_components)
160 , coordinate_system_offset(p)
170 const unsigned int component)
const
173 const std::array<double, dim>
sp =
175 return svalue(
sp, component);
183 const unsigned int )
const
196 constexpr int dim = 3;
198 const std::array<double, dim>
sp =
200 const std::array<double, dim>
sg = sgradient(
sp, component);
239 const unsigned int )
const
252 constexpr int dim = 3;
254 const std::array<double, dim>
sp =
256 const std::array<double, dim>
sg = sgradient(
sp, component);
257 const std::array<double, 6> sh = shessian(
sp, component);
264 const double r =
sp[0];
271 const double r2 =
r *
r;
283 const double c_ur2 = sh[0];
319 const unsigned int )
const
328 std::array<double, dim>
330 const unsigned int )
const
333 return std::array<double, dim>();
339 std::array<double, 6>
341 const unsigned int )
const
344 return std::array<double, 6>();
virtual SymmetricTensor< 2, dim > hessian(const Point< dim > &p, const unsigned int component=0) const override
Spherical(const Point< dim > ¢er=Point< dim >(), const unsigned int n_components=1)
virtual Tensor< 1, dim > gradient(const Point< dim > &p, const unsigned int component=0) const override
virtual std::array< double, dim > sgradient(const std::array< double, dim > &sp, const unsigned int component) const
virtual double value(const Point< dim > &point, const unsigned int component=0) const override
virtual std::size_t memory_consumption() const override
virtual std::array< double, 6 > shessian(const std::array< double, dim > &sp, const unsigned int component) const
virtual double svalue(const std::array< double, dim > &sp, const unsigned int component) const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define DEAL_II_NOT_IMPLEMENTED()
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
static ::ExceptionBase & ExcDivideByZero()
#define AssertThrow(cond, exc)
std::array< double, dim > to_spherical(const Point< dim > &point)
::VectorizedArray< Number, width > cos(const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > sin(const ::VectorizedArray< Number, width > &)