应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 专家会诊 » 分享一个计算孔径分布的程序同时求助该程序使用时碰到的问题
101/1返回列表
查看: 2184  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 shiyiwenren 的 18 个金币

shiyiwenren

至尊木虫 (著名写手)

[交流] 分享一个计算孔径分布的程序同时求助该程序使用时碰到的问题

以前做多孔材料的模拟的时候，一直想看下材料的孔径分布，但是不知道怎么进行这方面的模拟，今天在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 header for the header for C++ includes, or instead of the deprecated header . To disable this warning use -Wno-deprecated.
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
    • 2013-04-27 16:42:27, 13.99 K
  • 附件 2 : PSDsolvMPI.tar.gz
    • 2013-04-27 16:42:35, 14.72 K

» 本帖已获得的红花（最新10朵）

安德

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2013-04-27 16:57:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gzl9901

铁杆木虫 (文坛精英)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
shiyiwenren: 金币+1 2013-05-02 08:33:58
祝顺利，科研加油！
赞一下 回复此楼
2楼2013-04-27 17:50:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

安德

铁杆木虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
shiyiwenren: 金币+1 2013-05-02 08:34:02
虽然不是很容易看明白，但是为楼主这种认真的精神感动的泪流满面，送你一朵红花聊表我的敬意！
赞一下 回复此楼
3楼2013-04-27 21:13:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dying999

铁虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我也遇到楼主一样的问题。在网上查找过，找到了三种可以可以孔径计算的方法，目前好像只有楼主发的这种有现成软件，但是安装遇到一样的问题，完全不懂C++，不知道楼主找到了办法了吗
赞一下 回复此楼
6楼2015-03-11 12:39:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niujingxin

新虫 (初入文坛)
不错的
回复此楼
7楼2016-04-02 23:23:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clumsyinsect

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我想问一下，楼主，这个是可以作为ms的脚本使用吗？

发自小木虫Android客户端
赞一下 回复此楼
8楼2016-10-16 07:58:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

压路鸡

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
6楼: Originally posted by dying999 at 2015-03-11 12:39:43
我也遇到楼主一样的问题。在网上查找过，找到了三种可以可以孔径计算的方法，目前好像只有楼主发的这种有现成软件，但是安装遇到一样的问题，完全不懂C++，不知道楼主找到了办法了吗

求前辈分享经验
赞一下 回复此楼
9楼2017-03-31 11:02:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小猫真可爱

禁虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
10楼2019-06-10 19:59:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
亦岳98634楼
2013-04-27 21:32   回复  
rkl1135楼
2013-04-27 21:35   回复  
相关版块跳转 我要订阅楼主 shiyiwenren 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 关于论文第一通讯作者的问题？ +4 winsaint 2024-06-26 4/200 2024-06-27 11:12 by gancs7894
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +7 来了个 2024-06-26 7/350 2024-06-27 11:11 by wubiansiya
[硕博家园] 联培博士文章第一单位署名问题交流~ +10 橙成成c 2024-06-23 23/1150 2024-06-27 10:11 by 橙成成c
[基金申请] 国基在研影响申请结果吗 +13 WOWO159357 2024-06-26 21/1050 2024-06-27 09:44 by WOWO159357
[育儿交流] 暑假孩子一个人在家可以干什么 +8 唠叨的映 2024-06-25 9/450 2024-06-27 09:28 by 晓目崇
[有机交流] 做什么表征可以检测塑料中碳氢氧的含量 5+3 pzr的sci之路 2024-06-25 7/350 2024-06-26 19:32 by wrgeng
[硕博家园] 回家两天，不想呆了 +7 368ghnf 2024-06-22 13/650 2024-06-26 19:07 by ZZZemmm
[基金申请] 省基金对国自然申请是有利还是有弊 +6 zycumt04 2024-06-26 7/350 2024-06-26 15:02 by 潇湘之迷
[基金申请] 基金申请书名称有变化 +5 xuel2011 2024-06-25 7/350 2024-06-26 00:07 by 老虎当猫养
[基金申请] 博后面上到底什么时候出结果？？？ +7 爱学lsy 2024-06-24 7/350 2024-06-25 19:24 by 暴走的蒜泥
[有机交流] 对苯乙烯磺酰氯的合成机理 25+3 该死的科研 2024-06-24 5/250 2024-06-25 17:30 by 王学士
[硕博家园] 数据不好 +5 Hetai 2024-06-23 7/350 2024-06-25 12:37 by 1591099
[基金申请] 演员杨幂在核心期刊发表论文，中国广播电视学期刊 +9 babu2015 2024-06-22 12/600 2024-06-25 11:05 by 我是勇敢牛牛
[教师之家] 复旦夏同学提出高校成年人学生退学不应该让家长审核，大家认同吗？ +10 苏东坡二世 2024-06-22 17/850 2024-06-24 16:52 by wanghuawei
[基金申请] 国自然资助比率是不是要下降了？？ +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[基金申请] 青年和面上，哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[有机交流] 生成亚胺的反应怎么能进行完全 +3 1369836 2024-06-23 3/150 2024-06-23 18:44 by hwqMSE
[论文投稿] OSA期刊审稿逾期 +3 Thomas_Squid 2024-06-22 3/150 2024-06-23 15:20 by wspglt
[基金申请] 工材口青年基金大概什么样能上会？ +15 今晚推荐22 2024-06-20 21/1050 2024-06-22 23:04 by qbn0326
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
信息提示
关闭
请填处理意见
关闭