explore QudaDiracFieldOrder to see if we can reduce / remove reordering overheads when QUDA is used

[kostrzewa]

[kostrzewa]

This is somewhat of a pressing problem as with the MG on a dense machine for quark masses at or above the strange region (where the MG still outperforms double-half eo-CG), field remapping has a similar cost as the inversion itself...