core
Shared MPI block decomposition + asynchronous ghost-layer exchange (header-only C++20)
Loading...
Searching...
No Matches
face_csr.hpp File Reference
#include "peclet/core/common/types.hpp"
#include <morton/morton.hpp>

Go to the source code of this file.

Classes

struct  peclet::core::amr::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::amr::FaceCsrOpT< D, I >
 A backend-agnostic view of an assembled face-CSR operator. More...
 
struct  peclet::core::amr::FvCsrOpT< D, I >
 A backend-agnostic view of an assembled FV (weight-CSR) operator. More...
 

Namespaces

namespace  peclet
 
namespace  peclet::core
 
namespace  peclet::core::amr
 

Macros

#define MORTON_HD
 

Functions

template<class Op , class U >
MORTON_HD double peclet::core::amr::faceCsrApplyRow (const Op &op, Index i, const U &u)
 (A u)_i — one assembled-operator row.
 
template<class Op , class U >
MORTON_HD void peclet::core::amr::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).
 
MORTON_HD double peclet::core::amr::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 >
MORTON_HD double peclet::core::amr::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 >
MORTON_HD double peclet::core::amr::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.
 

Macro Definition Documentation

◆ MORTON_HD

#define MORTON_HD

Definition at line 32 of file face_csr.hpp.