cs_cdofb_vecteq.c File Reference

#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_bc.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_static_condensation.h"
#include "cs_cdofb_priv.h"
#include "cs_cdofb_vecteq.h"

Include dependency graph for cs_cdofb_vecteq.c:

Functions
void	cs_cdofb_vecteq_init_cell_system (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_vecteq_t *eqc, const cs_real_t dir_values[], const short int neu_tags[], const cs_real_t field_tn[], cs_real_t t_eval, cs_cell_sys_t *csys, cs_cell_builder_t *cb)
	Initialize the local structure for the current cell. More...
bool	cs_cdofb_vecteq_is_initialized (void)
	Check if the generic structures for building a CDO-Fb scheme are allocated. More...
void	cs_cdofb_vecteq_init_common (const cs_cdo_quantities_t *quant, const cs_cdo_connect_t *connect, const cs_time_step_t *time_step, const cs_matrix_structure_t *ms)
	Allocate work buffer and general structures related to CDO vector-valued face-based schemes. Set shared pointers from the main domain members. More...
const cs_matrix_structure_t *	cs_cdofb_vecteq_matrix_structure (void)
	Get the pointer to the related cs_matrix_structure_t. More...
void	cs_cdofb_vecteq_get (cs_cell_sys_t **csys, cs_cell_builder_t **cb)
	Retrieve work buffers used for building a CDO system cellwise. More...
void	cs_cdofb_vecteq_finalize_common (void)
	Free work buffer and general structure related to CDO face-based schemes. More...
void *	cs_cdofb_vecteq_init_context (const cs_equation_param_t *eqp, int var_id, int bflux_id, cs_equation_builder_t *eqb)
	Initialize a cs_cdofb_vecteq_t structure storing data useful for building and managing such a scheme. More...
void *	cs_cdofb_vecteq_free_context (void *data)
	Destroy a cs_cdofb_vecteq_t structure. More...
void	cs_cdofb_vecteq_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_vecteq_set_dir_bc (const cs_mesh_t *mesh, const cs_equation_param_t *eqp, cs_equation_builder_t *eqb, cs_real_t t_eval, cs_real_t field_val[])
	Set the boundary conditions known from the settings when the fields stem from a vector CDO face-based scheme. More...
void	cs_cdofb_vecteq_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 vector convection/diffusion equation with a CDO face-based scheme. One works cellwise and then process to the assembly. More...
void	cs_cdofb_vecteq_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_vecteq_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...
cs_real_t *	cs_cdofb_vecteq_get_cell_values (void *context)
	Get the computed values at mesh cells from the inverse operation w.r.t. the static condensation (DoF used in the linear system are located at primal faces) The lifecycle of this array is managed by the code. So one does not have to free the return pointer. More...
cs_real_t *	cs_cdofb_vecteq_get_face_values (void *context)
	Retrieve an array of values at mesh faces for the current context. The lifecycle of this array is managed by the code. So one does not have to free the return pointer. More...
void	cs_cdofb_vecteq_read_restart (cs_restart_t *restart, const char *eqname, void *scheme_context)
	Read additional arrays (not defined as fields) but useful for the checkpoint/restart process. More...
void	cs_cdofb_vecteq_write_restart (cs_restart_t *restart, const char *eqname, void *scheme_context)
	Write additional arrays (not defined as fields) but useful for the checkpoint/restart process. More...

Function Documentation

cs_cdofb_vecteq_build_system()

void cs_cdofb_vecteq_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 vector convection/diffusion equation with a CDO face-based 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_cdofb_vecteq_t structure
[in,out]	rhs	right-hand side
[in,out]	matrix	pointer to cs_matrix_t structure to compute

cs_cdofb_vecteq_extra_op()

void cs_cdofb_vecteq_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_cdofb_vecteq_t structure

cs_cdofb_vecteq_finalize_common()

void cs_cdofb_vecteq_finalize_common

(

void

)

Free work buffer and general structure related to CDO face-based schemes.

cs_cdofb_vecteq_free_context()

void* cs_cdofb_vecteq_free_context

(

void *

data

)

Destroy a cs_cdofb_vecteq_t structure.

Parameters

[in,out]

data

pointer to a cs_cdofb_vecteq_t structure

Returns

a NULL pointer

cs_cdofb_vecteq_get()

void cs_cdofb_vecteq_get	(	cs_cell_sys_t **	csys,
		cs_cell_builder_t **	cb
	)

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

cs_cdofb_vecteq_get_cell_values()

cs_real_t* cs_cdofb_vecteq_get_cell_values

(

void *

context

)

