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 inside th...
Definition cs_cdo_blas.h:79
#define BEGIN_C_DECLS
Definition cs_defs.h:542
double cs_real_t
Floating-point value.
Definition cs_defs.h:342
#define END_C_DECLS
Definition cs_defs.h:543
int cs_lnum_t
local mesh entity id
Definition cs_defs.h:335
@ x2
Definition cs_field_pointer.h:202
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.cpp:4125
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.cpp:2619
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.cpp:2712
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.cpp:3095
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.cpp:4663
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.cpp:3055
void cs_saddle_solver_free(cs_saddle_solver_t **p_solver)
Free a cs_saddle_solver_t structure.
Definition cs_saddle_solver.cpp:2208
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.cpp:2762
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.cpp:3868
void cs_saddle_solver_context_simple_free(cs_saddle_solver_context_simple_t **p_ctx)
Free the context structure associated to a Uzawa-CG algorithm.
Definition cs_saddle_solver.cpp:3427
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.cpp:2130
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.cpp:4352
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.cpp:2389
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.cpp:3455
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.cpp:3030
void cs_saddle_solver_finalize(void)
Free all remaining structures related to saddle-point solvers.
Definition cs_saddle_solver.cpp:2193
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.cpp:2736
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.cpp:3683
void cs_saddle_solver_context_simple_create(cs_lnum_t b22_max_size, cs_saddle_solver_t *solver)
Create and initialize the context structure for an algorithm related to the SIMPLE algorithm.
Definition cs_saddle_solver.cpp:3281
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.cpp:2443
void cs_saddle_solver_log_monitoring(void)
Log the monitoring performance of all saddle-point systems.
Definition cs_saddle_solver.cpp:2412
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.cpp:2917
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:73
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.cpp:4514
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.cpp:2649
void cs_saddle_solver_simple(cs_saddle_solver_t *solver, cs_real_t *x1, cs_real_t *x2)
Apply the SIMPLE-like algorithm to solve a saddle point problem.
Definition cs_saddle_solver.cpp:4908
void cs_saddle_solver_context_simple_clean(cs_saddle_solver_context_simple_t *ctx)
Free main memory consuming part of the context structure associated to a SIMPLE algorithm.
Definition cs_saddle_solver.cpp:3395
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.cpp:3219
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.cpp:2992
int cs_saddle_solver_get_n_systems(void)
Retrieve the number of saddle-point systems which have been added.
Definition cs_saddle_solver.cpp:2114
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.cpp:2864
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.cpp:2550
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.cpp:2503
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.cpp:2584
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.cpp:2305
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.cpp:3255
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.cpp:2893
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.cpp:2157
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:290
Structure storing all metadata related to the resolution of a saddle-point linear system....
Definition cs_param_saddle.h:273
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:143
cs_sles_t * init_sles
Definition cs_saddle_solver.h:157
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition cs_saddle_solver.h:163
cs_real_t * rhs
Definition cs_saddle_solver.h:153
cs_real_t * b1_tilda
Definition cs_saddle_solver.h:152
cs_real_t * inv_m22
Definition cs_saddle_solver.h:147
const cs_adjacency_t * m21_adj
Definition cs_saddle_solver.h:172
const cs_real_t * m21_val
Definition cs_saddle_solver.h:168
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition cs_saddle_solver.h:164
cs_real_t * m21x1
Definition cs_saddle_solver.h:150
cs_real_t * res2
Definition cs_saddle_solver.h:149
cs_cdo_blas_square_norm_t * square_norm_b11
Definition cs_saddle_solver.h:162
Definition cs_saddle_solver.h:190
cs_real_t * schur_diag
Definition cs_saddle_solver.h:212
cs_range_set_t * b11_range_set
Definition cs_saddle_solver.h:195
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition cs_saddle_solver.h:224
cs_real_t * m11_inv_diag
Definition cs_saddle_solver.h:217
cs_matrix_t * m11
Definition cs_saddle_solver.h:194
cs_sles_t * xtra_sles
Definition cs_saddle_solver.h:220
const cs_property_t * pty_22
Definition cs_saddle_solver.h:229
cs_matrix_t * schur_matrix
Definition cs_saddle_solver.h:206
const cs_adjacency_t * m21_adj
Definition cs_saddle_solver.h:231
const cs_real_t * m21_val
Definition cs_saddle_solver.h:230
double schur_scaling
Definition cs_saddle_solver.h:208
cs_sles_t * schur_sles
Definition cs_saddle_solver.h:207
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition cs_saddle_solver.h:225
cs_lnum_t b11_max_size
Definition cs_saddle_solver.h:200
cs_real_t * m22_mass_diag
Definition cs_saddle_solver.h:218
cs_real_t * schur_xtra
Definition cs_saddle_solver.h:213
cs_lnum_t b22_max_size
Definition cs_saddle_solver.h:201
Definition cs_saddle_solver.h:244
cs_real_t zeta_square_sum
Definition cs_saddle_solver.h:256
cs_sles_t * init_sles
Definition cs_saddle_solver.h:276
cs_real_t * rhs_tilda
Definition cs_saddle_solver.h:272
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition cs_saddle_solver.h:281
cs_real_t alpha
Definition cs_saddle_solver.h:248
cs_real_t * m21v
Definition cs_saddle_solver.h:262
cs_real_t zeta
Definition cs_saddle_solver.h:250
cs_real_t * m12q
Definition cs_saddle_solver.h:268
cs_real_t * inv_m22
Definition cs_saddle_solver.h:263
cs_real_t * zeta_array
Definition cs_saddle_solver.h:255
const cs_adjacency_t * m21_adj
Definition cs_saddle_solver.h:290
const cs_real_t * m22
Definition cs_saddle_solver.h:265
const cs_real_t * m21_val
Definition cs_saddle_solver.h:286
cs_real_t beta
Definition cs_saddle_solver.h:249
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition cs_saddle_solver.h:282
cs_real_t * q
Definition cs_saddle_solver.h:260
int zeta_size
Definition cs_saddle_solver.h:254
cs_real_t * x1_tilda
Definition cs_saddle_solver.h:270
cs_real_t * d
Definition cs_saddle_solver.h:261
cs_cdo_blas_square_norm_t * square_norm_b11
Definition cs_saddle_solver.h:280
cs_real_t * w
Definition cs_saddle_solver.h:267
cs_real_t * v
Definition cs_saddle_solver.h:269
Definition cs_saddle_solver.h:298
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition cs_saddle_solver.h:302
const cs_adjacency_t * m21_adj
Definition cs_saddle_solver.h:310
const cs_real_t * m21_val
Definition cs_saddle_solver.h:306
Definition cs_saddle_solver.h:383
cs_sles_t * init_sles
Definition cs_saddle_solver.h:429
cs_real_t * schur_diag
Definition cs_saddle_solver.h:421
cs_range_set_t * b11_range_set
Definition cs_saddle_solver.h:403
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition cs_saddle_solver.h:397
cs_real_t * m11_inv_diag
Definition cs_saddle_solver.h:426
cs_real_t * rhs
Definition cs_saddle_solver.h:392
cs_matrix_t * m11
Definition cs_saddle_solver.h:402
cs_real_t * b1_tilda
Definition cs_saddle_solver.h:391
cs_sles_t * xtra_sles
Definition cs_saddle_solver.h:428
const cs_property_t * pty_22
Definition cs_saddle_solver.h:433
cs_matrix_t * schur_matrix
Definition cs_saddle_solver.h:416
cs_real_t * inv_m22
Definition cs_saddle_solver.h:414
const cs_adjacency_t * m21_adj
Definition cs_saddle_solver.h:435
const cs_real_t * m21_val
Definition cs_saddle_solver.h:434
cs_sles_t * schur_sles
Definition cs_saddle_solver.h:417
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition cs_saddle_solver.h:398
cs_real_t * m21x1
Definition cs_saddle_solver.h:389
cs_real_t * res2
Definition cs_saddle_solver.h:388
cs_lnum_t b11_max_size
Definition cs_saddle_solver.h:408
cs_cdo_blas_square_norm_t * square_norm_b11
Definition cs_saddle_solver.h:396
cs_real_t * schur_xtra
Definition cs_saddle_solver.h:422
cs_lnum_t b22_max_size
Definition cs_saddle_solver.h:409
Definition cs_saddle_solver.h:318
cs_sles_t * init_sles
Definition cs_saddle_solver.h:369
cs_real_t * schur_diag
Definition cs_saddle_solver.h:361
cs_range_set_t * b11_range_set
Definition cs_saddle_solver.h:343
cs_saddle_solver_matvec_t * m12_vector_multiply
Definition cs_saddle_solver.h:337
cs_real_t * gk
Definition cs_saddle_solver.h:328
cs_real_t * m11_inv_diag
Definition cs_saddle_solver.h:366
cs_real_t * rhs
Definition cs_saddle_solver.h:332
cs_matrix_t * m11
Definition cs_saddle_solver.h:342
cs_real_t * b1_tilda
Definition cs_saddle_solver.h:330
cs_sles_t * xtra_sles
Definition cs_saddle_solver.h:368
cs_real_t alpha
Definition cs_saddle_solver.h:322
const cs_property_t * pty_22
Definition cs_saddle_solver.h:373
cs_matrix_t * schur_matrix
Definition cs_saddle_solver.h:356
cs_real_t * dzk
Definition cs_saddle_solver.h:331
cs_real_t * inv_m22
Definition cs_saddle_solver.h:354
const cs_adjacency_t * m21_adj
Definition cs_saddle_solver.h:375
const cs_real_t * m21_val
Definition cs_saddle_solver.h:374
cs_sles_t * schur_sles
Definition cs_saddle_solver.h:357
cs_saddle_solver_matvec_t * m21_vector_multiply
Definition cs_saddle_solver.h:338
cs_real_t * m21x1
Definition cs_saddle_solver.h:327
cs_real_t * res2
Definition cs_saddle_solver.h:326
cs_lnum_t b11_max_size
Definition cs_saddle_solver.h:348
cs_cdo_blas_square_norm_t * square_norm_b11
Definition cs_saddle_solver.h:336
cs_real_t * schur_xtra
Definition cs_saddle_solver.h:362
cs_lnum_t b22_max_size
Definition cs_saddle_solver.h:349
Definition cs_saddle_solver.h:89
cs_cdo_system_helper_t * system_helper
Definition cs_saddle_solver.h:101
cs_sles_t * main_sles
Definition cs_saddle_solver.h:109
cs_lnum_t n1_scatter_dofs
Definition cs_saddle_solver.h:118
bool do_setup
Definition cs_saddle_solver.h:103
cs_lnum_t n2_elts
Definition cs_saddle_solver.h:115
unsigned n_calls
Definition cs_saddle_solver.h:128
cs_lnum_t n1_elts
Definition cs_saddle_solver.h:113
cs_iter_algo_t * algo
Definition cs_saddle_solver.h:97
cs_lnum_t n2_scatter_dofs
Definition cs_saddle_solver.h:119
unsigned n_iter_max
Definition cs_saddle_solver.h:130
int n2_dofs_by_elt
Definition cs_saddle_solver.h:116
unsigned n_iter_tot
Definition cs_saddle_solver.h:131
unsigned n_iter_min
Definition cs_saddle_solver.h:129
const cs_param_saddle_t * param
Definition cs_saddle_solver.h:93
int n1_dofs_by_elt
Definition cs_saddle_solver.h:114
void * context
Definition cs_saddle_solver.h:124