#include

此方法来自武剑的论文，见 Jian Wu and Yuansheng Jiang, J. Comp. Chem. 2000, 21, 856.       此处的main函数只是用于测试算法是否正确，此算法最后会作为另一个程序的子函数。/* luozhen, 2015-05-20*///Available only for Sz = 0#include #include

因为没有给定参数，所以参数t和u全都用字符串代替，而且之后的解矩阵运算也没办法进行。如果给定参数而不是用字符串表示，那么整个程序的效率会更高，而且之后可以继续写具体的矩阵运算过程。#include #include #include #include using namespace std;const int atoms = 4;const int configuratio

只求完工而不求代码效率，现在终于写完这个程序了。接下来就要做代码的优化和并行化了——事实上需要优化的内容简直太多。        最后求解本征能量与本征函数的时候使用了 intel Math Kernel Library 提供的函数，即 lapack 中的 LAPACKE_dsyevd 函数。编译的时候需要正确调用 intel MKL 才可以编译完成。/* luozhen, 2015-

对两个 xyz 结构文件的线性插值程序。这是计算化学领域优化几何构型时最常见的获取初始猜测的策略之一。

针对 64 位处理器的基于二进制形式的占据数矢量（行列式）以及多组态波函数的构建与操作。

量子化学哈密顿矩阵元的计算，基于上一篇文章构建的占据数矢量。

哈密顿矩阵元计算的效率问题和存储问题。