Changes:
Non-EO Operators double+single degenerate+ND
Common gpu solver interface for all existing operators
OpenMP for temporalgauge
GPU blocks in main code wrapped with if(usegpu_flag) + #ifdef HAVE_GPU