itersolve.h

浏览该文件的文档.

#ifndef LMMC_ITERSOLVE_H
#define LMMC_ITERSOLVE_H
 
#include <stddef.h>
#include "lmmc/numeric.h"
#include "lmmc/diagnostic.h"
#include "lmmc/precond.h"
#include "lmmc/sparse.h"
 
#ifdef __cplusplus
extern "C" {
#endif

typedef struct {
    int converged;                          
    size_t num_iter;                        
    lmmc_real_t initial_residual_norm;      
    lmmc_real_t final_residual_norm;        
} lmmc_itersolve_result_t;

typedef lmmc_status_t (*lmmc_matvec_op_t)(
    const lmmc_vec_t* x, lmmc_vec_t* y, void* user_data);

typedef struct {
    lmmc_real_t abs_tol;                              
    lmmc_real_t rel_tol;                              
    size_t max_iter;                                  
    size_t restart;                                   
    lmmc_diagnostic_sink_t diagnostics;               
    lmmc_matvec_op_t apply_op;                        
    void* op_user_data;                               
    lmmc_matvec_op_t apply_transpose_op;              
} lmmc_itersolve_config_t;

lmmc_status_t lmmc_itersolve_default_config(size_t problem_size, lmmc_itersolve_config_t* out_cfg);

lmmc_status_t lmmc_cg_solve(
    const lmmc_sparse_mat_t* a,
    const lmmc_vec_t* b,
    const lmmc_precond_t* precond,
    const lmmc_itersolve_config_t* cfg,
    lmmc_vec_t* x,
    lmmc_itersolve_result_t* out_result
);

lmmc_status_t lmmc_bicgstab_solve(
    const lmmc_sparse_mat_t* a,
    const lmmc_vec_t* b,
    const lmmc_precond_t* precond,
    const lmmc_itersolve_config_t* cfg,
    lmmc_vec_t* x,
    lmmc_itersolve_result_t* out_result
);

lmmc_status_t lmmc_gmres_solve(
    const lmmc_sparse_mat_t* a,
    const lmmc_vec_t* b,
    const lmmc_precond_t* precond,
    const lmmc_itersolve_config_t* cfg,
    lmmc_vec_t* x,
    lmmc_itersolve_result_t* out_result
);

lmmc_status_t lmmc_minres_solve(
    const lmmc_sparse_mat_t* a,
    const lmmc_vec_t* b,
    const lmmc_precond_t* precond,
    const lmmc_itersolve_config_t* cfg,
    lmmc_vec_t* x,
    lmmc_itersolve_result_t* out_result
);

lmmc_status_t lmmc_lsqr_solve(
    const lmmc_sparse_mat_t* a,
    const lmmc_vec_t* b,
    const lmmc_itersolve_config_t* cfg,
    lmmc_vec_t* x,
    lmmc_itersolve_result_t* out_result
);
 
#ifdef __cplusplus
}
#endif
 
#endif