版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

yjcmwgk

禁蟲 (文壇精英)

密度泛函·小卒

QC強(qiáng)帖: 38
應(yīng)助: 165 (高中生)
貴賓: 16.999
金幣: 554.6
散金: 39880
紅花: 395
沙發(fā): 2
帖子: 10406
在線: 1528.1小時(shí)
蟲號(hào): 448368
注冊(cè): 2007-11-01
專業(yè): 理論和計(jì)算化學(xué)
管轄: 科研工具資源

[求助] 用Python讀Gaussian的cube文件，速度很慢，怎么破？謝謝�。�！已有4人參與

向各位python大俠求教一個(gè)問題，謝謝各位！
我正在用python讀一個(gè)由Gaussian生成的cube文件。但是讀文件讀得非常慢，讀一個(gè)6M多的cube文件，居然用了3分鐘！

具體問題如下：
用cubefile來打開某個(gè)cube文件。然后繼續(xù)

CODE:

    tmpcubedata = [] !先設(shè)一個(gè)空列表，待會(huì)兒用來存儲(chǔ)數(shù)據(jù)

    tmpcubedata = tmpcubedata + tmpline.split() !tmpline是用來存儲(chǔ)某一行的數(shù)據(jù)的。由于先前已經(jīng)經(jīng)過了多次tmpline.split()之后，實(shí)際上現(xiàn)在的tmpline存儲(chǔ)的就是我已經(jīng)想要的那些數(shù)據(jù)的第一行了。所以，先讀一下并存入tmpcubedata中。

    i = 1 !這是個(gè)計(jì)數(shù)器。每讀一行數(shù)據(jù)就加1

    tmplines=cubefile.readlines() !把剩下的行全讀進(jìn)來

    for i in range(1,linecount): !linecount就是應(yīng)該讀入的cube文件中數(shù)據(jù)部分的總行數(shù)，一般來說，Gaussian生成的cube文件中，如果生成過程全默認(rèn)，那么cube文件的數(shù)據(jù)部分應(yīng)該是8-9萬行的樣子。大部分行，每行有6個(gè)數(shù)字（大概95%的行是這樣的），但也不排除只有1-5個(gè)數(shù)字的情況（大概5%的行是這樣的）。

        tmpline = tmplines[i-1] !從tmplines中獲得某一行

        if tmpline == 'EOF' : break

        if tmpline == [] : break

        if tmpline == None : break

        if i == linecount: break !如果計(jì)數(shù)器i達(dá)到了總行數(shù)，就不讀了

        if i%10000 == 0: print "Converting", cubefilename, "......  ", int(100.0*float(i)/float(linecount)), "%" !每讀一萬行數(shù)據(jù)，就顯示一個(gè)進(jìn)度百分比

        if i == linecount-1 : print "Converting", cubefilename, "...... ", "100 %" !讀到最后一行數(shù)據(jù)的時(shí)候，就顯示一個(gè)百分之百

        i = i+1

        tmpcubedata = tmpcubedata + tmpline.split() !將讀進(jìn)來的數(shù)據(jù)存入tmpcubedata中

我懷疑tmpline = tmplines[i-1]這句話拖累的速度？應(yīng)該如何修改？非常感謝高手幫助��！
感謝��！

回復(fù)此樓

1樓 2014-10-04 11:04:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

rh1127

捐助貴賓 (著名寫手)

一名工程師@上海

應(yīng)助: 2 (幼兒園)
金幣: 21406.7
紅花: 5
帖子: 1616
在線: 139.5小時(shí)
蟲號(hào): 116848
注冊(cè): 2005-11-25
性別: GG
專業(yè): 膠體與界面化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
yjcmwgk: 金幣+20, ★有幫助, 非常感謝，我后來用tmpp = ((' '.join(cubefile.readlines())).replace("/n"," ")).split()解決了 2014-10-05 16:10:21
gmy1990: 金幣+2 2014-10-05 18:49:14

用numpy和pandas模塊兒來處理大型數(shù)組，比只用for循環(huán)速度快幾十倍不止。python自帶的循環(huán)未經(jīng)優(yōu)化，很忌諱來處理大型數(shù)組！
希望對(duì)你有幫助。

