mesh_gen::anonymous_namespace{annulusMesh3D.cpp} Namespace Reference

Functions
void	assertCuboidAnnulusValid (const geom::Cuboid &outer, const geom::Cuboid &inner)
void	assertSphereAnnulusValid (const geom::Sphere &outer, const geom::Sphere &inner)
static int	firstTagOfDim (const std::vector< std::pair< int, int > > &ov, int dim)
util::Point	embedPointInAnnulus3D (const geom::AnnulusGeomObject &a)
void	finish3DCutEmbedVolume (const std::vector< std::pair< int, int > > &ov, const geom::AnnulusGeomObject &a, double h)
void	buildCuboidAnnulus3D (const geom::AnnulusGeomObject &a, double h)
void	buildSphereAnnulus3D (const geom::AnnulusGeomObject &a, double h)

Function Documentation

assertCuboidAnnulusValid()

void mesh_gen::anonymous_namespace{annulusMesh3D.cpp}::assertCuboidAnnulusValid	(	const geom::Cuboid &	outer,
		const geom::Cuboid &	inner
	)

Definition at line 22 of file annulusMesh3D.cpp.

                                                                                {
  const auto &blo = outer.d_vertices[0];
  const auto &bhi = outer.d_vertices[6];
  const auto &bli = inner.d_vertices[0];
  const auto &bhi_i = inner.d_vertices[6];
 
  constexpr double eps = 1.0e-12;
  if (!(bli.d_x > blo.d_x + eps && bhi_i.d_x < bhi.d_x - eps && bli.d_y > blo.d_y + eps &&
        bhi_i.d_y < bhi.d_y - eps && bli.d_z > blo.d_z + eps && bhi_i.d_z < bhi.d_z - eps))
    throw std::runtime_error(
        "buildAnnulus3DInCurrentModel: inner cuboid must be strictly inside the outer cuboid.");
}

References assertCuboidAnnulusValid(), and geom::Cuboid::d_vertices.

Referenced by assertCuboidAnnulusValid(), and buildCuboidAnnulus3D().

Here is the call graph for this function:

Here is the caller graph for this function:

assertSphereAnnulusValid()

void mesh_gen::anonymous_namespace{annulusMesh3D.cpp}::assertSphereAnnulusValid	(	const geom::Sphere &	outer,
		const geom::Sphere &	inner
	)

Definition at line 35 of file annulusMesh3D.cpp.

                                                                                {
  constexpr double eps = 1.0e-10;
  if ((outer.d_x - inner.d_x).lengthSq() > eps)
    throw std::runtime_error(
        "buildAnnulus3DInCurrentModel: sphere annulus requires coincident centers (concentric spheres).");
  if (inner.d_r <= 0. || outer.d_r <= inner.d_r)
    throw std::runtime_error(
        "buildAnnulus3DInCurrentModel: require 0 < r_inner < r_outer for sphere − sphere.");
}

References assertSphereAnnulusValid(), geom::Sphere::d_r, and geom::Sphere::d_x.

Referenced by assertSphereAnnulusValid(), and buildSphereAnnulus3D().

Here is the call graph for this function:

Here is the caller graph for this function:

buildCuboidAnnulus3D()

void mesh_gen::anonymous_namespace{annulusMesh3D.cpp}::buildCuboidAnnulus3D	(	const geom::AnnulusGeomObject &	a,
		double	h
	)

Definition at line 94 of file annulusMesh3D.cpp.

                                                                    {
  const auto &outer = *static_cast<const geom::Cuboid *>(a.d_outObj_p);
  const auto &inner = *static_cast<const geom::Cuboid *>(a.d_inObj_p);
  assertCuboidAnnulusValid(outer, inner);
 
  const auto &lo = outer.d_vertices[0];
  const int out_vol =
      gmsh::model::occ::addBox(lo.d_x, lo.d_y, lo.d_z, outer.d_Lx, outer.d_Ly, outer.d_Lz);
  gmsh::model::occ::synchronize();
 
  const auto &li = inner.d_vertices[0];
  const int in_vol =
      gmsh::model::occ::addBox(li.d_x, li.d_y, li.d_z, inner.d_Lx, inner.d_Ly, inner.d_Lz);
  gmsh::model::occ::synchronize();
 
  std::vector<std::pair<int, int>> ov;
  std::vector<std::vector<std::pair<int, int>>> ovv;
  gmsh::model::occ::cut({{3, out_vol}}, {{3, in_vol}}, ov, ovv, -1, true, true);
  gmsh::model::occ::synchronize();
  gmsh::model::occ::removeAllDuplicates();
  gmsh::model::occ::synchronize();
 
  finish3DCutEmbedVolume(ov, a, h);
}

References assertCuboidAnnulusValid(), buildCuboidAnnulus3D(), geom::AnnulusGeomObject::d_inObj_p, geom::AnnulusGeomObject::d_outObj_p, and finish3DCutEmbedVolume().

Referenced by buildCuboidAnnulus3D().

Here is the call graph for this function:

Here is the caller graph for this function:

buildSphereAnnulus3D()

void mesh_gen::anonymous_namespace{annulusMesh3D.cpp}::buildSphereAnnulus3D	(	const geom::AnnulusGeomObject &	a,
		double	h
	)

