# driving forces (g stands for k^\perp.g/(2i\pi|k|) # a simple test driving force function g_test( kx::Int64, ky::Int64 ) if(kx==2 && ky==-1) return 0.5+sqrt(3)/2*1im elseif (kx==-2 && ky==1) return 0.5-sqrt(3)/2*1im else return 0. end end