core 1.5.0
Shared MPI block decomposition, asynchronous ghost-layer exchange, SDF geometry, load balancing (header-only C++20)
Loading...
Searching...
No Matches
peclet::core::geom::quad_detail Namespace Reference

Gauss-Legendre nodes/weights on [0,1] for n = 1..8. More...

Functions

template<class Real >
PECLET_HD void gaussLegendre01 (int n, int i, Real &x, Real &w)
 
PECLET_HD int gaussCount (int n)
 
template<class Real , class Phi >
PECLET_HD Real positiveLength (const Phi &phi, Vec3< Real > p0, Vec3< Real > d, int nseg)
 Length of {phi > 0} along the segment p0 + s*d, s in [0,1], by bracketed bisection.
 

Detailed Description

Gauss-Legendre nodes/weights on [0,1] for n = 1..8.

Returned by reference into a static constexpr table so the caller pays nothing; n outside the range clamps.

Function Documentation

◆ gaussLegendre01()

template<class Real >
PECLET_HD void peclet::core::geom::quad_detail::gaussLegendre01 ( int  n,
int  i,
Real &  x,
Real &  w 
)

Definition at line 66 of file quadrature.hpp.

◆ gaussCount()

PECLET_HD int peclet::core::geom::quad_detail::gaussCount ( int  n)

◆ positiveLength()

template<class Real , class Phi >
PECLET_HD Real peclet::core::geom::quad_detail::positiveLength ( const Phi &  phi,
Vec3< Real >  p0,
Vec3< Real >  d,
int  nseg 
)

Length of {phi > 0} along the segment p0 + s*d, s in [0,1], by bracketed bisection.

nseg initial subintervals: a feature thinner than 1/nseg of the segment is INVISIBLE here – failure mode 1 in the header comment.

Definition at line 137 of file quadrature.hpp.

Referenced by peclet::core::geom::cellVolumeFraction(), and peclet::core::geom::faceAperture().