9.2
general documentation
Enumerations | Functions
cs_convection_diffusion.h File Reference
#include "base/cs_defs.h"
#include "base/cs_base.h"
#include "base/cs_dispatch.h"
#include "base/cs_halo.h"
#include "base/cs_math.h"
#include "mesh/cs_mesh_quantities.h"
#include "cdo/cs_equation_param.h"
+ Include dependency graph for cs_convection_diffusion.h:

Go to the source code of this file.

Enumerations

enum  cs_nvd_type_t {
  CS_NVD_GAMMA = 0 , CS_NVD_SMART = 1 , CS_NVD_CUBISTA = 2 , CS_NVD_SUPERBEE = 3 ,
  CS_NVD_MUSCL = 4 , CS_NVD_MINMOD = 5 , CS_NVD_CLAM = 6 , CS_NVD_STOIC = 7 ,
  CS_NVD_OSHER = 8 , CS_NVD_WASEB = 9 , CS_NVD_VOF_HRIC = 10 , CS_NVD_VOF_CICSAM = 11 ,
  CS_NVD_VOF_STACS = 12 , CS_NVD_N_TYPES = 13
}
 

Functions

cs_real_t * cs_get_v_slope_test (int f_id, const cs_equation_param_t eqp)
 
void cs_beta_limiter_building (int f_id, int inc, const cs_real_t rovsdt[])
 Compute the beta blending coefficient of the beta limiter (ensuring preservation of a given min/max pair of values). More...
 
void cs_convection_diffusion_scalar (const cs_field_t *f, const cs_equation_param_t eqp, int icvflb, int inc, int imasac, const cs_real_t *pvar, const int icvfli[], 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[], const cs_real_t *c_weight, cs_real_t *rhs, cs_real_2_t i_flux[], cs_real_t b_flux[])
 
void cs_face_convection_scalar (int idtvar, int f_id, const cs_equation_param_t eqp, int icvflb, int inc, int imasac, cs_real_t *pvar, const cs_real_t *pvara, const int icvfli[], cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_massflux[], const cs_real_t b_massflux[], cs_real_2_t i_conv_flux[], cs_real_t b_conv_flux[])
 
void cs_convection_diffusion_vector (int idtvar, int f_id, const cs_equation_param_t eqp, int icvflb, int inc, int ivisep, int imasac, cs_real_3_t *pvar, const cs_real_3_t *pvara, const int icvfli[], 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[], const cs_real_t i_secvis[], const cs_real_t b_secvis[], cs_real_3_t *i_pvar, cs_real_3_t *b_pvar, cs_real_3_t *rhs)
 
void cs_convection_diffusion_tensor (int idtvar, int f_id, const cs_equation_param_t eqp, int icvflb, int inc, int imasac, cs_real_6_t *pvar, const cs_real_6_t *pvara, 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_6_t *rhs)
 
void cs_convection_diffusion_thermal (const cs_field_t *f, const cs_equation_param_t eqp, int inc, int imasac, const cs_real_t *pvar, 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[], const cs_real_t *c_weight, const cs_real_t xcpp[], cs_real_t *rhs)
 
void cs_anisotropic_diffusion_scalar (int idtvar, int f_id, const cs_equation_param_t eqp, int inc, const cs_real_t *pvar, const cs_real_t *pvara, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_6_t *viscel, const cs_real_2_t weighf[], const cs_real_t weighb[], cs_real_t *rhs)
 
void cs_anisotropic_left_diffusion_vector (int idtvar, int f_id, const cs_equation_param_t eqp, int inc, int ivisep, cs_real_3_t *pvar, const cs_real_3_t *pvara, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_33_t i_visc[], const cs_real_t b_visc[], const cs_real_t i_secvis[], cs_real_3_t *rhs)
 
