#include "base/cs_defs.h"
#include "base/cs_math.h"
#include "base/cs_parameters.h"
#include "base/cs_base.h"
#include "base/cs_halo.h"

Include dependency graph for cs_matrix_building.h:

Functions
void	cs_matrix_compute_coeffs (cs_matrix_t a, const cs_field_t f, int iconvp, int idiffp, int ndircp, double thetap, double relaxp, int imucpp, const cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t rovsdt[], const cs_real_t i_massflux[], const cs_real_t b_massflux[], const cs_real_t i_visc[], const cs_real_t b_visc[], const cs_real_t xcpp[])
	Build the diagonal of the advection/diffusion matrix for determining the variable time step, flow, Fourier. More...

template<cs_lnum_t stride>
void	cs_matrix_compute_coeffs (cs_matrix_t a, const cs_field_t f, int iconvp, int idiffp, int tensorial_diffusion, int ndircp, cs_lnum_t eb_size, double thetap, double relaxp, const cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t fimp[][stride][stride], const cs_real_t i_massflux[], const cs_real_t b_massflux[], const cs_real_t i_visc[], const cs_real_t b_visc[])
	Build the diagonal of the advection/diffusion matrix for determining the variable time step, flow, Fourier. More...

void	cs_matrix_wrapper (int iconvp, int idiffp, int ndircp, int isym, double thetap, int imucpp, const cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t rovsdt[], const cs_real_t i_massflux[], const cs_real_t b_massflux[], const cs_real_t i_visc[], const cs_real_t b_visc[], const cs_real_t xcpp[], cs_real_t da[], cs_real_t xa[])

template<cs_lnum_t stride>
void	cs_matrix_wrapper (int iconvp, int idiffp, int tensorial_diffusion, int ndircp, int isym, cs_lnum_t eb_size, double thetap, const cs_field_bc_coeffs_t *bc_coeffs_v, const cs_real_t fimp[][stride][stride], const cs_real_t i_massflux[], const cs_real_t b_massflux[], const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_t da[][stride][stride], cs_real_t xa[])
	Build the matrix for a vector or tensor field. More...

void	cs_matrix_time_step (const cs_mesh_t m, int iconvp, int idiffp, int isym, const cs_field_bc_coeffs_t bc_coeffs, const cs_real_t i_massflux[], const cs_real_t b_massflux[], const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_t *da)

Function Documentation

◆ cs_matrix_compute_coeffs() [1/2]

void cs_matrix_compute_coeffs	(	cs_matrix_t *	a,
		const cs_field_t *	f,
		int	iconvp,
		int	idiffp,
		int	ndircp,
		double	thetap,
		double	relaxp,
		int	imucpp,
		const cs_field_bc_coeffs_t *	bc_coeffs,
		const cs_real_t	rovsdt[],
		const cs_real_t	i_massflux[],
		const cs_real_t	b_massflux[],
		const cs_real_t	i_visc[],
		const cs_real_t	b_visc[],
		const cs_real_t	xcpp[]
	)

Build the diagonal of the advection/diffusion matrix for determining the variable time step, flow, Fourier.

Parameters

[in,out]	a	pointer to matrix structure
[in]	f	pointer to field, or null
[in]	iconvp	indicator 1 advection 0 otherwise
[in]	idiffp	indicator 1 diffusion 0 otherwise
[in]	ndircp	number of Dirichlet BCs
[in]	thetap	time scheme parameter
[in]	relaxp	relaxation coefficient (if < 1)
[in]	imucp	1 for temperature (with Cp), 0 otherwise
[in]	bc_coeffs	boundary condition structure
[in]	rovsdt	implicit terms (rho / dt)
[in]	i_massflux	mass flux at interior faces
[in]	b_massflux	mass flux at border faces
[in]	i_visc	$\mu_\fij \dfrac{S_\fij}{\ipf \jpf}$ at interior faces for the matrix
[in]	b_visc	$S_\fib$ at border faces for the matrix
[in]	xcpp	Cp per cell, or null

◆ cs_matrix_compute_coeffs() [2/2]

void cs_matrix_compute_coeffs	(	cs_matrix_t *	a,
		const cs_field_t *	f,
		int	iconvp,
		int	idiffp,
		int	tensorial_diffusion,
		int	ndircp,
		cs_lnum_t	eb_size,
		double	thetap,
		double	relaxp,
		const cs_field_bc_coeffs_t *	bc_coeffs,
		const cs_real_t	fimp[][stride][stride],
		const cs_real_t	i_massflux[],
		const cs_real_t	b_massflux[],
		const cs_real_t	i_visc[],
		const cs_real_t	b_visc[]
	)

