单格点的格点格林函数

单格点的格点格林函数是指在一个有限的格点网格中，某一格点与自身之间的格点格林函数。由于格点网格的离散性质，格点格林函数的计算通常需要用到矩阵运算的方法。

具体而言，设格点网格的大小为 $N\times N$，则格点格林函数 $G_{ii}$ 可以表示为矩阵形式：

$$ \textbf{G} = (\textbf{I} - \textbf{H})^{-1} $$

其中，$\textbf{I}$ 是 $N^2\times N^2$ 的单位矩阵，$\textbf{H}$ 是大小为 $N^2\times N^2$ 的对称矩阵，其元素 $H_{ij}$ 表示格点 $i$ 和格点 $j$ 之间的格点格林函数。具体而言，如果格点 $i$ 和格点 $j$ 相邻，则 $H_{ij}=-1$；否则 $H_{ij}=0$。注意到 $\textbf{H}$ 的对称性，因此 $\textbf{G}$ 也是对称矩阵。

在实际计算中，可以使用迭代法或者 LU 分解等方法求解矩阵 $(\textbf{I} - \textbf{H})$ 的逆。求解出 $\textbf{G}$ 后，格点格林函数 $G_{ii}$ 即为矩阵 $\textbf{G}$ 的第 $i$ 个对角线上的元素