|
| 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.
|
| |