24小时热门版块排行榜    

CyRhmU.jpeg
查看: 143  |  回复: 2
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 casjxm 的 5 个金币

casjxm

铜虫 (正式写手)

[交流] 【求助】做磁性计算时 gz2.py遇到的问题

在ABINIT的spin lesson of the tutorial里面,有一处要用到python gz2.py计算两个Fe原子的磁性,不知道我的为什么在计算时提示出错
Number of atoms =  2
Atomic coordinates
[['0.000000E+00', '0.000000E+00', '0.000000E+00'], ['1.650000E+00', '-1.650000E+00', '3.300000E+00']]
number of integration points: 32
Treating  atom # 0
Treating  atom # 1
0 0
Traceback (most recent call last):
  File "gz2.py", line 58, in ?
    integral.append(sum/npts_integral)
ZeroDivisionError: integer division or modulo by zero

我看了下源代码:
很可能是if(string.find(line,"Calculated">0)没有通过导致npts_integral为初始值0,故有上面的错误;不知道那位高人能解决这个问题,得到正确的结果。

源代码:http://www.abinit.org/documentat ... /lesson_spin/gz2.py
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenxianshifu

至尊木虫 (知名作家)

宗师

帮楼主顶一下,希望高手看到
2楼2009-09-06 21:51:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenjlu

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
freshgirl(金币+1,VIP+0):谢谢参与~ 9-16 22:05
casjxm(金币+10,VIP+0): 9-27 08:45
casjxm(金币+5,VIP+0): 10-16 14:46
我有修改过源代码的程序,并且已编译成windows版本了。
下载地址:
http://d.namipan.com/d/969e6f24b ... dba975a50d2f87a3500
3楼2009-09-16 16:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 casjxm 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见