Definition at line 119 of file annulusMesh3D.cpp.

                                                                    {
  const auto &outer = *static_cast<const geom::Sphere *>(a.d_outObj_p);
  const auto &inner = *static_cast<const geom::Sphere *>(a.d_inObj_p);
  assertSphereAnnulusValid(outer, inner);
 
  const int out_vol =
      gmsh::model::occ::addSphere(outer.d_x.d_x, outer.d_x.d_y, outer.d_x.d_z, outer.d_r);
  gmsh::model::occ::synchronize();
  const int in_vol =
      gmsh::model::occ::addSphere(inner.d_x.d_x, inner.d_x.d_y, inner.d_x.d_z, inner.d_r);
  gmsh::model::occ::synchronize();
 
  std::vector<std::pair<int, int>> ov;
  std::vector<std::vector<std::pair<int, int>>> ovv;
  gmsh::model::occ::cut({{3, out_vol}}, {{3, in_vol}}, ov, ovv, -1, true, true);
  gmsh::model::occ::synchronize();
  gmsh::model::occ::removeAllDuplicates();
  gmsh::model::occ::synchronize();
 
  finish3DCutEmbedVolume(ov, a, h);
}

References assertSphereAnnulusValid(), buildSphereAnnulus3D(), geom::AnnulusGeomObject::d_inObj_p, geom::AnnulusGeomObject::d_outObj_p, and finish3DCutEmbedVolume().

Referenced by buildSphereAnnulus3D().

Here is the call graph for this function:

Here is the caller graph for this function:

embedPointInAnnulus3D()

util::Point mesh_gen::anonymous_namespace{annulusMesh3D.cpp}::embedPointInAnnulus3D

(

const geom::AnnulusGeomObject &

a

)

See embedPointInAnnulus2D: centroid of shell can lie in the void; embed in solid material.

Definition at line 53 of file annulusMesh3D.cpp.

                                                                {
  if (a.d_inObj_p->d_name == "sphere" && a.d_outObj_p->d_name == "sphere") {
    const auto &outer = *static_cast<const geom::Sphere *>(a.d_outObj_p);
    const auto &inner = *static_cast<const geom::Sphere *>(a.d_inObj_p);
    const double r_mid = 0.5 * (outer.d_r + inner.d_r);
    return {outer.d_x.d_x + r_mid, outer.d_x.d_y, outer.d_x.d_z};
  }
  if (a.d_inObj_p->d_name == "cuboid" && a.d_outObj_p->d_name == "cuboid") {
    const auto &outer = *static_cast<const geom::Cuboid *>(a.d_outObj_p);
    const auto &inner = *static_cast<const geom::Cuboid *>(a.d_inObj_p);
    const auto &blo = outer.d_vertices[0];
    const auto &bli = inner.d_vertices[0];
    return {0.5 * (blo.d_x + bli.d_x), 0.5 * (blo.d_y + bli.d_y), 0.5 * (blo.d_z + bli.d_z)};
  }
  return a.center();
}

References geom::AnnulusGeomObject::center(), geom::AnnulusGeomObject::d_inObj_p, geom::GeomObject::d_name, geom::AnnulusGeomObject::d_outObj_p, geom::Cuboid::d_vertices, geom::Sphere::d_x, util::Point::d_x, and embedPointInAnnulus3D().

Referenced by embedPointInAnnulus3D(), and finish3DCutEmbedVolume().

Here is the call graph for this function:

Here is the caller graph for this function:

finish3DCutEmbedVolume()

void mesh_gen::anonymous_namespace{annulusMesh3D.cpp}::finish3DCutEmbedVolume	(	const std::vector< std::pair< int, int > > &	ov,
		const geom::AnnulusGeomObject &	a,
		double	h
	)

Definition at line 70 of file annulusMesh3D.cpp.

                                      {
 
  int vol_tag = firstTagOfDim(ov, 3);
  if (vol_tag < 0) {
    std::vector<std::pair<int, int>> ents;
    gmsh::model::getEntities(ents, 3);
    if (!ents.empty())
      vol_tag = ents.back().second;
  }
  if (vol_tag < 0) {
    std::ostringstream oss;
    oss << "buildAnnulus3DInCurrentModel: boolean cut did not produce a volume (ov size = "
        << ov.size() << ").";
    throw std::runtime_error(oss.str());
  }
 
  const util::Point c = embedPointInAnnulus3D(a);
  const int p = gmsh::model::occ::addPoint(c.d_x, c.d_y, c.d_z, h);
  gmsh::model::occ::synchronize();
  gmsh::model::mesh::embed(0, {p}, 3, vol_tag);
  gmsh::model::occ::synchronize();
}

References util::Point::d_x, util::Point::d_y, util::Point::d_z, embedPointInAnnulus3D(), finish3DCutEmbedVolume(), and firstTagOfDim().

Referenced by buildCuboidAnnulus3D(), buildSphereAnnulus3D(), and finish3DCutEmbedVolume().

Here is the call graph for this function:

Here is the caller graph for this function:

firstTagOfDim()

static int mesh_gen::anonymous_namespace{annulusMesh3D.cpp}::firstTagOfDim ( const std::vector< std::pair< int, int > > &  ov, int  dim  )

static

Definition at line 45 of file annulusMesh3D.cpp.

                                                                          {
  for (const auto &pr : ov)
    if (pr.first == dim)
      return pr.second;
  return -1;
}

References firstTagOfDim().

Referenced by finish3DCutEmbedVolume(), and firstTagOfDim().

Here is the call graph for this function:

Here is the caller graph for this function: