prec.h

#pragma once

#ifdef __cplusplus
extern "C" {
#endif

#include "bsr.h"

typedef
struct prec_t
{
    // lower triangular factor
    bsr_matrix *L;
    // diagonal factor
    bsr_matrix *D;
    // upper triangular factor
    bsr_matrix *U;
    // number of off-diagonal ilu0 updates
    int noffsets;
    // triplets uniquely identifying off-diagonal ilu0 updates
    int(*offsets)[3];
}
prec_t;

prec_t *prec_alloc();

void prec_free(prec_t *P);

void prec_init(prec_t *P, bsr_matrix const *A);

int  prec_analyze(bsr_matrix *M, int (*offsets)[3]);

void prec_dilu_factorize(prec_t *P, bsr_matrix *A);

void prec_ilu0_factorize(prec_t *P, bsr_matrix *A);

void prec_mapply3c(prec_t *P, double *x);

void prec_dapply3c(prec_t *P, double *x);

void prec_downcast(prec_t *P);

void prec_info(prec_t *P);

#ifdef __cplusplus
}
#endif