From: heltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Wed, 21 Aug 2013 15:59:11 +0000 (+0000)
Subject: Added manifold_lib.cc
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=16e1ea6421bf4213bee25554386cba04cbd169ec;p=dealii-svn.git

Added manifold_lib.cc

git-svn-id: https://svn.dealii.org/branches/branch_manifold_id@30375 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/examples/step-1/step-1.cc b/deal.II/examples/step-1/step-1.cc
index 09599c7935..6019451d59 100644
--- a/deal.II/examples/step-1/step-1.cc
+++ b/deal.II/examples/step-1/step-1.cc
@@ -73,9 +73,9 @@ void first_grid ()
   // Now we want to write a graphical representation of the mesh to an output
   // file. The GridOut class of deal.II can do that in a number of different
   // output formats; here, we choose encapsulated postscript (eps) format:
-  std::ofstream out ("grid-1.eps");
+  std::ofstream out ("grid-1.msh");
   GridOut grid_out;
-  grid_out.write_eps (triangulation, out);
+  grid_out.write_msh (triangulation, out);
 }
 
 
@@ -187,9 +187,9 @@ void second_grid ()
   // Finally, after these five iterations of refinement, we want to again
   // write the resulting mesh to a file, again in eps format. This works just
   // as above:
-  std::ofstream out ("grid-2.eps");
+  std::ofstream out ("grid-2.msh");
   GridOut grid_out;
-  grid_out.write_eps (triangulation, out);
+  grid_out.write_msh (triangulation, out);
 
 
   // At this point, all objects created in this function will be destroyed in
diff --git a/deal.II/include/deal.II/grid/manifold.h b/deal.II/include/deal.II/grid/manifold.h
index 8354314640..922d7d1cab 100644
--- a/deal.II/include/deal.II/grid/manifold.h
+++ b/deal.II/include/deal.II/grid/manifold.h
@@ -160,6 +160,57 @@ public:
 };
 
 
+/**
+ *   A chart of dimension dim, which is part of a Manifold<spacedim>.
+ *   This object specializes a Manifold of dimension dim embedded in a
+ *   manifold of dimension spacedim, for which you have an explicit
+ *   pull_back and push_forward transformations.
+ *
+ *   @ingroup manifold
+ *
+ *   @author Luca Heltai, 2013
+ */
+template <int dim, int spacedim>
+class ManifoldChart: public Manifold<spacedim>
+{
+public:
+    
+  /**
+   * Destructor. Does nothing here, but
+   * needs to be declared to make it
+   * virtual.
+   */
+  virtual ~ManifoldChart ();
+
+
+  /**
+   * Let the new point be the
+   * arithmetic mean of the two
+   * vertices of the line.
+   *
+   * Refer to the general
+   * documentation of this class
+   * and the documentation of the
+   * base class for more
+   * information.
+   */
+    virtual Point<spacedim>
+    get_new_point(const std::vector<Point<spacedim> > &surrounding_points,
+		  const std::vector<double> &weights) const;
+
+    virtual Point<dim>
+    pull_back(const Point<spacedim> space_point) const = 0;
+
+    virtual Point<spacedim>
+    push_forward(const Point<dim> chart_point) const = 0;
+    
+
+  private:
+    FlatManifold<dim> sub_manifold;
+    std::vector<Point<dim> > chart_points;
+};
+
+
 
 /* -------------- declaration of explicit specializations ------------- */
 
diff --git a/deal.II/include/deal.II/grid/manifold_lib.h b/deal.II/include/deal.II/grid/manifold_lib.h
index 7bd1cf5e29..2aab365c68 100644
--- a/deal.II/include/deal.II/grid/manifold_lib.h
+++ b/deal.II/include/deal.II/grid/manifold_lib.h
@@ -22,6 +22,28 @@
 #include <deal.II/grid/manifold.h>
 
 DEAL_II_NAMESPACE_OPEN
+
+
+/**
+ * Manifold description for a polar space coordinate system. Given any
+ * two points in space, tranform them in polar coordinate centered at
+ * center, compute intermediate points in polar coordinate systems,
+ * and returns its space representation.
+ */
+template <int spacedim>
+class PolarManifold : public FlatManifold<spacedim> 
+{
+  public:
+    PolarManifold(const Point<spacedim> center = Point<spacedim>());
+
+    virtual Point<spacedim>
+    get_new_point (const std::vector<Point<spacedim> > &surrounding_points,
+		   const std::vector<double> &weights) const;
+
+  private:
+    Point<spacedim> center;
+}
+
 //  
 //  /**
 //   * Boundary object for the hull of a cylinder.  In three dimensions,
diff --git a/deal.II/source/grid/manifold.cc b/deal.II/source/grid/manifold.cc
index 5135baf965..eff55fa9ce 100644
--- a/deal.II/source/grid/manifold.cc
+++ b/deal.II/source/grid/manifold.cc
@@ -65,6 +65,31 @@ get_new_point (const std::vector<Point<spacedim> > &surrounding_points,
   return p;
 }
 
+
+/* -------------------------- ManifoldChart --------------------- */
+
+template <int spacedim>
+ManifoldChart<spacedim>::~ManifoldChart ()
+{}
+
+
+template <int spacedim>
+Point<spacedim>
+FlatManifold<spacedim>::
+get_new_point (const std::vector<Point<spacedim> > &surrounding_points,
+	       const std::vector<double> &weights) const
+{
+  Assert(surrounding_points.size() == weights.size(),
+	 ExcDimensionMismatch(surrounding_points.size(), weights.size()));
+  
+  Point<spacedim> p;
+  for(unsigned int i=0; i<surrounding_points.size(); ++i)
+    p += surrounding_points[i]*weights[i];
+  
+  return p;
+}
+
+
 // explicit instantiations
 #include "manifold.inst"
 
diff --git a/deal.II/source/grid/manifold_lib.cc b/deal.II/source/grid/manifold_lib.cc
new file mode 100644
index 0000000000..16c8f61769
--- /dev/null
+++ b/deal.II/source/grid/manifold_lib.cc
@@ -0,0 +1,58 @@
+// ---------------------------------------------------------------------
+// $Id: manifold_lib.cc 30130 2013-07-23 13:01:18Z heltai $
+//
+// Copyright (C) 2013 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/grid/manifold_lib.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/base/tensor.h>
+#include <cmath>
+
+template <int spacedim>
+Point<spacedim>
+PolarManifold<spacedim>::PolarManifold(const Point<spacedim> center):
+  center(center)
+{}
+
+
+template <int spacedim>
+Point<spacedim>
+PolarManifold<spacedim>::get_new_point (const std::vector<Point<spacedim> > &sp,
+					const std::vector<double> &sw) const
+{
+  AssertDimension(sp.size(), sw.size());
+  
+  std::vector<Point<spacedim> > sp_(sp.size());
+  std::vector<double> sw_(sp.size());
+  
+  
+  Point<spacedim> middle = StraightBoundary<dim,spacedim>::get_new_point_on_line (line);
+
+  middle -= center;
+
+  double r=0;
+  if (compute_radius_automatically)
+    {
+      const Point<spacedim> vertex_relative = line->vertex(0) - center;
+      r = std::sqrt(vertex_relative.square());
+    }
+  else
+    r=radius;
+  // project to boundary
+  middle *= r / std::sqrt(middle.square());
+  middle += center;
+  return middle;
+}