gmres.h

Go to the documentation of this file.

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

#ifndef MYRAMATH_ITERATIVE_GMRES_H
#define MYRAMATH_ITERATIVE_GMRES_H

#include <myramath/MYRAMATH_EXPORT.h>
#include <myramath/utility/Number.h>

#include <vector>

namespace myra {

// Forward declarations.
template<class Number> class Action;
template<class Number> class CVectorRange;
template<class Number> class  VectorRange;
template<class Number> class CMatrixRange;
template<class Number> class  MatrixRange;

template<class Number> class MYRAMATH_EXPORT gmres_output
  {
  public:

    typedef typename ReflectPrecision<Number>::type Precision;

    // Did gmres() converge?
    bool converged;
  
    // History of relative residual, |A*x-b|/|b|
    std::vector<Precision> history;

    // Constructs from components.
    gmres_output(bool in_converged, const std::vector<Precision>& in_history)
      {
      converged = in_converged;
      history = in_history;
      }
    
  };

MYRAMATH_EXPORT gmres_output<NumberS> gmres(const Action<NumberS>& M1, const Action<NumberS>& A, const Action<NumberS>& M2, const CVectorRange<NumberS>& b, const VectorRange<NumberS>& x, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberD> gmres(const Action<NumberD>& M1, const Action<NumberD>& A, const Action<NumberD>& M2, const CVectorRange<NumberD>& b, const VectorRange<NumberD>& x, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberC> gmres(const Action<NumberC>& M1, const Action<NumberC>& A, const Action<NumberC>& M2, const CVectorRange<NumberC>& b, const VectorRange<NumberC>& x, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberZ> gmres(const Action<NumberZ>& M1, const Action<NumberZ>& A, const Action<NumberZ>& M2, const CVectorRange<NumberZ>& b, const VectorRange<NumberZ>& x, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);

MYRAMATH_EXPORT gmres_output<NumberS> gmres(const Action<NumberS>& M1, const Action<NumberS>& A, const CVectorRange<NumberS>& b, const VectorRange<NumberS>& x, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberD> gmres(const Action<NumberD>& M1, const Action<NumberD>& A, const CVectorRange<NumberD>& b, const VectorRange<NumberD>& x, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberC> gmres(const Action<NumberC>& M1, const Action<NumberC>& A, const CVectorRange<NumberC>& b, const VectorRange<NumberC>& x, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberZ> gmres(const Action<NumberZ>& M1, const Action<NumberZ>& A, const CVectorRange<NumberZ>& b, const VectorRange<NumberZ>& x, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);

MYRAMATH_EXPORT gmres_output<NumberS> gmres(const Action<NumberS>& M1, const Action<NumberS>& A, const Action<NumberS>& M2, const CMatrixRange<NumberS>& B, const MatrixRange<NumberS>& X, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberD> gmres(const Action<NumberD>& M1, const Action<NumberD>& A, const Action<NumberD>& M2, const CMatrixRange<NumberD>& B, const MatrixRange<NumberD>& X, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberC> gmres(const Action<NumberC>& M1, const Action<NumberC>& A, const Action<NumberC>& M2, const CMatrixRange<NumberC>& B, const MatrixRange<NumberC>& X, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberZ> gmres(const Action<NumberZ>& M1, const Action<NumberZ>& A, const Action<NumberZ>& M2, const CMatrixRange<NumberZ>& B, const MatrixRange<NumberZ>& X, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);

MYRAMATH_EXPORT gmres_output<NumberS> gmres(const Action<NumberS>& M1, const Action<NumberS>& A, const CMatrixRange<NumberS>& B, const MatrixRange<NumberS>& X, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberD> gmres(const Action<NumberD>& M1, const Action<NumberD>& A, const CMatrixRange<NumberD>& B, const MatrixRange<NumberD>& X, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberC> gmres(const Action<NumberC>& M1, const Action<NumberC>& A, const CMatrixRange<NumberC>& B, const MatrixRange<NumberC>& X, NumberS tolerance = 1.0e-4, int restart = 100, int iterations = 1);
MYRAMATH_EXPORT gmres_output<NumberZ> gmres(const Action<NumberZ>& M1, const Action<NumberZ>& A, const CMatrixRange<NumberZ>& B, const MatrixRange<NumberZ>& X, NumberD tolerance = 1.0e-4, int restart = 100, int iterations = 1);

} // namespace myra

#endif