asolver_SAP.h

Go to the documentation of this file.

 
#ifndef ASOLVER_SAP_H
#define ASOLVER_SAP_H
 
#include <cstdio>
#include <cstdlib>
#include <vector>
using std::vector;
#include <string>
using std::string;
 
#include "lib_alt/Solver/asolver.h"
#include "lib_alt/Fopr/afopr_dd.h"
#include "asolver_SAP_MINRES.h"
 
template<typename AFIELD>
class ASolver_SAP : public ASolver<AFIELD>
{
 public:
  typedef typename AFIELD::real_t                  real_t;
  typedef AIndex_block_lex<real_t, AFIELD::IMPL>   block_index_t;
  using ASolver<AFIELD>::m_vl;
  static const std::string class_name;
 
 protected:
 
  const block_index_t *m_block_index;
  AFopr_dd<AFIELD> *m_fopr;   // given from outside
  unique_ptr<ASolver_SAP_MINRES<AFIELD> > m_sap_minres;
 
  int m_Niter;           
  real_t m_Stop_cond;    
  int m_Nconv;           
 
  const int m_min_res_iter = 6;
  //const int m_min_res_iter = 20;
 
  int m_nconv;
 
  ASolver_SAP() { }
 
  AFIELD m_x, m_r, m_p;
 
 
 public:
  //  ASolver_SAP(AFopr<AFIELD>* fopr, const block_index_t *block_index)
  ASolver_SAP(AFopr_dd<AFIELD> *fopr, const block_index_t *block_index)
    : m_Niter(0), m_Stop_cond(0.0L)
  {
    m_fopr        = fopr;
    m_block_index = block_index;
    this->init();
  }
 
  ~ASolver_SAP() { 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_H
//============================================================END=====