julia/driving.jl


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

# 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