void cs_anisotropic_right_diffusion_vector (int idtvar, int f_id, const cs_equation_param_t eqp, int inc, cs_real_3_t *pvar, const cs_real_3_t *pvara, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_6_t *viscel, const cs_real_2_t weighf[], const cs_real_t weighb[], cs_real_3_t *rhs)
 
void cs_anisotropic_diffusion_tensor (int idtvar, int f_id, const cs_equation_param_t eqp, int inc, cs_real_6_t *pvar, const cs_real_6_t *pvara, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_6_t *viscel, const cs_real_2_t weighf[], const cs_real_t weighb[], cs_real_6_t *rhs)
 
void cs_face_diffusion_potential (const cs_field_t *f, const cs_equation_param_t *eqp, const cs_mesh_t *m, cs_mesh_quantities_t *fvq, int init, int inc, int iphydp, cs_real_3_t *frcxt, cs_real_t *pvar, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_t *visel, cs_real_t *i_massflux, cs_real_t *b_massflux)
 
void cs_face_anisotropic_diffusion_potential (const cs_field_t *f, const cs_equation_param_t *eqp, const cs_mesh_t *m, cs_mesh_quantities_t *fvq, int init, int inc, int iphydp, cs_real_3_t *frcxt, cs_real_t *pvar, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_6_t *viscel, const cs_real_2_t weighf[], const cs_real_t weighb[], cs_real_t *i_massflux, cs_real_t *b_massflux)
 
void cs_diffusion_potential (const cs_field_t *f, const cs_equation_param_t *eqp, const cs_mesh_t *m, cs_mesh_quantities_t *fvq, int init, int inc, int iphydp, cs_real_3_t *frcxt, cs_real_t *pvar, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_t visel[], cs_real_t *diverg)
 
void cs_anisotropic_diffusion_potential (const cs_field_t *f, const cs_equation_param_t *eqp, const cs_mesh_t *m, cs_mesh_quantities_t *fvq, int init, int inc, int iphydp, cs_real_3_t *restrict frcxt, cs_real_t *restrict pvar, cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_6_t *restrict viscel, const cs_real_2_t weighf[], const cs_real_t weighb[], cs_real_t *restrict diverg)
 Add the explicit part of the divergence of the mass flux due to the pressure gradient (analog to cs_anisotropic_diffusion_scalar). More...
 
template<typename T >
void cs_slope_test_gradient (int f_id, cs_dispatch_context &ctx, const T(*grad)[3], T(*grdpa)[3], const cs_real_t *pvar, const cs_real_t val_f[], const cs_real_t *i_massflux)
 Compute the upwind gradient used in the slope tests. More...
 
template<cs_lnum_t stride, typename T >
void cs_slope_test_gradient_strided (cs_dispatch_context &ctx, const T grad[][stride][3], T(*restrict grdpa)[stride][3], const cs_real_t pvar[][stride], const cs_real_t val_f[][stride], const cs_real_t *i_massflux)
 Compute the upwind gradient used in the slope tests. More...
 
template<typename T >
void cs_upwind_gradient (cs_dispatch_context &ctx, const int inc, const cs_halo_type_t halo_type, const cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_massflux[], const cs_real_t b_massflux[], const cs_real_t *pvar, T(*grdpa)[3])
 
template<cs_lnum_t stride, typename T >
void cs_upwind_gradient_strided (cs_dispatch_context &ctx, const int inc, const cs_halo_type_t halo_type, const cs_field_bc_coeffs_t *bc_coeffs, const cs_real_t i_massflux[], const cs_real_t b_massflux[], const cs_real_t *restrict pvar[stride], T(*restrict grdpa)[stride][3])
 
void cs_cell_courant_number (const cs_field_t *f, cs_dispatch_context &ctx, cs_real_t *courant)
 
void cs_convection_diffusion_secvis (cs_dispatch_context &ctx, const cs_mesh_t *m, const cs_mesh_quantities_t *mq, cs_real_t thetap, const cs_real_t i_visc[], const cs_real_t i_secvis[], const cs_real_t b_secvis[], const cs_rreal_t gradv[][3][3], cs_real_3_t *restrict rhs)
 Compute balance contribution of the transpose grad(vel) term and grad(-2/3 div(vel)) More...
 
