![]() |
programmer's documentation
|
#include "cs_defs.h"
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>
#include <assert.h>
#include <string.h>
#include <bft_mem.h>
#include "cs_cdo_advection.h"
#include "cs_cdo_bc.h"
#include "cs_cdo_diffusion.h"
#include "cs_equation_common.h"
#include "cs_hodge.h"
#include "cs_log.h"
#include "cs_math.h"
#include "cs_mesh_location.h"
#include "cs_post.h"
#include "cs_quadrature.h"
#include "cs_reco.h"
#include "cs_search.h"
#include "cs_source_term.h"
#include "cs_cdofb_scaleq.h"
Macros | |
#define | CS_CDOFB_SCALEQ_DBG 0 |
#define | CS_CDOFB_SCALEQ_MODULO 10 |
Functions | |
static void | _init_cell_structures (const cs_flag_t cell_flag, const cs_cell_mesh_t *cm, const cs_equation_param_t *eqp, const cs_equation_builder_t *eqb, const cs_cdofb_scaleq_t *eqd, const cs_real_t dir_values[], const short int neu_tags[], const cs_real_t field_tn[], cs_cell_sys_t *csys, cs_cell_bc_t *cbc, cs_cell_builder_t *cb) |
Initialize the local structure for the current cell. More... | |
static void | _condense_and_store (const cs_adjacency_t *c2f, cs_cdofb_scaleq_t *eqd, cs_cell_builder_t *cb, cs_cell_sys_t *csys) |
Proceed to a static condensation of the local system and keep information inside the builder to be able to compute the values at cell centers. More... | |
void | cs_cdofb_scaleq_set_shared_pointers (const cs_cdo_quantities_t *quant, const cs_cdo_connect_t *connect, const cs_time_step_t *time_step) |
Set shared pointers from the main domain members. More... | |
void | cs_cdofb_scaleq_initialize (void) |
Allocate work buffer and general structures related to CDO face-based schemes. More... | |
void | cs_cdofb_scaleq_finalize (void) |
Free work buffer and general structure related to CDO face-based schemes. More... | |
void * | cs_cdofb_scaleq_init_data (const cs_equation_param_t *eqp, cs_equation_builder_t *eqb) |
Initialize a cs_cdofb_scaleq_t structure storing data useful for managing such a scheme. More... | |
void * | cs_cdofb_scaleq_free_data (void *data) |
Destroy a cs_cdofb_scaleq_t structure. More... | |
void | cs_cdofb_scaleq_compute_source (const cs_equation_param_t *eqp, cs_equation_builder_t *eqb, void *data) |
Compute the contributions of source terms (store inside data) More... | |
void | cs_cdofb_scaleq_initialize_system (const cs_equation_param_t *eqp, cs_equation_builder_t *eqb, void *data, cs_matrix_t **system_matrix, cs_real_t **system_rhs) |
Create the matrix of the current algebraic system. Allocate and initialize the right-hand side associated to the given data structure. More... | |
void | cs_cdofb_scaleq_build_system (const cs_mesh_t *mesh, const cs_real_t *field_val, double dt_cur, const cs_equation_param_t *eqp, cs_equation_builder_t *eqb, void *data, cs_real_t *rhs, cs_matrix_t *matrix) |
Build the linear system arising from a scalar convection/diffusion equation with a CDO face-based scheme. One works cellwise and then process to the assembly. More... | |
void | cs_cdofb_scaleq_update_field (const cs_real_t *solu, const cs_real_t *rhs, const cs_equation_param_t *eqp, cs_equation_builder_t *eqb, void *data, cs_real_t *field_val) |
Store solution(s) of the linear system into a field structure Update extra-field values if required (for hybrid discretization) More... | |
void | cs_cdofb_scaleq_extra_op (const char *eqname, const cs_field_t *field, const cs_equation_param_t *eqp, cs_equation_builder_t *eqb, void *data) |
Predefined extra-operations related to this equation. More... | |
double * | cs_cdofb_scaleq_get_face_values (const void *data) |
Get the computed values at each face. More... | |
Variables | |
static cs_cell_sys_t ** | cs_cdofb_cell_sys = NULL |
static cs_cell_bc_t ** | cs_cdofb_cell_bc = NULL |
static cs_cell_builder_t ** | cs_cdofb_cell_bld = NULL |
static const cs_cdo_quantities_t * | cs_shared_quant |
static const cs_cdo_connect_t * | cs_shared_connect |
static const cs_time_step_t * | cs_shared_time_step |
#define CS_CDOFB_SCALEQ_DBG 0 |
#define CS_CDOFB_SCALEQ_MODULO 10 |
|
static |
Proceed to a static condensation of the local system and keep information inside the builder to be able to compute the values at cell centers.
[in] | c2f | pointer to a cs_adjacency_t structure |
[in,out] | eqd | pointer to a cs_cdovcb_scaleq_t structure |
[in,out] | cb | pointer to a cs_cell_builder_t structure |
[in,out] | csys | pointer to a cs_cell_sys_t structure |
|
static |
Initialize the local structure for the current cell.
[in] | cell_flag | flag related to the current cell |
[in] | cm | pointer to a cellwise view of the mesh |
[in] | eqp | pointer to a cs_equation_param_t structure |
[in] | eqb | pointer to a cs_equation_builder_t structure |
[in] | eqd | pointer to a cs_cdofb_scaleq_t structure |
[in] | dir_values | Dirichlet values associated to each face |
[in] | neu_tags | Definition id related to each Neumann face |
[in] | field_tn | values of the field at the last computed time |
[in,out] | csys | pointer to a cellwise view of the system |
[in,out] | cbc | pointer to a cellwise view of the BCs |
[in,out] | cb | pointer to a cellwise builder |
void cs_cdofb_scaleq_build_system | ( | const cs_mesh_t * | mesh, |
const cs_real_t * | field_val, | ||
double | dt_cur, | ||
const cs_equation_param_t * | eqp, | ||
cs_equation_builder_t * | eqb, | ||
void * | data, | ||
cs_real_t * | rhs, | ||
cs_matrix_t * | matrix | ||
) |
Build the linear system arising from a scalar convection/diffusion equation with a CDO face-based scheme. One works cellwise and then process to the assembly.
[in] | mesh | pointer to a cs_mesh_t structure |
[in] | field_val | pointer to the current value of the vertex field |
[in] | dt_cur | current value of the time step |
[in] | eqp | pointer to a cs_equation_param_t structure |
[in,out] | eqb | pointer to a cs_equation_builder_t structure |
[in,out] | data | pointer to cs_cdofb_scaleq_t structure |
[in,out] | rhs | right-hand side |
[in,out] | matrix | pointer to cs_matrix_t structure to compute |
void cs_cdofb_scaleq_compute_source | ( | const cs_equation_param_t * | eqp, |
cs_equation_builder_t * | eqb, | ||
void * | data | ||
) |
Compute the contributions of source terms (store inside data)
[in] | eqp | pointer to a cs_equation_param_t structure |
[in,out] | eqb | pointer to a cs_equation_builder_t structure |
[in,out] | data | pointer to a cs_cdofb_scaleq_t structure |
void cs_cdofb_scaleq_extra_op | ( | const char * | eqname, |
const cs_field_t * | field, | ||
const cs_equation_param_t * | eqp, | ||
cs_equation_builder_t * | eqb, | ||
void * | data | ||
) |
Predefined extra-operations related to this equation.
[in] | eqname | name of the equation |
[in] | field | pointer to a field structure |
[in] | eqp | pointer to a cs_equation_param_t structure |
[in,out] | eqb | pointer to a cs_equation_builder_t structure |
[in,out] | data | pointer to cs_cdofb_scaleq_t structure |
void cs_cdofb_scaleq_finalize | ( | void | ) |
Free work buffer and general structure related to CDO face-based schemes.
void* cs_cdofb_scaleq_free_data | ( | void * | data | ) |
Destroy a cs_cdofb_scaleq_t structure.
[in,out] | data | pointer to a cs_cdofb_scaleq_t structure |
double* cs_cdofb_scaleq_get_face_values | ( | const void * | data | ) |
Get the computed values at each face.
[in] | data | pointer to cs_cdofb_scaleq_t structure |
void* cs_cdofb_scaleq_init_data | ( | const cs_equation_param_t * | eqp, |
cs_equation_builder_t * | eqb | ||
) |
Initialize a cs_cdofb_scaleq_t structure storing data useful for managing such a scheme.
[in] | eqp | pointer to a cs_equation_param_t structure |
[in,out] | eqb | pointer to a cs_equation_builder_t structure |
void cs_cdofb_scaleq_initialize | ( | void | ) |
Allocate work buffer and general structures related to CDO face-based schemes.
void cs_cdofb_scaleq_initialize_system | ( | const cs_equation_param_t * | eqp, |
cs_equation_builder_t * | eqb, | ||
void * | data, | ||
cs_matrix_t ** | system_matrix, | ||
cs_real_t ** | system_rhs | ||
) |
Create the matrix of the current algebraic system. Allocate and initialize the right-hand side associated to the given data structure.
[in] | eqp | pointer to a cs_equation_param_t structure |
[in,out] | eqb | pointer to a cs_equation_builder_t structure |
[in,out] | data | pointer to cs_cdofb_scaleq_t structure |
[in,out] | system_matrix | pointer of pointer to a cs_matrix_t struct. |
[in,out] | system_rhs | pointer of pointer to an array of double |
void cs_cdofb_scaleq_set_shared_pointers | ( | const cs_cdo_quantities_t * | quant, |
const cs_cdo_connect_t * | connect, | ||
const cs_time_step_t * | time_step | ||
) |
Set shared pointers from the main domain members.
[in] | quant | additional mesh quantities struct. |
[in] | connect | pointer to a cs_cdo_connect_t struct. |
[in] | time_step | pointer to a time step structure |
void cs_cdofb_scaleq_update_field | ( | const cs_real_t * | solu, |
const cs_real_t * | rhs, | ||
const cs_equation_param_t * | eqp, | ||
cs_equation_builder_t * | eqb, | ||
void * | data, | ||
cs_real_t * | field_val | ||
) |
Store solution(s) of the linear system into a field structure Update extra-field values if required (for hybrid discretization)
[in] | solu | solution array |
[in] | rhs | rhs associated to this solution array |
[in] | eqp | pointer to a cs_equation_param_t structure |
[in,out] | eqb | pointer to a cs_equation_builder_t structure |
[in,out] | data | pointer to cs_cdofb_scaleq_t structure |
[in,out] | field_val | pointer to the current value of the field |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |