Default solver function for solving shifted systems and factoring B. More...

#include <string.h>
#include "def.h"
#include "struct.h"
#include "internal_proto.h"
#include "cholmod.h"
#include "umfpack.h"
#include "evsl_direct.h"

Data Structures
struct	_BSolDataDirect

Macros
#define	CHOLMOD(name) cholmod_ ## name

Typedefs
typedef struct _BSolDataDirect	BSolDataDirect

Functions
cholmod_sparse *	csrMat_to_cholmod_sparse (csrMat *A, int stype)
	Create cholmod_sparse matrix just as a wrapper of a csrMat This will be useful for the B matrix, since it will be factored. More...

void	arr_to_cholmod_dense (int m, int n, double A, cholmod_dense B)
	wrap an array in cholmod dense (no copying) More...

void	BSolDirect (double b, double x, void *data)
	Solver function of B. More...

int	SetupBSolDirect (csrMat B, void *data)
	Setup the B-sol by computing the Cholesky factorization of B. More...

void	FreeBSolDirectData (void *vdata)
	Free solver data. More...

void	LTSolDirect (double b, double x, void *data)
	Solver function of L^{T} x = L^{-T}*b. More...

void	ASIGMABSolDirect (int n, double br, double bi, double xr, double xz, void *data)
	default complex linear solver routine passed to evsl More...

int	SetupASIGMABSolDirect (csrMat A, csrMat BB, int num, complex double zk, void *data)
	setup the default solver for A - SIGMA B More...

void	FreeASIGMABSolDirect (int num, void **data)
	free the data needed by UMFpack More...

Variables
cholmod_dense	cholmod_X

cholmod_dense	cholmod_B

cholmod_dense *	cholmod_Y =NULL

cholmod_dense *	cholmod_E =NULL

cholmod_dense *	cholmod_W =NULL

Detailed Description

Default solver function for solving shifted systems and factoring B.

This file contains the default solver function for solving shifted systems with A - SIGMA*I or A - SIGMA*B, and for factoring the B matrix. The default solvers are UMFPACK and CHOLMOD from SuiteSparse

Definition in file evsl_suitesparse.c.

Macro Definition Documentation

#define CHOLMOD ( name ) cholmod_ ## name

Definition at line 13 of file evsl_suitesparse.c.

Referenced by BSolDirect(), FreeBSolDirectData(), LTSolDirect(), and SetupBSolDirect().

Typedef Documentation

typedef struct _BSolDataDirect BSolDataDirect

Function Documentation

void arr_to_cholmod_dense	(	int	m,
		int	n,
		double *	A,
		cholmod_dense *	B
	)

wrap an array in cholmod dense (no copying)

Definition at line 79 of file evsl_suitesparse.c.

Referenced by SetupBSolDirect().

Here is the caller graph for this function:

void ASIGMABSolDirect	(	int	n,
		double *	br,
		double *	bi,
		double *	xr,
		double *	xz,
		void *	data
	)

default complex linear solver routine passed to evsl

complex linear solver routine passed to evsl

Parameters

	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)

Warning: : This function MUST be of this prototype

Definition at line 249 of file evsl_suitesparse.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.

Here is the call graph for this function:

void BSolDirect	(	double *	b,
		double *	x,
		void *	data
	)

Solver function of B.

Definition at line 118 of file evsl_suitesparse.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.

Here is the call graph for this function:

cholmod_sparse* csrMat_to_cholmod_sparse	(	csrMat *	A,
		int	stype
	)

Create cholmod_sparse matrix just as a wrapper of a csrMat This will be useful for the B matrix, since it will be factored.

Warning: cholmod_sparse is a CSC format. But since B is symmetric, it is the same

Definition at line 34 of file evsl_suitesparse.c.

References _csrMat::a, _csrMat::ia, _csrMat::ja, Malloc, _csrMat::ncols, _csrMat::nrows, and SuiteSparse_long.

Referenced by SetupBSolDirect().

Here is the caller graph for this function:

void FreeASIGMABSolDirect	(	int	num,
		void **	data
	)

free the data needed by UMFpack

free the data needed by CXSparse

Definition at line 380 of file evsl_suitesparse.c.

References _ASBSolDataDirect::b, cs_ci_nfree(), cs_ci_sfree(), _ASBSolDataDirect::N, _ASBSolDataDirect::S, and _ASBSolDataDirect::x.

Here is the call graph for this function:

void FreeBSolDirectData ( void * vdata )

Free solver data.

Definition at line 198 of file evsl_suitesparse.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.

Here is the call graph for this function:

void LTSolDirect	(	double *	b,
		double *	x,
		void *	data
	)

Solver function of L^{T} x = L^{-T}*b.

Definition at line 215 of file evsl_suitesparse.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.

Here is the call graph for this function:

int SetupASIGMABSolDirect	(	csrMat *	A,
		csrMat *	BB,
		int	num,
		complex double *	zk,
		void **	data
	)

setup the default solver for A - SIGMA B

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

Parameters

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