+++ /dev/null
-//-----------------------------------------------------------
-//
-// Copyright (C) 2014 - 2015 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.
-//
-//-----------------------------------------------------------
-
-#include "../tests.h"
-#include <fstream>
-#include <deal.II/base/logstream.h>
-
-#include <deal.II/opencascade/utilities.h>
-#include <deal.II/opencascade/boundary_lib.h>
-#include <deal.II/grid/tria.h>
-#include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/grid_generator.h>
-
-#include <TopTools.hxx>
-#include <TopoDS_Shape.hxx>
-#include <Standard_Stream.hxx>
-
-// Create a Triangulation, interpolate its boundary points to a smooth
-// BSpline, and use that as an arlength Boundary Descriptor.
-
-using namespace OpenCASCADE;
-
-int main ()
-{
- std::ofstream logfile("output");
-
- // Create a bspline passign through the points
- std::vector<Point<3> > pts;
- pts.push_back(Point<3>(0,0,0));
- pts.push_back(Point<3>(0,1,0));
- pts.push_back(Point<3>(1,1,0));
- pts.push_back(Point<3>(1,0,0));
-
- TopoDS_Edge edge = interpolation_curve(pts);
- ArclengthProjectionLineManifold<1,3> manifold(edge);
-
- Triangulation<1,3> tria;
- GridGenerator::hyper_cube(tria);
-
- tria.begin_active()->set_all_manifold_ids(1);
- tria.set_manifold(1, manifold);
-
- tria.refine_global(4);
- GridOut gridout;
- gridout.write_msh(tria, logfile);
-
- return 0;
-}
-
+++ /dev/null
-$NOD
-17
-1 0 0 0
-2 1 0 0
-3 0.5 1.125 0
-4 -0.213053 0.770812 0
-5 1.21305 0.770812 0
-6 -0.209733 0.361552 0
-7 0.090632 1.04613 0
-8 0.909368 1.04613 0
-9 1.20973 0.361552 0
-10 -0.118931 0.172725 0
-11 -0.251438 0.566341 0
-12 -0.0875601 0.936891 0
-13 0.291313 1.106 0
-14 0.708687 1.106 0
-15 1.08756 0.936891 0
-16 1.25144 0.566341 0
-17 1.11893 0.172725 0
-$ENDNOD
-$ELM
-16
-1 1 0 0 2 1 10
-2 1 0 0 2 10 6
-3 1 0 0 2 6 11
-4 1 0 0 2 11 4
-5 1 0 0 2 4 12
-6 1 0 0 2 12 7
-7 1 0 0 2 7 13
-8 1 0 0 2 13 3
-9 1 0 0 2 3 14
-10 1 0 0 2 14 8
-11 1 0 0 2 8 15
-12 1 0 0 2 15 5
-13 1 0 0 2 5 16
-14 1 0 0 2 16 9
-15 1 0 0 2 9 17
-16 1 0 0 2 17 2
-$ENDELM
+++ /dev/null
-//-----------------------------------------------------------
-//
-// Copyright (C) 2014 - 2015 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 sphere, intersect it with a plane, and use the resulting
-// object as a boundary.
-
-#include "../tests.h"
-
-#include <deal.II/opencascade/utilities.h>
-#include <deal.II/opencascade/boundary_lib.h>
-
-#include <fstream>
-#include <deal.II/base/logstream.h>
-#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_accessor.h>
-#include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
-
-#include <gp_Pnt.hxx>
-#include <gp_Dir.hxx>
-#include <gp_Ax2.hxx>
-#include <GC_MakeCircle.hxx>
-#include <BRepBuilderAPI_MakeEdge.hxx>
-#include <BRepBuilderAPI_MakeFace.hxx>
-#include <BRepBuilderAPI_MakeWire.hxx>
-#include <BRepPrimAPI_MakeSphere.hxx>
-#include <TopExp_Explorer.hxx>
-#include <TopoDS_Wire.hxx>
-#include <TopoDS_Edge.hxx>
-#include <TopoDS_Face.hxx>
-#include <TopoDS.hxx>
-using namespace OpenCASCADE;
-
-int main()
-{
- std::ofstream logfile("output");
- deallog.attach(logfile);
-
- // Create a sphere
- gp_Pnt center(.5,.5,.5);
- Standard_Real radius(Point<3>().distance(point(center)));
-
- TopoDS_Face face = BRepPrimAPI_MakeSphere(center, radius);
-
- // Intersect the sphere with the plane z = 0
- TopoDS_Shape intersect = intersect_plane(face, 0, 0, 1, 0);
- // This shape should be the standard circle. Only one edge should be
- // in intersect: our edge.
- TopExp_Explorer exp(intersect, TopAbs_EDGE);
- TopoDS_Edge edge = TopoDS::Edge(exp.Current());
-
- ArclengthProjectionLineManifold<2,3> boundary_line(edge);
-
- // The unit square.
- Triangulation<2,3> tria;
- GridGenerator::hyper_cube(tria);
-
- // Set the exterior boundary
- tria.set_manifold(1, boundary_line);
-
- // This is here to ignore the points created in the interior of the
- // face.
- tria.begin()->set_all_manifold_ids(1);
- tria.begin()->set_manifold_id(0);
-
- // We refine twice, and expect the outer points to end up on a
- // smooth curve interpolating the square vertices.
- tria.refine_global(2);
-
-
- // You can open the generated file with gmsh.
- GridOut gridout;
- gridout.write_msh (tria, logfile);
- return 0;
-}
+++ /dev/null
-
-$NOD
-25
-1 0.00000 0.00000 0.00000
-2 1.00000 0.00000 0.00000
-3 0.00000 1.00000 0.00000
-4 1.00000 1.00000 0.00000
-5 0.500000 -0.207107 0.00000
-6 -0.207107 0.500000 0.00000
-7 1.20711 0.500000 0.00000
-8 0.500000 1.20711 0.00000
-9 0.500000 0.500000 0.00000
-10 0.229402 -0.153281 0.00000
-11 0.770598 -0.153281 0.00000
-12 -0.153281 0.229402 0.00000
-13 -0.153281 0.770598 0.00000
-14 1.15328 0.229402 0.00000
-15 1.15328 0.770598 0.00000
-16 0.229402 1.15328 0.00000
-17 0.770598 1.15328 0.00000
-18 0.500000 0.146447 0.00000
-19 0.500000 0.853553 0.00000
-20 0.146447 0.500000 0.00000
-21 0.853553 0.500000 0.00000
-22 0.189433 0.189433 0.00000
-23 0.810567 0.189433 0.00000
-24 0.189433 0.810567 0.00000
-25 0.810567 0.810567 0.00000
-$ENDNOD
-$ELM
-16
-1 3 0 0 4 1 10 22 12
-2 3 0 0 4 10 5 18 22
-3 3 0 0 4 12 22 20 6
-4 3 0 0 4 22 18 9 20
-5 3 0 0 4 5 11 23 18
-6 3 0 0 4 11 2 14 23
-7 3 0 0 4 18 23 21 9
-8 3 0 0 4 23 14 7 21
-9 3 0 0 4 6 20 24 13
-10 3 0 0 4 20 9 19 24
-11 3 0 0 4 13 24 16 3
-12 3 0 0 4 24 19 8 16
-13 3 0 0 4 9 21 25 19
-14 3 0 0 4 21 7 15 25
-15 3 0 0 4 19 25 17 8
-16 3 0 0 4 25 15 4 17
-$ENDELM