IDR(s) is an efficient method for solving large nonsymmetric systems of linear equations.
More...
|
std::unique_ptr< LinOp > | transpose () const override |
| Returns a LinOp representing the transpose of the Transposable object.
|
std::unique_ptr< LinOp > | conj_transpose () const override |
| Returns a LinOp representing the conjugate transpose of the Transposable object.
|
bool | apply_uses_initial_guess () const override |
| Return true as iterative solvers use the data in x as an initial guess.
|
size_type | get_subspace_dim () const |
| Gets the subspace dimension of the solver.
|
void | set_subspace_dim (const size_type other) |
| Sets the subspace dimension of the solver.
|
remove_complex< ValueType > | get_kappa () const |
| Gets the kappa parameter of the solver.
|
void | set_kappa (const remove_complex< ValueType > other) |
| Sets the kappa parameter of the solver.
|
bool | get_deterministic () const |
| Gets the deterministic parameter of the solver.
|
void | set_deterministic (const bool other) |
| Sets the deterministic parameter of the solver.
|
bool | get_complex_subspace () const |
| Gets the complex_subspace parameter of the solver.
|
void | set_complex_subpsace (const bool other) |
| Sets the complex_subspace parameter of the solver.
|
void | set_complex_subspace (const bool other) |
| Sets the complex_subspace parameter of the solver.
|
const parameters_type & | get_parameters () const |
const Idr< default_precision > * | apply (ptr_param< const LinOp > b, ptr_param< LinOp > x) const |
std::unique_ptr< Idr< default_precision > > | create_default (std::shared_ptr< const Executor > exec) const |
std::unique_ptr< Idr< default_precision > > | clone (std::shared_ptr< const Executor > exec) const |
Idr< default_precision > * | copy_from (const PolymorphicObject *other) |
Idr< default_precision > * | move_from (ptr_param< PolymorphicObject > other) |
Idr< default_precision > * | clear () |
LinOp * | apply (ptr_param< const LinOp > b, ptr_param< LinOp > x) |
| Applies a linear operator to a vector (or a sequence of vectors).
|
const LinOp * | apply (ptr_param< const LinOp > b, ptr_param< LinOp > x) const |
LinOp * | apply (ptr_param< const LinOp > alpha, ptr_param< const LinOp > b, ptr_param< const LinOp > beta, ptr_param< LinOp > x) |
| Performs the operation x = alpha * op(b) + beta * x.
|
const LinOp * | apply (ptr_param< const LinOp > alpha, ptr_param< const LinOp > b, ptr_param< const LinOp > beta, ptr_param< LinOp > x) const |
const dim< 2 > & | get_size () const noexcept |
| Returns the size of the operator.
|
LinOp & | operator= (const LinOp &)=default |
| Copy-assigns a LinOp.
|
LinOp & | operator= (LinOp &&other) |
| Move-assigns a LinOp.
|
| LinOp (const LinOp &)=default |
| Copy-constructs a LinOp.
|
| LinOp (LinOp &&other) |
| Move-constructs a LinOp.
|
std::unique_ptr< LinOp > | create_default (std::shared_ptr< const Executor > exec) const |
std::unique_ptr< LinOp > | clone (std::shared_ptr< const Executor > exec) const |
LinOp * | copy_from (const PolymorphicObject *other) |
LinOp * | move_from (ptr_param< PolymorphicObject > other) |
LinOp * | clear () |
PolymorphicObject & | operator= (const PolymorphicObject &) |
std::unique_ptr< PolymorphicObject > | create_default (std::shared_ptr< const Executor > exec) const |
| Creates a new "default" object of the same dynamic type as this object.
|
std::unique_ptr< PolymorphicObject > | create_default () const |
| Creates a new "default" object of the same dynamic type as this object.
|
std::unique_ptr< PolymorphicObject > | clone (std::shared_ptr< const Executor > exec) const |
| Creates a clone of the object.
|
std::unique_ptr< PolymorphicObject > | clone () const |
| Creates a clone of the object.
|
PolymorphicObject * | copy_from (const PolymorphicObject *other) |
| Copies another object into this object.
|
template<typename Derived, typename Deleter> |
std::enable_if_t< std::is_base_of< PolymorphicObject, std::decay_t< Derived > >::value, PolymorphicObject > * | copy_from (std::unique_ptr< Derived, Deleter > &&other) |
| Moves another object into this object.
|
template<typename Derived, typename Deleter> |
std::enable_if_t< std::is_base_of< PolymorphicObject, std::decay_t< Derived > >::value, PolymorphicObject > * | copy_from (const std::unique_ptr< Derived, Deleter > &other) |
| Copies another object into this object.
|
PolymorphicObject * | copy_from (const std::shared_ptr< const PolymorphicObject > &other) |
| Copies another object into this object.
|
PolymorphicObject * | move_from (ptr_param< PolymorphicObject > other) |
| Moves another object into this object.
|
PolymorphicObject * | clear () |
| Transforms the object into its default state.
|
std::shared_ptr< const Executor > | get_executor () const noexcept |
| Returns the Executor of the object.
|
void | add_logger (std::shared_ptr< const Logger > logger) override |
| Adds a new logger to the list of subscribed loggers.
|
void | remove_logger (const Logger *logger) override |
| Removes a logger from the list of subscribed loggers.
|
const std::vector< std::shared_ptr< const Logger > > & | get_loggers () const override |
| Returns the vector containing all loggers registered at this object.
|
void | clear_loggers () override |
| Remove all loggers registered at this object.
|
void | remove_logger (ptr_param< const Logger > logger) |
void | convert_to (result_type *result) const override |
void | move_to (result_type *result) override |
EnableSolverBase & | operator= (const EnableSolverBase &other) |
| Creates a shallow copy of the provided system matrix, clones it onto this executor if executors don't match.
|
int | get_num_workspace_ops () const override |
std::vector< std::string > | get_workspace_op_names () const override |
std::vector< int > | get_workspace_scalars () const override |
| Returns the IDs of all scalars (workspace vectors with system dimension-independent size, usually 1 x num_rhs).
|
std::vector< int > | get_workspace_vectors () const override |
| Returns the IDs of all vectors (workspace vectors with system dimension-dependent size, usually system_matrix_size x num_rhs).
|
std::shared_ptr< const LinOp > | get_system_matrix () const |
| Returns the system matrix, with its concrete type, used by the solver.
|
EnableIterativeBase & | operator= (const EnableIterativeBase &other) |
| Creates a shallow copy of the provided stopping criterion, clones it onto this executor if executors don't match.
|
void | set_stop_criterion_factory (std::shared_ptr< const stop::CriterionFactory > new_stop_factory) override |
| Sets the stopping criterion of the solver.
|
std::shared_ptr< const stop::CriterionFactory > | get_stop_criterion_factory () const |
| Gets the stopping criterion factory of the solver.
|
void | set_preconditioner (std::shared_ptr< const LinOp > new_precond) override |
| Sets the preconditioner operator used by the Preconditionable.
|
EnablePreconditionable & | operator= (const EnablePreconditionable &other) |
| Creates a shallow copy of the provided preconditioner, clones it onto this executor if executors don't match.
|
virtual std::shared_ptr< const LinOp > | get_preconditioner () const |
| Returns the preconditioner operator used by the Preconditionable.
|
template<typename ValueType = default_precision>
class gko::solver::Idr< ValueType >
IDR(s) is an efficient method for solving large nonsymmetric systems of linear equations.
The implemented version is the one presented in the paper "Algorithm 913: An elegant IDR(s) variant that efficiently exploits
biorthogonality properties" by M. B. Van Gijzen and P. Sonneveld.
The method is based on the induced dimension reduction theorem which provides a way to construct subsequent residuals that lie in a sequence of shrinking subspaces. These subspaces are spanned by s vectors which are first generated randomly and then orthonormalized. They are stored in a dense matrix.
- Template Parameters
-
ValueType | precision of the elements of the system matrix. |