贊一下

回復(fù)此樓

Actionnoworachievenever!

4樓2014-10-05 09:21:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 10 個(gè)回答

yjcmwgk

禁蟲 (文壇精英)

密度泛函·小卒

QC強(qiáng)帖: 38
應(yīng)助: 165 (高中生)
貴賓: 16.999
金幣: 554.6
散金: 39880
紅花: 395
沙發(fā): 2
帖子: 10406
在線: 1528.1小時(shí)
蟲號(hào): 448368
注冊(cè): 2007-11-01
專業(yè): 理論和計(jì)算化學(xué)
管轄: 科研工具資源

cube文件的基本結(jié)構(gòu)舉例如下

CODE:

 Title Card Required mo=1

 MO coefficients

   -1   -2.916225   -2.916225   -2.916225

   80    0.073828    0.000000    0.000000

   80    0.000000    0.073828    0.000000

   80    0.000000    0.000000    0.073828

    2    2.000000    0.000000    0.000000    0.000000

    1    1

  7.36784E-05  8.36290E-05  9.46155E-05  1.06698E-04  1.19934E-04  1.34374E-04

  1.50064E-04  1.67044E-04  1.85341E-04  2.04976E-04  2.25955E-04  2.48274E-04

  2.71913E-04  2.96838E-04  3.22995E-04  3.50318E-04  3.78721E-04  4.08098E-04

  4.38328E-04  4.69271E-04  5.00768E-04  5.32647E-04  5.64718E-04  5.96779E-04

  6.28614E-04  6.60000E-04  6.90707E-04  7.20497E-04  7.49136E-04  7.76386E-04

  8.02019E-04  8.25811E-04  8.47551E-04  8.67043E-04  8.84107E-04  8.98583E-04

  9.10334E-04  9.19248E-04  9.25239E-04  9.28250E-04  9.28250E-04  9.25239E-04

  9.19248E-04  9.10334E-04  8.98583E-04  8.84107E-04  8.67043E-04  8.47551E-04

  8.25811E-04  8.02019E-04  7.76386E-04  7.49136E-04  7.20497E-04  6.90707E-04

  6.60000E-04  6.28614E-04  5.96779E-04  5.64718E-04  5.32647E-04  5.00768E-04

  4.69271E-04  4.38328E-04  4.08098E-04  3.78721E-04  3.50318E-04  3.22995E-04

  2.96838E-04  2.71913E-04  2.48274E-04  2.25955E-04  2.04976E-04  1.85341E-04

  1.67044E-04  1.50064E-04  1.34374E-04  1.19934E-04  1.06698E-04  9.46155E-05

  8.36290E-05  7.36784E-05

  8.36290E-05  9.49233E-05  1.07394E-04  1.21108E-04  1.36131E-04  1.52521E-04

  1.70331E-04  1.89603E-04  2.10372E-04  2.32658E-04  2.56471E-04  2.81805E-04

  3.08636E-04  3.36926E-04  3.66617E-04  3.97630E-04  4.29868E-04  4.63213E-04

下面還有八萬多行，略了

對(duì)于這個(gè)cube文件來說，從第9行開始是第一個(gè)數(shù)據(jù)行，直到文件結(jié)束。值得注意的是，絕大多數(shù)數(shù)據(jù)行，每行有6個(gè)數(shù)字，但第22行就只有2個(gè)數(shù)字了。

贊一下

回復(fù)此樓

2樓2014-10-04 11:08:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

frk2008

木蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 3434.4
紅花: 1
帖子: 489
在線: 440.9小時(shí)
蟲號(hào): 841330
注冊(cè): 2009-09-05
專業(yè): 電化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
yjcmwgk: 金幣+20, ★有幫助, 非常感謝，我后來用tmpp = ((' '.join(cubefile.readlines())).replace("/n"," ")).split()解決了 2014-10-05 16:10:10
gmy1990: 金幣+3 2014-10-05 18:49:03

試試這個(gè)，看看效果如何：

[ code]

