bslv.h

#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