1 #ifndef __CS_SADDLE_SOLVER_H__
2 #define __CS_SADDLE_SOLVER_H__
cs_real_t() cs_cdo_blas_square_norm_t(const cs_real_t *array)
Generic function pointer for computing a square norm. Parallel synchronization is performed.
Definition: cs_cdo_blas.h:79
#define BEGIN_C_DECLS
Definition: cs_defs.h:528
double cs_real_t
Floating-point value.
Definition: cs_defs.h:332
#define END_C_DECLS
Definition: cs_defs.h:529
int cs_lnum_t
local mesh entity id
Definition: cs_defs.h:325
@ x2
Definition: cs_field_pointer.h:224
struct _cs_matrix_t cs_matrix_t
Definition: cs_matrix.h:110
Handle the settings of saddle-point systems. These systems arise from the monolithic coupling of the ...
void cs_saddle_solver_gcr(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply the GCR algorithm to a saddle point problem (the system is stored in a hybrid way)....
Definition: cs_saddle_solver.c:3827
cs_saddle_solver_t * cs_saddle_solver_add(cs_lnum_t n1_elts, int n1_dofs_by_elt, cs_lnum_t n2_elts, int n2_dofs_by_elt, const cs_param_saddle_t *saddlep, cs_cdo_system_helper_t *sh, cs_sles_t *main_sles)
Add a new solver for solving a saddle-point problem.
Definition: cs_saddle_solver.c:2083
void cs_saddle_solver_m21_multiply_scalar(cs_lnum_t n2_dofs, const cs_real_t *x1, const cs_adjacency_t *m21_adj, const cs_real_t *m21_val, cs_real_t *m21x1)
Compute the resulting vector of the operation m21*x1 The stride is equal to 1 for the operator m21 op...
Definition: cs_saddle_solver.c:2512
void cs_saddle_solver_context_alu_clean(cs_saddle_solver_context_alu_t *ctx)
Free main memory consuming part of the context structure associated to an ALU algorithm.
Definition: cs_saddle_solver.c:2604
void cs_saddle_solver_context_uzawa_cg_create(cs_lnum_t b22_max_size, cs_saddle_solver_t *solver)
Create and initialize the context structure for an algorithm related to the Uzawa-CG algorithm.
Definition: cs_saddle_solver.c:2986
void cs_saddle_solver_uzawa_cg(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply the Uzawa-CG algorithm to solve a saddle point problem (the system is stored in a hybrid way)....
Definition: cs_saddle_solver.c:4359
void cs_saddle_solver_context_notay_create(cs_saddle_solver_t *solver)
Create and initialize the context structure for the algorithm relying on the Notay's algebraic transf...
Definition: cs_saddle_solver.c:2947
void cs_saddle_solver_free(cs_saddle_solver_t **p_solver)
Free a cs_saddle_solver_t structure.
Definition: cs_saddle_solver.c:2134
void cs_saddle_solver_context_block_pcd_create(cs_lnum_t b22_max_size, cs_saddle_solver_t *solver)
Create and initialize the context structure for a block preconditioner used in combination with a Kry...
Definition: cs_saddle_solver.c:2656
void cs_saddle_solver_minres(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply the MINRES algorithm to a saddle point problem (the system is stored in a hybrid way)....
Definition: cs_saddle_solver.c:3571
cs_real_t * cs_saddle_solver_m11_inv_lumped(cs_saddle_solver_t *solver, const cs_matrix_t *m11, const cs_range_set_t *b11_rset, cs_sles_t *xtra_sles, int *n_iter)
Retrieve the lumped matrix the inverse of the diagonal of the (1,1)-block matrix. The storage of a ma...
Definition: cs_saddle_solver.c:2340
void cs_saddle_solver_gkb_inhouse(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply the GKB algorithm to a saddle point problem (the system is stored in a hybrid way)....
Definition: cs_saddle_solver.c:4052
void cs_saddle_solver_update_monitoring(cs_saddle_solver_t *solver, unsigned n_iter)
Update the current monitoring state with n_iter.
Definition: cs_saddle_solver.c:2286
void cs_saddle_solver_alu_incr(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply the Augmented Lagrangian-Uzawa algorithm to a saddle point problem (the system is stored in a h...
Definition: cs_saddle_solver.c:3172
void cs_saddle_solver_context_gkb_free(cs_saddle_solver_context_gkb_t **p_ctx)
Free the context structure associated to the GKB algorithm.
Definition: cs_saddle_solver.c:2922
void cs_saddle_solver_finalize(void)
Free all remaining structures related to saddle-point solvers.
Definition: cs_saddle_solver.c:2119
void cs_saddle_solver_context_alu_free(cs_saddle_solver_context_alu_t **p_ctx)
Free the context structure associated to an ALU algorithm.
Definition: cs_saddle_solver.c:2630
void cs_saddle_solver_notay(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply the Notay's transformation algorithm to solve a saddle point problem (the system is stored in a...
Definition: cs_saddle_solver.c:3398
void cs_saddle_solver_log_monitoring(void)
Log the monitoring performance of all saddle-point systems.
Definition: cs_saddle_solver.c:2309
void cs_saddle_solver_context_gkb_create(cs_saddle_solver_t *solver)
Create and initialize the context structure for the GKB algorithm.
Definition: cs_saddle_solver.c:2811
void cs_saddle_solver_sles_full_system(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply an (external) solver to solve a saddle point problem (the system is stored in a monolithic way)...
Definition: cs_saddle_solver.c:4212
void cs_saddle_solver_context_alu_create(cs_saddle_solver_t *solver)
Create and initialize the context structure for an algorithm related to the ALU algorithm.
Definition: cs_saddle_solver.c:2542
void cs_saddle_solver_context_uzawa_cg_clean(cs_saddle_solver_context_uzawa_cg_t *ctx)
Free main memory consuming part of the context structure associated to a Uzawa-CG algorithm.
Definition: cs_saddle_solver.c:3107
void cs_saddle_solver_context_gkb_clean(cs_saddle_solver_context_gkb_t *ctx)
Free the main memory consuming part of the context structure associated to the GKB algorithm.
Definition: cs_saddle_solver.c:2884
cs_saddle_solver_t * cs_saddle_solver_by_id(int id)
Get a pointer to saddle-point solver from its id.
Definition: cs_saddle_solver.c:2056
int cs_saddle_solver_get_n_systems(void)
Retrieve the number of saddle-point systems which have been added.
Definition: cs_saddle_solver.c:2040
void cs_saddle_solver_context_block_pcd_clean(cs_saddle_solver_context_block_pcd_t *ctx)
Free the context structure for a block preconditioner used in combination with a Krylov solver such a...
Definition: cs_saddle_solver.c:2757
void cs_saddle_solver_m12_multiply_scalar(cs_lnum_t n2_elts, const cs_real_t *x2, const cs_adjacency_t *m21_adj, const cs_real_t *m21_val, cs_real_t *m12x2)
Compute the resulting vector of the operation m12*x2 The stride is equal to 1 for the operator m21 (u...
Definition: cs_saddle_solver.c:2443
void cs_saddle_solver_m12_multiply_vector(cs_lnum_t n2_dofs, const cs_real_t *x2, const cs_adjacency_t *m21_adj, const cs_real_t *m21_val, cs_real_t *m12x2)
Compute the resulting vector of the operation m12*x2 The stride is equal to 3 for the operator m21 (u...
Definition: cs_saddle_solver.c:2396
void cs_saddle_solver_m21_multiply_vector(cs_lnum_t n2_dofs, const cs_real_t *x1, const cs_adjacency_t *m21_adj, const cs_real_t *m21_val, cs_real_t *m21x1)
Compute the resulting vector of the operation m21*x1 The stride is equal to 3 for the operator m21 op...
Definition: cs_saddle_solver.c:2477
void cs_saddle_solver_clean(cs_saddle_solver_t *solver)
Free/reset only a part of a cs_saddle_solver_t structure.
Definition: cs_saddle_solver.c:2214
void cs_saddle_solver_context_uzawa_cg_free(cs_saddle_solver_context_uzawa_cg_t **p_ctx)
Free the context structure associated to a Uzawa-CG algorithm.
Definition: cs_saddle_solver.c:3144
void cs_saddle_solver_context_block_pcd_free(cs_saddle_solver_context_block_pcd_t **p_ctx)
Free the context structure for a block preconditioner used in combination with a Krylov solver such a...
Definition: cs_saddle_solver.c:2787
void() cs_saddle_solver_matvec_t(cs_lnum_t n2_dofs, const cs_real_t *vec, const cs_adjacency_t *mat_adj, const cs_real_t *mat_op, cs_real_t *matvec)
Generic function prototype to perform a matrix vector operation This operation takes place between an...
Definition: cs_saddle_solver.h:72
struct _cs_sles_t cs_sles_t
Definition: cs_sles.h:68
Definition: cs_mesh_adjacencies.h:68
Definition: cs_cdo_system.h:377
Structure to handle the convergence of an iterative algorithm.
Definition: cs_iter_algo.h:289
Structure storing all metadata related to the resolution of a saddle-point linear system....
Definition: cs_param_saddle.h:272
Structure associated to the definition of a property relying on the cs_xdef_t structure.
Definition: cs_range_set.h:57
Definition: cs_saddle_solver.h:142
cs_sles_t * init_sles
Definition: cs_saddle_solver.h:156
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition: cs_saddle_solver.h:162
cs_real_t * rhs
Definition: cs_saddle_solver.h:152
cs_real_t * b1_tilda
Definition: cs_saddle_solver.h:151
cs_real_t * inv_m22
Definition: cs_saddle_solver.h:146
const cs_adjacency_t * m21_adj
Definition: cs_saddle_solver.h:171
const cs_real_t * m21_val
Definition: cs_saddle_solver.h:167
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition: cs_saddle_solver.h:163
cs_real_t * m21x1
Definition: cs_saddle_solver.h:149
cs_real_t * res2
Definition: cs_saddle_solver.h:148
cs_cdo_blas_square_norm_t * square_norm_b11
Definition: cs_saddle_solver.h:161
Definition: cs_saddle_solver.h:189
cs_real_t * schur_diag
Definition: cs_saddle_solver.h:211
cs_range_set_t * b11_range_set
Definition: cs_saddle_solver.h:194
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition: cs_saddle_solver.h:223
cs_real_t * m11_inv_diag
Definition: cs_saddle_solver.h:216
cs_matrix_t * m11
Definition: cs_saddle_solver.h:193
cs_sles_t * xtra_sles
Definition: cs_saddle_solver.h:219
const cs_property_t * pty_22
Definition: cs_saddle_solver.h:228
cs_matrix_t * schur_matrix
Definition: cs_saddle_solver.h:205
const cs_adjacency_t * m21_adj
Definition: cs_saddle_solver.h:230
const cs_real_t * m21_val
Definition: cs_saddle_solver.h:229
double schur_scaling
Definition: cs_saddle_solver.h:207
cs_sles_t * schur_sles
Definition: cs_saddle_solver.h:206
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition: cs_saddle_solver.h:224
cs_lnum_t b11_max_size
Definition: cs_saddle_solver.h:199
cs_real_t * m22_mass_diag
Definition: cs_saddle_solver.h:217
cs_real_t * schur_xtra
Definition: cs_saddle_solver.h:212
cs_lnum_t b22_max_size
Definition: cs_saddle_solver.h:200
Definition: cs_saddle_solver.h:243
cs_real_t zeta_square_sum
Definition: cs_saddle_solver.h:255
cs_sles_t * init_sles
Definition: cs_saddle_solver.h:275
cs_real_t * rhs_tilda
Definition: cs_saddle_solver.h:271
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition: cs_saddle_solver.h:280
cs_real_t alpha
Definition: cs_saddle_solver.h:247
cs_real_t * m21v
Definition: cs_saddle_solver.h:261
cs_real_t zeta
Definition: cs_saddle_solver.h:249
cs_real_t * m12q
Definition: cs_saddle_solver.h:267
cs_real_t * inv_m22
Definition: cs_saddle_solver.h:262
cs_real_t * zeta_array
Definition: cs_saddle_solver.h:254
const cs_adjacency_t * m21_adj
Definition: cs_saddle_solver.h:289
const cs_real_t * m22
Definition: cs_saddle_solver.h:264
const cs_real_t * m21_val
Definition: cs_saddle_solver.h:285
cs_real_t beta
Definition: cs_saddle_solver.h:248
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition: cs_saddle_solver.h:281
cs_real_t * q
Definition: cs_saddle_solver.h:259
int zeta_size
Definition: cs_saddle_solver.h:253
cs_real_t * x1_tilda
Definition: cs_saddle_solver.h:269
cs_real_t * d
Definition: cs_saddle_solver.h:260
cs_cdo_blas_square_norm_t * square_norm_b11
Definition: cs_saddle_solver.h:279
cs_real_t * w
Definition: cs_saddle_solver.h:266
cs_real_t * v
Definition: cs_saddle_solver.h:268
Definition: cs_saddle_solver.h:297
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition: cs_saddle_solver.h:301
const cs_adjacency_t * m21_adj
Definition: cs_saddle_solver.h:309
const cs_real_t * m21_val
Definition: cs_saddle_solver.h:305
Definition: cs_saddle_solver.h:317
cs_sles_t * init_sles
Definition: cs_saddle_solver.h:368
cs_real_t * schur_diag
Definition: cs_saddle_solver.h:360
cs_range_set_t * b11_range_set
Definition: cs_saddle_solver.h:342
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition: cs_saddle_solver.h:336
cs_real_t * gk
Definition: cs_saddle_solver.h:327
cs_real_t * m11_inv_diag
Definition: cs_saddle_solver.h:365
cs_real_t * rhs
Definition: cs_saddle_solver.h:331
cs_matrix_t * m11
Definition: cs_saddle_solver.h:341
cs_real_t * b1_tilda
Definition: cs_saddle_solver.h:329
cs_sles_t * xtra_sles
Definition: cs_saddle_solver.h:367
cs_real_t alpha
Definition: cs_saddle_solver.h:321
const cs_property_t * pty_22
Definition: cs_saddle_solver.h:372
cs_matrix_t * schur_matrix
Definition: cs_saddle_solver.h:355
cs_real_t * dzk
Definition: cs_saddle_solver.h:330
cs_real_t * inv_m22
Definition: cs_saddle_solver.h:353
const cs_adjacency_t * m21_adj
Definition: cs_saddle_solver.h:374
const cs_real_t * m21_val
Definition: cs_saddle_solver.h:373
cs_sles_t * schur_sles
Definition: cs_saddle_solver.h:356
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition: cs_saddle_solver.h:337
cs_real_t * m21x1
Definition: cs_saddle_solver.h:326
cs_real_t * res2
Definition: cs_saddle_solver.h:325
cs_lnum_t b11_max_size
Definition: cs_saddle_solver.h:347
cs_cdo_blas_square_norm_t * square_norm_b11
Definition: cs_saddle_solver.h:335
cs_real_t * schur_xtra
Definition: cs_saddle_solver.h:361
cs_lnum_t b22_max_size
Definition: cs_saddle_solver.h:348
Definition: cs_saddle_solver.h:88
cs_cdo_system_helper_t * system_helper
Definition: cs_saddle_solver.h:100
cs_sles_t * main_sles
Definition: cs_saddle_solver.h:108
cs_lnum_t n1_scatter_dofs
Definition: cs_saddle_solver.h:117
bool do_setup
Definition: cs_saddle_solver.h:102
cs_lnum_t n2_elts
Definition: cs_saddle_solver.h:114
unsigned n_calls
Definition: cs_saddle_solver.h:127
cs_lnum_t n1_elts
Definition: cs_saddle_solver.h:112
cs_iter_algo_t * algo
Definition: cs_saddle_solver.h:96
cs_lnum_t n2_scatter_dofs
Definition: cs_saddle_solver.h:118
unsigned n_iter_max
Definition: cs_saddle_solver.h:129
int n2_dofs_by_elt
Definition: cs_saddle_solver.h:115
unsigned n_iter_tot
Definition: cs_saddle_solver.h:130
unsigned n_iter_min
Definition: cs_saddle_solver.h:128
const cs_param_saddle_t * param
Definition: cs_saddle_solver.h:92
int n1_dofs_by_elt
Definition: cs_saddle_solver.h:113
void * context
Definition: cs_saddle_solver.h:123