--- /dev/null
+//-----------------------------------------------------------
+//
+// Copyright (C) 2017 by the deal.II authors
+//
+// This file is subject to LGPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//-----------------------------------------------------------
+
+// Create a NURBS with a quarter of a disc, make a Manifold out of it, create a
+// Triangulation with it, attach the manifold, and make a few refinements
+
+#include "../tests.h"
+
+#include <deal.II/opencascade/boundary_lib.h>
+#include <deal.II/opencascade/utilities.h>
+#include <BRepPrimAPI_MakeRevol.hxx>
+
+#include <deal.II/grid/grid_out.h>
+
+
+using namespace OpenCASCADE;
+
+int main ()
+{
+ std::ofstream logfile("output");
+
+ // Create an edge, and make a face out of it by revolving along the z
+ // axis.
+ std::vector<Point<2> > ps;
+ ps.push_back(Point<2>(0, 1));
+ ps.push_back(Point<2>(0, 2));
+ auto edge = interpolation_curve(ps);
+
+ BRepPrimAPI_MakeRevol::BRepPrimAPI_MakeRevol revol(edge, gp_Ax1(gp_Pnt(0,0,0), gp_Dir(0,0,1)), -numbers::PI/2);
+ revol.Build();
+ auto sh = revol.Shape();
+
+ std::vector<TopoDS_Face> faces;
+ std::vector<TopoDS_Edge> edges;
+ std::vector<TopoDS_Vertex> vertices;
+
+ extract_geometrical_shapes(sh, faces, edges, vertices);
+
+ // Create a boundary projector on the first face.
+ NURBSPatchManifold<2,2> manifold(faces[0]);
+
+ // Create a Triangulation with a single cell
+ Triangulation<2,2> tria;
+ create_triangulation(faces[0], tria);
+
+ // Set the manifold
+ tria.set_manifold(1, manifold);
+ tria.set_all_manifold_ids(1);
+
+ tria.refine_global(3);
+
+ // You can open the generated file with gmsh
+ GridOut gridout;
+ gridout.write_msh (tria, logfile);
+}
+
projects / dealii.git / commitdiff