20 #ifdef DEAL_II_WITH_OPENCASCADE
26 # include <IGESControl_Controller.hxx>
27 # include <IGESControl_Reader.hxx>
28 # include <IGESControl_Writer.hxx>
29 # include <STEPControl_Controller.hxx>
30 # include <STEPControl_Reader.hxx>
31 # include <STEPControl_Writer.hxx>
32 # include <TopExp_Explorer.hxx>
33 # include <TopoDS.hxx>
34 # include <TopoDS_Edge.hxx>
35 # include <TopoDS_Face.hxx>
36 # include <TopoDS_Shape.hxx>
40 # if DEAL_II_OPENCASCADE_VERSION_GTE(7, 0, 0)
41 # include <Standard_Transient.hxx>
43 # include <Handle_Standard_Transient.hxx>
46 # include <BRepAdaptor_Curve.hxx>
47 # if DEAL_II_OPENCASCADE_VERSION_GTE(7, 6, 0)
48 # include <BRepAlgoAPI_Section.hxx>
50 # include <BRepAdaptor_HCompCurve.hxx>
51 # include <BRepAdaptor_HCurve.hxx>
52 # include <BRepAlgo_Section.hxx>
54 # include <BRepAdaptor_Surface.hxx>
55 # include <BRepBndLib.hxx>
56 # include <BRepBuilderAPI_MakeEdge.hxx>
57 # include <BRepBuilderAPI_Sewing.hxx>
58 # include <BRepBuilderAPI_Transform.hxx>
59 # include <BRepMesh_IncrementalMesh.hxx>
60 # include <BRepTools.hxx>
61 # include <BRep_Builder.hxx>
62 # include <GCPnts_AbscissaPoint.hxx>
63 # include <GeomAPI_Interpolate.hxx>
64 # include <GeomAPI_ProjectPointOnCurve.hxx>
65 # include <GeomAPI_ProjectPointOnSurf.hxx>
66 # include <GeomConvert_CompCurveToBSplineCurve.hxx>
67 # include <GeomLProp_SLProps.hxx>
68 # include <Geom_BoundedCurve.hxx>
69 # include <Geom_Plane.hxx>
70 # include <IntCurvesFace_ShapeIntersector.hxx>
71 # include <Poly_Triangulation.hxx>
72 # include <ShapeAnalysis_Surface.hxx>
73 # include <StlAPI_Reader.hxx>
74 # include <StlAPI_Writer.hxx>
75 # include <TColStd_HSequenceOfTransient.hxx>
76 # include <TColStd_SequenceOfTransient.hxx>
77 # include <TColgp_HArray1OfPnt.hxx>
78 # include <TopLoc_Location.hxx>
79 # include <gp_Lin.hxx>
80 # include <gp_Pnt.hxx>
81 # include <gp_Vec.hxx>
91 std::tuple<unsigned int, unsigned int, unsigned int>
95 unsigned int n_faces = 0, n_edges = 0, n_vertices = 0;
96 for (exp.Init(shape, TopAbs_FACE); exp.More(); exp.Next(), ++n_faces)
99 for (exp.Init(shape, TopAbs_EDGE); exp.More(); exp.Next(), ++n_edges)
102 for (exp.Init(shape, TopAbs_VERTEX); exp.More(); exp.Next(), ++n_vertices)
105 return std::tuple<unsigned int, unsigned int, unsigned int>(n_faces,
112 std::vector<TopoDS_Face> &faces,
113 std::vector<TopoDS_Edge> &edges,
114 std::vector<TopoDS_Vertex> &
vertices)
121 for (exp.Init(shape, TopAbs_FACE); exp.More(); exp.Next())
123 faces.push_back(TopoDS::Face(exp.Current()));
125 for (exp.Init(shape, TopAbs_EDGE); exp.More(); exp.Next())
127 edges.push_back(TopoDS::Edge(exp.Current()));
129 for (exp.Init(shape, TopAbs_VERTEX); exp.More(); exp.Next())
138 std::vector<TopoDS_Compound> &compounds,
139 std::vector<TopoDS_CompSolid> &compsolids,
140 std::vector<TopoDS_Solid> &solids,
141 std::vector<TopoDS_Shell> &shells,
142 std::vector<TopoDS_Wire> &wires)
145 compsolids.resize(0);
151 for (exp.Init(shape, TopAbs_COMPOUND); exp.More(); exp.Next())
153 compounds.push_back(TopoDS::Compound(exp.Current()));
155 for (exp.Init(shape, TopAbs_COMPSOLID); exp.More(); exp.Next())
157 compsolids.push_back(TopoDS::CompSolid(exp.Current()));
159 for (exp.Init(shape, TopAbs_SOLID); exp.More(); exp.Next())
161 solids.push_back(TopoDS::Solid(exp.Current()));
163 for (exp.Init(shape, TopAbs_SHELL); exp.More(); exp.Next())
165 shells.push_back(TopoDS::Shell(exp.Current()));
167 for (exp.Init(shape, TopAbs_WIRE); exp.More(); exp.Next())
169 wires.push_back(TopoDS::Wire(exp.Current()));
173 template <
int spacedim>
180 return gp_Pnt(p[0], 0, 0);
182 return gp_Pnt(p[0], p[1], 0);
184 return gp_Pnt(p[0], p[1], p[2]);
190 template <
int spacedim>
192 point(
const gp_Pnt &p,
const double tolerance)
198 Assert(std::abs(p.Y()) <= tolerance,
200 "Cannot convert OpenCASCADE point to 1d if p.Y() != 0."));
201 Assert(std::abs(p.Z()) <= tolerance,
203 "Cannot convert OpenCASCADE point to 1d if p.Z() != 0."));
206 Assert(std::abs(p.Z()) <= tolerance,
208 "Cannot convert OpenCASCADE point to 2d if p.Z() != 0."));
222 const double tolerance)
224 const double rel_tol =
226 if (direction.
norm() > 0.0)
227 return (p1 * direction < p2 * direction - rel_tol);
229 for (
int d = dim;
d >= 0; --
d)
230 if (p1[
d] < p2[
d] - rel_tol)
232 else if (p2[
d] < p1[
d] - rel_tol)
242 read_IGES(
const std::string &filename,
const double scale_factor)
244 IGESControl_Reader reader;
245 IFSelect_ReturnStatus stat;
246 stat = reader.ReadFile(filename.c_str());
249 Standard_Boolean failsonly = Standard_False;
250 IFSelect_PrintCount mode = IFSelect_ItemsByEntity;
251 reader.PrintCheckLoad(failsonly, mode);
253 Standard_Integer nRoots = reader.TransferRoots();
262 scale.SetScale(Origin, scale_factor);
264 TopoDS_Shape sh = reader.OneShape();
265 BRepBuilderAPI_Transform trans(sh,
scale);
267 return trans.Shape();
271 write_IGES(
const TopoDS_Shape &shape,
const std::string &filename)
273 IGESControl_Controller::Init();
274 IGESControl_Writer ICW(
"MM", 0);
275 Standard_Boolean ok = ICW.AddShape(shape);
278 Standard_Boolean OK = ICW.Write(filename.c_str());
286 StlAPI_Reader reader;
288 reader.Read(shape, filename.c_str());
295 const std::string &filename,
296 const double deflection,
297 const bool sew_different_faces,
298 const double sewer_tolerance,
299 const bool is_relative,
300 const double angular_deflection,
301 const bool in_parallel)
304 std::vector<TopoDS_Vertex>
vertices;
305 std::vector<TopoDS_Edge> edges;
306 std::vector<TopoDS_Face> faces;
308 const bool mesh_is_present =
309 std::none_of(faces.begin(), faces.end(), [&Loc](
const TopoDS_Face &face) {
310 Handle(Poly_Triangulation) theTriangulation =
311 BRep_Tool::Triangulation(face, Loc);
312 return theTriangulation.IsNull();
314 TopoDS_Shape shape_to_be_written = shape;
315 if (!mesh_is_present)
317 if (sew_different_faces)
319 BRepBuilderAPI_Sewing sewer(sewer_tolerance);
320 sewer.Add(shape_to_be_written);
322 shape_to_be_written = sewer.SewedShape();
325 shape_to_be_written = shape;
329 BRepMesh_IncrementalMesh mesh_im(shape_to_be_written,
336 StlAPI_Writer writer;
338 # if DEAL_II_OPENCASCADE_VERSION_GTE(6, 9, 0)
340 const auto error = writer.Write(shape_to_be_written, filename.c_str());
343 # if !DEAL_II_OPENCASCADE_VERSION_GTE(7, 2, 0)
354 writer.Write(shape_to_be_written, filename.c_str());
359 read_STEP(
const std::string &filename,
const double scale_factor)
361 STEPControl_Reader reader;
362 IFSelect_ReturnStatus stat;
363 stat = reader.ReadFile(filename.c_str());
366 Standard_Boolean failsonly = Standard_False;
367 IFSelect_PrintCount mode = IFSelect_ItemsByEntity;
368 reader.PrintCheckLoad(failsonly, mode);
370 Standard_Integer nRoots = reader.TransferRoots();
379 scale.SetScale(Origin, scale_factor);
381 TopoDS_Shape sh = reader.OneShape();
382 BRepBuilderAPI_Transform trans(sh,
scale);
384 return trans.Shape();
388 write_STEP(
const TopoDS_Shape &shape,
const std::string &filename)
390 STEPControl_Controller::Init();
391 STEPControl_Writer SCW;
392 IFSelect_ReturnStatus status;
393 status = SCW.Transfer(shape, STEPControl_AsIs);
395 ExcMessage(
"Failed to add shape to STEP controller."));
397 status = SCW.Write(filename.c_str());
400 ExcMessage(
"Failed to write translated shape to STEP file."));
406 double tolerance = 0.0;
408 std::vector<TopoDS_Face> faces;
409 std::vector<TopoDS_Edge> edges;
410 std::vector<TopoDS_Vertex>
vertices;
415 tolerance = std::fmax(tolerance, BRep_Tool::Tolerance(vertex));
417 for (
const auto &edge : edges)
418 tolerance = std::fmax(tolerance, BRep_Tool::Tolerance(edge));
420 for (
const auto &face : faces)
421 tolerance = std::fmax(tolerance, BRep_Tool::Tolerance(face));
435 Handle(Geom_Plane) plane =
new Geom_Plane(c_x, c_y, c_z, c);
436 # if DEAL_II_OPENCASCADE_VERSION_GTE(7, 6, 0)
437 BRepAlgoAPI_Section section(in_shape, plane);
439 BRepAlgo_Section section(in_shape, plane);
441 TopoDS_Shape edges = section.Shape();
446 join_edges(
const TopoDS_Shape &in_shape,
const double tolerance)
448 TopoDS_Edge out_shape;
449 const TopoDS_Shape &edges = in_shape;
450 std::vector<Handle_Geom_BoundedCurve> intersections;
454 gp_Pnt PIn(0.0, 0.0, 0.0);
455 gp_Pnt PFin(0.0, 0.0, 0.0);
456 gp_Pnt PMid(0.0, 0.0, 0.0);
457 TopExp_Explorer edgeExplorer(edges, TopAbs_EDGE);
459 while (edgeExplorer.More())
461 edge = TopoDS::Edge(edgeExplorer.Current());
462 Handle(Geom_Curve) curve = BRep_Tool::Curve(edge,
L, First, Last);
463 intersections.push_back(Handle(Geom_BoundedCurve)::DownCast(curve));
469 unsigned int numIntersEdges = intersections.size();
472 GeomConvert_CompCurveToBSplineCurve convert_bspline(intersections[0]);
474 bool check =
false, one_added =
true, one_failed =
true;
475 std::vector<bool> added(numIntersEdges,
false);
477 while (one_added ==
true)
481 for (
unsigned int i = 1; i < numIntersEdges; ++i)
482 if (added[i] ==
false)
484 Handle(Geom_Curve) curve = intersections[i];
485 Handle(Geom_BoundedCurve) bcurve =
486 Handle(Geom_BoundedCurve)::DownCast(curve);
487 check = convert_bspline.Add(bcurve, tolerance,
false,
true, 0);
492 Handle(Geom_BoundedCurve) bcurve =
493 Handle(Geom_BoundedCurve)::DownCast(curve);
495 convert_bspline.Add(bcurve, tolerance,
false,
true, 0);
497 one_failed = one_failed || (
check ==
false);
498 one_added = one_added || (
check ==
true);
503 Assert(one_failed ==
false,
504 ExcMessage(
"Joining some of the Edges failed."));
506 Handle(Geom_Curve) bspline = convert_bspline.BSplineCurve();
508 out_shape = BRepBuilderAPI_MakeEdge(bspline);
517 const double tolerance)
521 gp_Pnt P0 =
point(origin);
524 dim > 1 ? direction[1] : 0,
525 dim > 2 ? direction[2] : 0));
529 gp_Pnt Pproj(0.0, 0.0, 0.0);
533 IntCurvesFace_ShapeIntersector Inters;
534 Inters.Load(in_shape, tolerance);
538 Inters.Perform(line, -RealLast(), +RealLast());
541 double minDistance = 1e7;
543 for (
int i = 0; i < Inters.NbPnt(); ++i)
545 const double distance =
point(origin).Distance(Inters.Pnt(i + 1));
548 if (distance < minDistance)
550 minDistance = distance;
551 result = point<dim>(Inters.Pnt(i + 1));
563 const double tolerance)
565 unsigned int n_vertices = curve_points.size();
567 if (direction * direction > 0)
569 std::sort(curve_points.begin(),
572 return OpenCASCADE::point_compare(p1,
580 Handle(TColgp_HArray1OfPnt)
vertices =
581 new TColgp_HArray1OfPnt(1, n_vertices);
582 for (
unsigned int vertex = 0; vertex < n_vertices; ++vertex)
588 GeomAPI_Interpolate bspline_generator(
vertices, closed, tolerance);
589 bspline_generator.Perform();
590 Assert((bspline_generator.IsDone()),
591 ExcMessage(
"Interpolated bspline generation failed"));
593 Handle(Geom_BSplineCurve) bspline = bspline_generator.Curve();
594 TopoDS_Edge out_shape = BRepBuilderAPI_MakeEdge(bspline);
595 out_shape.Closed(closed);
601 template <
int spacedim>
602 std::vector<TopoDS_Edge>
610 std::map<unsigned int, std::pair<unsigned int, unsigned int>> vert_to_faces;
611 std::map<unsigned int, std::pair<unsigned int, unsigned int>> face_to_verts;
612 std::map<unsigned int, bool> visited_faces;
613 std::map<unsigned int, Point<spacedim>> vert_to_point;
615 unsigned int face_index;
617 for (
const auto &cell :
triangulation.active_cell_iterators())
619 if (cell->face(f)->at_boundary())
622 face_index = cell->face(f)->index();
623 const unsigned int v0 = cell->face(f)->vertex_index(0);
624 const unsigned int v1 = cell->face(f)->vertex_index(1);
625 face_to_verts[face_index].first =
v0;
626 face_to_verts[face_index].second =
v1;
627 visited_faces[face_index] =
false;
632 f, 0,
true,
false,
false)];
634 f, 1,
true,
false,
false)];
637 if (vert_to_faces.find(
v0) == vert_to_faces.end())
639 vert_to_faces[
v0].first = face_index;
643 vert_to_faces[
v0].second = face_index;
645 if (vert_to_faces.find(
v1) == vert_to_faces.end())
647 vert_to_faces[
v1].first = face_index;
651 vert_to_faces[
v1].second = face_index;
657 std::vector<TopoDS_Edge> interpolation_curves;
658 bool finished = (face_to_verts.empty());
659 face_index = finished ? 0 : face_to_verts.begin()->first;
661 while (finished ==
false)
663 const unsigned int start_point_index = face_to_verts[face_index].first;
664 unsigned int point_index = start_point_index;
667 std::vector<Point<spacedim>> pointlist;
670 visited_faces[face_index] =
true;
671 auto current_point = vert_to_point[point_index];
672 pointlist.push_back(current_point);
675 if (face_to_verts[face_index].
first != point_index)
676 point_index = face_to_verts[face_index].first;
678 point_index = face_to_verts[face_index].second;
681 if (vert_to_faces[point_index].
first != face_index)
682 face_index = vert_to_faces[point_index].first;
684 face_index = vert_to_faces[point_index].second;
686 while (point_index != start_point_index);
688 interpolation_curves.push_back(
692 for (
const auto &f : visited_faces)
693 if (f.second ==
false)
695 face_index = f.first;
700 return interpolation_curves;
705 std::tuple<Point<dim>, TopoDS_Shape,
double,
double>
708 const double tolerance)
711 gp_Pnt Pproj =
point(origin);
713 double minDistance = 1e7;
714 gp_Pnt tmp_proj(0.0, 0.0, 0.0);
716 [[maybe_unused]]
unsigned int counter = 0;
717 unsigned int face_counter = 0;
719 TopoDS_Shape out_shape;
723 for (exp.Init(in_shape, TopAbs_FACE); exp.More(); exp.Next())
725 TopoDS_Face face = TopoDS::Face(exp.Current());
729 Handle(Geom_Surface) SurfToProj = BRep_Tool::Surface(face);
731 ShapeAnalysis_Surface projector(SurfToProj);
732 gp_Pnt2d proj_params = projector.ValueOfUV(
point(origin), tolerance);
734 SurfToProj->D0(proj_params.X(), proj_params.Y(), tmp_proj);
736 double distance = point<dim>(tmp_proj).distance(origin);
737 if (distance < minDistance)
739 minDistance = distance;
753 if (face_counter == 0)
754 for (exp.Init(in_shape, TopAbs_EDGE); exp.More(); exp.Next())
756 TopoDS_Edge edge = TopoDS::Edge(exp.Current());
757 if (!BRep_Tool::Degenerated(edge))
765 Handle(Geom_Curve) CurveToProj =
766 BRep_Tool::Curve(edge,
L, First, Last);
768 GeomAPI_ProjectPointOnCurve Proj(
point(origin), CurveToProj);
769 unsigned int num_proj_points = Proj.NbPoints();
770 if ((num_proj_points > 0) && (Proj.LowerDistance() < minDistance))
772 minDistance = Proj.LowerDistance();
773 Pproj = Proj.NearestPoint();
775 u = Proj.LowerDistanceParameter();
782 return std::tuple<Point<dim>, TopoDS_Shape,
double,
double>(
783 point<dim>(Pproj), out_shape, u, v);
791 const double tolerance)
793 std::tuple<Point<dim>, TopoDS_Shape,
double,
double> ref =
795 return std::get<0>(ref);
801 const double tolerance)
804 std::tuple<Point<3>, TopoDS_Shape,
double,
double> shape_and_params =
807 TopoDS_Shape &out_shape = std::get<1>(shape_and_params);
808 double &u = std::get<2>(shape_and_params);
809 double &v = std::get<3>(shape_and_params);
813 std::tuple<unsigned int, unsigned int, unsigned int>
numbers =
820 "Could not find normal: the shape containing the closest point has 0 faces."));
824 "Could not find normal: the shape containing the closest point has more than 1 face."));
828 exp.Init(out_shape, TopAbs_FACE);
829 TopoDS_Face face = TopoDS::Face(exp.Current());
835 push_forward(
const TopoDS_Shape &in_shape,
const double u,
const double v)
837 switch (in_shape.ShapeType())
841 BRepAdaptor_Surface surf(TopoDS::Face(in_shape));
842 return point<dim>(surf.Value(u, v));
846 BRepAdaptor_Curve curve(TopoDS::Edge(in_shape));
847 return point<dim>(curve.Value(u));
861 Handle(Geom_Surface) SurfToProj = BRep_Tool::Surface(face);
862 GeomLProp_SLProps props(SurfToProj, u, v, 1, 1
e-7);
863 gp_Pnt Value = props.Value();
865 gp_Dir Normal = props.Normal();
866 Assert(props.IsCurvatureDefined(),
867 ExcMessage(
"Curvature is not well defined!"));
868 Standard_Real Min_Curvature = props.MinCurvature();
869 Standard_Real Max_Curvature = props.MaxCurvature();
876 if (face.Orientation() == TopAbs_REVERSED)
891 template <
int spacedim>
896 BRepAdaptor_Surface surf(face);
897 const double u0 = surf.FirstUParameter();
898 const double u1 = surf.LastUParameter();
899 const double v0 = surf.FirstVParameter();
900 const double v1 = surf.LastVParameter();
902 std::vector<CellData<2>> cells;
903 std::vector<Point<spacedim>>
vertices;
906 vertices.push_back(point<spacedim>(surf.Value(u0,
v0)));
907 vertices.push_back(point<spacedim>(surf.Value(u1,
v0)));
908 vertices.push_back(point<spacedim>(surf.Value(u0,
v1)));
909 vertices.push_back(point<spacedim>(surf.Value(u1,
v1)));
912 for (
unsigned int i = 0; i < 4; ++i)
915 cells.push_back(cell);
919 # include "utilities.inst"
Abstract base class for mapping classes.
virtual boost::container::small_vector< Point< spacedim >, GeometryInfo< dim >::vertices_per_cell > get_vertices(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
numbers::NumberTraits< Number >::real_type norm() const
virtual void create_triangulation(const std::vector< Point< spacedim >> &vertices, const std::vector< CellData< dim >> &cells, const SubCellData &subcelldata)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
static ::ExceptionBase & ExcNotImplemented()
static ::ExceptionBase & ExcUnsupportedShape()
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
graph_traits< Graph >::vertex_descriptor Vertex
std::tuple< Point< 3 >, Tensor< 1, 3 >, double, double > push_forward_and_differential_forms(const TopoDS_Face &face, const double u, const double v, const double tolerance=1e-7)
void extract_compound_shapes(const TopoDS_Shape &shape, std::vector< TopoDS_Compound > &compounds, std::vector< TopoDS_CompSolid > &compsolids, std::vector< TopoDS_Solid > &solids, std::vector< TopoDS_Shell > &shells, std::vector< TopoDS_Wire > &wires)
Point< dim > push_forward(const TopoDS_Shape &in_shape, const double u, const double v)
std::tuple< unsigned int, unsigned int, unsigned int > count_elements(const TopoDS_Shape &shape)
void write_STEP(const TopoDS_Shape &shape, const std::string &filename)
std::tuple< Point< 3 >, Tensor< 1, 3 >, double, double > closest_point_and_differential_forms(const TopoDS_Shape &in_shape, const Point< 3 > &origin, const double tolerance=1e-7)
void extract_geometrical_shapes(const TopoDS_Shape &shape, std::vector< TopoDS_Face > &faces, std::vector< TopoDS_Edge > &edges, std::vector< TopoDS_Vertex > &vertices)
TopoDS_Edge join_edges(const TopoDS_Shape &in_shape, const double tolerance=1e-7)
TopoDS_Shape read_STEP(const std::string &filename, const double scale_factor=1e-3)
bool point_compare(const Point< dim > &p1, const Point< dim > &p2, const Tensor< 1, dim > &direction=Tensor< 1, dim >(), const double tolerance=1e-10)
Point< spacedim > point(const gp_Pnt &p, const double tolerance=1e-10)
TopoDS_Edge interpolation_curve(std::vector< Point< dim >> &curve_points, const Tensor< 1, dim > &direction=Tensor< 1, dim >(), const bool closed=false, const double tolerance=1e-7)
Point< dim > closest_point(const TopoDS_Shape &in_shape, const Point< dim > &origin, const double tolerance=1e-7)
std::tuple< Point< dim >, TopoDS_Shape, double, double > project_point_and_pull_back(const TopoDS_Shape &in_shape, const Point< dim > &origin, const double tolerance=1e-7)
Point< dim > line_intersection(const TopoDS_Shape &in_shape, const Point< dim > &origin, const Tensor< 1, dim > &direction, const double tolerance=1e-7)
void write_STL(const TopoDS_Shape &shape, const std::string &filename, const double deflection, const bool sew_different_faces=false, const double sewer_tolerance=1e-6, const bool is_relative=false, const double angular_deflection=0.5, const bool in_parallel=false)
void create_triangulation(const TopoDS_Face &face, Triangulation< 2, spacedim > &tria)
TopoDS_Shape read_STL(const std::string &filename)
TopoDS_Shape intersect_plane(const TopoDS_Shape &in_shape, const double c_x, const double c_y, const double c_z, const double c, const double tolerance=1e-7)
std::vector< TopoDS_Edge > create_curves_from_triangulation_boundary(const Triangulation< 2, spacedim > &triangulation, const Mapping< 2, spacedim > &mapping=StaticMappingQ1< 2, spacedim >::mapping)
double get_shape_tolerance(const TopoDS_Shape &shape)
void write_IGES(const TopoDS_Shape &shape, const std::string &filename)
TopoDS_Shape read_IGES(const std::string &filename, const double scale_factor=1e-3)
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
SymmetricTensor< 2, dim, Number > e(const Tensor< 2, dim, Number > &F)
bool check(const ConstraintKinds kind_in, const unsigned int dim)
const ::parallel::distributed::Triangulation< dim, spacedim > * triangulation
std::vector< unsigned int > vertices
static unsigned int face_to_cell_vertices(const unsigned int face, const unsigned int vertex, const bool face_orientation=true, const bool face_flip=false, const bool face_rotation=false)
const ::Triangulation< dim, spacedim > & tria