In the electronk.jl example defining the "Adagdn" operator for SiteType "ElecK", the current implementation is:
Op[s' => 3, s => 1] = 1.0
return Op[s' => 2, s => 4] = 1.0
It seems the second matrix element should instead be:
Op[s' => 4, s => 2] = 1.0
Could you please confirm whether this is a typo?
Best