bslv.h

Go to the documentation of this file.

#pragma once
 
#ifdef __cplusplus
extern "C" {
#endif
 
 
#include "bsr.h"
#include "prec.h"
 
#include <stdbool.h>
 
typedef
struct bslv_memory
{
    // perform dilu (ilu0) of true (False)
    bool use_dilu;
 
    // tolerance
    double tol;
    // maximum number of iterations
    int    max_iter;
    // convergencey history
    double *e;
 
    // size of linear system
    int n;
    // pointer to temporary arrays
    double **dtmp;
 
    // pointer to preconditioner
    prec_t *P;
 
    // spmv function pointer
    void (*bsr_spmv)(bsr_matrix *A, const double *x, double *y);
 
    // ILU0/DILU factorization function pointer
    void (*prec_factorize)(prec_t *P, bsr_matrix *A, bool use_dilu);
 
    // ILU0/DILU apply function pointer
    void (*prec_apply)(prec_t *P, double *x);
}
bslv_memory;
 
bslv_memory *bslv_alloc();
 
void bslv_free(bslv_memory *mem);
 
void bslv_info(bslv_memory *mem, int count);
 
void bslv_init(bslv_memory *mem, double tol, int max_iter, bsr_matrix const *A, bool use_dilu);
 
int  bslv_pbicgstabm(bslv_memory *mem, bsr_matrix *A, const double *b, double *x);
 
int  bslv_pbicgstab3d(bslv_memory *mem, bsr_matrix *A, const double *b, double *x);
 
#ifdef __cplusplus
}
#endif