From 9e8cb24a229281bde1a7204bcb153ea777ea6c07 Mon Sep 17 00:00:00 2001
From: Stefano Dominici <sfndmn@gmail.com>
Date: Sun, 21 Oct 2018 10:10:33 +0200
Subject: [PATCH] Data members of TensorProductManifold converted to
 std::unique_ptr<...>

Mimor modifications to:
 - Remove unnecessary parts.
---
 include/deal.II/grid/manifold_lib.h           | 15 ---------
 .../deal.II/grid/tensor_product_manifold.h    | 29 ++++++-----------
 tests/manifold/tensor_product_manifold_03.cc  | 32 +++++++++++++++----
 .../tensor_product_manifold_03.output         |  2 ++
 4 files changed, 37 insertions(+), 41 deletions(-)

diff --git a/include/deal.II/grid/manifold_lib.h b/include/deal.II/grid/manifold_lib.h
index a4498acaf3..f6da483202 100644
--- a/include/deal.II/grid/manifold_lib.h
+++ b/include/deal.II/grid/manifold_lib.h
@@ -73,11 +73,6 @@ public:
    */
   PolarManifold(const Point<spacedim> center = Point<spacedim>());
 
-  /**
-   * Virtual destructor
-   */
-  virtual ~PolarManifold() = default;
-
   /**
    * Make a clone of this Manifold object.
    */
@@ -404,10 +399,6 @@ public:
                       const Point<spacedim> &    point_on_axis,
                       const double               tolerance = 1e-10);
 
-  /**
-   * Virtual destructor.
-   */
-  virtual ~CylindricalManifold() = default;
   /**
    * Make a clone of this Manifold object.
    */
@@ -507,7 +498,6 @@ public:
     const Tensor<1, chartdim> &periodicity = Tensor<1, chartdim>(),
     const double               tolerance   = 1e-10);
 
-
   /**
    * Expressions constructor. Takes the expressions of the push_forward
    * function of spacedim components, and of the pull_back function of @p
@@ -677,11 +667,6 @@ public:
    */
   TorusManifold(const double R, const double r);
 
-  /**
-   * Virtual destructor
-   */
-  virtual ~TorusManifold() = default;
-
   /**
    * Make a clone of this Manifold object.
    */
diff --git a/include/deal.II/grid/tensor_product_manifold.h b/include/deal.II/grid/tensor_product_manifold.h
index 73c6561132..5f371c3437 100644
--- a/include/deal.II/grid/tensor_product_manifold.h
+++ b/include/deal.II/grid/tensor_product_manifold.h
@@ -21,6 +21,7 @@
 #include <deal.II/base/point.h>
 #include <deal.II/base/std_cxx14/memory.h>
 #include <deal.II/base/subscriptor.h>
+#include <deal.II/base/utilities.h>
 
 #include <deal.II/grid/manifold.h>
 
@@ -89,11 +90,6 @@ public:
     const ChartManifold<dim_A, spacedim_A, chartdim_A> &manifold_A,
     const ChartManifold<dim_B, spacedim_B, chartdim_B> &manifold_B);
 
-  /**
-   * Virtual destructor
-   */
-  virtual ~TensorProductManifold() = default;
-
   /**
    * Clone this manifold.
    */
@@ -119,10 +115,10 @@ public:
   push_forward_gradient(const Point<chartdim> &chart_point) const override;
 
 private:
-  std::shared_ptr<const ChartManifold<dim_A, spacedim_A, chartdim_A>>
+  std::unique_ptr<const ChartManifold<dim_A, spacedim_A, chartdim_A>>
     manifold_A;
 
-  std::shared_ptr<const ChartManifold<dim_B, spacedim_B, chartdim_B>>
+  std::unique_ptr<const ChartManifold<dim_B, spacedim_B, chartdim_B>>
     manifold_B;
 };
 
@@ -196,18 +192,13 @@ TensorProductManifold<dim,
       internal::TensorProductManifoldImplementation::concat(
         manifold_A.get_periodicity(),
         manifold_B.get_periodicity()))
-{
-  std::shared_ptr<Manifold<dim_A, spacedim_A>> tmpA(
-    std::move(manifold_A.clone()));
-  std::shared_ptr<Manifold<dim_B, spacedim_B>> tmpB(
-    std::move(manifold_B.clone()));
-  this->manifold_A =
-    std::static_pointer_cast<ChartManifold<dim_A, spacedim_A, chartdim_A>>(
-      tmpA);
-  this->manifold_B =
-    std::static_pointer_cast<ChartManifold<dim_B, spacedim_B, chartdim_B>>(
-      tmpB);
-}
+  , manifold_A(Utilities::dynamic_unique_cast<
+               ChartManifold<dim_A, spacedim_A, chartdim_A>,
+               Manifold<dim_A, spacedim_A>>(manifold_A.clone()))
+  , manifold_B(Utilities::dynamic_unique_cast<
+               ChartManifold<dim_B, spacedim_B, chartdim_B>,
+               Manifold<dim_B, spacedim_B>>(manifold_B.clone()))
+{}
 
 template <int dim,
           int dim_A,
diff --git a/tests/manifold/tensor_product_manifold_03.cc b/tests/manifold/tensor_product_manifold_03.cc
index ecfea9cccb..6d96b0b4af 100644
--- a/tests/manifold/tensor_product_manifold_03.cc
+++ b/tests/manifold/tensor_product_manifold_03.cc
@@ -26,19 +26,28 @@
 #include "../tests.h"
 
 
-void
+bool
 test()
 {
   Triangulation<3, 3> tria;
   GridGenerator::hyper_cube(tria);
   {
-    FunctionManifold<1, 1>                     F("x", "x");
-    PolarManifold<2, 2>                        G(Point<2>(0.5, 0.5));
-    TensorProductManifold<3, 2, 2, 2, 1, 1, 1> manifold(G, F);
-    tria.set_all_manifold_ids(0);
-    tria.set_manifold(0, manifold);
+    FunctionManifold<1, 1> F("x", "x");
+    PolarManifold<2, 2>    G(Point<2>(0.5, 0.5));
+    try
+      {
+        TensorProductManifold<3, 2, 2, 2, 1, 1, 1> manifold(G, F);
+        tria.set_all_manifold_ids(0);
+        tria.set_manifold(0, manifold);
+      }
+    catch (std::bad_cast &e)
+      {
+        (void)e;
+        return false;
+      }
   }
   tria.refine_global(1);
+  return true;
 }
 
 
@@ -46,7 +55,16 @@ test()
 int
 main()
 {
-  test();
+  initlog();
+
+  if (test())
+    {
+      deallog << "OK" << std::endl;
+    }
+  else
+    {
+      deallog << "FAILED" << std::endl;
+    }
 
   return 0;
 }
diff --git a/tests/manifold/tensor_product_manifold_03.output b/tests/manifold/tensor_product_manifold_03.output
index e69de29bb2..0fd8fc12f0 100644
--- a/tests/manifold/tensor_product_manifold_03.output
+++ b/tests/manifold/tensor_product_manifold_03.output
@@ -0,0 +1,2 @@
+
+DEAL::OK
-- 
2.39.5