Added some orthogonalization routines
Added two variants of the CG3 solver
Tweaked the hisq tests to eliminate some warnings
Modified make.inc.in so that the code only links to the pthread library
when POSIX_THREADS = yes
QUDA residual type is now set to include one or multiple residual types to reflect behavior of MILC
update MILC interface to allow L2 and heavy quark residual to be specified
Added prototype s-step CG. Unfortunately, still not working properly.
Merge remote branch 'origin/quda-0.7' into sstep
fixed check whether heavy-quark residual is set in SetInvertParams for MILC multi shift invert
Matrix-Powers CG (MPCG) seems to be working for a single outer iteration.
Merge branch 'quda-0.7' into mw/milc_invert
The matrix powers kernel with Nstep = 2 now works, but it's still
a complete mess. The formulation presented in Hoemenn's PhD thesis
is incorrect.
Got Matrix-powers CG working.
Beginning to clean-up matrix-powers CG
Merged the j=0 and j!=0 loops of the Matrix-powers CG. Substantial
cleanup.
Now use the Gram matrix to compute r2 and rAr in each iteration.
The Gram matrix is currently stored as a complex 2d array, but
only the real components are actually used.
Now store the Gram matrix in a 2d array of doubles
Removed the MPCG3 solver, which was for testing purposes.
The Matrix-powers CG solver now returns the correct value of x.
Need to refine the stopping criterion a bit.
Removed MPCG3 solver declaration in invert_quda.h.
This was only needed for testing purposes.