用Fortran 90改写Modern Quantum Chemistry附录B中的HF自洽场源码

zhaixiaoyi001 · 发表于 Post on 2026-1-29 08:47:37

本帖最后由 zhaixiaoyi001 于 2026-1-30 09:20 编辑

作为练习，用Fortran 90改写了HeH+自洽场代码，主要改动如下：
1、去掉了跳转，改写了do循环和if语句。
2、保留了原代码中的变量和电子积分计算方式。
3、用lapack做对角化，有部分变量是冗余的没有删除。

SeptemberMy · 发表于 Post on 2026-1-29 19:39:15

如果你愿意的话可以把这个程序添加到一个Modern Quantum Chemistry汉化的项目中，直接提PR即可，https://github.com/nominhanggai/szaboqc

slxc920113 · 发表于 Post on 2026-1-29 23:46:47

代码的327行有语法错误，编译后运行会报错，改成

write(*,'(/,4X,A,F10.6,/)') 'DELTA(CONVERGENCE OF DENSITY MATRIX) = ', delta

复制代码

我重写了一个python的版本，逻辑和结构上会清晰很多。

HeH .py (6.41 KB, 下载次数 Times of downloads: 7)

zhaixiaoyi001 · 发表于 Post on 2026-1-30 08:46:00

slxc920113 发表于 2026-1-29 23:46
代码的327行有语法错误，编译后运行会报错，改成

应该是那个39H字符串报错了，这块的Fortran 77语法我没修改，有的编译器识别不了，之前用ifx和gfortran编译倒是没报错。我再修改一下

zhaixiaoyi001 · 发表于 Post on 2026-1-30 10:06:38

SeptemberMy 发表于 2026-1-29 19:39
如果你愿意的话可以把这个程序添加到一个Modern Quantum Chemistry汉化的项目中，直接提PR即可，https://gi ...

又修改了一下代码，已提PR。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[算法与编程] 用Fortran 90改写Modern Quantum Chemistry附录B中的HF自洽场源码

评分 Rate

浏览过的版块