#include "Teuchos_LAPACK.hpp"
#include "ROL_StdVector.hpp"

Functions
template<class Real >
void	trisolve (ROL::Ptr< Teuchos::LAPACK< int, Real > > lapack, const std::vector< Real > &a, const std::vector< Real > &b, const std::vector< Real > &c, const std::vector< Real > &r, std::vector< Real > &x)
	Solve a tridiagonal system. More...

template<class Real >
void	lusolve (ROL::Ptr< Teuchos::LAPACK< int, Real > > lapack, const std::vector< Real > &A, const std::vector< Real > &B, std::vector< Real > &X)

template<class Real >
void	cholsolve (ROL::Ptr< Teuchos::LAPACK< int, Real > > lapack, const std::vector< Real > &A, const std::vector< Real > &B, std::vector< Real > &X)

Function Documentation

template<class Real >

void trisolve	(	ROL::Ptr< Teuchos::LAPACK< int, Real > >	lapack,
		const std::vector< Real > &	a,
		const std::vector< Real > &	b,
		const std::vector< Real > &	c,
		const std::vector< Real > &	r,
		std::vector< Real > &	x
	)

Solve a tridiagonal system.

Parameters

[in]	lapack	pointer to LAPACK interface
[in]	a	subdiagonal band ( length N-1 )
[in]	b	diagonal band ( length N )
[in]	c	superidiagonal band ( length N-1 )
[in]	r	right-hand-side ( length N )
[out]	x	solution vector ( length N )

Definition at line 15 of file LinearAlgebra.hpp.

Referenced by rec_lobatto().

template<class Real >

void lusolve	(	ROL::Ptr< Teuchos::LAPACK< int, Real > >	lapack,
		const std::vector< Real > &	A,
		const std::vector< Real > &	B,
		std::vector< Real > &	X
	)

Solve a general system AX=B, with optional multiple (M) right-hand-sides using the LU factorization

Parameters

[in]	lapack	pointer to LAPACK interface
[in]	A	vector of column-stacked LHS matrix elements ( length N*N )
[in]	B	vector of column-stacked RHS matrix elements ( length N*M )
[in]	X	vector of column-stacked solution matrix elements ( length N*M )

Definition at line 50 of file LinearAlgebra.hpp.

template<class Real >

void cholsolve	(	ROL::Ptr< Teuchos::LAPACK< int, Real > >	lapack,
		const std::vector< Real > &	A,
		const std::vector< Real > &	B,
		std::vector< Real > &	X
	)

Solve a general system AX=B, with optional multiple (M) right-hand-sides using the Cholesky factorization

Parameters

[in]	lapack	pointer to LAPACK interface
[in]	A	vector of column-stacked LHS matrix elements ( length N*N )
[in]	B	vector of column-stacked RHS matrix elements ( length N*M )
[in]	X	vector of column-stacked solution matrix elements ( length N*M )

Definition at line 81 of file LinearAlgebra.hpp.

Functions