libMesh::PeriodicBoundaries类 参考

We're using a class instead of a typedef to allow forward declarations and future flexibility. 更多...

#include <periodic_boundaries.h>

类 libMesh::PeriodicBoundaries 继承关系图:

Public 成员函数
PeriodicBoundaryBase *	boundary (boundary_id_type id)
const PeriodicBoundaryBase *	boundary (boundary_id_type id) const
	PeriodicBoundaries ()=default
	~PeriodicBoundaries ()
const Elem *	neighbor (boundary_id_type boundary_id, const PointLocatorBase &point_locator, const Elem *e, unsigned int side, unsigned int *neigh_side=nullptr) const

详细描述

We're using a class instead of a typedef to allow forward declarations and future flexibility.

作者

Roy Stogner

日期

2010 Maps between boundary ids and PeriodicBoundaryBase objects.

在文件 periodic_boundaries.h 第 49 行定义.

构造及析构函数说明

libMesh::PeriodicBoundaries::PeriodicBoundaries ( )

default

libMesh::PeriodicBoundaries::~PeriodicBoundaries ( )

default

成员函数说明

PeriodicBoundaryBase * libMesh::PeriodicBoundaries::boundary

(

boundary_id_type

id

)

在文件 periodic_boundaries.C 第 38 行定义.

参考自 neighbor().

{
  iterator i = this->find(id);
  if (i == this->end())
    return nullptr;
  return i->second.get();
}

const PeriodicBoundaryBase * libMesh::PeriodicBoundaries::boundary

(

boundary_id_type

id

)

const

在文件 periodic_boundaries.C 第 48 行定义.

{
  const_iterator i = this->find(id);
  if (i == this->end())
    return nullptr;
  return i->second.get();
}

const Elem * libMesh::PeriodicBoundaries::neighbor	(	boundary_id_type	boundary_id,
		const PointLocatorBase &	point_locator,
		const Elem *	e,
		unsigned int	side,
		unsigned int *	neigh_side = nullptr
	)		const

在文件 periodic_boundaries.C 第 59 行定义.

参考 boundary(), libMesh::PeriodicBoundaryBase::get_corresponding_pos(), libMesh::invalid_uint , 以及 libMesh::PeriodicBoundaryBase::pairedboundary.

{
  std::unique_ptr<const Elem> neigh_side_proxy;

  // Find a point on that side (and only that side)
  Point p = e->build_side_ptr(side)->vertex_average();

  const PeriodicBoundaryBase * b = this->boundary(boundary_id);
  libmesh_assert (b);
  p = b->get_corresponding_pos(p);

  std::set<const Elem *> candidate_elements;
  point_locator.operator()(p, candidate_elements);

  // We might have found multiple elements, e.g. if two distinct periodic
  // boundaries are overlapping (see systems_of_equations_ex9, for example).
  // As a result, we need to search for the element that has boundary_id.
  const MeshBase & mesh = point_locator.get_mesh();
  for(const Elem * elem_it : candidate_elements)
    {
      std::vector<unsigned int> neigh_sides =
        mesh.get_boundary_info().sides_with_boundary_id(elem_it, b->pairedboundary);

      for (auto ns : neigh_sides)
        {
          if (neigh_side)
            {
              elem_it->build_side_ptr(neigh_side_proxy, ns);
              if (neigh_side_proxy->contains_point(p))
                {
                  *neigh_side = ns;
                  return elem_it;
                }
            }
          else
            // checking contains_point is too expensive if we don't
            // definitely need it to find neigh_side
            return elem_it;
        }
    }

  // If we should have found a periodic neighbor but didn't then
  // either we're on a ghosted element with a remote periodic neighbor
  // or we're on a mesh with an inconsistent periodic boundary.
  libmesh_error_msg_if(mesh.is_serial() ||
                       (e->processor_id() == mesh.processor_id()),
                       "Periodic boundary neighbor not found");

  if (neigh_side)
    *neigh_side = libMesh::invalid_uint;
  return remote_elem;
}

---

该类的文档由以下文件生成:

• /home/lwz/libmesh/include/base/periodic_boundaries.h
• /home/lwz/libmesh/src/base/periodic_boundaries.C