void cs_convection_anisotropic_leff_diffusion_secvis (cs_dispatch_context &ctx, const cs_mesh_t *m, const cs_mesh_quantities_t *mq, const cs_real_33_t i_visc[], const cs_real_t i_secvis[], const cs_real_t gradv[][3][3], cs_real_3_t *restrict rhs)
 Compute balance contribution of the transpose grad(vel) term and grad(-2/3 div(vel)) with anisotropic. More...
 
int cs_convection_diffusion_get_scheme_version (void)
 Query convection-diffusion scheme variants. More...
 
void cs_convection_diffusion_set_scheme_version (int version)
 Allow reverting to older convection-diffusion scheme variants. More...
 

Enumeration Type Documentation

◆ cs_nvd_type_t

enum cs_nvd_type_t
Enumerator
CS_NVD_GAMMA 
CS_NVD_SMART 
CS_NVD_CUBISTA 
CS_NVD_SUPERBEE 
CS_NVD_MUSCL 
CS_NVD_MINMOD 
CS_NVD_CLAM 
CS_NVD_STOIC 
CS_NVD_OSHER 
CS_NVD_WASEB 
CS_NVD_VOF_HRIC 
CS_NVD_VOF_CICSAM 
CS_NVD_VOF_STACS 
CS_NVD_N_TYPES 

Function Documentation

◆ cs_anisotropic_diffusion_potential()

void cs_anisotropic_diffusion_potential ( const cs_field_t *  f,
const cs_equation_param_t *  eqp,
const cs_mesh_t *  m,
cs_mesh_quantities_t *  fvq,
int  init,
int  inc,
int  iphydp,
cs_real_3_t *restrict  frcxt,
cs_real_t *restrict  pvar,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_visc[],
const cs_real_t  b_visc[],
cs_real_6_t *restrict  viscel,
const cs_real_2_t  weighf[],
const cs_real_t  weighb[],
cs_real_t *restrict  diverg 
)

Add the explicit part of the divergence of the mass flux due to the pressure gradient (analog to cs_anisotropic_diffusion_scalar).

More precisely, the divergence of the mass flux side $ \sum_{\fij \in \Facei{\celli}} \dot{m}_\fij $ is updated as follows:

\[ \sum_{\fij \in \Facei{\celli}} \dot{m}_\fij = \sum_{\fij \in \Facei{\celli}} \dot{m}_\fij - \sum_{\fij \in \Facei{\celli}} \left( \tens{\mu}_\fij \gradv_\fij P \cdot \vect{S}_\ij \right) \]

Parameters
[in]fpointer to field or nullptr
[in]eqpequation parameters
[in]mpointer to mesh
[in]fvqpointer to finite volume quantities
[in]initindicator
  • 1 initialize the mass flux to 0
  • 0 otherwise
[in]incindicator
  • 0 when solving an increment
  • 1 otherwise
[in]iphydpindicator
  • 1 hydrostatic pressure taken into account
  • 0 otherwise
[in]frcxtbody force creating the hydrostatic pressure
[in]pvarsolved variable (pressure)
[in]bc_coeffsboundary condition structure for the variable
[in]i_visc$ \mu_\fij \dfrac{S_\fij}{\ipf \jpf} $ at interior faces for the r.h.s.
[in]b_visc$ \mu_\fib \dfrac{S_\fib}{\ipf \centf} $ at border faces for the r.h.s.
[in]viscelsymmetric cell tensor $ \tens{\mu}_\celli $
[in]weighfinternal face weight between cells i j in case of tensor diffusion
[in]weighbboundary face weight for cells i in case of tensor diffusion
[in,out]divergdivergence of the mass flux

◆ cs_anisotropic_diffusion_scalar()

