Incomplete LU solver class. More...

#include <monolish_equation.hpp>

Inheritance diagram for monolish::equation::ILU< MATRIX, Float >:

Collaboration diagram for monolish::equation::ILU< MATRIX, Float >:

Public Member Functions
void	apply_precond (const vector< Float > &r, vector< Float > &z)

void	create_precond (MATRIX &A)

std::string	name () const
	get solver name "monolish::equation::ILU" More...

int	solve (MATRIX &A, vector< Float > &x, vector< Float > &b)
	solve with incomplete LU factorization More...

std::string	solver_name () const
	get solver name "ILU" More...

	~ILU ()

Public Member Functions inherited from monolish::solver::solver< MATRIX, Float >
size_t	get_final_iter ()

double	get_final_residual ()

initvec_scheme	get_initvec_scheme () const
	get handling scheme of initial vector handling More...

int	get_lib () const
	get library option More...

size_t	get_maxiter () const
	get maxiter More...

size_t	get_miniter () const
	get miniter More...

Float	get_omega ()
	get the relaxation coefficient omega for SOR method (Default: 1.9) More...

bool	get_print_rhistory () const
	get print rhistory status More...

int	get_reorder ()
	0: no ordering 1: symrcm, 2: symamd, 3: csrmetisnd is used to reduce zero fill-in. More...

size_t	get_residual_method () const
	get residual method(default=0) More...

int	get_singularity ()
	-1 if A is symmetric positive definite. default reorder algorithm is csrmetisnd More...

double	get_tol () const
	get tolerance More...

template<class PRECOND >
void	set_apply_precond (PRECOND &p)
	set precondition apply function More...

template<class PRECOND >
void	set_create_precond (PRECOND &p)
	set precondition create function More...

void	set_initvec_scheme (initvec_scheme scheme)
	set how to handle initial vector More...

void	set_lib (int l)
	set library option (to change library, monolish, cusolver, etc.) More...

void	set_maxiter (size_t max)
	set max iter. (default = SIZE_MAX) More...

void	set_miniter (size_t min)
	set min iter. (default = 0) More...

void	set_omega (Float w)
	set the relaxation coefficient omega for SOR method (0 < w < 2, Default: 1.9) More...

void	set_print_rhistory (bool flag)
	print rhistory to standart out true/false. (default = false) More...

void	set_reorder (int r)
	0: no ordering 1: symrcm, 2: symamd, 3: csrmetisnd is used to reduce zero fill-in. More...

void	set_residual_method (size_t r)
	set residual method (default=0) More...

void	set_rhistory_filename (std::string file)
	rhistory filename More...

void	set_tol (double t)
	set tolerance (default:1.0e-8) More...

	solver ()
	create solver class More...

	~solver ()
	delete solver class More...

Private Member Functions
int	cusparse_ILU (MATRIX &A, vector< Float > &x, vector< Float > &b)

Private Attributes
monolish::vector< double >	buf

int	bufsize

void *	cusparse_handle = nullptr

void *	infoL = 0

void *	infoM = 0

void *	infoU = 0

void *	matL = 0

void *	matM = 0

void *	matU = 0

monolish::vector< Float >	zbuf

Additional Inherited Members
Protected Member Functions inherited from monolish::solver::solver< MATRIX, Float >
Float	get_residual (vector< Float > &x)

Protected Attributes inherited from monolish::solver::solver< MATRIX, Float >
size_t	final_iter = 0

double	final_resid = 0

initvec_scheme	initvecscheme = initvec_scheme::RANDOM

int	lib = 0

size_t	maxiter = SIZE_MAX

size_t	miniter = 0

Float	omega = 1.9

precondition< MATRIX, Float >	precond

bool	print_rhistory = false

int	reorder = 3

size_t	resid_method = 0

std::string	rhistory_file

std::ostream *	rhistory_stream

int	singularity

double	tol = 1.0e-8

Detailed Description

template<typename MATRIX, typename Float>
class monolish::equation::ILU< MATRIX, Float >

Incomplete LU solver class.

Note

attribute:

solver : true
preconditioner : false

input / architecture

Dense / Intel : false
Dense / NVIDIA : false
Dense / OSS : false
Sparse / Intel : false
Sparse / NVIDIA : true
Sparse / OSS : false

Definition at line 411 of file monolish_equation.hpp.

Constructor & Destructor Documentation

◆ ~ILU()

