| 查看: 2188 | 回复: 9 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 shiyiwenren 的 18 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
分享一个计算孔径分布的程序同时求助该程序使用时碰到的问题
|
|||
|
以前做多孔材料的模拟的时候，一直想看下材料的孔径分布，但是不知道怎么进行这方面的模拟，今天在google上搜索pore size distribution的时候，惊喜地发现有人提供了这方面计算的程序，网址如下： http://supriyo.net/research/psd/psd.htm 我将上面的程序下载了下来，并上传到小木虫上，大家有兴趣可以看下，相互交流一下。网上提供了两个版本的程序：普通的和可以并行的。 我将这两个程序都尝试安装到自己组里的机群上了，但是在用的时候发现有很多问题。 首先我安装了普通版本的，我将PSDsolv.tar.gz文件解压到自己新建的目录下面，然后按照readme里的说明执行操作的时候，在第一步就碰到了错误，readme文件内容如下： ******************* PSDsolv v.1.0 ********************** Author of psd.cpp: Supriyo Bhattacharya solvopt.cpp and solvopt.h were included upon permission from Prof. Franz Kappel. ran_nr.cpp is taken from "Numerical Recipes in C++" Feel free to modify / distribute the original code for nonprofit / academic purposes. I give permission for commercial use of my own code, however solvopt is forbidden from being used commercially. Therefore commercial users can use the code as long as they provide their own nonlinear optimization routine and give proper credits to the author and cite the reference given at the end of this document. Also, commercial establishments may not prohibit others from modifying / distributing the original codes either for nonprofit or for commercial purposes. GENERAL DESCRIPTION This program calculates the pore size distribution from a structure file. The program does not need the connect records, only the locations and atom types, so the xyz format is the simplest to work with. If your structure is in a pdb or mol2 format, use babel to convert: i.e. openbabel -ipdb struc.pdb -oxyz struc.xyz The atoms are modeled as hard spheres (vdW radii). The program uses a nonlinear optimization routine (SOLVOPT) INSTALL Run the install script (./install) Then place the executable (psd) somewhere in your path or create a link in /usr/local/bin. Place radii_list.dat in the same directory where you run the program. This file contains the vdW radii of the common atom types used by the program. New user defined atom types can be added to radii_list.dat. RUN psd A bin width of 0.25A seems to be the optimum. The output produces three files: psd_diff.dat (the actual distribution), psd_cumm.dat (cummulative PSD, see paper for explanation) and error.dat (shows the average error over time / MC steps). Once you start the program, a rough profile can be seen within an hour. However, it may take upto 24 hours to get a smooth profile (depending on the system size). Please continue to monitor the distribution and the estimated error. Once you are satisfied with the results, just kill the process, simple!! OTHER INFO Please cite the following: psd calculation: S. Bhattacharya and K. E. Gubbins, fast Method for Computing Pore Size Distributions of Model Materials, Langmuir 2006, 22, 7726. solvopt: F. Kappel and A. V. Kuntsevich, Comput. Optim. Appl. 2000, 15, 193. Send all comments to sup27606@yahoo.com --------------------------------------------------------------------------------------------------- 我执行了./install后，跳出如下说明： cc1plus: error: unrecognized command line option "-fast" cc1plus: error: unrecognized command line option "-fast" 我对install进行了编译，删除了原来里面包含的-fast后，再重新执行./install后，跳出如下说明： solvopt.cpp: In function 鈊ouble solvopt(short unsigned int, double*, double (*)(double*), double (*)(double*, double*), double*, double (*)(double*), double (*)(double*, double*))? solvopt.cpp:186: warning: converting to 鈙hort unsigned int?from 鈊ouble? solvopt.cpp:205: warning: converting to 鈙hort int?from 鈊ouble? In file included from /usr/lib/gcc/x86_64-redhat-linux/4.1.2/../../../../include/c++/4.1.2/backward/iostream.h:31, from psd.cpp:1: /usr/lib/gcc/x86_64-redhat-linux/4.1.2/../../../../include/c++/4.1.2/backward/backward_warning.h:32:2: warning: #warning This file includes at least one deprecated or antiquated header. Please consider using one of the 32 headers found in section 17.4.1.2 of the C++ standard. Examples include substituting the In file included from psd.cpp:7: ran_nr.cpp:58:2: warning: no newline at end of file gettime.c: In function 鈏nt gettime(double*)? gettime.c:8: error: aggregate 鈚imezone tz?has incomplete type and cannot be defined gettime.c:11: error: 鈍ettimeofday?was not declared in this scope psd.cpp: In function 鈒ong int nearest(double)? psd.cpp:17: warning: converting to 鈒ong int?from 鈊ouble 所以程序在进行install的时候就执行不下去了，也就没有能用这个程序进行测试。我用C++对程序里包含的.cpp和.c文件进行了compile，发现都有提示程序有错误。但是因为我自己不会C++编程，所以也不知道怎么解决。希望有用过这个程序的人能给与帮助吧。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : PSDsolv.tar.gz
- 附件 2 : PSDsolvMPI.tar.gz
2013-04-27 16:42:27, 13.99 K
2013-04-27 16:42:35, 14.72 K
» 本帖已获得的红花（最新10朵）
安德» 猜你喜欢
科研之友最近访客那出现来自北京的未知用户，什么情况？
已经有20人回复
-
75批博士后面上项目公布了
已经有5人回复
-
名单中没有，评审结果反馈也没有信息，有没有相同的小伙伴指教
已经有6人回复
-
放榜是先短信还是先系统？
已经有9人回复
-
提交了只要等着就好，这样的评审会有吗？
已经有4人回复
-
博后基金意难平……
已经有20人回复
-
2024国自然会评日期
已经有4人回复
-
no news is good news
已经有9人回复
-
两个月二审还没完成
已经有7人回复
-
收到科研之友的邮件
已经有18人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助氮气吸脱曲线给出的t-plot数据怎么解释？
已经有11人回复
- 求助：氮气吸脱附曲线及BJH孔径分布的的作图（急）
已经有27人回复
- 求助，氮气吸附脱附曲线及孔径分布图
已经有5人回复
- BET孔径分布图求助
已经有22人回复
- 【求助】孔径分布图 纵坐标代表什么意思 急！！！！！
已经有22人回复
- 【求助】有关BET孔径分布
已经有6人回复
- 【求助】请教关于氮气吸附等温线的问题？
已经有13人回复
- 【求助】孔径分布解析
已经有11人回复
- 【求助】氮吸附吸附线与孔径大小的关系
已经有16人回复
- 【求助】ASAP2020做孔径分布
已经有7人回复
- 【求助】请教一下美国康塔或麦克的物理化学吸附仪什么型号比较好，价格怎么样？
已经有19人回复
- 【求助】N2-bet的结果中，用哪个数据做孔径分布图呀
已经有8人回复
- 【请教】求助有关BET的问题
已经有8人回复
- 【求助】比表面积，孔容，孔径分布求助 50金币
已经有26人回复
» 抢金币啦！回帖就可以得到:
-
#重庆本地女生，诚征男友# 最后一次发帖了
+5/232
-
博后基金散金求中
+5/112
-
2024暑期夏令营招生【西北工业大学翱翔领军团队倪骏康特任研究员招生】
+1/85
-
浙大博后招聘、访学招聘---课题方向不限
+1/84
-
牧原实验室（郑州）招聘研究人员
+2/44
-
南京大学现代工程与应用科学学院“有机半导体材料及其器件应用”方向招聘博士后
+1/36
-
哈工大（深圳）材料学院张嘉恒课题组（青千）—— 招聘科研助理
+1/33
-
中国农业科学院研究团队招收机器学习算法人工智能数据分析等专业暑期兼职
+1/33
-
厦门市某高中急需招聘高中老师若干名待遇丰厚
+1/32
-
大湾区大学（筹）张国强课题组招聘光催化方向博士后及特聘研究员
+1/30
-
省基金结题问题
+1/25
-
（2024年7月5日截止报名）深圳大学大湾区国际创新学院“2025年全国优秀大学生夏令营”
+1/11
-
中科院青岛生物能源与过程研究所氢能装备先进制造中心诚招科研助理，有读博意愿者优先
+1/9
-
中国科学院重庆绿色智能技术研究院招聘科研助理（无机非金属复合材料）
+1/7
-
Postdoc Positions in the University College Dublin
+1/7
-
西北工业大学 力学与土木建筑学院 断裂力学计算仿真方向 招收硕士研究生
+1/5
-
申请2025年博士，可提前开展实验（或进组做科研助理）
+1/3
-
求大佬告知，怎么通过cif看晶体是否为同构，需要看哪些参数？
+1/2
-
知名外资仪器公司急招Internal Sales Engineer
+1/1
-
安利一个老师
+1/1
+5/232
6楼2015-03-11 12:39:43
祝顺利，科研加油！
2楼2013-04-27 17:50:16
3楼2013-04-27 21:13:41
8楼2016-10-16 07:58:54
