7.0
general documentation
cs_cdofb_priv.h
Go to the documentation of this file.
1 #ifndef __CS_CDOFB_PRIV_H__
2 #define __CS_CDOFB_PRIV_H__
3 
4 /*============================================================================
5  * Definition of cs_cdofb_scaleq_t and cs_cdofb_vecteq structures
6  *============================================================================*/
7 
8 /*
9  This file is part of Code_Saturne, a general-purpose CFD tool.
10 
11  Copyright (C) 1998-2021 EDF S.A.
12 
13  This program is free software; you can redistribute it and/or modify it under
14  the terms of the GNU General Public License as published by the Free Software
15  Foundation; either version 2 of the License, or (at your option) any later
16  version.
17 
18  This program is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
20  FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
21  details.
22 
23  You should have received a copy of the GNU General Public License along with
24  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
25  Street, Fifth Floor, Boston, MA 02110-1301, USA.
26 */
27 
28 /*----------------------------------------------------------------------------
29  * Local headers
30  *----------------------------------------------------------------------------*/
31 
32 #include "cs_defs.h"
33 #include "cs_hodge.h"
34 #include "cs_cdo_advection.h"
35 #include "cs_equation_assemble.h"
36 #include "cs_equation_bc.h"
37 #include "cs_equation_common.h"
38 
39 /*----------------------------------------------------------------------------*/
40 
42 
43 /*============================================================================
44  * Macro definitions
45  *============================================================================*/
46 
47 /*============================================================================
48  * Type definitions
49  *============================================================================*/
50 
51 /* Algebraic system for CDO face-based discretization */
52 struct _cs_cdofb_t {
53 
54  /* Ids related to the variable field and to the boundary flux field */
57 
58  /* System size (n_faces + n_cells) */
61 
62  /* Solution of the algebraic system DoF unknowns (x) + BCs */
63  cs_real_t *face_values; /* At the last iteration */
64  cs_real_t *face_values_pre; /* At the previous iteration */
65 
66  /* Assembly process */
68 
69  /* Members related to the static condensation */
70  cs_real_t *rc_tilda; /* Acc^-1 * RHS_cell */
71  cs_real_t *acf_tilda; /* Acc^-1 * Acf
72  Cell-faces lower-left block of the full matrix
73  Access to the values thanks to the c2f
74  connectivity */
75 
76  /* Array storing the value arising from the contribution of all source
77  terms (only allocated to n_cells) */
79 
80  /* Pointer of function to build the diffusion term */
86 
87  /* Pointer of functions to define the advection term:
88  * advection_open is called first, then advection_build which calls
89  * advection_scheme and after the build step, advection is called last
90  */
95 
97 
98  /* If one needs to build a local hodge op. for time and reaction */
102 };
103 
104 typedef struct _cs_cdofb_t cs_cdofb_priv_t;
105 
106 /*============================================================================
107  * Public function prototypes
108  *============================================================================*/
109 
110 /*----------------------------------------------------------------------------*/
119 /*----------------------------------------------------------------------------*/
120 
121 void
124  cs_cdofb_priv_t *eqc);
125 
126 /*----------------------------------------------------------------------------*/
127 
129 
130 #endif /* __CS_CDOFB_PRIV_H__ */
cs_cdo_enforce_bc_t * enforce_dirichlet
Definition: cs_cdofb_priv.h:83
Store common elements used when building an algebraic system related to an equation.
cs_equation_assembly_t * assemble
Definition: cs_cdofb_priv.h:67
cs_cdofb_adv_build_t * advection_build
Definition: cs_cdofb_priv.h:92
cs_lnum_t n_faces
Definition: cs_cdofb_priv.h:59
void() cs_cdofb_adv_close_hook_t(const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, cs_cell_sys_t *csys, cs_cell_builder_t *cb, cs_sdm_t *adv)
Operation done after the matrix related to the advection term has been defined.
Definition: cs_cdo_advection.h:136
void() cs_cdofb_adv_build_t(const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_cell_sys_t *csys, cs_cdofb_adv_scheme_t *scheme_func, cs_cell_builder_t *cb)
Build the cellwise advection operator for CDO-Fb schemes The local matrix related to this operator is...
Definition: cs_cdo_advection.h:116
Set of parameters to handle an unsteady convection-diffusion-reaction equation with term sources...
Definition: cs_equation_param.h:202
cs_real_t * rc_tilda
Definition: cs_cdofb_priv.h:70
cs_hodge_t ** diffusion_hodge
Definition: cs_cdofb_priv.h:81
void() cs_cdo_enforce_bc_t(const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, cs_face_mesh_t *fm, cs_hodge_t *hodge, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Enforcement of a boundary condition (Dirichlet, Robin, sliding...)
Definition: cs_equation_bc.h:89
int var_field_id
Definition: cs_cdofb_priv.h:55
cs_hodge_compute_t * get_mass_matrix
Definition: cs_cdofb_priv.h:101
#define BEGIN_C_DECLS
Definition: cs_defs.h:495
void() cs_equation_assembly_t(const cs_sdm_t *m, const cs_lnum_t *dof_ids, const cs_range_set_t *rset, cs_equation_assemble_t *eqa, cs_matrix_assembler_values_t *mav)
Assemble a cellwise matrix into the global matrix Block or no block versions are handled.
Definition: cs_equation_assemble.h:72
cs_cdofb_adv_scheme_t * advection_scheme
Definition: cs_cdofb_priv.h:93
cs_hodge_compute_t * get_stiffness_matrix
Definition: cs_cdofb_priv.h:82
cs_hodge_param_t mass_hodgep
Definition: cs_cdofb_priv.h:99
Structure associated to a discrete Hodge operator.
Definition: cs_hodge.h:189
cs_cdofb_adv_close_hook_t * advection_close
Definition: cs_cdofb_priv.h:94
cs_hodge_t ** mass_hodge
Definition: cs_cdofb_priv.h:100
double cs_real_t
Floating-point value.
Definition: cs_defs.h:307
void() cs_hodge_compute_t(const cs_cell_mesh_t *cm, cs_hodge_t *hodge, cs_cell_builder_t *cb)
Build a discrete Hodge operator or a related operator (such as the stiffmess matrix) for a given cell...
Definition: cs_hodge.h:213
int bflux_field_id
Definition: cs_cdofb_priv.h:56
cs_cdo_enforce_bc_t * enforce_sliding
Definition: cs_cdofb_priv.h:85
cs_real_t * face_values_pre
Definition: cs_cdofb_priv.h:64
void * advection_input
Definition: cs_cdofb_priv.h:96
cs_real_t * source_terms
Definition: cs_cdofb_priv.h:78
cs_lnum_t n_dofs
Definition: cs_cdofb_priv.h:60
cs_real_t * acf_tilda
Definition: cs_cdofb_priv.h:71
cs_cdo_enforce_bc_t * enforce_robin_bc
Definition: cs_cdofb_priv.h:84
void() cs_cdofb_adv_open_hook_t(const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_cell_sys_t *csys, void *input, cs_cell_builder_t *cb)
Perform preprocessing such as the computation of the advection flux at the expected location in order...
Definition: cs_cdo_advection.h:71
int cs_lnum_t
local mesh entity id
Definition: cs_defs.h:301
#define END_C_DECLS
Definition: cs_defs.h:496
void() cs_cdofb_adv_scheme_t(int dim, const cs_cell_mesh_t *cm, const cs_cell_sys_t *csys, cs_cell_builder_t *cb, cs_sdm_t *adv)
Define the local convection operator in CDO-Fb schemes. Case of an operator that should be assemble i...
Definition: cs_cdo_advection.h:92
cs_cdofb_adv_open_hook_t * advection_open
Definition: cs_cdofb_priv.h:91
Definition: cs_cdofb_priv.h:52
void cs_cdofb_set_advection_function(const cs_equation_param_t *eqp, cs_equation_builder_t *eqb, cs_cdofb_priv_t *eqc)
Set the advection-related parameters in the context structure of CDO face-based schemes.
Definition: cs_cdofb_priv.c:90
cs_real_t * face_values
Definition: cs_cdofb_priv.h:63
Structure storing all metadata/parameters related to the usage of a discrete Hodge operator...
Definition: cs_hodge.h:165