ERF_EBCutCell.H File Reference

#include <AMReX_REAL.H>
#include <AMReX_Array.H>
#include <AMReX_RealBox.H>
#include <AMReX_RealVect.H>
#include <AMReX_Algorithm.H>
#include <AMReX_EBCellFlag.H>
#include "ERF_EBPolygon.H"

Include dependency graph for ERF_EBCutCell.H:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	eb_cut_cell_
struct	eb_cut_cell_::path_data

Functions
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE int	intersect_plane_edge (amrex::RealVect const &a_plane_point, amrex::RealVect const &a_plane_normal, amrex::RealVect const &a_edge_point0, amrex::RealVect const &a_edge_point1, amrex::RealVect &a_intersection_point, amrex::Real &a_intersection_dist)

Function Documentation

intersect_plane_edge()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE int intersect_plane_edge	(	amrex::RealVect const &	a_plane_point,
		amrex::RealVect const &	a_plane_normal,
		amrex::RealVect const &	a_edge_point0,
		amrex::RealVect const &	a_edge_point1,
		amrex::RealVect &	a_intersection_point,
		amrex::Real &	a_intersection_dist
	)

{
  amrex::RealVect const edge(a_edge_point1 - a_edge_point0);
  amrex::Real const edge_length = edge.vectorLength();
 
  AMREX_ALWAYS_ASSERT(edge_length > 0.);
 
  amrex::RealVect edge_normal = edge / edge_length;
 
  amrex::Real np_dot_ne = a_plane_normal.dotProduct(edge_normal);
 
  // if ( amrex::Math::abs(np_dot_ne) < std::numeric_limits<amrex::Real>::min() )
  if ( amrex::Math::abs(np_dot_ne) < 10.0 * std::numeric_limits<amrex::Real>::epsilon() )
  { return 0; }
 
  a_intersection_dist = a_plane_normal.dotProduct(a_plane_point)
                      - a_plane_normal.dotProduct(a_edge_point0);
  a_intersection_dist /= np_dot_ne;
 
  a_intersection_point = a_edge_point0 + a_intersection_dist*edge_normal;
 
  if (0. <= a_intersection_dist && a_intersection_dist <= edge_length) { return 1;}
 
  return 0;
}

Referenced by eb_cut_cell_::calc_edge_intersections(), and eb_cut_cell_::path_data::set().

Here is the caller graph for this function: