cs_cdofb_monolithic_sles.h File Reference

#include "cs_defs.h"
#include "cs_navsto_param.h"

Include dependency graph for cs_cdofb_monolithic_sles.h:

Go to the source code of this file.

Data Structures
struct	cs_cdofb_monolithic_sles_t

Functions
cs_cdofb_monolithic_sles_t *	cs_cdofb_monolithic_sles_create (void)
	Create an empty cs_cdofb_monolithic_sles_t structure. More...
void	cs_cdofb_monolithic_sles_init (cs_lnum_t n_cells, cs_lnum_t n_faces, cs_cdofb_monolithic_sles_t *msles)
	Allocate and initialize the rhs. More...
void	cs_cdofb_monolithic_sles_reset (cs_cdofb_monolithic_sles_t *msles)
	Reset to zero rhs and clean the cs_sles_t structure. More...
void	cs_cdofb_monolithic_sles_clean (cs_cdofb_monolithic_sles_t *msles)
	Free a part of the structure. More...
void	cs_cdofb_monolithic_sles_free (cs_cdofb_monolithic_sles_t **p_msles)
	Free memory related to cs_cdofb_monolithic_sles_t structure. More...
void	cs_cdofb_monolithic_sles_set_shared (const cs_cdo_connect_t *connect, const cs_cdo_quantities_t *quant, const cs_range_set_t *rset)
	Set pointers to shared structures. More...
void	cs_cdofb_monolithic_set_sles (cs_navsto_param_t *nsp, void *context)
	Start setting-up the Navier-Stokes equations when a monolithic algorithm is used to couple the system. No mesh information is available at this stage. nsp is not declared as cnost to avoid compilation warnings but it should be modified at this stage. More...
int	cs_cdofb_monolithic_solve (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Solve a linear system arising from the discretization of the Navier-Stokes equation with a CDO face-based approach. The full system is treated as one block and solved as it is. In this situation, PETSc or MUMPS are usually considered. More...
int	cs_cdofb_monolithic_krylov_block_precond (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Solve a linear system arising from the discretization of the Navier-Stokes equation with a CDO face-based approach. The system is split into blocks to enable more efficient preconditioning techniques. The main iterative solver is a Krylov solver such as GCR, GMRES or MINRES. More...
int	cs_cdofb_monolithic_by_blocks_solve (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Solve a linear system arising from the discretization of the Navier-Stokes equation with a CDO face-based approach. The system is split into blocks to enable more efficient preconditioning techniques. More...
int	cs_cdofb_monolithic_gkb_solve (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Use the GKB algorithm to solve the saddle-point problem arising from CDO-Fb schemes for Stokes and Navier-Stokes with a monolithic coupling. More...
int	cs_cdofb_monolithic_uzawa_cg_solve (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Use the preconditioned Uzawa-CG algorithm to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling This algorithm is based on Koko's paper "Uzawa conjugate gradient method for the Stokes problem: Matlab implementation with P1-iso-P2/ P1 finite element". More...
int	cs_cdofb_monolithic_uzawa_n3s_solve (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Use the preconditioned Uzawa-CG algorithm to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling This algorithm is based on the EDF report H-E40-1991-03299-FR devoted the numerical algorithms used in the code N3S. Specifically a Cahout-Chabard preconditioning is used to approximate the Schur complement. More...
int	cs_cdofb_monolithic_uzawa_al_solve (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Use the Uzawa algorithm with an Augmented Lagrangian technique to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling. More...
int	cs_cdofb_monolithic_uzawa_al_incr_solve (const cs_navsto_param_t *nsp, const cs_equation_param_t *eqp, cs_cdofb_monolithic_sles_t *msles)
	Use the Uzawa algorithm with an Augmented Lagrangian technique in an incremental way to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling. More...

Function Documentation

cs_cdofb_monolithic_by_blocks_solve()

int cs_cdofb_monolithic_by_blocks_solve	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Solve a linear system arising from the discretization of the Navier-Stokes equation with a CDO face-based approach. The system is split into blocks to enable more efficient preconditioning techniques.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the (cumulated) number of iterations of the solver

cs_cdofb_monolithic_gkb_solve()

int cs_cdofb_monolithic_gkb_solve	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Use the GKB algorithm to solve the saddle-point problem arising from CDO-Fb schemes for Stokes and Navier-Stokes with a monolithic coupling.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the cumulated number of iterations of the solver

cs_cdofb_monolithic_krylov_block_precond()

int cs_cdofb_monolithic_krylov_block_precond	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Solve a linear system arising from the discretization of the Navier-Stokes equation with a CDO face-based approach. The system is split into blocks to enable more efficient preconditioning techniques. The main iterative solver is a Krylov solver such as GCR, GMRES or MINRES.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the (cumulated) number of iterations of the solver

cs_cdofb_monolithic_set_sles()

void cs_cdofb_monolithic_set_sles	(	cs_navsto_param_t *	nsp,
		void *	context
	)

Start setting-up the Navier-Stokes equations when a monolithic algorithm is used to couple the system. No mesh information is available at this stage. nsp is not declared as cnost to avoid compilation warnings but it should be modified at this stage.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in,out]	context	pointer to a context structure cast on-the-fly

cs_cdofb_monolithic_sles_clean()

void cs_cdofb_monolithic_sles_clean

(

cs_cdofb_monolithic_sles_t *

msles

)

Free a part of the structure.

Parameters

[in,out]

msles

pointer to the structure to clean

cs_cdofb_monolithic_sles_create()

cs_cdofb_monolithic_sles_t* cs_cdofb_monolithic_sles_create

(

void

)

Create an empty cs_cdofb_monolithic_sles_t structure.

Returns

a pointer to a newly allocated structure

cs_cdofb_monolithic_sles_free()

void cs_cdofb_monolithic_sles_free

(

cs_cdofb_monolithic_sles_t **

p_msles

)

Free memory related to cs_cdofb_monolithic_sles_t structure.

Parameters

[in,out]

p_msles

double pointer to the structure to free

cs_cdofb_monolithic_sles_init()

void cs_cdofb_monolithic_sles_init	(	cs_lnum_t	n_cells,
		cs_lnum_t	n_faces,
		cs_cdofb_monolithic_sles_t *	msles
	)

Allocate and initialize the rhs.

Parameters

[in]	n_cells	local number of cells
[in]	n_faces	local number of faces
[in,out]	msles	pointer to the structure to reset

cs_cdofb_monolithic_sles_reset()

void cs_cdofb_monolithic_sles_reset

(

cs_cdofb_monolithic_sles_t *

msles

)

Reset to zero rhs and clean the cs_sles_t structure.

Parameters

[in,out]

msles

pointer to the structure to reset

cs_cdofb_monolithic_sles_set_shared()

void cs_cdofb_monolithic_sles_set_shared	(	const cs_cdo_connect_t *	connect,
		const cs_cdo_quantities_t *	quant,
		const cs_range_set_t *	rset
	)

Set pointers to shared structures.

Parameters

[in]	connect	pointer to cdo connectivities
[in]	quant	pointer to additional mesh quantities
[in]	rset	pointer to a cs_range_set_t structure

cs_cdofb_monolithic_solve()

int cs_cdofb_monolithic_solve	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Solve a linear system arising from the discretization of the Navier-Stokes equation with a CDO face-based approach. The full system is treated as one block and solved as it is. In this situation, PETSc or MUMPS are usually considered.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the (cumulated) number of iterations of the solver

cs_cdofb_monolithic_uzawa_al_incr_solve()

int cs_cdofb_monolithic_uzawa_al_incr_solve	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Use the Uzawa algorithm with an Augmented Lagrangian technique in an incremental way to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the cumulated number of iterations of the solver

cs_cdofb_monolithic_uzawa_al_solve()

int cs_cdofb_monolithic_uzawa_al_solve	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Use the Uzawa algorithm with an Augmented Lagrangian technique to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the cumulated number of iterations of the solver

cs_cdofb_monolithic_uzawa_cg_solve()

int cs_cdofb_monolithic_uzawa_cg_solve	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Use the preconditioned Uzawa-CG algorithm to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling This algorithm is based on Koko's paper "Uzawa conjugate gradient method for the Stokes problem: Matlab implementation with P1-iso-P2/ P1 finite element".

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the cumulated number of iterations of the solver

cs_cdofb_monolithic_uzawa_n3s_solve()

int cs_cdofb_monolithic_uzawa_n3s_solve	(	const cs_navsto_param_t *	nsp,
		const cs_equation_param_t *	eqp,
		cs_cdofb_monolithic_sles_t *	msles
	)

Use the preconditioned Uzawa-CG algorithm to solve the saddle-point problem arising from CDO-Fb schemes for Stokes, Oseen and Navier-Stokes with a monolithic coupling This algorithm is based on the EDF report H-E40-1991-03299-FR devoted the numerical algorithms used in the code N3S. Specifically a Cahout-Chabard preconditioning is used to approximate the Schur complement.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	msles	pointer to a cs_cdofb_monolithic_sles_t structure

Returns

the cumulated number of iterations of the solver