tmplines = cubefile.readlines()
tempstr = "".join(tmplines)          #把list轉(zhuǎn)化為str，形式應(yīng)該是'num  num  ...num\nnum  ...'
tempstr = tempstr.replace('\n','  ') #把str中的'\n'替換為'  '
tmpcubedata = tempstr.split('  ')    #把str轉(zhuǎn)為list
n=0
numcount=len(tmpcubedata)             #list的數(shù)據(jù)總數(shù)
for i in tmpcubedata:
n+=1
if i%60000==0:
      print "Converting", cubefilename, "......  ", int(100.0*float(i)/float(numcount)), "%"
#if i=="":
      #tmpcubedata.remove(i)       #如果cubefile 文件中中間有空行，則需要這個(gè)功能
if tmpcubedata[-1]=="":
del tmpcubedata[-1]                #如果cubefile 文件以‘num\n‘結(jié)尾
print "Converting", cubefilename, "...... ", "100 %"

[ /code]

贊一下

回復(fù)此樓

3樓2014-10-04 15:23:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

magicmonk

至尊木蟲 (著名寫手)

應(yīng)助: 172 (高中生)
金幣: 11451.3
紅花: 24
帖子: 1124
在線: 3108.8小時(shí)
蟲號(hào): 1191315
注冊(cè): 2011-01-17
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

個(gè)人呢建議用fortran或MATLAB

贊一下

回復(fù)此樓

我們最深的恐懼不是我們能力不夠，我們最深的恐懼是能力超越了極限。

5樓2014-10-05 09:50:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 10 個(gè)回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 化學(xué)工程專碩324分，一志愿中國(guó)礦業(yè)大學(xué)求調(diào)劑 +6	耿耿1314 2026-04-01	6/300	2026-04-01 23:32 by chaolymer
[考研] 311求調(diào)劑 +13	zchqwer 2026-04-01	15/750	2026-04-01 22:30 by 2026材料調(diào)劑
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-30	16/800	2026-04-01 17:58 by 邱gl
[考研] 362求調(diào)劑 +13	西南交材料專碩3 2026-03-31	13/650	2026-04-01 17:38 by JYD2011
[考研] 311求調(diào)劑 +11	藍(lán)月亮亮 2026-03-30	11/550	2026-04-01 16:33 by 七度不信任
[考研] 調(diào)劑推薦 +9	清酒714 2026-03-26	10/500	2026-04-01 16:25 by yanflower7133
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +6	篤行致遠(yuǎn). 2026-03-28	6/300	2026-04-01 09:28 by oooqiao
[考研] 080200學(xué)碩，機(jī)械工程專業(yè)277分，求帶走！ +4	瓶子PZ 2026-03-31	4/200	2026-03-31 20:16 by vgtyfty
[考研] 299求調(diào)劑 +8	嗯嗯嗯嗯2 2026-03-27	8/400	2026-03-31 18:23 by lizhi8172
[考研] 343求調(diào)劑 +8	愛羈絆 2026-03-28	8/400	2026-03-31 16:12 by 不吃魚的貓
[考研] 085601 329分調(diào)劑 +6	yzsa12 2026-03-31	6/300	2026-03-31 15:23 by yanflower7133
[考研] 282求調(diào)劑不挑專業(yè) 求收留 +4	Yam. 2026-03-30	5/250	2026-03-31 14:41 by 王亮_大連醫(yī)科大
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-29	3/150	2026-03-31 10:06 by cal0306
[考研] 323分食品與營(yíng)養(yǎng)調(diào)劑 +3	嘿ooo 2026-03-31	3/150	2026-03-31 09:38 by longlotian
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3	塔比烏斯 2026-03-30	3/150	2026-03-30 22:55 by 無際的草原
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +5	天賦帶到THU 2026-03-29	6/300	2026-03-30 20:49 by dick_runner
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3	Lyy930824@ 2026-03-29	3/150	2026-03-30 20:15 by DHUSHUAI
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

yjcmwgk

[求助] 用Python讀Gaussian的cube文件，速度很慢，怎么破？謝謝�。�！ 已有4人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

rh1127

【答案】應(yīng)助回帖

yjcmwgk

frk2008

【答案】應(yīng)助回帖

magicmonk

【答案】應(yīng)助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 用Python讀Gaussian的cube文件，速度很慢，怎么破？謝謝�。�！已有4人參與

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: