LSQRInterface.h

#ifndef OPTIMIZATION_LSQR_INTERFACE_H
#define OPTIMIZATION_LSQR_INTERFACE_H

#include <KrisLibrary/Logger.h>
#include <KrisLibrary/math/sparsematrix.h>

namespace Optimization
{
using namespace Math;

struct LSQRInterface
{
  LSQRInterface();
  bool Solve(const SparseMatrix& A,const Vector& b);

  //input quantities
  Vector x0;       
  Real dampValue;  
  Real relError;   
  Real condLimit;  
  int maxIters;    
  int verbose;     

  //output quantities
  Vector x;        
  Vector stdErr;  
  int numIters;    
  Real condEstA;   
  Real residualNorm; //<the estimate of the final residual's norm
};

} //namespace Optimization

#endif