#include <cstddef>
#include "peclet/core/common/types.hpp"
#include "peclet/core/common/view.hpp"
#include "peclet/core/solver/coloring.hpp"
#include "peclet/core/solver/face_csr.hpp"

Go to the source code of this file.

Classes
struct	peclet::core::solver::MomentumOp
	Assembled momentum operator on the device: (A u)_i = diag_i u_i + Σ coef·u[nbr], with an optional implicit-FOU advection part (rebuilt each step from the lagged velocity): a per-cell outflow diagonal `advDiag` + per-face inflow coefficients over a second (face-geometry) CSR. More...

Namespaces
namespace	peclet

namespace	peclet::core

namespace	peclet::core::solver

Functions
FaceCsrOpT< View< const double >, View< const Index > >	peclet::core::solver::momView (const MomentumOp &op)
	View the assembled momentum operator through the shared, backend-agnostic FaceCsrOpT, so the device kernels and the host serial solver (cut_cell.hpp) run the same row arithmetic (face_csr.hpp) and cannot drift.

void	peclet::core::solver::applyMom (const MomentumOp &op, View< const double > u, View< double > Au)
	Au = A u (cut-cell operator + optional implicit-FOU advection).

void	peclet::core::solver::residualMom (const MomentumOp &op, View< const double > u, View< const double > b, View< double > res)
	res = b − A u.

void	peclet::core::solver::jacobiMom (const MomentumOp &op, View< double > u, View< const double > b, View< double > tmp, double omega)
	One weighted-Jacobi sweep of A u = b (in place).

void	peclet::core::solver::multicolorGSMom (const MomentumOp &op, View< double > u, View< const double > b, const Coloring &col, double omega)
	One symmetric multicolour Gauss–Seidel sweep of A u = b in place (momentum operator: diag + face CSR + optional implicit-FOU advection): a forward pass over colours 0…C-1 followed by a reverse pass C-1…0.

Classes

Namespaces

Functions