37 CGAL::Polyhedron_incremental_builder_3<HDS>
B(hds,
true);
39 B.begin_surface(s_.nPoints(), s_.size());
41 typedef typename HDS::Vertex
Vertex;
42 typedef typename Vertex::Point
Point;
44 for (
const auto& pt : s_.points())
46 B.add_vertex(
Point(pt.x(), pt.y(), pt.z()));
49 for (
const auto&
f : s_)
53 for (
const label verti :
f)
55 B.add_vertex_to_facet(verti);
69 Build_triangle<HalfedgeDS> triangle(
s);
73 Foam::label nTris = 0;
83 fi->region =
s[nTris].region();
static const Foam::dimensionedScalar B("", Foam::dimless, 18.678)
Polyhedron::Facet_iterator Facet_iterator
CGAL::Polyhedron_3< K, My_items > Polyhedron
graph_traits< Graph >::vertex_descriptor Vertex
Ostream & operator()() const
Output stream (master only).
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))