asolver_SAP_QWS.h

Go to the documentation of this file.

 
#ifndef ASOLVER_SAP_QWS_H
#define ASOLVER_SAP_QWS_H
 
#include <cstdio>
#include <cstdlib>
#include <vector>
using std::vector;
#include <string>
using std::string;
 
#include "lib_alt/Solver/asolver.h"
#include "lib_alt_QXS/Fopr/afopr_Clover_QWS_dd.h"
 
#ifdef USE_QWSLIB
#include "lib_alt_QXS/extra/qws_lib.h"
#endif
 
 
template<typename AFIELD>
class ASolver_SAP_QWS : public ASolver<AFIELD>
{
 public:
  typedef typename AFIELD::real_t real_t;
  using ASolver<AFIELD>::m_vl;
  static const std::string class_name;
 
 protected:
 
  //  const block_index_t *m_block_index;
  AFopr_Clover_QWS_dd<AFIELD> *m_fopr; // given from outside
 
  int m_Niter;                         
  real_t m_Stop_cond;                  
  int m_Nconv;                         
  double m_flop;                       
  const int m_nm = 10;                 
 
  int m_nconv;
 
  ASolver_SAP_QWS() { }
 
  AFIELD m_x, m_r, m_b;
 
#ifdef USE_QWSLIB
  scs_t *m_s, *m_q;
#endif
 
 public:
  ASolver_SAP_QWS(AFopr_dd<AFIELD> *fopr)
    : m_Niter(0), m_Stop_cond(0.0L)
  {
#ifndef USE_QWSLIB
    vout.crucial("%s: USE_QWSLIB is not defined\n", class_name.c_str());
    exit(EXIT_FAILURE);
#endif
 
    m_fopr = dynamic_cast<AFopr_Clover_QWS_dd<AFIELD> *>(fopr);
    if (m_fopr == nullptr) {
      vout.crucial(m_vl, "asolver_SAP_QWS: bad fopr is given, must be AFopr_Clover_QWS_dd\n");
      exit(EXIT_FAILURE);
    }
    this->init();
  }
 
  ~ASolver_SAP_QWS() { this->tidyup(); }
 
  void set_parameters(const Parameters& params);
 
  void set_parameters(const int Niter, const real_t Stop_cond);
 
  void solve(AFIELD& xq, const AFIELD& b, int& nconv, real_t& diff);
 
  AFopr<AFIELD> *get_fopr() { return m_fopr; }
 
  double flop_count();
 
 protected:
 
  void init(void);
 
  void tidyup(void);
};
 
#endif // ASOLVER_SAP_QWS_H
//============================================================END=====