Grid 0.7.0
Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type > Class Template Reference

#include <MixedPrecisionFlexibleGeneralisedMinimalResidual.h>

Inheritance diagram for Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >:
OperatorFunction< FieldD >

Public Member Functions

 MixedPrecisionFlexibleGeneralisedMinimalResidual (RealD tol, Integer maxit, GridBase *sp_grid, LinearFunction< FieldF > &Prec, Integer restart_length, bool err_on_no_conv=true)
void operator() (LinearOperatorBase< FieldD > &LinOp, const FieldD &src, FieldD &psi)
RealD outerLoopBody (LinearOperatorBase< FieldD > &LinOp, const FieldD &src, FieldD &psi, RealD rsq)
void arnoldiStep (LinearOperatorBase< FieldD > &LinOp, std::vector< FieldD > &v, std::vector< FieldD > &z, FieldD &w, int iter)
void qrUpdate (int iter)
void computeSolution (std::vector< FieldD > const &z, FieldD &psi, int iter)

Data Fields

bool ErrorOnNoConverge
RealD Tolerance
Integer MaxIterations
Integer RestartLength
Integer MaxNumberOfRestarts
Integer IterationCount
GridStopWatch MatrixTimer
GridStopWatch PrecTimer
GridStopWatch LinalgTimer
GridStopWatch QrTimer
GridStopWatch CompSolutionTimer
GridStopWatch ChangePrecTimer
Eigen::MatrixXcd H
std::vector< ComplexD > y
std::vector< ComplexD > gamma
std::vector< ComplexD > c
std::vector< ComplexD > s
GridBase * SinglePrecGrid
LinearFunction< FieldF > & Preconditioner

Detailed Description

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
class Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >

Definition at line 35 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Constructor & Destructor Documentation

◆ MixedPrecisionFlexibleGeneralisedMinimalResidual()

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::MixedPrecisionFlexibleGeneralisedMinimalResidual ( RealD tol,
Integer maxit,
GridBase * sp_grid,
LinearFunction< FieldF > & Prec,
Integer restart_length,
bool err_on_no_conv = true )
inline

Definition at line 69 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

References c, ErrorOnNoConverge, gamma, H, MaxIterations, MaxNumberOfRestarts, Preconditioner, RestartLength, s, SinglePrecGrid, Tolerance, and y.

Member Function Documentation

◆ operator()()

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
void Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::operator() ( LinearOperatorBase< FieldD > & LinOp,
const FieldD & src,
FieldD & psi )
inlinevirtual

Implements OperatorFunction< FieldD >.

Definition at line 88 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

References axpy(), ChangePrecTimer, CompSolutionTimer, conformable(), GridStopWatch::Elapsed(), ErrorOnNoConverge, GridLogIterative, GridLogMessage, IterationCount, LinalgTimer, MatrixTimer, MaxNumberOfRestarts, norm2(), LinearOperatorBase< Field >::Op(), outerLoopBody(), PrecTimer, QrTimer, sqrt(), GridStopWatch::Start(), GridStopWatch::Stop(), and Tolerance.

◆ outerLoopBody()

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
RealD Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::outerLoopBody ( LinearOperatorBase< FieldD > & LinOp,
const FieldD & src,
FieldD & psi,
RealD rsq )
inline

Definition at line 156 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

References arnoldiStep(), computeSolution(), gamma, GridLogIterative, IterationCount, LinalgTimer, MatrixTimer, MaxIterations, norm2(), LinearOperatorBase< Field >::Op(), qrUpdate(), RestartLength, and sqrt().

Referenced by operator()().

◆ arnoldiStep()

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
void Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::arnoldiStep ( LinearOperatorBase< FieldD > & LinOp,
std::vector< FieldD > & v,
std::vector< FieldD > & z,
FieldD & w,
int iter )
inline

Definition at line 204 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

References ChangePrecTimer, H, innerProduct(), LinalgTimer, MatrixTimer, norm2(), LinearOperatorBase< Field >::Op(), precisionChange(), Preconditioner, PrecTimer, SinglePrecGrid, and sqrt().

Referenced by outerLoopBody().

