EVSL
1.1.0
EigenValues Slicing Library
|
Definitions used for direct solver interface. More...
Go to the source code of this file.
Functions | |
int | SetupBSolDirect (csrMat *B, void **data) |
Setup the B-sol by computing the Cholesky factorization of B. More... | |
void | BSolDirect (double *b, double *x, void *data) |
Solver function of B. More... | |
void | LTSolDirect (double *b, double *x, void *data) |
Solver function of L^{T} x = L^{-T}*b. More... | |
void | FreeBSolDirectData (void *data) |
Free solver data. More... | |
int | SetupASIGMABSolDirect (csrMat *A, csrMat *BB, int num, complex double *zk, void **data) |
setup CXsparse solver for A - SIGMA B More... | |
void | ASIGMABSolDirect (int n, double *br, double *bi, double *xr, double *xz, void *data) |
complex linear solver routine passed to evsl More... | |
void | FreeASIGMABSolDirect (int num, void **data) |
free the data needed by CXSparse More... | |
Definitions used for direct solver interface.
Note that this file is meant to be the header file of both evsl_suitesparse.c, evsl_pardiso.c, and evsl_cxsparse.c. If more direct solver options will be added later, this should serve them as well.
Definition in file evsl_direct.h.
void ASIGMABSolDirect | ( | int | n, |
double * | br, | ||
double * | bi, | ||
double * | xr, | ||
double * | xz, | ||
void * | data | ||
) |
complex linear solver routine passed to evsl
n | size of the system | |
br,bi | vectors of length n, complex right-hand side (real and imaginary) | |
data | all data that are needed for solving the system | |
[out] | xr,xz | vectors of length n, complex solution (real and imaginary) |
complex linear solver routine passed to evsl
n | size of the system | |
br,bi | vectors of length n, complex right-hand side (real and imaginary) | |
data | all data that are needed for solving the system | |
[out] | xr,xz | vectors of length n, complex solution (real and imaginary) |
Definition at line 324 of file evsl_cxsparse.c.
References _ASBSolDataDirect::b, CHKERR, cs_ci_ipvec(), cs_ci_lsolve(), cs_ci_usolve(), cs_complex_t, cs_ci_numeric::L, _ASBSolDataDirect::n, _ASBSolDataDirect::N, cs_ci_numeric::pinv, cs_ci_symbolic::q, _ASBSolDataDirect::S, cs_ci_numeric::U, and _ASBSolDataDirect::x.
Referenced by main(), and set_asigmabsol_direct().
void BSolDirect | ( | double * | b, |
double * | x, | ||
void * | data | ||
) |
Solver function of B.
Definition at line 135 of file evsl_cxsparse.c.
References CHOLMOD, cholmod_B, cholmod_E, cholmod_X, cholmod_Y, _BSolDataDirect::cm, cs_di_ipvec(), cs_di_lsolve(), cs_di_ltsolve(), cs_di_pvec(), cs_di_numeric::L, _BSolDataDirect::LB, _BSolDataDirect::n, _BSolDataDirect::N, cs_di_symbolic::pinv, _BSolDataDirect::S, and _BSolDataDirect::w.
Referenced by main(), and setup_bsol_direct().
void FreeASIGMABSolDirect | ( | int | num, |
void ** | data | ||
) |
free the data needed by CXSparse
free the data needed by CXSparse
Definition at line 352 of file evsl_cxsparse.c.
References _ASBSolDataDirect::b, cs_ci_nfree(), cs_ci_sfree(), _ASBSolDataDirect::N, _ASBSolDataDirect::S, and _ASBSolDataDirect::x.
Referenced by free_asigmabsol_direct(), and main().
void FreeBSolDirectData | ( | void * | data | ) |
Free solver data.
Definition at line 165 of file evsl_cxsparse.c.
References CHOLMOD, cholmod_E, cholmod_W, cholmod_Y, _BSolDataDirect::cm, cs_di_nfree(), cs_di_sfree(), _BSolDataDirect::LB, _BSolDataDirect::N, _BSolDataDirect::S, and _BSolDataDirect::w.
Referenced by free_bsol_direct(), and main().
void LTSolDirect | ( | double * | b, |
double * | x, | ||
void * | data | ||
) |
Solver function of L^{T} x = L^{-T}*b.
Definition at line 151 of file evsl_cxsparse.c.
References CHOLMOD, cholmod_B, cholmod_E, cholmod_W, cholmod_X, cholmod_Y, _BSolDataDirect::cm, cs_di_ltsolve(), cs_di_pvec(), cs_di_numeric::L, _BSolDataDirect::LB, _BSolDataDirect::n, _BSolDataDirect::N, cs_di_symbolic::pinv, _BSolDataDirect::S, and _BSolDataDirect::w.
Referenced by main(), and setup_bsol_direct().
setup CXsparse solver for A - SIGMA B
The setup invovles shifting the matrix A - SIGMA B and factorizing the shifted matrix The solver function and the data will be saved data Generally speaking, each pole can have a different solver
A | matrix A |
BB | matrix B, if NULL, it means B is identity |
num | the number of SIGMA's |
zk | array of SIGMA's of length num |
data | all data that are needed for solving the system |
setup CXsparse solver for A - SIGMA B
The setup invovles shifting the matrix A - SIGMA B and factorizing the shifted matrix The solver function and the data will be saved data Generally speaking, each pole can have a different solver
A | matrix A |
BB | matrix B, if NULL, it means B is identity |
num | the number of SIGMA's |
zk | array of SIGMA's of length num |
data | all data that are needed for solving the system |
Definition at line 190 of file evsl_cxsparse.c.
References _csrMat::a, _ASBSolDataDirect::b, Calloc, CHKERR, cs_ci_lu(), cs_ci_sqr(), cs_complex_t, evsl_timer(), evslstat, free_csr(), cs_ci_sparse::i, _csrMat::ia, _csrMat::ja, cs_ci_sparse::m, Malloc, matadd(), _ASBSolDataDirect::n, _ASBSolDataDirect::N, cs_ci_sparse::n, _csrMat::ncols, _csrMat::nrows, cs_ci_sparse::nz, cs_ci_sparse::nzmax, cs_ci_sparse::p, _ASBSolDataDirect::S, speye(), SuiteSparse_long, _evslstat::t_setASigBsv, _ASBSolDataDirect::x, and cs_ci_sparse::x.
Referenced by main(), and setup_asigmabsol_direct().
int SetupBSolDirect | ( | csrMat * | B, |
void ** | data | ||
) |
Setup the B-sol by computing the Cholesky factorization of B.
B | matrix B |
data | Struct which will be initialized |
Definition at line 57 of file evsl_cxsparse.c.
References _csrMat::a, arr_to_cholmod_dense(), CHKERR, CHOLMOD, cholmod_B, cholmod_X, _BSolDataDirect::cm, cs_di_chol(), cs_di_schol(), cs_di_spfree(), csrMat_to_cholmod_sparse(), evsl_timer(), evslstat, cs_di_sparse::i, _csrMat::ia, _csrMat::ja, _BSolDataDirect::LB, cs_di_sparse::m, Malloc, _BSolDataDirect::n, _BSolDataDirect::N, cs_di_sparse::n, _csrMat::nrows, cs_di_sparse::nz, cs_di_sparse::nzmax, cs_di_sparse::p, _BSolDataDirect::S, _evslstat::t_setBsv, _BSolDataDirect::w, and cs_di_sparse::x.
Referenced by main(), and setup_bsol_direct().