lu_solve.h

Go to the documentation of this file.

// -*- C++ -*-
/*! \file
 * \brief LU solver
 *
 * Predictors for HMC
 */
 
#ifndef LU_SOLVE_H
#define LU_SOLVE_H
 
#include "chromabase.h"
 
using namespace QDP;
 
namespace Chroma 
{
 
  //! Solve M a = b by LU decomposition with partial pivoting
  /*! @ingroup predictor */
  void LUSolve( multi1d<DComplex>& a, 
                const multi2d<DComplex>& M, 
                const multi1d<DComplex>& b );
 
}
 
#endif