◆ qrUpdate()

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
void Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::qrUpdate ( int iter)
inline

Definition at line 237 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

References c, conjugate(), gamma, H, QrTimer, s, and sqrt().

Referenced by outerLoopBody().

◆ computeSolution()

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
void Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::computeSolution ( std::vector< FieldD > const & z,
FieldD & psi,
int iter )
inline

Definition at line 260 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

References CompSolutionTimer, gamma, H, and y.

Referenced by outerLoopBody().

Field Documentation

◆ ErrorOnNoConverge

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
bool Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::ErrorOnNoConverge

Definition at line 40 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by MixedPrecisionFlexibleGeneralisedMinimalResidual(), and operator()().

◆ Tolerance

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
RealD Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::Tolerance

Definition at line 43 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by MixedPrecisionFlexibleGeneralisedMinimalResidual(), and operator()().

◆ MaxIterations

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
Integer Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::MaxIterations

Definition at line 45 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by MixedPrecisionFlexibleGeneralisedMinimalResidual(), and outerLoopBody().

◆ RestartLength

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
Integer Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::RestartLength

Definition at line 46 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by MixedPrecisionFlexibleGeneralisedMinimalResidual(), and outerLoopBody().

◆ MaxNumberOfRestarts

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
Integer Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::MaxNumberOfRestarts

Definition at line 47 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by MixedPrecisionFlexibleGeneralisedMinimalResidual(), and operator()().

◆ IterationCount

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
Integer Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::IterationCount

Definition at line 48 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by operator()(), and outerLoopBody().

◆ MatrixTimer

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
GridStopWatch Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::MatrixTimer

Definition at line 51 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by arnoldiStep(), operator()(), and outerLoopBody().

◆ PrecTimer

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
GridStopWatch Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::PrecTimer

Definition at line 52 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by arnoldiStep(), and operator()().

◆ LinalgTimer

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
GridStopWatch Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::LinalgTimer

Definition at line 53 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by arnoldiStep(), operator()(), and outerLoopBody().

◆ QrTimer

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
GridStopWatch Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::QrTimer

Definition at line 54 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by operator()(), and qrUpdate().

◆ CompSolutionTimer

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
GridStopWatch Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::CompSolutionTimer

Definition at line 55 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by computeSolution(), and operator()().

◆ ChangePrecTimer

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
GridStopWatch Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::ChangePrecTimer

Definition at line 56 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by arnoldiStep(), and operator()().

◆ H

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
Eigen::MatrixXcd Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::H

Definition at line 58 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by arnoldiStep(), computeSolution(), MixedPrecisionFlexibleGeneralisedMinimalResidual(), and qrUpdate().

◆ y

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
std::vector<ComplexD> Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::y

Definition at line 60 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by computeSolution(), and MixedPrecisionFlexibleGeneralisedMinimalResidual().

◆ gamma

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
std::vector<ComplexD> Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::gamma

Definition at line 61 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by computeSolution(), MixedPrecisionFlexibleGeneralisedMinimalResidual(), outerLoopBody(), and qrUpdate().

◆ c

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
std::vector<ComplexD> Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::c

Definition at line 62 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by MixedPrecisionFlexibleGeneralisedMinimalResidual(), and qrUpdate().

◆ s

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
std::vector<ComplexD> Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::s

Definition at line 63 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by MixedPrecisionFlexibleGeneralisedMinimalResidual(), and qrUpdate().

◆ SinglePrecGrid

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
GridBase* Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::SinglePrecGrid

Definition at line 65 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by arnoldiStep(), and MixedPrecisionFlexibleGeneralisedMinimalResidual().

◆ Preconditioner

template<class FieldD, class FieldF, typename std::enable_if< getPrecision< FieldD >::value==2, int >::type = 0, typename std::enable_if< getPrecision< FieldF >::value==1, int >::type = 0>
LinearFunction<FieldF>& Grid::MixedPrecisionFlexibleGeneralisedMinimalResidual< FieldD, FieldF, type, type >::Preconditioner

Definition at line 67 of file MixedPrecisionFlexibleGeneralisedMinimalResidual.h.

Referenced by arnoldiStep(), and MixedPrecisionFlexibleGeneralisedMinimalResidual().

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