tensorium::solver::Gauss< K > Class Template Reference

Direct Gaussian elimination solver with SIMD acceleration. More...

#include <LinearSolver.hpp>

Collaboration diagram for tensorium::solver::Gauss< K >:

Public Types
using	SimdT = simd::SimdTraits<K, DefaultISA>
using	regT = typename SimdT::reg

Public Member Functions
size_t	rows () const
	__attribute__ ((always_inline, hot, flatten)) static inline Vector< K > solve(const Matrix< K > &A_in
	Solve the linear system \( Ax = b \).
	assert (n==A_in.cols &&n==b_in.size())
Vector< K >	x (n)
aligned_vector< K >	rowi (n)
aligned_vector< K >	rowj (n)
	for (size_t i=0;i< n;++i)
	for (size_t ii=n;ii-- > 0;)

Static Public Member Functions
static void	raw_row_echelon (Matrix< K > &A, Vector< K > *b=nullptr, K eps=1e-12)

Public Attributes
size_t	block_size
aligned_vector< K >	data
const Vector< K > &	b_in
const size_t	n = A_in.rows
if(n >=1024) return Jacobi< K > Matrix< K >	M = A_in
Vector< K >	B = b_in
const size_t	W = SimdT::width
return	x

Detailed Description

template<typename K>
class tensorium::solver::Gauss< K >

Direct Gaussian elimination solver with SIMD acceleration.

This solver uses LU-style elimination with partial pivoting. For large systems ( \( n \geq 1024 \)), it redirects to Jacobi.

Template Parameters

K	Scalar type (must be floating-point)

Member Typedef Documentation

regT

template<typename K>

using tensorium::solver::Gauss< K >::regT = typename SimdT::reg

SimdT

template<typename K>

using tensorium::solver::Gauss< K >::SimdT = simd::SimdTraits<K, DefaultISA>

Member Function Documentation

__attribute__()

template<typename K>

tensorium::solver::Gauss< K >::__attribute__

(

(always_inline, hot, flatten)

)

const &

Solve the linear system \( Ax = b \).

Performs Gaussian elimination followed by back-substitution. Uses unrolled and SIMD-optimized loops for performance.

Parameters

A_in	Input matrix \( A \in \mathbb{R}^{n \times n} \)
b_in	Right-hand side vector \( b \in \mathbb{R}^n \)

Returns

Solution vector \( x \in \mathbb{R}^n \)

Exceptions

std::runtime_error

if matrix is singular or ill-conditioned

assert()

template<typename K>

tensorium::solver::Gauss< K >::assert

(

n

= =A_in.cols &&n==b_in.size()

)

References b_in, and n.

Referenced by raw_row_echelon().

Here is the caller graph for this function:

for() [1/2]

template<typename K>

tensorium::solver::Gauss< K >::for ( )

inline

References MathsUtils::_abs(), MathsUtils::_swap(), B, M, n, rowi(), rowj(), and W.

Here is the call graph for this function:

for() [2/2]

template<typename K>

tensorium::solver::Gauss< K >::for ( size_t ii = n; ii--, 0;  )

inline

References B, M, n, rowi(), W, and x.

Here is the call graph for this function:

raw_row_echelon()

template<typename K>

void tensorium::solver::Gauss< K >::raw_row_echelon ( Matrix< K > & A, Vector< K > * b = nullptr, K eps = 1e-12 )

inlinestatic

References MathsUtils::_abs(), MathsUtils::_swap(), assert(), tensorium::Matrix< K, RowMajor >::cols, n, tensorium::Matrix< K, RowMajor >::rows, and tensorium::Matrix< K, RowMajor >::swap_rows().

Referenced by tensorium::row_echelon().

Here is the call graph for this function:

Here is the caller graph for this function:

rowi()

template<typename K>

aligned_vector< K > tensorium::solver::Gauss< K >::rowi

(

n

)

References n.

Referenced by for(), and for().

Here is the caller graph for this function:

rowj()

template<typename K>

aligned_vector< K > tensorium::solver::Gauss< K >::rowj

(

n

)

Referenced by for().

Here is the caller graph for this function:

rows()

template<typename K>

size_t tensorium::solver::Gauss< K >::rows

(

)

const

x()

template<typename K>

Vector< K > tensorium::solver::Gauss< K >::x

(

n

)

References n.

Member Data Documentation

B

template<typename K>

Vector<K> tensorium::solver::Gauss< K >::B = b_in

Referenced by for(), and for().

b_in

template<typename K>

const Vector<K>& tensorium::solver::Gauss< K >::b_in

Initial value:

{
        static_assert(std::is_floating_point<K>::value, "")

Referenced by assert().

block_size

template<typename K>

size_t tensorium::solver::Gauss< K >::block_size

data

template<typename K>

aligned_vector<K> tensorium::solver::Gauss< K >::data

M

template<typename K>

if (n >= 1024) return Jacobi<K> Matrix<K> tensorium::solver::Gauss< K >::M = A_in

Referenced by for(), and for().

n

template<typename K>

const size_t tensorium::solver::Gauss< K >::n = A_in.rows

Referenced by assert(), for(), for(), raw_row_echelon(), rowi(), and x().

W

template<typename K>

const size_t tensorium::solver::Gauss< K >::W = SimdT::width

Referenced by for(), and for().

x

template<typename K>

return tensorium::solver::Gauss< K >::x

Referenced by for().

---

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

• includes/Tensorium/Core/LinearSolver.hpp