The per-cut-cell-scalar state (ScalarField::cut). More...

#include <scalar_cutcell_operator.hpp>

Public Attributes
ScalarWallSpec	wallDefault
	every body without an instance override

std::vector< std::pair< int, ScalarWallSpec > >	wallInstance
	instance overrides, last wins

double	rtol = 1e-10
	max-norm relative residual tolerance (§5.1)

int	maxit = 200
	BiCGStab iteration cap (§5.1)

double	sourceConst = 0.0
	volumetric source S (c/T), uniform part

CCField	sourceField
	or per cell (physical), when allocated

bool	sourceIsField = false
	set_scalar_source was given an array (sourceField)

bool	hasProfile [6] = {false, false, false, false, false, false}

std::vector< double >	profile [6]
	Dirichlet domain-face profile (t1 fastest)

int	profileN [6][2] = {}
	its tangential extents

double	meanGrad [3] = {0.0, 0.0, 0.0}
	Mean-gradient (closure) mode (§1.5, WO-6): the field holds theta, c = G.x + theta, in the frame moving with the fluid's mean velocity.

CCField	SAC
	level-0 surrogate diagonal (bands + lumped wall)

Kokkos::View< double *, CCMem >	cw

Kokkos::View< double *, CCMem >	rw
	per facet (§4.1)

Kokkos::View< double *, CCMem >	gFace [6]
	Dirichlet domain-face value per face cell.

bool	dirFace [6] = {false, false, false, false, false, false}

double	lam = 0.0
	internal Lam = D tRef/hRef^2 of that operator

double	idt = 0.0
	1/dt' (0 steady)

double	dt = 0.0
	dt' (internal)

double	srcFactor = 0.0
	S' = S * srcFactor (tRef)

std::vector< int >	wtType
	the resolved per-body wall table of that operator

std::vector< double >	wtK

std::vector< double >	wtG

std::vector< double >	wtQ
	(internal k', g, q')

long	numUnknowns = 0
	global fluid unknowns (the singular mean)

bool	steady = false

bool	built = false

bool	singular = false

bool	meanGradOn = false
	that operator carried the mean-gradient terms (§1.5)

double	gPhys [3] = {0.0, 0.0, 0.0}
	its G (physical, c/L)

double	gInt [3] = {0.0, 0.0, 0.0}
	its G' = G hRef (c per internal length)

double	ubar [3] = {0.0, 0.0, 0.0}
	its moving-frame velocity U'bar = sum U_i / sum kappa_i V

CCField	kr

CCField	krh

CCField	kp

CCField	kv

CCField	kt

CCField	kz

CCField	kz2

Kokkos::View< double *, CCMem >	ks
	WO-9b: the resident Krylov's scalar slots.

Kokkos::View< double *, Kokkos::HostSpace >	pk
	and its host packet

CCField	kb

CCField	kq
	singular case: the projected rhs, and the preconditioner's rhs copy

int	iterations = 0

double	residual = 0.0
	final TRUE residual max\|b - A c\| / ref, after any gauge shift

bool	converged = true

bool	warnedNoConv = false

long	warnedGeometryVersion = -1
	The geometry version whose §9 resolution warnings were issued (WO-8: once per scalar and geometry; a rebuilt geometry warns again, review finding 7).

bool	warnedKoren = false
	the koren-above-bulk-Courant-1/2 warning was issued (once)

Kokkos::View< int *, CCMem >	wtTypeV

Kokkos::View< int *, CCMem >	mtConjV

Kokkos::View< double *, CCMem >	wtKV

Kokkos::View< double *, CCMem >	wtGV

Kokkos::View< double *, CCMem >	wtQV

Kokkos::View< double *, CCMem >	mtLamV

Kokkos::View< double *, CCMem >	mtCKV

Kokkos::View< double *, CCMem >	mtKV

Kokkos::View< double *, CCMem >	mtRcV

double	incompatibility = 0.0
	steady singular case: \|sum b\| / sum \|b\| before projection

std::shared_ptr< ScalarMG >	mg

long	mgVersion = -1
	the geometry version its level table was built for

std::size_t	mgN = 0
	the block size it was built for

int	mgLevels = 1
	levels the last solve's V-cycle used (census `mg_levels`)

