本版本的代码专注于能量流的计算,不涉及任何热力学量。实现过程较为缓慢,教育中理解NRG方法。例如,代码未考虑状态间仅差一个占据数或S_z差1/2的矩阵元素,需对耗时部分进行矢量化以实现隐式并行化。