Build the diagonal of the advection/diffusion matrix for determining the variable time step, flow, Fourier.

Template Parameters

stride 3 for vectors, 6 for tensors

Parameters

[in,out]	a	pointer to matrix structure
[in]	f	pointer to field, or null
[in]	iconvp	indicator 1 advection 0 otherwise
[in]	idiffp	indicator 1 diffusion 0 otherwise
[in]	tensorial_diffusion	indicator
[in]	ndircp	number of Dirichlet BCs
[in]	thetap	time scheme parameter
[in]	relaxp	relaxation coefficient (if < 1)
[in]	eb_size	extra-diagonal block size (1 or 3 for stride 3, 1 for stride 6)
[in]	bc_coeffs	boundary conditions structure
[in]	fimp	implicit terms, or null
[in]	i_massflux	mass flux at interior faces
[in]	b_massflux	mass flux at border faces
[in]	i_visc	$\mu_\fij \dfrac{S_\fij}{\ipf \jpf}$ at interior faces for the matrix
[in]	b_visc	$S_\fib$ at boundary faces for the matrix

◆ cs_matrix_time_step()

void cs_matrix_time_step	(	const cs_mesh_t *	m,
		int	iconvp,
		int	idiffp,
		int	isym,
		const cs_field_bc_coeffs_t *	bc_coeffs,
		const cs_real_t	i_massflux[],
		const cs_real_t	b_massflux[],
		const cs_real_t	i_visc[],
		const cs_real_t	b_visc[],
		cs_real_t *	da
	)

◆ cs_matrix_wrapper() [1/2]

void cs_matrix_wrapper	(	int	iconvp,
		int	idiffp,
		int	ndircp,
		int	isym,
		double	thetap,
		int	imucpp,
		const cs_field_bc_coeffs_t *	bc_coeffs,
		const cs_real_t	rovsdt[],
		const cs_real_t	i_massflux[],
		const cs_real_t	b_massflux[],
		const cs_real_t	i_visc[],
		const cs_real_t	b_visc[],
		const cs_real_t	xcpp[],
		cs_real_t	da[],
		cs_real_t	xa[]
	)

◆ cs_matrix_wrapper() [2/2]

void cs_matrix_wrapper	(	int	iconvp,
		int	idiffp,
		int	tensorial_diffusion,
		int	ndircp,
		int	isym,
		cs_lnum_t	eb_size,
		double	thetap,
		const cs_field_bc_coeffs_t *	bc_coeffs_v,
		const cs_real_t	fimp[][stride][stride],
		const cs_real_t	i_massflux[],
		const cs_real_t	b_massflux[],
		const cs_real_t	i_visc[],
		const cs_real_t	b_visc[],
		cs_real_t	da[][stride][stride],
		cs_real_t	xa[]
	)

Build the matrix for a vector or tensor field.

The advection (if present) is upwind. The diffusion is not reconstructed. The matrix is split into a diagonal part (stride*stride blocks) and an extra diagonal part.

Template Parameters

stride 3 for vectors, 6 for tensors

Parameters

[in]	m	pointer to mesh structure
[in]	idiffp	indicator 1 diffusion 0 otherwise
[in]	thetap	weighting coefficient for the theta-scheme, thetap = 0: explicit scheme thetap = 0.5: time-centered scheme (mix between Crank-Nicolson and Adams-Bashforth) thetap = 1: implicit scheme
[in]	bc_coeffs_v	boundary condition structure for the variable
[in]	fimp	$\tens{f_s}^{imp}$ , or null
[in]	i_visc	$\mu_\fij \dfrac{S_\fij}{\ipf \jpf}$ at interior faces for the matrix
[in]	b_visc	$\mu_\fib \dfrac{S_\fib}{\ipf \centf}$ at border faces for the matrix
[out]	da	diagonal part of the matrix
[out]	xa	extra interleaved diagonal part of the matrix

Functions

Function Documentation

◆ cs_matrix_compute_coeffs() [1/2]

◆ cs_matrix_compute_coeffs() [2/2]

◆ cs_matrix_time_step()

◆ cs_matrix_wrapper() [1/2]

◆ cs_matrix_wrapper() [2/2]