54 Commits

Author	SHA1	Message	Date
Peter Boyle	b92060f511	5d OpDir direction interface refers to the 5d dims, not 4d to present a ... sensible and consistent external interface.	2015-06-09 22:41:59 +01:00
Peter Boyle	1e5b015ee3	Some unary ops and coarse grid support	2015-06-09 10:26:19 +01:00
Peter Boyle	d6f1ddf99c	Conjugate residual algorithm; some more unary functions	2015-06-08 12:04:59 +01:00
Peter Boyle	3845f267cb	Domain wall fermions now invert ; have the basis set up for ... Tanh/Zolo * (Cayley/PartFrac/ContFrac) * (Mobius/Shamir/Wilson) Approx Representation Kernel. All are done with space-time taking part in checkerboarding, Ls uncheckerboarded Have only so far tested the Domain Wall limit of mobius, and at that only checked that it i) Inverts ii) 5dim DW == Ls copies of 4dim D2 iii) MeeInv Mee == 1 iv) Meo+Mee+Moe+Moo == M unprec. v) MpcDagMpc is hermitan vi) Mdag is the adjoint of M between stochastic vectors. That said, the RB schur solve, RB MpcDagMpc solve, Unprec solve all converge and the true residual becomes small; so pretty good tests.	2015-06-02 16:57:12 +01:00