ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Public Member Functions | Private Attributes | Static Private Attributes | List of all members
polygon_ Class Reference

#include <ERF_EBPolygon.H>

Collaboration diagram for polygon_:
Collaboration graph
[legend]

Public Member Functions

AMREX_GPU_HOST_DEVICE polygon_ (amrex::RealVect a_point, amrex::RealVect a_normal)
 
AMREX_GPU_HOST_DEVICE polygon_ ()
 
AMREX_GPU_HOST_DEVICE void add_vertex (amrex::RealVect const &a_v)
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void set_area (amrex::Real const &a_area)
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void define ()
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE int ok () const noexcept
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real area () const noexcept
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real distance (amrex::RealVect const &a_point) const noexcept
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::RealVect get_centroid () const noexcept
 
void report (int const a_id, amrex::RealVect a_v0)
 

Private Attributes

int const m_cell_face
 
amrex::RealVect const m_eb_point
 
amrex::RealVect const m_eb_normal
 
int m_defined
 
int m_vertices
 
amrex::Real m_area
 
int m_sorted
 
amrex::Array< amrex::RealVect, m_max_vertices > m_vertex
 
amrex::GpuArray< amrex::Real, m_max_vertices > m_theta
 
amrex::RealVect m_zdir
 

Static Private Attributes

static constexpr int m_max_vertices = 6
 

Constructor & Destructor Documentation

◆ polygon_() [1/2]

AMREX_GPU_HOST_DEVICE polygon_::polygon_ ( amrex::RealVect  a_point,
amrex::RealVect  a_normal 
)
inline
16  : m_cell_face(1)
17  , m_eb_point(a_point)
18  , m_eb_normal(a_normal)
19  , m_defined(0)
20  , m_vertices(0)
21  , m_area(0.)
22  , m_sorted(0)
23  , m_vertex({amrex::RealVect(0.), amrex::RealVect(0.), amrex::RealVect(0.),
24  amrex::RealVect(0.), amrex::RealVect(0.), amrex::RealVect(0.)})
25  , m_zdir(0.)
26  {}
polygon_::m_defined
int m_defined
Definition: ERF_EBPolygon.H:180
polygon_::m_area
amrex::Real m_area
Definition: ERF_EBPolygon.H:184
polygon_::m_cell_face
int const m_cell_face
Definition: ERF_EBPolygon.H:175
polygon_::m_sorted
int m_sorted
Definition: ERF_EBPolygon.H:186
polygon_::m_vertices
int m_vertices
Definition: ERF_EBPolygon.H:182
polygon_::m_vertex
amrex::Array< amrex::RealVect, m_max_vertices > m_vertex
Definition: ERF_EBPolygon.H:188
polygon_::m_eb_normal
amrex::RealVect const m_eb_normal
Definition: ERF_EBPolygon.H:178
polygon_::m_eb_point
amrex::RealVect const m_eb_point
Definition: ERF_EBPolygon.H:177
polygon_::m_zdir
amrex::RealVect m_zdir
Definition: ERF_EBPolygon.H:192

◆ polygon_() [2/2]

AMREX_GPU_HOST_DEVICE polygon_::polygon_ ( )
inline
30  : m_cell_face(0)
31  , m_eb_point(0.)
32  , m_eb_normal(0.)
33  , m_defined(0)
34  , m_vertices(0)
35  , m_area(0.)
36  , m_sorted(0)
37  , m_vertex({amrex::RealVect(0.), amrex::RealVect(0.), amrex::RealVect(0.),
38  amrex::RealVect(0.), amrex::RealVect(0.), amrex::RealVect(0.)})
39  , m_zdir(0.)
40  {}

Member Function Documentation

◆ add_vertex()

AMREX_GPU_HOST_DEVICE void polygon_::add_vertex ( amrex::RealVect const &  a_v)
inline
43  {
44  AMREX_ASSERT( m_vertices < m_max_vertices );
45  m_vertex[m_vertices] = a_v;
46  m_vertices++;
47  }
polygon_::m_max_vertices
static constexpr int m_max_vertices
Definition: ERF_EBPolygon.H:173

