过冷水最近一段时间接触了一个非常牛X的程序包——parsec,花了九牛二虎之力才弄懂它,它能干什么呢?
使用局域密度近似(LDA)的实空间网格电子结构计算代码, 计算使用赝势. 代码可以用含时密度泛函理论计算局部体系的光学特性. 对于团簇, 可以进行从头分子动力学计算. 代码目前还不支持并行计算.
懂的人自然懂,不懂的人吗!*看星星一灿明。我自然会讲一点你能听懂的知识,过冷水从网上下载下来这个程序包,程序包的使用很麻烦,对编程不是很精通的根本无法成功启动该程序包,本着独乐乐不如众乐乐的精神给有缘人分享一下正确使用该程序包的方法。程序包的原下载地址如上,不过过冷书运行不了这个程序包,对里面的相关代码有做更改,过冷书的程序包分享文末有附。
parsec程序包是fortran语言基于Linux写的,Windows系统运行肯定是需要移植的。所以解决办法是:
1.安装cygwin,并安装以下包:
gcc、gcc-g 、gcc-gfortran、make、libopenblas、liblapack-devel、libarpack-devel、libfftw3-devel
如下图所示:
2.将源码包parsec解压到cygwin中的home目录下.
3.打开cygwin,并解压
4.进入源码下的src文件夹
5.打开config/make.ices文件
6.将原有的ifort编译器改成gfortran,因为intel的编译器要钱,使用免费的
修改前
修改后
7.保存并退出该文件,并输入make命令则开始编译(-j2是使用两个进程同时编译)
8.编译完成后,找到parsec.ser这个文件,这是一个可执行文件
9.提前将 parse.dat、parsec.in这两个文件放在/parsec.se所在的文件夹下,输入./parsec.ser便能运行。
10. 运行结束后会生成两个输出文件parsec.out、wfn.dat,该程序运行结束。得到我们想要的数据。
我们得到parsec.out、wfn.dat两个文件后就万事大吉了吗?不!
这样的文件正常人能够看的懂吗?How do?我们要进行重新编码,对其可视化。可视化这部分涉及到python和matlab的使用,你既不会python也不会matlab那还玩毛线?回家洗洗睡吧!可见公众号平常推广的python课程,过冷水分享Matlab知识有多重要。PVOX可视化工具包文末附上。
1.解压PVOX,打开pvox.m文件如图:
2.运行程序导入parsec.out、wfn.dat文件
3.如果没有出错的话,你应该会看到这幅图:
4.但实际过冷水大半个月看到的都是:
Why?因为2.4.1版本的python码、7.0 的MATLAB 码的不兼容以及其它错误坑死人。大家在下载网上共享资源时,很多时候的资源都是有缺陷或者是需要灵活处理的,在灵活处理过程中就是掌握知识的过程中,接下来的很长一段时间,过冷时将会和大家一起学习如何将Error 转变成Success!。这就是我之前说的Big project。
原创文章,过冷水,未经授权,禁止私自转载,转载请联系作者,如果您希望加入Matlab仿真秀官方交流群和资料下载可加群:927550334