manifold_lib.h

// ---------------------------------------------------------------------
//
// Copyright (C) 1999 - 2016 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 at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------

#ifndef dealii__manifold_lib_h
#define dealii__manifold_lib_h


#include <deal.II/base/config.h>
#include <deal.II/grid/manifold.h>
#include <deal.II/base/function.h>
#include <deal.II/base/function_parser.h>

DEAL_II_NAMESPACE_OPEN

template <int dim, int spacedim = dim>
class SphericalManifold : public ChartManifold<dim, spacedim, spacedim>
{
public:
  SphericalManifold(const Point<spacedim> center = Point<spacedim>());

  virtual Point<spacedim>
  pull_back(const Point<spacedim> &space_point) const;

  virtual Point<spacedim>
  push_forward(const Point<spacedim> &chart_point) const;

  virtual Point<spacedim>
  get_new_point(const Quadrature<spacedim> &quad) const;

  const Point<spacedim> center;
private:

  static Point<spacedim> get_periodicity();
};


template <int dim, int spacedim = dim>
class CylindricalManifold : public Manifold<dim,spacedim>
{
public:
  CylindricalManifold (const unsigned int axis = 0,
                       const double tolerance = 1e-10);

  CylindricalManifold (const Point<spacedim> &direction,
                       const Point<spacedim> &point_on_axis,
                       const double tolerance = 1e-10);

  virtual Point<spacedim>
  get_new_point(const Quadrature<spacedim> &quad) const;

protected:
  const Point<spacedim> direction;

  const Point<spacedim> point_on_axis;

private:
  FlatManifold<dim,spacedim> flat_manifold;

  double tolerance;
};


template <int dim, int spacedim=dim, int chartdim=dim>
class FunctionManifold : public ChartManifold<dim, spacedim, chartdim>
{
public:
  FunctionManifold(const Function<chartdim> &push_forward_function,
                   const Function<spacedim> &pull_back_function,
                   const Point<chartdim> periodicity=Point<chartdim>(),
                   const double tolerance=1e-10);

  FunctionManifold(const std::string push_forward_expression,
                   const std::string pull_back_expression,
                   const Point<chartdim> periodicity=Point<chartdim>(),
                   const typename FunctionParser<spacedim>::ConstMap = typename FunctionParser<spacedim>::ConstMap(),
                   const std::string chart_vars=FunctionParser<chartdim>::default_variable_names(),
                   const std::string space_vars=FunctionParser<spacedim>::default_variable_names(),
                   const double tolerance=1e-10);

  ~FunctionManifold();

  virtual Point<spacedim>
  push_forward(const Point<chartdim> &chart_point) const;

  virtual Point<chartdim>
  pull_back(const Point<spacedim> &space_point) const;

private:
  const typename FunctionParser<spacedim>::ConstMap const_map;

  SmartPointer<const Function<chartdim>,
               FunctionManifold<dim,spacedim,chartdim> > push_forward_function;

  SmartPointer<const Function<spacedim>,
               FunctionManifold<dim,spacedim,chartdim> > pull_back_function;

  const double tolerance;

  const bool owns_pointers;
};


DEAL_II_NAMESPACE_CLOSE

#endif