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
Added prototype s-step CG. Unfortunately, still not working properly.
Merge remote branch 'origin/quda-0.7' into sstep
Matrix-Powers CG (MPCG) seems to be working for a single outer iteration.
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.
Have SimpleBiCGstab solver working. Note that this solver will only
be used to test Communication-avoiding BiCGstab.
Added a prototype matrix-powers bicgstab routine, which is still buggy.
Modified SimpleBiCGstab so that updates to r and p do not use the
intermediate vector s. This is less efficient than the formulation
that uses s, but will be used for debugging purposes.
MPBiCGstab seems to be working.