“…The ability to compute LGFs can be useful, for example, for simulations of Hubbard models 1 and for non-perturbative renormalization group studies of scalar boson models. 2,3 For i = j, G ij (ω) can be expressed as closed forms in terms of named special functions (mainly elliptic integrals or generalized hypergeometric functions) for square, bcc, honeycomb, diamond, cubic, hypercubic, triangular, and fcc lattices. [4][5][6][7][8][9][10][11][12][13][14][15][16] For i = j, G ij (ω) can be expressed as closed forms for square, bcc, 14 honeycomb, 14 triangular, 14 kagome, 14 diced, 14 and cubic 10 lattices; spatial recurrence relations exist but are often numerically unstable.…”