CCField	phi [3]
	F/V through the LOW a-face of cell i (internal 1/T), guarded (sco::faceFluxes)

CCField	Phi [3]
	explicit face fluxes phi c*(c^n), 0 on implicit faces (§6.2)

CCField	small
	1 on a small cell (§6.3), inner cells + exchanged ghosts

CCField	outflow
	lumped implicit outflow per cell (inside SAC; ScalarMG coarse mass, §5.2)

CCField	omegaOpen
	A2 (§6.7, D-WO5-3): the implicit outflow the open-face rows add to AC, summed over the cell's open faces and clamped at 0 (implicit backflow), 0 elsewhere — the advective path's coarse mass.

double	maxCellPeclet = 0.0
	census `max_cell_peclet` (A2): max \|phi_a\| / (Lam' w_a)

bool	openFace [6] = {false, false, false, false, false, false}
	open face rows built (WO-5b)

bool	openInflow [6] = {false, false, false, false, false, false}
	... and it is an inflow

bool	advecting = false
	some face carried flux: else every advection kernel was skipped

long	numSmall = 0
	census `num_small_cells`

long	numImplicitFaces = 0
	census `num_implicit_faces` (faces carrying flux, implicit)

long	numFluxFaces = 0
	faces carrying flux (the implicit fraction's denominator)

long	numGuardedFaces = 0
	faces whose projection flux the guard zeroed (a no-op reading)

long	numBackflowFaces = 0
	census `num_backflow_faces`: outflow-face rows with F_out < 0

bool	warnedBackflow = false
	the steady backflow warning was issued (once per scalar)

bool	warnedBottom = false
	the ScalarMG shallow-bottom warning was issued (once)

double	bulkCourant = 0.0
	census `bulk_courant`: C_bulk of §6.3 (0 steady)

CCField	solid
	the registered field `<name>_solid` (psi_s), when allocated

bool	conj = false
	the last operator had a conjugate body (two fields)

long	solidVersion = -1
	geometry version of the record

std::vector< char >	solidKey
	the conjugate flag per body it was built for

CCField	mat
	material (body id) per cell, haloed; -1: no solid

CCField	sunk
	solid-unknown flags (haloed, 0 beyond non-periodic faces)

Kokkos::View< double *, CCMem >	sS
	solid probe distance per facet (conjugate facets)

Kokkos::View< int *[8], Kokkos::LayoutLeft, CCMem >	offS

Kokkos::View< double *[8], Kokkos::LayoutLeft, CCMem >	wS

Kokkos::View< std::uint8_t *, CCMem >	rungS

Kokkos::View< std::uint8_t *, CCMem >	conjF
	1 on a facet of a conjugate body

long	numSolidUnknowns = 0
	census `num_solid_unknowns` (global)

double	solidVolume = 0.0
	census `solid_volume`: conjugate solid, physical

double	sealedSolidVolume = 0.0
	census `sealed_solid_volume` (D-WO7-1), physical

long	solidRungs [4] = {0, 0, 0, 0}
	census probe_rungs['solid'] (global)

CCField	spx

CCField	spy

CCField	spz
	level-0 solid face products Lam_face a_s (guarded), low faces

CCField	ms
	solid mass C_s K kappa_s idt on solid unknowns

CCField	ASC

CCField	ASW

CCField	ASE

CCField	ASS

CCField	ASN

CCField	ASB

CCField	AST

CCField	bS

CCField	sOld

CCField	SACs

CCField	W0
	solid surrogate diagonal; the same-cell coupling sum cc

Kokkos::View< double *, CCMem >	cc
	per facet: alpha G_c (conjugate facets, else 0)

std::vector< double >	mtLam

std::vector< double >	mtCK

std::vector< double >	mtK

std::vector< double >	mtRc
	internal per-body material (conjugate bodies)

double	maxSolidD = 0.0
	max internal D_s' = Lam_s'/(C_s K) (the level rule)

CCField	krS

CCField	krhS

CCField	kpS

CCField	kvS

CCField	ktS

CCField	kzS

CCField	kz2S

CCField	kbS

CCField	kqS
	Krylov scratch, solid phase.

Detailed Description

The per-cut-cell-scalar state (ScalarField::cut).

Settings are stored verbatim (physical); the operator members are rebuilt on every advance or steady solve (ruling D-WO3-3).

Definition at line 56 of file scalar_cutcell_operator.hpp.

Member Data Documentation

◆ wallDefault

ScalarWallSpec peclet::flow::ScalarCutState::wallDefault

every body without an instance override

Definition at line 58 of file scalar_cutcell_operator.hpp.

◆ wallInstance

std::vector<std::pair<int, ScalarWallSpec> > peclet::flow::ScalarCutState::wallInstance

instance overrides, last wins

Definition at line 59 of file scalar_cutcell_operator.hpp.

◆ rtol

double peclet::flow::ScalarCutState::rtol = 1e-10

max-norm relative residual tolerance (§5.1)

Definition at line 60 of file scalar_cutcell_operator.hpp.

◆ maxit

int peclet::flow::ScalarCutState::maxit = 200

BiCGStab iteration cap (§5.1)

Definition at line 61 of file scalar_cutcell_operator.hpp.

◆ sourceConst

double peclet::flow::ScalarCutState::sourceConst = 0.0

volumetric source S (c/T), uniform part

Definition at line 62 of file scalar_cutcell_operator.hpp.

◆ sourceField

CCField peclet::flow::ScalarCutState::sourceField

or per cell (physical), when allocated

Definition at line 63 of file scalar_cutcell_operator.hpp.

◆ sourceIsField

bool peclet::flow::ScalarCutState::sourceIsField = false

set_scalar_source was given an array (sourceField)

Definition at line 64 of file scalar_cutcell_operator.hpp.

◆ hasProfile

bool peclet::flow::ScalarCutState::hasProfile[6] = {false, false, false, false, false, false}

Definition at line 65 of file scalar_cutcell_operator.hpp.

◆ profile

std::vector<double> peclet::flow::ScalarCutState::profile[6]

Dirichlet domain-face profile (t1 fastest)

Definition at line 66 of file scalar_cutcell_operator.hpp.

◆ profileN

int peclet::flow::ScalarCutState::profileN[6][2] = {}

its tangential extents

Definition at line 67 of file scalar_cutcell_operator.hpp.

◆ meanGrad

double peclet::flow::ScalarCutState::meanGrad[3] = {0.0, 0.0, 0.0}

Mean-gradient (closure) mode (§1.5, WO-6): the field holds theta, c = G.x + theta, in the frame moving with the fluid's mean velocity.

G physical (c/L); (0, 0, 0) is off.

Definition at line 70 of file scalar_cutcell_operator.hpp.

◆ SAC

CCField peclet::flow::ScalarCutState::SAC

level-0 surrogate diagonal (bands + lumped wall)

Definition at line 72 of file scalar_cutcell_operator.hpp.

◆ cw

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::cw

Definition at line 73 of file scalar_cutcell_operator.hpp.

◆ rw

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::rw

per facet (§4.1)

Definition at line 73 of file scalar_cutcell_operator.hpp.

◆ gFace

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::gFace[6]

Dirichlet domain-face value per face cell.

Definition at line 74 of file scalar_cutcell_operator.hpp.

◆ dirFace

bool peclet::flow::ScalarCutState::dirFace[6] = {false, false, false, false, false, false}

Definition at line 75 of file scalar_cutcell_operator.hpp.

◆ lam

double peclet::flow::ScalarCutState::lam = 0.0

internal Lam = D tRef/hRef^2 of that operator

Definition at line 76 of file scalar_cutcell_operator.hpp.

Referenced by anderson_synthetic_maps.LinearMapSolver::step().

◆ idt

double peclet::flow::ScalarCutState::idt = 0.0

1/dt' (0 steady)

Definition at line 77 of file scalar_cutcell_operator.hpp.

◆ dt

double peclet::flow::ScalarCutState::dt = 0.0

dt' (internal)

Definition at line 78 of file scalar_cutcell_operator.hpp.

Referenced by steady_acceleration_gates.Case::configure().

◆ srcFactor

double peclet::flow::ScalarCutState::srcFactor = 0.0

S' = S * srcFactor (tRef)

Definition at line 79 of file scalar_cutcell_operator.hpp.

◆ wtType

std::vector<int> peclet::flow::ScalarCutState::wtType

the resolved per-body wall table of that operator

Definition at line 80 of file scalar_cutcell_operator.hpp.

◆ wtK

std::vector<double> peclet::flow::ScalarCutState::wtK

Definition at line 81 of file scalar_cutcell_operator.hpp.

◆ wtG

std::vector<double> peclet::flow::ScalarCutState::wtG

Definition at line 81 of file scalar_cutcell_operator.hpp.

◆ wtQ

std::vector<double> peclet::flow::ScalarCutState::wtQ

(internal k', g, q')

Definition at line 81 of file scalar_cutcell_operator.hpp.

◆ numUnknowns

long peclet::flow::ScalarCutState::numUnknowns = 0

global fluid unknowns (the singular mean)

Definition at line 82 of file scalar_cutcell_operator.hpp.

◆ steady

bool peclet::flow::ScalarCutState::steady = false

Definition at line 83 of file scalar_cutcell_operator.hpp.

◆ built

bool peclet::flow::ScalarCutState::built = false

Definition at line 83 of file scalar_cutcell_operator.hpp.

◆ singular

bool peclet::flow::ScalarCutState::singular = false

Definition at line 83 of file scalar_cutcell_operator.hpp.

◆ meanGradOn

bool peclet::flow::ScalarCutState::meanGradOn = false

that operator carried the mean-gradient terms (§1.5)

Definition at line 84 of file scalar_cutcell_operator.hpp.

◆ gPhys

double peclet::flow::ScalarCutState::gPhys[3] = {0.0, 0.0, 0.0}

its G (physical, c/L)

Definition at line 85 of file scalar_cutcell_operator.hpp.

◆ gInt

double peclet::flow::ScalarCutState::gInt[3] = {0.0, 0.0, 0.0}

its G' = G hRef (c per internal length)

Definition at line 86 of file scalar_cutcell_operator.hpp.

◆ ubar

double peclet::flow::ScalarCutState::ubar[3] = {0.0, 0.0, 0.0}

its moving-frame velocity U'bar = sum U_i / sum kappa_i V

Definition at line 87 of file scalar_cutcell_operator.hpp.

◆ kr

CCField peclet::flow::ScalarCutState::kr

Definition at line 89 of file scalar_cutcell_operator.hpp.

◆ krh

CCField peclet::flow::ScalarCutState::krh

Definition at line 89 of file scalar_cutcell_operator.hpp.

◆ kp

CCField peclet::flow::ScalarCutState::kp

Definition at line 89 of file scalar_cutcell_operator.hpp.

◆ kv

CCField peclet::flow::ScalarCutState::kv

Definition at line 89 of file scalar_cutcell_operator.hpp.

◆ kt

CCField peclet::flow::ScalarCutState::kt

Definition at line 89 of file scalar_cutcell_operator.hpp.

◆ kz

CCField peclet::flow::ScalarCutState::kz

Definition at line 89 of file scalar_cutcell_operator.hpp.

◆ kz2

CCField peclet::flow::ScalarCutState::kz2

Definition at line 89 of file scalar_cutcell_operator.hpp.

◆ ks

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::ks

WO-9b: the resident Krylov's scalar slots.

Definition at line 90 of file scalar_cutcell_operator.hpp.

◆ pk

Kokkos::View<double*, Kokkos::HostSpace> peclet::flow::ScalarCutState::pk

and its host packet

Definition at line 91 of file scalar_cutcell_operator.hpp.

◆ kb

CCField peclet::flow::ScalarCutState::kb

Definition at line 92 of file scalar_cutcell_operator.hpp.

◆ kq

CCField peclet::flow::ScalarCutState::kq

singular case: the projected rhs, and the preconditioner's rhs copy

Definition at line 92 of file scalar_cutcell_operator.hpp.

◆ iterations

int peclet::flow::ScalarCutState::iterations = 0

Definition at line 93 of file scalar_cutcell_operator.hpp.

◆ residual

double peclet::flow::ScalarCutState::residual = 0.0

final TRUE residual max|b - A c| / ref, after any gauge shift

Definition at line 94 of file scalar_cutcell_operator.hpp.

Referenced by anderson_oracle.AndersonOracle::step().

◆ converged

bool peclet::flow::ScalarCutState::converged = true

Definition at line 95 of file scalar_cutcell_operator.hpp.

◆ warnedNoConv

bool peclet::flow::ScalarCutState::warnedNoConv = false

Definition at line 96 of file scalar_cutcell_operator.hpp.

◆ warnedGeometryVersion

long peclet::flow::ScalarCutState::warnedGeometryVersion = -1

The geometry version whose §9 resolution warnings were issued (WO-8: once per scalar and geometry; a rebuilt geometry warns again, review finding 7).

-1: none yet.

Definition at line 99 of file scalar_cutcell_operator.hpp.

◆ warnedKoren

bool peclet::flow::ScalarCutState::warnedKoren = false

the koren-above-bulk-Courant-1/2 warning was issued (once)

Definition at line 100 of file scalar_cutcell_operator.hpp.

◆ wtTypeV

Kokkos::View<int*, CCMem> peclet::flow::ScalarCutState::wtTypeV

Definition at line 103 of file scalar_cutcell_operator.hpp.

◆ mtConjV

Kokkos::View<int*, CCMem> peclet::flow::ScalarCutState::mtConjV

Definition at line 103 of file scalar_cutcell_operator.hpp.

◆ wtKV

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::wtKV

Definition at line 104 of file scalar_cutcell_operator.hpp.

◆ wtGV

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::wtGV

Definition at line 104 of file scalar_cutcell_operator.hpp.

◆ wtQV

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::wtQV

Definition at line 104 of file scalar_cutcell_operator.hpp.

◆ mtLamV

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::mtLamV

Definition at line 104 of file scalar_cutcell_operator.hpp.

◆ mtCKV

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::mtCKV

Definition at line 104 of file scalar_cutcell_operator.hpp.

◆ mtKV

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::mtKV

Definition at line 104 of file scalar_cutcell_operator.hpp.

◆ mtRcV

Kokkos::View<double*, CCMem> peclet::flow::ScalarCutState::mtRcV

Definition at line 104 of file scalar_cutcell_operator.hpp.

◆ incompatibility

double peclet::flow::ScalarCutState::incompatibility = 0.0

steady singular case: |sum b| / sum |b| before projection

Definition at line 105 of file scalar_cutcell_operator.hpp.

◆ mg

std::shared_ptr<ScalarMG> peclet::flow::ScalarCutState::mg

Definition at line 107 of file scalar_cutcell_operator.hpp.

◆ mgVersion

long peclet::flow::ScalarCutState::mgVersion = -1

the geometry version its level table was built for

Definition at line 108 of file scalar_cutcell_operator.hpp.

◆ mgN

std::size_t peclet::flow::ScalarCutState::mgN = 0

the block size it was built for

Definition at line 109 of file scalar_cutcell_operator.hpp.

◆ mgLevels

int peclet::flow::ScalarCutState::mgLevels = 1

levels the last solve's V-cycle used (census mg_levels)

Definition at line 110 of file scalar_cutcell_operator.hpp.

◆ phi

CCField peclet::flow::ScalarCutState::phi[3]

F/V through the LOW a-face of cell i (internal 1/T), guarded (sco::faceFluxes)

Definition at line 113 of file scalar_cutcell_operator.hpp.

◆ Phi

CCField peclet::flow::ScalarCutState::Phi[3]

explicit face fluxes phi c*(c^n), 0 on implicit faces (§6.2)

Definition at line 114 of file scalar_cutcell_operator.hpp.

◆ small

CCField peclet::flow::ScalarCutState::small

1 on a small cell (§6.3), inner cells + exchanged ghosts

Definition at line 115 of file scalar_cutcell_operator.hpp.

◆ outflow

CCField peclet::flow::ScalarCutState::outflow

lumped implicit outflow per cell (inside SAC; ScalarMG coarse mass, §5.2)

Definition at line 116 of file scalar_cutcell_operator.hpp.

◆ omegaOpen

CCField peclet::flow::ScalarCutState::omegaOpen

A2 (§6.7, D-WO5-3): the implicit outflow the open-face rows add to AC, summed over the cell's open faces and clamped at 0 (implicit backflow), 0 elsewhere — the advective path's coarse mass.

Written by steady advecting solves only (outflow stays the transient path's input).

Definition at line 120 of file scalar_cutcell_operator.hpp.