Referenced by eb_cut_cell_::calc_edge_intersections().

Here is the caller graph for this function:

◆ area()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real polygon_::area ( ) const
inlinenoexcept
123  {
124  AMREX_ALWAYS_ASSERT( ok() ); // <-------------------------------------- TODO remove ALWAYS
125  return m_area;
126  }
polygon_::ok
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE int ok() const noexcept
Definition: ERF_EBPolygon.H:119

Referenced by eb_cut_cell_::volume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ define()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void polygon_::define ( )
inline
53  {
54 
55  AMREX_ALWAYS_ASSERT( m_defined == 0 ); // TODO ---------------------------- remove ALWAYS
56 
57  m_defined = 1;
58 
59  // We need at least 3 faces.
60  if (m_vertices < 3) { return; }
61 
62  // Check to see if the vertices of the face are inside or outside of the plane.
63  // If they are outside, this face doesn't belong to the volume.
64  if ( m_cell_face ) {
65  for ( int i(0); i<m_vertices; ++i) {
66  if ((m_eb_normal.dotProduct(m_vertex[i] - m_eb_point)) >= 0.) {
67  } else {
68  }
69  }
70  }
71 
72  // Calculate the centroid of the polygon.
73  amrex::RealVect centroid = get_centroid();
74 
75  // Shift the vertices relative to the centroid
76  amrex::Array<amrex::RealVect,m_max_vertices> vertex_cent;
77  for ( int i(0); i<m_vertices; ++i) { vertex_cent[i] = m_vertex[i] - centroid; }
78 
79  // Find a vertex that isn't inline with v[0] and the centroid so we can
80  // compute a vector normal to the face.
81  for ( int i(1); i<m_vertices; ++i ) {
82  m_zdir = vertex_cent[0].crossProduct(vertex_cent[i]);
83  if ( !amrex::almostEqual(m_zdir.vectorLength(),0.) ) { break; }
84  }
85  AMREX_ALWAYS_ASSERT(!amrex::almostEqual(m_zdir.vectorLength(), 0.));
86  m_zdir /= m_zdir.vectorLength();
87 
88  m_theta[0] = 0.0;
89  for ( int i(1); i<m_vertices; ++i) {
90 
91  m_theta[i] = std::atan2(m_zdir.dotProduct( vertex_cent[0].crossProduct(vertex_cent[i]) ),
92  vertex_cent[0].dotProduct(vertex_cent[i]));
93 
94  m_theta[i] += ((m_theta[i] >= 0.) ? 0.0 : 2.0*PI);
95  }
96 
97  // Sort counter clockwise based on theta.
98  for (int i(0); i<m_vertices; ++i) {
99  for (int j(0); j < m_vertices-i-1; ++j) {
100  if ( m_theta[j] > m_theta[j+1] ) {
101  amrex::Swap(m_theta[j], m_theta[j+1]);
102  amrex::Swap(m_vertex[j], m_vertex[j+1]);
103  amrex::Swap(vertex_cent[j], vertex_cent[j+1]);
104  }
105  } // j-loop
106  } // i-loop
107  m_sorted = 1;
108 
109  // Compute areas of triangles
110  for (int i(0); i<m_vertices; ++i) {
111  int const j( (i+1 == m_vertices) ? 0 : i+1 );
112  amrex::RealVect vi_cross_vj = vertex_cent[i].crossProduct(vertex_cent[j]);
113  m_area += 0.5*vi_cross_vj.vectorLength();
114  }
115  AMREX_ALWAYS_ASSERT( m_area > 0. ); // <------------------------------- TODO remove ALWAYS
116  }
PI
constexpr amrex::Real PI
Definition: ERF_Constants.H:6
polygon_::get_centroid
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::RealVect get_centroid() const noexcept
Definition: ERF_EBPolygon.H:138
polygon_::m_theta
amrex::GpuArray< amrex::Real, m_max_vertices > m_theta
Definition: ERF_EBPolygon.H:190

