23 #include <cusp/detail/config.h>
94 template <
typename ValueType,
typename MemorySpace>
100 ValueType default_omega;
118 template <
typename MatrixType>
119 jacobi(
const MatrixType& A, ValueType omega=1.0);
127 template<
typename MemorySpace2>
129 : default_omega(A.default_omega), diagonal(A.diagonal), temp(A.temp){}
142 template <
typename MatrixType,
typename VectorType1,
typename VectorType2>
143 void operator()(
const MatrixType& A,
const VectorType1& b, VectorType2& x);
156 template <
typename MatrixType,
typename VectorType1,
typename VectorType2>
157 void operator()(
const MatrixType& A,
const VectorType1& b, VectorType2& x,
const ValueType omega);
165 #include <cusp/relaxation/detail/jacobi.inl>
jacobi(const jacobi< ValueType, MemorySpace2 > &A)
Represents a Jacobi relaxation scheme.
void operator()(const MatrixType &A, const VectorType1 &b, VectorType2 &x)
Abstract representation of a linear operator.
Abstract interface for iterative solvers.