void cs_anisotropic_diffusion_scalar ( int  idtvar,
int  f_id,
const cs_equation_param_t  eqp,
int  inc,
const cs_real_t *  pvar,
const cs_real_t *  pvara,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_visc[],
const cs_real_t  b_visc[],
cs_real_6_t *  viscel,
const cs_real_2_t  weighf[],
const cs_real_t  weighb[],
cs_real_t *  rhs 
)

◆ cs_anisotropic_diffusion_tensor()

void cs_anisotropic_diffusion_tensor ( int  idtvar,
int  f_id,
const cs_equation_param_t  eqp,
int  inc,
cs_real_6_t *  pvar,
const cs_real_6_t *  pvara,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_visc[],
const cs_real_t  b_visc[],
cs_real_6_t *  viscel,
const cs_real_2_t  weighf[],
const cs_real_t  weighb[],
cs_real_6_t *  rhs 
)

◆ cs_anisotropic_left_diffusion_vector()

void cs_anisotropic_left_diffusion_vector ( int  idtvar,
int  f_id,
const cs_equation_param_t  eqp,
int  inc,
int  ivisep,
cs_real_3_t *  pvar,
const cs_real_3_t *  pvara,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_33_t  i_visc[],
const cs_real_t  b_visc[],
const cs_real_t  i_secvis[],
cs_real_3_t *  rhs 
)

◆ cs_anisotropic_right_diffusion_vector()

void cs_anisotropic_right_diffusion_vector ( int  idtvar,
int  f_id,
const cs_equation_param_t  eqp,
int  inc,
cs_real_3_t *  pvar,
const cs_real_3_t *  pvara,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_visc[],
const cs_real_t  b_visc[],
cs_real_6_t *  viscel,
const cs_real_2_t  weighf[],
const cs_real_t  weighb[],
cs_real_3_t *  rhs 
)

◆ cs_beta_limiter_building()

void cs_beta_limiter_building ( int  f_id,
int  inc,
const cs_real_t  rovsdt[] 
)

Compute the beta blending coefficient of the beta limiter (ensuring preservation of a given min/max pair of values).

Parameters
[in]f_idfield id
[in]inc"not an increment" flag
[in]rovsdtrho * volume / dt

◆ cs_cell_courant_number()

void cs_cell_courant_number ( const cs_field_t *  f,
cs_dispatch_context &  ctx,
cs_real_t *  courant 
)

◆ cs_convection_anisotropic_leff_diffusion_secvis()

void cs_convection_anisotropic_leff_diffusion_secvis ( cs_dispatch_context &  ctx,
const cs_mesh_t *  m,
const cs_mesh_quantities_t *  mq,
const cs_real_33_t  i_visc[],
const cs_real_t  i_secvis[],
const cs_rreal_t  gradv[][3][3],
cs_real_3_t *restrict  rhs 
)

Compute balance contribution of the transpose grad(vel) term and grad(-2/3 div(vel)) with anisotropic.

Compute $ \mu \transpose{\gradt\vect{\varia}} + \lambda \trace{\gradt\vect{\varia}} $, where $ \lambda $ is the secondary viscosity, i.e. usually $ -\frac{2}{3} \mu $.

Warning
  • $ \vect{Rhs} $ must already have been initialized.
  • mind the minus sign
Parameters
[in,out]ctxreference to dispatch context
[in]mpointer to mesh structure
[in]mqpointer to mesh quantities
[in]i_visc$ \mu_\fij \dfrac{S_\fij}{\ipf \jpf} $ at interior faces for the r.h.s.
[in]i_secvissecondary viscosity at interior faces
[in]b_secvissecondary viscosity at boundary faces
[in]gradvelocity gradient
[in,out]rhsright hand side $ \vect{Rhs} $

◆ cs_convection_diffusion_get_scheme_version()

int cs_convection_diffusion_get_scheme_version ( void  )

Query convection-diffusion scheme variants.