Referenced by eb_cut_cell_::eb_cut_cell_().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ distance()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real polygon_::distance ( amrex::RealVect const &  a_point) const
inlinenoexcept
130  {
131  AMREX_ALWAYS_ASSERT( m_defined == 1 ); // <--------------------------- TODO remove ALWAYS
132  amrex::RealVect x0 = a_point - m_vertex[0];
133  return amrex::Math::abs(x0.dotProduct(m_zdir));
134  }

Referenced by report(), and eb_cut_cell_::volume().

Here is the caller graph for this function:

◆ get_centroid()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::RealVect polygon_::get_centroid ( ) const
inlinenoexcept
138  {
139  amrex::RealVect cent(0.);
140  for ( int i(0); i<m_vertices; ++i) { cent += m_vertex[i]; }
141  return cent / static_cast<amrex::Real>(m_vertices);
142  }

Referenced by define(), and report().

Here is the caller graph for this function:

◆ ok()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE int polygon_::ok ( ) const
inlinenoexcept
120  { return ((m_area > 0. || (m_area == 0. && m_defined == 1)) ? 1 : 0); }

Referenced by area(), report(), and eb_cut_cell_::volume().

Here is the caller graph for this function:

◆ report()

void polygon_::report ( int const  a_id,
amrex::RealVect  a_v0 
)
inline
149  {
150 
151 #ifndef AMREX_USE_GPU
152  amrex::RealVect centroid = get_centroid();
153  amrex::Print() << "Face " << a_id
154  << " -------------------------------------------\n"
155  << "\nok? " << (ok() ? "yes" : "no") << "\n\n";
156  for (int i(0); i<m_vertices; ++i) {
157  amrex::Print() << "v" << i << ": " << m_vertex[i] << '\n';
158  }
159  amrex::Print() << "\nvc: " << " " << centroid << "\n";
160 
161  amrex::Real const dist = distance( a_v0 );
162  // amrex::Real const vol = ok() ? m_area * dist : 0.;
163 
164  amrex::Print() << "\narea: " << m_area
165  << "\ndistance: " << dist
166  << "\nvolume: " << dist*m_area
167  << "\n==================================================\n\n";
168 #endif
169  }
polygon_::distance
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real distance(amrex::RealVect const &a_point) const noexcept
Definition: ERF_EBPolygon.H:130
Here is the call graph for this function:

◆ set_area()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void polygon_::set_area ( amrex::Real const &  a_area)
inline
50 { m_area = a_area; }

Referenced by eb_cut_cell_::set_regular().

Here is the caller graph for this function:

Member Data Documentation

◆ m_area

amrex::Real polygon_::m_area
private

Referenced by area(), define(), ok(), report(), and set_area().

◆ m_cell_face

int const polygon_::m_cell_face
private

Referenced by define().

◆ m_defined

int polygon_::m_defined
private

Referenced by define(), distance(), and ok().

◆ m_eb_normal

amrex::RealVect const polygon_::m_eb_normal
private

Referenced by define().

◆ m_eb_point

amrex::RealVect const polygon_::m_eb_point
private

Referenced by define().

◆ m_max_vertices

constexpr int polygon_::m_max_vertices = 6
staticconstexprprivate

Referenced by add_vertex().

◆ m_sorted

int polygon_::m_sorted
private

Referenced by define().

◆ m_theta

amrex::GpuArray<amrex::Real,m_max_vertices> polygon_::m_theta
private

Referenced by define().

◆ m_vertex

amrex::Array<amrex::RealVect,m_max_vertices> polygon_::m_vertex
private

Referenced by add_vertex(), define(), distance(), get_centroid(), and report().

◆ m_vertices

int polygon_::m_vertices
private

Referenced by add_vertex(), define(), get_centroid(), and report().

◆ m_zdir

amrex::RealVect polygon_::m_zdir
private

Referenced by define(), and distance().

The documentation for this class was generated from the following file: