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;
}
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_pbicgstab3m(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