44 #ifndef ROL_CONJUGATERESIDUALS_H
45 #define ROL_CONJUGATERESIDUALS_H
63 ROL::Ptr<Vector<Real> >
r_;
64 ROL::Ptr<Vector<Real> >
v_;
65 ROL::Ptr<Vector<Real> >
p_;
66 ROL::Ptr<Vector<Real> >
Ap_;
67 ROL::Ptr<Vector<Real> >
MAp_;
70 ConjugateResiduals( Real absTol = 1.e-4, Real relTol = 1.e-2,
int maxit = 100,
bool useInexact =
false )
75 int &iter,
int &flag ) {
86 Real rnorm = b.
norm();
88 Real itol = std::sqrt(ROL_EPSILON<Real>());
112 Real kappa(0), beta(0), alpha(0), tmp(0);
113 Real gHg =
r_->dot(
v_->dual());
115 for (iter = 0; iter < (int)Krylov<Real>::getMaximumIteration(); iter++) {
116 itol = std::sqrt(ROL_EPSILON<Real>());
118 kappa =
MAp_->dot(
Ap_->dual());
129 if ( rnorm < rtol ) {
138 gHg =
r_->dot(
v_->dual());
virtual ROL::Ptr< Vector > clone() const =0
Clone to make a new (uninitialized) vector.
Provides definition of the Conjugate Residual solver.
virtual void axpy(const Real alpha, const Vector &x)
Compute where .
Contains definitions of custom data types in ROL.
virtual void apply(Vector< Real > &Hv, const Vector< Real > &v, Real &tol) const =0
Apply linear operator.
virtual void zero()
Set to zero vector.
Defines the linear algebra or vector space interface.
ROL::Ptr< Vector< Real > > r_
ROL::Ptr< Vector< Real > > Ap_
virtual void applyInverse(Vector< Real > &Hv, const Vector< Real > &v, Real &tol) const
Apply inverse of linear operator.
Real run(Vector< Real > &x, LinearOperator< Real > &A, const Vector< Real > &b, LinearOperator< Real > &M, int &iter, int &flag)
Provides definitions for Krylov solvers.
Provides the interface to apply a linear operator.
ROL::Ptr< Vector< Real > > v_
ConjugateResiduals(Real absTol=1.e-4, Real relTol=1.e-2, int maxit=100, bool useInexact=false)
ROL::Ptr< Vector< Real > > p_
ROL::Ptr< Vector< Real > > MAp_
virtual Real norm() const =0
Returns where .