function_level_set.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2019 - 2021 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, 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 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------
 
#ifndef dealii_function_level_set_h
#define dealii_function_level_set_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/function.h>
 
DEAL_II_NAMESPACE_OPEN
 
namespace Functions
{
  namespace LevelSet
  {
    template <int dim>
    class Sphere : public Function<dim>
    {
    public:
      Sphere(const Point<dim> &center = Point<dim>(), const double radius = 1);
 
      double
      value(const Point<dim> & point,
            const unsigned int component = 0) const override;
 
      Tensor<1, dim>
      gradient(const Point<dim> & point,
               const unsigned int component = 0) const override;
 
      SymmetricTensor<2, dim>
      hessian(const Point<dim> & point,
              const unsigned int component = 0) const override;
 
    private:
      const Point<dim> center;
      const double     radius;
    };
 
 
    template <int dim>
    class Plane : public Function<dim>
    {
    public:
      Plane(const Point<dim> &point, const Tensor<1, dim> &normal);
 
      double
      value(const Point<dim> & point,
            const unsigned int component = 0) const override;
 
      Tensor<1, dim>
      gradient(const Point<dim> &,
               const unsigned int component = 0) const override;
 
      SymmetricTensor<2, dim>
      hessian(const Point<dim> &,
              const unsigned int component = 0) const override;
 
    private:
      const Point<dim>     point_in_plane;
      const Tensor<1, dim> normal;
    };
 
  } // namespace LevelSet
} // namespace Functions
 
DEAL_II_NAMESPACE_CLOSE
 
#endif /* CE8BBB3E_B726_40A7_B963_561AE7B84973 */