function_spherical.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later
// Copyright (C) 2017 - 2026 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Detailed license information governing the source code and contributions
// can be found in LICENSE.md and CONTRIBUTING.md at the top level directory.
//
// -----------------------------------------------------------------------------
 
#ifndef dealii_function_spherical_h
#define dealii_function_spherical_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/function.h>
#include <deal.II/base/point.h>
#include <deal.II/base/symmetric_tensor.h>
#include <deal.II/base/tensor.h>
 
#include <array>
#include <cstddef>
 
DEAL_II_NAMESPACE_OPEN
 
namespace Functions
{
  template <int dim>
  class Spherical : public Function<dim>
  {
  public:
    Spherical(const Point<dim>  &center       = Point<dim>(),
              const unsigned int n_components = 1);
 
    virtual double
    value(const Point<dim>  &point,
          const unsigned int component = 0) const override;
 
    virtual Tensor<1, dim>
    gradient(const Point<dim>  &p,
             const unsigned int component = 0) const override;
 
    virtual SymmetricTensor<2, dim>
    hessian(const Point<dim>  &p,
            const unsigned int component = 0) const override;
 
    virtual std::size_t
    memory_consumption() const override;
 
  private:
    virtual double
    svalue(const std::array<double, dim> &sp,
           const unsigned int             component) const;
 
    virtual std::array<double, dim>
    sgradient(const std::array<double, dim> &sp,
              const unsigned int             component) const;
 
    virtual std::array<double, 6>
    shessian(const std::array<double, dim> &sp,
             const unsigned int             component) const;
 
    const Tensor<1, dim> coordinate_system_offset;
  };
} // namespace Functions
 
DEAL_II_NAMESPACE_CLOSE
 
#endif