avp_altivecf_solver.h

Go to the documentation of this file.

// -*- C++ -*-
// $Id: avp_altivecf_solver.h,v 3.3 2007/01/17 02:39:27 edwards Exp $
/*! \file
 *  \brief DWF/Bluegene altivec solver
 */

#ifndef AVP_ALTIVEC_SOLVER_H
#define AVP_ALTIVEC_SOLVER_H

#include "chromabase.h"
#include "avp_inverter_interface.h"

extern "C" 
{
  struct MIT_altivecf_DWF_Gauge;
  struct MIT_altivecf_DWF_Fermion;
};

namespace Chroma 
{ 
  //! Bluegene single-prec solver
  /*!
   * \ingroup qprop
   *
   * @{
   */
  namespace AVPSolver 
  { 
    //! Bluegene altivec solver
    class AltiVecDWFSolverF : public AVPSolverInterface< MIT_altivecf_DWF_Gauge, MIT_altivecf_DWF_Fermion > 
    {
    public:
    protected:
      MIT_altivecf_DWF_Fermion* loadFermionRHS(const void* OuterFermion) const;
      MIT_altivecf_DWF_Fermion* loadFermionGuess(const void *OuterFermion) const;
      MIT_altivecf_DWF_Fermion* allocateFermion(void) const ;
      void saveFermionSolver(void *OuterFermion, 
                             MIT_altivecf_DWF_Fermion* CGFermion) const;

      void saveFermionOperator(void *OuterFermion, 
                               MIT_altivecf_DWF_Fermion* CGFermion) const;

      void deleteFermion(MIT_altivecf_DWF_Fermion* ptr) const;
      int cgInternal(MIT_altivecf_DWF_Fermion       *psi,
                     double        *out_eps,
                     int           *out_iter,
                     double        M,
                     double        m_f,
                     const MIT_altivecf_DWF_Fermion *x0,
                     const MIT_altivecf_DWF_Fermion *eta,
                     double        eps,
                     int           min_iter,
                     int           max_iter)  const;
    public:
      void loadGauge(const void *u,
                     const void *v);
     
      void deleteGauge(void);
     
      // Init the system -- Constructor call?
      int init(const int lattice[5],
               void *(*allocator)(size_t size),
               void (*deallocator)(void *));
     
      // Finalize - destructor call
      void fini(void);
    private:
      MIT_altivecf_DWF_Gauge *g;
    };
  }

  /*! @} */   // end of group qprop
}

#endif

---