libcint的使用求助

薛定谔的老猫 · 发表于 Post on 2019-7-25 23:21:13

sobereva 发表于 2019-7-25 23:16
PS：过几个月，两年前写了一半的《量子化学程序编写入门》会写完然后发布在网上。里面的程序就是基于libcin ...

哇咔咔，感觉运气很好啊，刚接触计算就上了初级班，基础班，想学量化程序编写，社长的大作就出来了，非常感谢！期待！！！

niobium · 发表于 Post on 2019-7-27 07:38:57

sobereva 发表于 2019-7-25 23:16
PS：过几个月，两年前写了一半的《量子化学程序编写入门》会写完然后发布在网上。里面的程序就是基于libcin ...

期待啊！

薛定谔的老猫 · 发表于 Post on 2019-9-11 21:00:12

sobereva 发表于 2019-7-25 23:16
PS：过几个月，两年前写了一半的《量子化学程序编写入门》会写完然后发布在网上。里面的程序就是基于libcin ...

社长大人，请问您的书有着落了么，啥时候发布呀，期待

sobereva · 发表于 Post on 2019-9-12 11:06:05

薛定谔的老猫发表于 2019-9-11 21:00
社长大人，请问您的书有着落了么，啥时候发布呀，期待

不知道，也许半年

薛定谔的老猫 · 发表于 Post on 2019-9-12 17:53:04

好的，多谢

ajz34 · 发表于 Post on 2019-9-22 21:27:10

本帖最后由 ajz34 于 2019-9-22 21:44 编辑

我现在在尝试用 Python 写后自洽场程序，以双杂化 hessian 为目标
我没试过 libcint，我猜用 libcint 的话还要处理 primitive 到 CGTO 的转换问题，这个挺烦的，特别角动量一高
只是代码的效率和内存优化我玩不来……

但是如果是为了学习后自洽场用，分子足够小的话，调 PySCF 的 gto 模块函数 (intor) 我认为足够了
我觉得一定程度上 PySCF 当成 libcint 的外包，想调什么积分问 PySCF 其实就行了，所以现在其实还完全不会 libcint……

以及重要的是活用 numpy.einsum(**kargs, optimize=True)，这个真的很强
以及 numpy 的 Boardcasting 功能
如果打算写一个 RMP2 相关能的程序，电子积分全部存在内存里
核心的 ERI AO->MO 和 RMP2 相关能量计算可以在 5 行代码搞定

现在也有一些英文教程，我还没仔细看过，但质量应该都不差，是 Psi4 系列的教程
C++: https://github.com/CrawfordGroup/ProgrammingProjects
Python + Jupyter: https://github.com/psi4/psi4numpy
以后多半也会学一下
(但 Jupyter 调用 Psi4 的程序，代码写得不好很容易 kernel died，于是就放弃 Psi4 投入几乎纯 Python 的 PySCF 了……)

ajz34 · 发表于 Post on 2019-9-22 21:58:29

说实话我自己最近也在尝试写教程文档，但现在处于未完成缓更状态233
https://py-xdh.readthedocs.io/zh_CN/latest/qcbasic/index.html
尽管现在的目标单纯地是用 Python 程序解决 Hessian 的计算，但写到现在为止的内容希望多少也能对大家有所帮助

薛定谔的老猫 · 发表于 Post on 2019-9-24 16:00:45

ajz34 发表于 2019-9-22 21:58
说实话我自己最近也在尝试写教程文档，但现在处于未完成缓更状态233
https://py-xdh.readthedocs.io/zh_CN ...

好的，感谢

wangxubo · 发表于 Post on 2019-9-25 07:14:11

就用pyscf.gto.moleintor()就非常方便，对于写个教学程序足够了，而且还能和pyscf对照，非常方便。

薛定谔的老猫 · 发表于 Post on 2019-9-25 11:09:00

wangxubo 发表于 2019-9-25 07:14
就用pyscf.gto.moleintor()就非常方便，对于写个教学程序足够了，而且还能和pyscf对照，非常方便。

感谢回复，好的

海棠 · 发表于 Post on 2022-6-29 13:28:24

niobium 发表于 2019-7-22 08:25
写了个简单的使用教程，希望对你有用

https://mp.weixin.qq.com/s?__biz=MzU5NjMxNjkzMw==&mid=22474845 ...

你好，自己写的Fortran程序通过什么编译选项可以链接libcint的库？

niobium · 发表于 Post on 2022-6-29 15:16:09

海棠发表于 2022-6-29 13:28
你好，自己写的Fortran程序通过什么编译选项可以链接libcint的库？

我是先把Libcint编译成了.so文件，然后在编译自己的程序的时候把这个库文件链接进来

zjxitcc · 发表于 Post on 2022-6-29 20:31:32

海棠发表于 2022-6-29 13:28
你好，自己写的Fortran程序通过什么编译选项可以链接libcint的库？

示例：ifort -c a.f90 -o aaa -L/home/libcint/lib -lcint

海棠 · 发表于 Post on 2022-7-2 12:41:36

zjxitcc 发表于 2022-6-29 20:31
示例：ifort -c a.f90 -o aaa -L/home/libcint/lib -lcint

谢谢。但我这边用-L选项不管用（我不太懂），只能把.so文件放到/usr/lib/里。

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

[算法与编程] libcint的使用求助

浏览过的版块