Returns
specific scheme (i.e. 90 for v9.0) or -1 for current.

◆ cs_convection_diffusion_scalar()

void cs_convection_diffusion_scalar ( const cs_field_t *  f,
const cs_equation_param_t  eqp,
int  icvflb,
int  inc,
int  imasac,
const cs_real_t *  pvar,
const int  icvfli[],
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[],
const cs_real_t *  c_weight,
cs_real_t *  rhs,
cs_real_2_t  i_flux[],
cs_real_t  b_flux[] 
)

◆ cs_convection_diffusion_secvis()

void cs_convection_diffusion_secvis ( cs_dispatch_context &  ctx,
const cs_mesh_t *  m,
const cs_mesh_quantities_t *  mq,
cs_real_t  thetap,
const cs_real_t  i_visc[],
const cs_real_t  i_secvis[],
const cs_real_t  b_secvis[],
const cs_rreal_t  gradv[][3][3],
cs_real_3_t *restrict  rhs 
)

Compute balance contribution of the transpose grad(vel) term and grad(-2/3 div(vel))

Compute $ \mu \transpose{\gradt\vect{\varia}} + \lambda \trace{\gradt\vect{\varia}} $, where $ \lambda $ is the secondary viscosity, i.e. usually $ -\frac{2}{3} \mu $.

Warning
  • $ \vect{Rhs} $ must already have been initialized.
  • mind the minus sign
Parameters
[in,out]ctxreference to dispatch context
[in]mpointer to mesh structure
[in]mqpointer to mesh quantities
[in]thetaptheta-scheme value
[in]i_visc$ \mu_\fij \dfrac{S_\fij}{\ipf \jpf} $ at interior faces for the r.h.s.
[in]i_secvissecondary viscosity at interior faces
[in]b_secvissecondary viscosity at boundary faces
[in]gradvvelocity gradient
[in,out]rhsright hand side $ \vect{Rhs} $

◆ cs_convection_diffusion_set_scheme_version()

void cs_convection_diffusion_set_scheme_version ( int  version)

Allow reverting to older convection-diffusion scheme variants.

Parameters
[in]versionScheme version (90 for v9.0) -1 for current)

◆ cs_convection_diffusion_tensor()

void cs_convection_diffusion_tensor ( int  idtvar,
int  f_id,
const cs_equation_param_t  eqp,
int  icvflb,
int  inc,
int  imasac,
cs_real_6_t *  pvar,
const cs_real_6_t *  pvara,
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_6_t *  rhs 
)

◆ cs_convection_diffusion_thermal()

void cs_convection_diffusion_thermal ( const cs_field_t *  f,
const cs_equation_param_t  eqp,
int  inc,
int  imasac,
const cs_real_t *  pvar,
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[],
const cs_real_t *  c_weight,
const cs_real_t  xcpp[],
cs_real_t *  rhs 
)

◆ cs_convection_diffusion_vector()

void cs_convection_diffusion_vector ( int  idtvar,
int  f_id,
const cs_equation_param_t  eqp,
int  icvflb,
int  inc,
int  ivisep,
int  imasac,
cs_real_3_t *  pvar,
const cs_real_3_t *  pvara,
const int  icvfli[],
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[],
const cs_real_t  i_secvis[],
const cs_real_t  b_secvis[],
cs_real_3_t *  i_pvar,
cs_real_3_t *  b_pvar,
cs_real_3_t *  rhs 
)

◆ cs_diffusion_potential()

void cs_diffusion_potential ( const cs_field_t *  f,
const cs_equation_param_t *  eqp,
const cs_mesh_t *  m,
cs_mesh_quantities_t *  fvq,
int  init,
int  inc,
int  iphydp,
cs_real_3_t *  frcxt,
cs_real_t *  pvar,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_visc[],
const cs_real_t  b_visc[],
cs_real_t  visel[],
cs_real_t *  diverg 
)

