krylov.h

Go to the documentation of this file.

 
#ifndef KRYLOV_H
#define KRYLOV_H
#include "mkl.h"
#include "ucfd_types.h"
#include "config.h"
 
#ifndef EPS
    #define eps 2.22e-16 // Machine epsilon of double precision
#endif
 
 
ucfd_status_t serial_gmres(sparse_matrix_t op, ucfd_precon_type_t precon_type, int bn, int block, int m, int *iter, double tol,
                           int *row_ptr, int *col_ind, int *diag_ind, double *precon_nnz_data,
                           double *x, double *b, double *H, double *V, double *g, double *y, double *w, double *r);
 
 
ucfd_status_t step_gmres(sparse_matrix_t op, ucfd_precon_solve psolve, const struct matrix_descr descr,
                         int bn, int m, int *flag,
                         int *row_ptr, int *col_ind, int *diag_ind, double *precon_nnz_data,
                         double *x, double *b, double *H, double *V, double *g, double *y, double *w, double *r);
 
 
ucfd_status_t serial_bicgstab(sparse_matrix_t op, ucfd_precon_type_t precon_type, int bn, int *iter, double tol,
                              int *row_ptr, int *col_ind, int *diag_ind, double *precon_nnz_data,
                              double *x, double *b, double *r, double *p, double *v, double *s, double *t);
 
#endif // KRYLOV_H