prec.h

Go to the documentation of this file.

#pragma once
 
#ifdef __cplusplus
extern "C" {
#endif
 
#include "bsr.h"
#include <stdbool.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_ilu0_factorize2(prec_t *P, bsr_matrix *A, bool use_dilu);
void prec_ilu0_factorize3(prec_t *P, bsr_matrix *A, bool use_dilu);
void prec_ilu0_factorize4(prec_t *P, bsr_matrix *A, bool use_dilu);
 
void prec_mapply2c(prec_t *P, double *x);
void prec_mapply3c(prec_t *P, double *x);
void prec_mapply4c(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