◆ cs_face_anisotropic_diffusion_potential()

void cs_face_anisotropic_diffusion_potential ( const cs_field_t *  f,
const cs_equation_param_t *  eqp,
const cs_mesh_t *  m,
cs_mesh_quantities_t *  fvq,
int  init,
int  inc,
int  iphydp,
cs_real_3_t *  frcxt,
cs_real_t *  pvar,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_visc[],
const cs_real_t  b_visc[],
cs_real_6_t *  viscel,
const cs_real_2_t  weighf[],
const cs_real_t  weighb[],
cs_real_t *  i_massflux,
cs_real_t *  b_massflux 
)

◆ cs_face_convection_scalar()

void cs_face_convection_scalar ( int  idtvar,
int  f_id,
const cs_equation_param_t  eqp,
int  icvflb,
int  inc,
int  imasac,
cs_real_t *  pvar,
const cs_real_t *  pvara,
const int  icvfli[],
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_massflux[],
const cs_real_t  b_massflux[],
cs_real_2_t  i_conv_flux[],
cs_real_t  b_conv_flux[] 
)

◆ cs_face_diffusion_potential()

void cs_face_diffusion_potential ( const cs_field_t *  f,
const cs_equation_param_t *  eqp,
const cs_mesh_t *  m,
cs_mesh_quantities_t *  fvq,
int  init,
int  inc,
int  iphydp,
cs_real_3_t *  frcxt,
cs_real_t *  pvar,
cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_visc[],
const cs_real_t  b_visc[],
cs_real_t *  visel,
cs_real_t *  i_massflux,
cs_real_t *  b_massflux 
)

◆ cs_get_v_slope_test()

cs_real_t * cs_get_v_slope_test ( int  f_id,
const cs_equation_param_t  eqp 
)

◆ cs_slope_test_gradient()

void cs_slope_test_gradient ( int  f_id,
cs_dispatch_context &  ctx,
const T(*)  grad[3],
T(*)  grdpa[3],
const cs_real_t *  pvar,
const cs_real_t  val_f[],
const cs_real_t *  i_massflux 
)

Compute the upwind gradient used in the slope tests.

This function assumes the input gradient and pvar values have already been synchronized.

Parameters
[in]f_idfield id
[in]ctxReference to dispatch context
[in]gradstandard gradient
[out]grdpaupwind gradient
[in]pvarvalues
[in]val_fface values for gradient
[in]i_massfluxmass flux at interior faces

◆ cs_slope_test_gradient_strided()

void cs_slope_test_gradient_strided ( cs_dispatch_context &  ctx,
const T  grad[][stride][3],
T(*)  grdpa[stride][3],
const cs_real_t  pvar[][stride],
const cs_real_t  val_f[][stride],
const cs_real_t *  i_massflux 
)

Compute the upwind gradient used in the slope tests.

template parameters: stride 1 for scalars, 3 for vectors, 6 for symmetric tensors

This function assumes the input gradient and pvar values have already been synchronized.

Parameters
[in]ctxReference to dispatch context
[in]gradstandard gradient
[out]grdpaupwind gradient
[in]pvarvalues
[in]val_fface values for gradient
[in]i_massfluxmass flux at interior faces

◆ cs_upwind_gradient()

void cs_upwind_gradient ( cs_dispatch_context &  ctx,
const int  inc,
const cs_halo_type_t  halo_type,
const cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_massflux[],
const cs_real_t  b_massflux[],
const cs_real_t *  pvar,
T(*)  grdpa[3] 
)

◆ cs_upwind_gradient_strided()

void cs_upwind_gradient_strided ( cs_dispatch_context &  ctx,
const int  inc,
const cs_halo_type_t  halo_type,
const cs_field_bc_coeffs_t *  bc_coeffs,
const cs_real_t  i_massflux[],
const cs_real_t  b_massflux[],
const cs_real_t *restrict  pvar[stride],
T(*)  grdpa[stride][3] 
)