--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 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__composition_manifold_h
+#define dealii__composition_manifold_h
+
+
+/*---------------------------- composition_manifold.h ------------*/
+
+#include <deal.II/base/config.h>
+#include <deal.II/base/subscriptor.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/base/thread_management.h>
+#include <deal.II/base/point.h>
+#include <deal.II/base/derivative_form.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/manifold.h>
+
+DEAL_II_NAMESPACE_OPEN
+
+/**
+ * CompositionManifold. Take two ChartManifold objects, and make
+ * their composition. The CompositionManifold object is a
+ * ChartManifold going from the chart of the first ChartManifold to
+ * the embedding space of the second ChartManifold. If the first
+ * ChartManifold is periodic, so is the resulting ChartManifold, with
+ * the same periodicity. Periodicity on the second ChartManifold is not
+ * allowed, and the constructor will throw an axception if the second
+ * Manifold is periodic.
+ *
+ * This class only works for dim <= chartdim <= intermediate_spacedim
+ * <= spacedim. If you try to instantiate anything different, an
+ * Exception will be thrown in one of the ChartManifold classes that
+ * violates this condition.
+ *
+ * Given the ChartManifold F and the ChartManifold G, this class
+ * represents the composition of G after F.
+ *
+ * The template parameters have the following meaning:
+ *
+ * @tparam dim The dimension of the resulting ChartManifold
+ * @tparam spacedim The space dimension of the resulting ChartManifold
+ * @tparam chartdim The chart dimension of the resulting ChartManifold
+ * @tparam intermediate_dim The space dimension of the first ChartManifold
+ * @tparam dim1 The dimension of the first ChartManifold, which coincides also
+ * with the chart dimension of the second ChartManifold
+ * @tparam dim2 The dimension of the second ChartManifold
+ *
+ * @ingroup manifold
+ * @author Luca Heltai, Timo Heister, 2016
+ */
+template <int dim, int spacedim=dim, int chartdim=dim,
+ int intermediate_dim=dim, int dim1=dim, int dim2=dim>
+class CompositionManifold : public ChartManifold<dim, spacedim, chartdim>
+{
+public:
+
+ /**
+ * Construct the composition of the two given manifolds.
+ */
+ CompositionManifold(const ChartManifold<dim1, intermediate_dim, chartdim> &F,
+ const ChartManifold<dim2, spacedim, intermediate_dim> &G);
+
+ /**
+ * Pull back the given point in spacedim to the Euclidean chartdim
+ * dimensional space. This function calls the pull_back() function
+ * of G, and then the pull_back() function of F.
+ */
+ virtual
+ Point<chartdim>
+ pull_back(const Point<spacedim> &space_point) const;
+
+
+ /**
+ * Push forward the chartdim dimensional point to a spacedim
+ * Euclidean point. The function calls first the push_forward() of
+ * F, and then the push_foward() of G.
+ */
+ virtual
+ Point<spacedim>
+ push_forward(const Point<chartdim> &chart_point) const;
+
+ /**
+ * Return the derivative of the composition of G after F.
+ */
+ virtual
+ DerivativeForm<1,chartdim,spacedim>
+ push_forward_gradient(const Point<chartdim> &chart_point) const;
+
+private:
+ /**
+ * The first ChartManifold.
+ */
+ SmartPointer<const ChartManifold<dim1, intermediate_dim, chartdim>,
+ CompositionManifold<dim,spacedim,chartdim,dim1,dim2,intermediate_dim> > F;
+
+
+ /**
+ * The second ChartManifold.
+ */
+ SmartPointer<const ChartManifold<dim2, spacedim, intermediate_dim>,
+ CompositionManifold<dim,spacedim,chartdim,dim1,dim2,intermediate_dim> > G;
+};
+
+
+/*------------------Template Implementations------------------------*/
+
+template <int dim, int spacedim, int chartdim, int intermediate_dim,
+ int dim1, int dim2>
+CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::CompositionManifold
+(const ChartManifold<dim1, intermediate_dim, chartdim> &F,
+ const ChartManifold<dim2, spacedim, intermediate_dim> &G) :
+ ChartManifold<dim,spacedim,chartdim>(F.get_periodicity()),
+ F(&F),
+ G(&G)
+{
+ // We don't know what to do with a periodicity in the second manifold, so
+ // throw an assertion if the second manifold is periodic
+ Assert(G.get_periodicity().norm() == 0.0,
+ ExcMessage("The second manifold cannot be periodic."));
+}
+
+
+
+template <int dim, int spacedim, int chartdim, int intermediate_dim,
+ int dim1, int dim2>
+Point<chartdim>
+CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::pull_back(const Point<spacedim> &space_point) const
+{
+ return F->pull_back(G->pull_back(space_point));
+}
+
+
+
+template <int dim, int spacedim, int chartdim, int intermediate_dim,
+ int dim1, int dim2>
+Point<spacedim>
+CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::push_forward(const Point<chartdim> &chart_point) const
+{
+ return G->push_forward(F->push_forward(chart_point));
+}
+
+
+
+template <int dim, int spacedim, int chartdim, int intermediate_dim,
+ int dim1, int dim2>
+DerivativeForm<1,chartdim,spacedim>
+CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::push_forward_gradient(const Point<chartdim> &chart_point) const
+{
+ DerivativeForm<1,chartdim,intermediate_dim> DF =
+ F->push_forward_gradient(chart_point);
+
+ DerivativeForm<1,intermediate_dim,spacedim> DG =
+ G->push_forward_gradient(F->push_forward(chart_point));
+
+ DerivativeForm<1,chartdim,spacedim> DF_DG;
+
+ for (unsigned int d=0; d<spacedim; ++d)
+ for (unsigned int c=0; c<chartdim; ++c)
+ for (unsigned int s=0; s<intermediate_dim; ++s)
+ DF_DG[d][c] += DG[d][s]*DF[s][c];
+
+ return DF_DG;
+}
+
+
+DEAL_II_NAMESPACE_CLOSE
+
+#endif
+++ /dev/null
-// ---------------------------------------------------------------------
-//
-// Copyright (C) 1998 - 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.
-//
-// ---------------------------------------------------------------------
-
-#include <deal.II/base/tensor.h>
-#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_accessor.h>
-#include <deal.II/grid/manifold_tools.h>
-#include <deal.II/fe/fe_q.h>
-#include <cmath>
-
-DEAL_II_NAMESPACE_OPEN
-
-
-/* -------------------------- CompositionManifold --------------------- */
-
-template <int dim, int spacedim, int chartdim, int intermediate_dim,
- int dim1, int dim2>
-CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::CompositionManifold
-(const ChartManifold<dim1, intermediate_dim, chartdim> &F,
- const ChartManifold<dim2, spacedim, intermediate_dim> &G) :
- ChartManifold<dim,spacedim,chartdim>(F.get_periodicity()),
- F(&F),
- G(&G)
-{
- // We don't know what to do with a periodicity in the second manifold, so
- // throw an assertion if the second manifold is periodic
- Assert(G.get_periodicity().norm() == 0.0,
- ExcMessage("The second manifold cannot be periodic."));
-}
-
-
-
-template <int dim, int spacedim, int chartdim, int intermediate_dim,
- int dim1, int dim2>
-Point<chartdim>
-CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::pull_back(const Point<spacedim> &space_point) const
-{
- return F->pull_back(G->pull_back(space_point));
-}
-
-
-
-template <int dim, int spacedim, int chartdim, int intermediate_dim,
- int dim1, int dim2>
-Point<spacedim>
-CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::push_forward(const Point<chartdim> &chart_point) const
-{
- return G->push_forward(F->push_forward(chart_point));
-}
-
-
-
-template <int dim, int spacedim, int chartdim, int intermediate_dim,
- int dim1, int dim2>
-DerivativeForm<1,chartdim,spacedim>
-CompositionManifold<dim,spacedim,chartdim,intermediate_dim,dim1,dim2>::push_forward_gradient(const Point<chartdim> &chart_point) const
-{
- DerivativeForm<1,chartdim,intermediate_dim> DF =
- F->push_forward_gradient(chart_point);
-
- DerivativeForm<1,intermediate_dim,spacedim> DG =
- G->push_forward_gradient(F->push_forward(chart_point));
-
- DerivativeForm<1,chartdim,spacedim> DF_DG;
-
- for (unsigned int d=0; d<spacedim; ++d)
- for (unsigned int c=0; c<chartdim; ++c)
- for (unsigned int s=0; s<intermediate_dim; ++s)
- DF_DG[d][c] += DG[d][s]*DF[s][c];
-
- return DF_DG;
-}
-
-
-
-// explicit instantiations
-#include "manifold_tools.inst"
-
-DEAL_II_NAMESPACE_CLOSE
-
projects / dealii.git / commitdiff