Get the computed values at mesh cells from the inverse operation w.r.t. the static condensation (DoF used in the linear system are located at primal faces) The lifecycle of this array is managed by the code. So one does not have to free the return pointer.

Parameters

[in,out]

context

pointer to a data structure cast on-the-fly

Returns

a pointer to an array of cs_real_t

cs_cdofb_vecteq_get_face_values()

cs_real_t* cs_cdofb_vecteq_get_face_values

(

void *

context

)

Retrieve an array of values at mesh faces for the current context. The lifecycle of this array is managed by the code. So one does not have to free the return pointer.

Parameters

[in,out]

context

pointer to a data structure cast on-the-fly

Returns

a pointer to an array of cs_real_t (size n_faces)

cs_cdofb_vecteq_init_cell_system()

void cs_cdofb_vecteq_init_cell_system	(	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_vecteq_t *	eqc,
		const cs_real_t	dir_values[],
		const short int	neu_tags[],
		const cs_real_t	field_tn[],
		cs_real_t	t_eval,
		cs_cell_sys_t *	csys,
		cs_cell_builder_t *	cb
	)

Initialize the local structure for the current cell.

DOXYGEN_SHOULD_SKIP_THIS

Parameters

[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]	eqc	pointer to a cs_cdofb_vecteq_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]	t_eval	time at which one performs the evaluation
[in,out]	csys	pointer to a cellwise view of the system
[in,out]	cb	pointer to a cellwise builder

cs_cdofb_vecteq_init_common()

void cs_cdofb_vecteq_init_common	(	const cs_cdo_quantities_t *	quant,
		const cs_cdo_connect_t *	connect,
		const cs_time_step_t *	time_step,
		const cs_matrix_structure_t *	ms
	)

Allocate work buffer and general structures related to CDO vector-valued face-based schemes. Set shared pointers from the main domain members.

Parameters

[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]	ms	pointer to a cs_matrix_structure_t structure

cs_cdofb_vecteq_init_context()

void* cs_cdofb_vecteq_init_context	(	const cs_equation_param_t *	eqp,
		int	var_id,
		int	bflux_id,
		cs_equation_builder_t *	eqb
	)

Initialize a cs_cdofb_vecteq_t structure storing data useful for building and managing such a scheme.

Parameters

[in]	eqp	pointer to a cs_equation_param_t structure
[in]	var_id	id of the variable field
[in]	bflux_id	id of the boundary flux field
[in,out]	eqb	pointer to a cs_equation_builder_t structure

Returns

a pointer to a new allocated cs_cdofb_vecteq_t structure

cs_cdofb_vecteq_initialize_system()

void cs_cdofb_vecteq_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 cs_cdofb_vecteq_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

cs_cdofb_vecteq_is_initialized()

bool cs_cdofb_vecteq_is_initialized

(

void

)

Check if the generic structures for building a CDO-Fb scheme are allocated.

Returns

true or false

cs_cdofb_vecteq_matrix_structure()

const cs_matrix_structure_t* cs_cdofb_vecteq_matrix_structure

(

void

)

Get the pointer to the related cs_matrix_structure_t.

Returns

a pointer to a cs_matrix_structure_t structure

cs_cdofb_vecteq_read_restart()

void cs_cdofb_vecteq_read_restart	(	cs_restart_t *	restart,
		const char *	eqname,
		void *	scheme_context
	)

Read additional arrays (not defined as fields) but useful for the checkpoint/restart process.

Parameters

[in,out]	restart	pointer to cs_restart_t structure
[in]	eqname	name of the related equation
[in]	scheme_context	pointer to a data structure cast on-the-fly

cs_cdofb_vecteq_set_dir_bc()

void cs_cdofb_vecteq_set_dir_bc	(	const cs_mesh_t *	mesh,
		const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		cs_real_t	t_eval,
		cs_real_t	field_val[]
	)

Set the boundary conditions known from the settings when the fields stem from a vector CDO face-based scheme.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in]	t_eval	time at which one evaluates BCs
[in,out]	field_val	pointer to the values of the variable field

cs_cdofb_vecteq_update_field()

void cs_cdofb_vecteq_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)

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

cs_cdofb_vecteq_write_restart()

void cs_cdofb_vecteq_write_restart	(	cs_restart_t *	restart,
		const char *	eqname,
		void *	scheme_context
	)

Write additional arrays (not defined as fields) but useful for the checkpoint/restart process.

Parameters

[in,out]	restart	pointer to cs_restart_t structure
[in]	eqname	name of the related equation
[in]	scheme_context	pointer to a data structure cast on-the-fly