template<typename MATRIX , typename Float >

monolish::equation::ILU< MATRIX, Float >::~ILU ( )

Member Function Documentation

◆ apply_precond()

template<typename MATRIX , typename Float >

void monolish::equation::ILU< MATRIX, Float >::apply_precond	(	const vector< Float > &	r,
		vector< Float > &	z
	)

◆ create_precond()

template<typename MATRIX , typename Float >

void monolish::equation::ILU< MATRIX, Float >::create_precond ( MATRIX & A )

◆ cusparse_ILU()

template<typename MATRIX , typename Float >

int monolish::equation::ILU< MATRIX, Float >::cusparse_ILU	(	MATRIX &	A,
		vector< Float > &	x,
		vector< Float > &	b
	)

private

◆ name()

template<typename MATRIX , typename Float >

std::string monolish::equation::ILU< MATRIX, Float >::name ( ) const

inline

get solver name "monolish::equation::ILU"

Note

# of computation: 1
Multi-threading: false
GPU acceleration: false

Definition at line 439 of file monolish_equation.hpp.

◆ solve()

template<typename MATRIX , typename Float >

int monolish::equation::ILU< MATRIX, Float >::solve	(	MATRIX &	A,
		vector< Float > &	x,
		vector< Float > &	b
	)

solve with incomplete LU factorization

Warning: This solves Ax = b incompletely. In many cases the answer is wrong.

◆ solver_name()

template<typename MATRIX , typename Float >

std::string monolish::equation::ILU< MATRIX, Float >::solver_name ( ) const

inline

get solver name "ILU"

Note

# of computation: 1
Multi-threading: false
GPU acceleration: false

Definition at line 448 of file monolish_equation.hpp.

Member Data Documentation

◆ buf

template<typename MATRIX , typename Float >

monolish::vector<double> monolish::equation::ILU< MATRIX, Float >::buf

private

Definition at line 418 of file monolish_equation.hpp.

◆ bufsize

template<typename MATRIX , typename Float >

int monolish::equation::ILU< MATRIX, Float >::bufsize

private

Definition at line 417 of file monolish_equation.hpp.

◆ cusparse_handle

template<typename MATRIX , typename Float >

void* monolish::equation::ILU< MATRIX, Float >::cusparse_handle = nullptr

private

Definition at line 416 of file monolish_equation.hpp.

◆ infoL

template<typename MATRIX , typename Float >

void * monolish::equation::ILU< MATRIX, Float >::infoL = 0

private

Definition at line 415 of file monolish_equation.hpp.

◆ infoM

template<typename MATRIX , typename Float >

void* monolish::equation::ILU< MATRIX, Float >::infoM = 0

private

Definition at line 415 of file monolish_equation.hpp.

◆ infoU

template<typename MATRIX , typename Float >

void * monolish::equation::ILU< MATRIX, Float >::infoU = 0

private

Definition at line 415 of file monolish_equation.hpp.

◆ matL

template<typename MATRIX , typename Float >

void * monolish::equation::ILU< MATRIX, Float >::matL = 0

private

Definition at line 414 of file monolish_equation.hpp.

◆ matM

template<typename MATRIX , typename Float >

void* monolish::equation::ILU< MATRIX, Float >::matM = 0

private

Definition at line 414 of file monolish_equation.hpp.

◆ matU

template<typename MATRIX , typename Float >

void * monolish::equation::ILU< MATRIX, Float >::matU = 0

private

Definition at line 414 of file monolish_equation.hpp.

◆ zbuf

template<typename MATRIX , typename Float >

monolish::vector<Float> monolish::equation::ILU< MATRIX, Float >::zbuf

private

Definition at line 419 of file monolish_equation.hpp.

The documentation for this class was generated from the following file:

monolish_equation.hpp

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<typename MATRIX, typename Float> class monolish::equation::ILU< MATRIX, Float >

Constructor & Destructor Documentation

◆ ~ILU()

Member Function Documentation

◆ apply_precond()

◆ create_precond()

◆ cusparse_ILU()

◆ name()

◆ solve()

◆ solver_name()

Member Data Documentation

◆ buf

◆ bufsize

◆ cusparse_handle

◆ infoL

◆ infoM

◆ infoU

◆ matL

◆ matM

◆ matU

◆ zbuf

template<typename MATRIX, typename Float>
class monolish::equation::ILU< MATRIX, Float >