#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_boundary_zone.h"
#include "cs_equation_common.h"
#include "cs_equation_priv.h"
#include "cs_hho_builder.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_scheme_geometry.h"
#include "cs_search.h"
#include "cs_sdm.h"
#include "cs_source_term.h"
#include "cs_hho_stokes.h"

Include dependency graph for cs_hho_stokes.c:

Macros
#define	CS_HHO_STOKES_DBG 1

#define	CS_HHO_STOKES_MODULO 4

#define	_dp3 cs_math_3_dot_product

Functions
static cs_cell_builder_t *	_cell_builder_create (cs_param_space_scheme_t space_scheme, const cs_cdo_connect_t *connect)
	Initialize the local builder structure used for building the system cellwise. More...

void	cs_hho_stokes_initialize (cs_flag_t scheme_flag, const cs_cdo_quantities_t quant, const cs_cdo_connect_t connect, const cs_time_step_t time_step, const cs_matrix_assembler_t ma0, const cs_matrix_assembler_t ma1, const cs_matrix_assembler_t ma2, const cs_matrix_structure_t ms0, const cs_matrix_structure_t ms1, const cs_matrix_structure_t *ms2)
	Allocate work buffer and general structures related to HHO schemes Set shared pointers. More...

void	cs_hho_stokes_get (cs_cell_sys_t csys, cs_cell_builder_t cb, cs_hho_builder_t **hhob)
	Retrieve work buffers used for building a CDO system cellwise. More...

void	cs_hho_stokes_finalize (void)
	Free buffers and generic structures related to HHO schemes. More...

void *	cs_hho_stokes_init_context (const cs_equation_param_t eqp, cs_equation_builder_t eqb)
	Initialize a cs_hho_stokes_t structure storing data useful for managing such a scheme. More...

void *	cs_hho_stokes_free_context (void *data)
	Destroy a cs_hho_stokes_t structure. More...

void	cs_hho_stokes_compute_source (const cs_equation_param_t eqp, cs_equation_builder_t eqb, void *data)
	Compute the contributions of source terms (store inside builder) More...

void	cs_hho_stokes_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_hho_stokes_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 HHO scheme. One works cellwise and then process to the assembly. More...

void	cs_hho_stokes_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 required for hybrid discretization. More...

double *	cs_hho_stokes_get_face_values (const void *data)
	Get the computed values at faces (DoF used in the linear system are located at primal faces) More...

double *	cs_hho_stokes_get_cell_values (const void *data)
	Get the computed values at cells (DoF used in the linear system are located at primal faces) More...

void	cs_hho_stokes_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...

Variables
static cs_cell_sys_t **	cs_hho_cell_sys = NULL

static cs_cell_builder_t **	cs_hho_cell_bld = NULL

static cs_hho_builder_t **	cs_hho_builders = 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

static const cs_matrix_assembler_t *	cs_shared_ma0

static const cs_matrix_structure_t *	cs_shared_ms0

static const cs_matrix_assembler_t *	cs_shared_ma1

static const cs_matrix_structure_t *	cs_shared_ms1

static const cs_matrix_assembler_t *	cs_shared_ma2

static const cs_matrix_structure_t *	cs_shared_ms2

Macro Definition Documentation

◆ _dp3

#define _dp3 cs_math_3_dot_product

◆ CS_HHO_STOKES_DBG

#define CS_HHO_STOKES_DBG 1

◆ CS_HHO_STOKES_MODULO

#define CS_HHO_STOKES_MODULO 4

Function Documentation

◆ _cell_builder_create()

static cs_cell_builder_t* _cell_builder_create	(	cs_param_space_scheme_t	space_scheme,
		const cs_cdo_connect_t *	connect
	)

static

Initialize the local builder structure used for building the system cellwise.

Parameters

[in]	space_scheme	discretization scheme
[in]	connect	pointer to a cs_cdo_connect_t structure

Returns: a pointer to a new allocated cs_cell_builder_t structure

◆ cs_hho_stokes_build_system()

void cs_hho_stokes_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 HHO scheme. One works cellwise and then process to the assembly.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	field_val	pointer to the current value of the 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_hho_stokes_t structure
[in,out]	rhs	right-hand side
[in,out]	matrix	pointer to cs_matrix_t structure to compute

◆ cs_hho_stokes_compute_source()

void cs_hho_stokes_compute_source	(	const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		void *	data
	)

Compute the contributions of source terms (store inside builder)

Parameters

[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_hho_stokes_t structure

◆ cs_hho_stokes_extra_op()

void cs_hho_stokes_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.

Parameters

[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_hho_stokes_t structure

◆ cs_hho_stokes_finalize()

void cs_hho_stokes_finalize ( void )

Free buffers and generic structures related to HHO schemes.

◆ cs_hho_stokes_free_context()

void* cs_hho_stokes_free_context ( void * data )

Destroy a cs_hho_stokes_t structure.

Parameters

[in,out] data pointer to a cs_hho_stokes_t structure

Returns: a NULL pointer

◆ cs_hho_stokes_get()

void cs_hho_stokes_get	(	cs_cell_sys_t **	csys,
		cs_cell_builder_t **	cb,
		cs_hho_builder_t **	hhob
	)

Retrieve work buffers used for building a CDO system cellwise.

Parameters

[out]	csys	pointer to a pointer on a cs_cell_sys_t structure
[out]	cb	pointer to a pointer on a cs_cell_builder_t structure
[out]	hhob	pointer to a pointer on a cs_hho_builder_t structure

◆ cs_hho_stokes_get_cell_values()

double* cs_hho_stokes_get_cell_values ( const void * data )

Get the computed values at cells (DoF used in the linear system are located at primal faces)

Parameters

[in] data pointer to a data structure

Returns: a pointer to an array of double

◆ cs_hho_stokes_get_face_values()

double* cs_hho_stokes_get_face_values ( const void * data )

Get the computed values at faces (DoF used in the linear system are located at primal faces)

Parameters

[in] data pointer to a data structure

Returns: a pointer to an array of double

◆ cs_hho_stokes_init_context()

void* cs_hho_stokes_init_context	(	const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb
	)

Initialize a cs_hho_stokes_t structure storing data useful for managing such a scheme.

Parameters

[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure

Returns: a pointer to a new allocated cs_hho_stokes_t structure

◆ cs_hho_stokes_initialize()

void cs_hho_stokes_initialize	(	cs_flag_t	scheme_flag,
		const cs_cdo_quantities_t *	quant,
		const cs_cdo_connect_t *	connect,
		const cs_time_step_t *	time_step,
		const cs_matrix_assembler_t *	ma0,
		const cs_matrix_assembler_t *	ma1,
		const cs_matrix_assembler_t *	ma2,
		const cs_matrix_structure_t *	ms0,
		const cs_matrix_structure_t *	ms1,
		const cs_matrix_structure_t *	ms2
	)

Allocate work buffer and general structures related to HHO schemes Set shared pointers.

Parameters

[in]	scheme_flag	flag to identify which kind of numerical scheme is requested to solve the computational domain
[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
[in]	ma0	pointer to a cs_matrix_assembler_t structure (P0)
[in]	ma1	pointer to a cs_matrix_assembler_t structure (P1)
[in]	ma2	pointer to a cs_matrix_assembler_t structure (P2)
[in]	ms0	pointer to a cs_matrix_structure_t structure (P0)
[in]	ms1	pointer to a cs_matrix_structure_t structure (P1)
[in]	ms2	pointer to a cs_matrix_structure_t structure (P2)

◆ cs_hho_stokes_initialize_system()

void cs_hho_stokes_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.

Parameters

[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 generic data 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

◆ cs_hho_stokes_update_field()

void cs_hho_stokes_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 required for hybrid discretization.

Parameters

[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_hho_stokes_t structure
[in,out]	field_val	pointer to the current value of the field

Variable Documentation

◆ cs_hho_builders

cs_hho_builder_t** cs_hho_builders = NULL

static

◆ cs_hho_cell_bld

cs_cell_builder_t** cs_hho_cell_bld = NULL

static

◆ cs_hho_cell_sys

cs_cell_sys_t** cs_hho_cell_sys = NULL

static

◆ cs_shared_connect

const cs_cdo_connect_t* cs_shared_connect

static

◆ cs_shared_ma0

const cs_matrix_assembler_t* cs_shared_ma0

static

◆ cs_shared_ma1

const cs_matrix_assembler_t* cs_shared_ma1

static

◆ cs_shared_ma2

const cs_matrix_assembler_t* cs_shared_ma2

static

◆ cs_shared_ms0

const cs_matrix_structure_t* cs_shared_ms0

static

◆ cs_shared_ms1

const cs_matrix_structure_t* cs_shared_ms1

static

◆ cs_shared_ms2

const cs_matrix_structure_t* cs_shared_ms2

static

◆ cs_shared_quant

const cs_cdo_quantities_t* cs_shared_quant

static

◆ cs_shared_time_step

const cs_time_step_t* cs_shared_time_step

static

Macros

Functions

Variables

Macro Definition Documentation

◆ _dp3

◆ CS_HHO_STOKES_DBG

◆ CS_HHO_STOKES_MODULO

Function Documentation

◆ _cell_builder_create()

◆ cs_hho_stokes_build_system()

◆ cs_hho_stokes_compute_source()

◆ cs_hho_stokes_extra_op()

◆ cs_hho_stokes_finalize()

◆ cs_hho_stokes_free_context()

◆ cs_hho_stokes_get()

◆ cs_hho_stokes_get_cell_values()

◆ cs_hho_stokes_get_face_values()

◆ cs_hho_stokes_init_context()

◆ cs_hho_stokes_initialize()

◆ cs_hho_stokes_initialize_system()

◆ cs_hho_stokes_update_field()

Variable Documentation

◆ cs_hho_builders

◆ cs_hho_cell_bld

◆ cs_hho_cell_sys

◆ cs_shared_connect

◆ cs_shared_ma0

◆ cs_shared_ma1

◆ cs_shared_ma2

◆ cs_shared_ms0

◆ cs_shared_ms1

◆ cs_shared_ms2

◆ cs_shared_quant

◆ cs_shared_time_step