face_csr.hpp File Reference

#include "peclet/core/common/types.hpp"
#include "peclet/core/common/portable.hpp"

Go to the source code of this file.

Classes
struct	peclet::core::solver::HostArr< T >
	A uniform accessor over a raw host array, giving it the operator()(i) that Kokkos::View has, so the row kernels can be written once for both. More...
struct	peclet::core::solver::FaceCsrOpT< D, I >
	A backend-agnostic view of an assembled face-CSR operator. More...
struct	peclet::core::solver::FvCsrOpT< D, I >
	A backend-agnostic view of an assembled FV (weight-CSR) operator. More...

Namespaces
namespace	peclet
namespace	peclet::core
namespace	peclet::core::solver

Functions
template<class Op , class U >
PECLET_HD double	peclet::core::solver::faceCsrApplyRow (const Op &op, Index i, const U &u)
	(A u)_i — one assembled-operator row.
template<class Op , class U >
PECLET_HD void	peclet::core::solver::faceCsrOffDiag (const Op &op, Index i, const U &u, double &off, double &d)
	Off-diagonal sum and the (advection-inclusive) diagonal for the point smoothers: out off = Σ coef·u[nbr] (+ advection inflow), d = diag (+ advection outflow).
PECLET_HD double	peclet::core::solver::faceCsrPointUpdate (double b_i, double off, double d, double uOld, double omega)
	The damped point update used by both Jacobi and (multicolour) Gauss–Seidel: returns the new u_i given the right-hand side b_i, the off-diagonal sum off, the diagonal d, the old value uOld, and the relaxation omega.
template<class Op , class U >
PECLET_HD double	peclet::core::solver::fvApplyRow (const Op &op, Index i, const U &u)
	(H u)_i = c0·u_i + cD·( invVol_i·( Σ w·(u_nbr − u_i) − bcDiag_i·u_i ) ).
template<class Op , class U >
PECLET_HD double	peclet::core::solver::fvPointSolve (const Op &op, Index i, const U &u, double rhs_i, double uOld)
	The point solve of H u = rhs for one row (raw value, before damping), matching the host point solve and device jacobiFv: the pure-L path (c0=0,cD=1) keeps the exact original expression.