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

版塊導航: 正在加載中...

客戶端APP下載

調(diào)劑小程序

登錄注冊

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

24小時熱門版塊排行榜

hakuna

木蟲 (知名作家)

1ST強帖: 3
應助: 82 (初中生)
金幣: 8735.9
散金: 158
紅花: 126
帖子: 5599
在線: 1367.2小時
蟲號: 167622
注冊: 2006-01-15
專業(yè): 考古理論

[交流] 【分享】一個繪制band的python小程序已有4人參與

網(wǎng)上淘來的，向作者致敬！
不做能帶，不懂python，未經(jīng)測試

原文地址：http://people.vanderbilt.edu/~brandon.g.cook/vaspbandplot/vband.py

CODE:

!/usr/bin/env python

#make a simple plot of band structure from vasp calculation

# -Brandon Cook April 7, 2010

from numpy import *

import pylab

class OUTCAR():

    def __init__(self):

        f = open('OUTCAR', 'r')

        s = f.read()

        f.close()

        i = s.find('E-fermi')

        self.ef = float(s[i:i+19].split()[-1])

        i = s.find('NBANDS')

        self.nbands = int(s[i:i+15].split()[-1])

        i = s.find('NKPTS')

        self.nkpts = int(s[i:i+15].split()[-1])

class KPOINTS():

    def __init__(self):

        # kpts file should be in line mode

        # only take 1 line for now

        f = open('KPOINTS', 'r')

        s = f.readlines()

        f.close()

        l1 = s[-2].split()

        l2 = s[-1].split()

        self.p1 = array([ l1[0], l1[1], l1[2] ], float)

        self.p2 = array([ l2[0], l2[1], l2[2] ], float)

class BandStructure():

    # Just take all the points and make a scatter plot

    # todo: connect the lines

    def __init__(self):

        f = open('EIGENVAL', 'r')

        s = f.readlines()

        f.close()

        oc = OUTCAR()

        kp = KPOINTS()

        axis = kp.p2 - kp.p1

        axis /= linalg.norm(axis)

        nbands = oc.nbands

        nkpts  = oc.nkpts

        kpts = zeros( nkpts*nbands )

        en   = zeros( nkpts*nbands )

        istart = 7

        for k in range(0, nkpts):

            ik = istart + k*(nbands+2)

            kpt = s[ik].split()

            kpt = array([kpt[0], kpt[1], kpt[2]], float)

            x = dot(kpt, axis)

            i1 = k*nbands

            i2 = (k+1)*nbands

            kpts[i1:i2] = x

            for ib in range(0, nbands):

                j1 = ik + 1 + ib

                j2 = i1 + ib

                en[j2] = float(s[j1].split()[-1])

        pylab.scatter(kpts, en, s=1)

        # add a line for the fermi energy

        xmin = 0.0

        xmax = linalg.norm(kp.p2 - kp.p1)

        x = linspace(xmin, xmax, num=100, endpoint=True)

        y = oc.ef * ones(100)

        pylab.plot(x,y)

        pylab.show()

o = OUTCAR()

print 'fermi energy, nbands, nkpts'

print o.ef,o.nbands,o.nkpts

b = BandStructure()

» 猜你喜歡

北科281學碩材料求調(diào)劑已經(jīng)有13人回復
289求調(diào)劑已經(jīng)有11人回復
網(wǎng)絡空間安全0839招調(diào)劑已經(jīng)有5人回復
材料277求調(diào)劑已經(jīng)有4人回復
274求調(diào)劑已經(jīng)有7人回復
085600 材料與化工 329分求調(diào)劑已經(jīng)有3人回復
302求調(diào)劑已經(jīng)有3人回復
312求調(diào)劑已經(jīng)有5人回復
329求調(diào)劑已經(jīng)有3人回復
【2026考研調(diào)劑】制藥工程 284分求相關專業(yè)調(diào)劑名額已經(jīng)有8人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

推薦幾個vasp處理計算結果的python腳本已經(jīng)有66人回復
Python寫的arcgis中點坐標txt轉面狀shp的小程序已經(jīng)有90人回復
求一個簡單的讀寫文件的小程序。。。已經(jīng)有4人回復
用phonopy計算聲子譜已得到FORCE_SETS，之后我怎么得到聲子譜？已經(jīng)有7人回復
《python科學計算》已經(jīng)出版了已經(jīng)有36人回復
運行一個python程序，出現(xiàn)問題，請高手賜教已經(jīng)有9人回復
【原創(chuàng)】處理elk能帶的程序已經(jīng)有16人回復
【求助】求助安裝abinitbandstructuremaker.py? 已經(jīng)有3人回復

1樓 2010-12-26 17:32:54

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

sg18408926

至尊木蟲 (著名寫手)

應助: 0 (幼兒園)
金幣: 10945.9
散金: 246
紅花: 1
帖子: 1075
在線: 431.4小時
蟲號: 727956
注冊: 2009-03-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

不錯啊

2樓2010-12-26 20:43:52

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

gavinliu7390

木蟲 (著名寫手)

葉落鷹飛

1ST強帖: 3
應助: 51 (初中生)
貴賓: 0.481
金幣: 1564.7
散金: 2975
紅花: 34
帖子: 1278
在線: 241.4小時
蟲號: 675997
注冊: 2008-12-17
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

執(zhí)行不了，錯誤很多。

贊一下(1人)

真理是一點點接近的！

3樓2010-12-27 10:44:18

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

westmonster

金蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 5637.3
帖子: 401
在線: 236.4小時
蟲號: 456224
注冊: 2007-11-11
專業(yè): 凝聚態(tài)物性 II ：電子結構

pylab是嘛？孤陋寡聞了，改改還是能用的，謝謝樓主分享。

蟬，于其脫殼展翅，則蟄居地下，似無聲無息，實則以備有所為。

4樓2010-12-27 12:54:53

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

akakcolin

金蟲 (著名寫手)

應助: 3 (幼兒園)
金幣: 685.9
散金: 51
紅花: 2
帖子: 1010
在線: 524.6小時
蟲號: 321213
注冊: 2007-03-10
專業(yè): 力學

python 是個好語言啊

為人應當誠實正直不能心懷叵測心眼明亮才能迎來幸福避開災禍盲目者若有人指點才不會迷失道路若單獨上路就難以保證不誤入岐途

5樓2010-12-27 23:09:27

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

相關版塊跳轉我要訂閱樓主 hakuna 的主題更新

普通表情龍兔虎貓高級回復 (可上傳附件)

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

	[考研] 085602 289分求調(diào)劑 +7	WWW西西弗斯 2026-03-24	7/350	2026-03-25 14:28 by 3Strings
	[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
	[考研] 282求調(diào)劑 +3	wcq131415 2026-03-24	3/150	2026-03-25 12:16 by userper
	[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper
	[考研] 085600材料與化工調(diào)劑 +9	A-哆啦Z夢 2026-03-23	15/750	2026-03-25 11:18 by Ainin_
	[考研] 299求調(diào)劑 +7	shxchem 2026-03-20	9/450	2026-03-25 10:41 by lbsjt
	[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-
	[考研] 085601求調(diào)劑總分293英一數(shù)二 +3	鋼鐵大炮 2026-03-24	3/150	2026-03-24 22:03 by bingxueer79
	[考研] 招08考數(shù)學 +7	laoshidan 2026-03-20	16/800	2026-03-24 17:52 by 烏拉兒山脈
	[考研] 0854 考研調(diào)劑招生了！AI 方向 +5	pk3725069 2026-03-19	17/850	2026-03-24 17:30 by zhouxuan..
	[考研] 一志愿華東理工大學081700，初試分數(shù)271 +5	kotoko_ik 2026-03-23	6/300	2026-03-24 10:29 by 學術搬磚er
	[考研] 341求調(diào)劑(一志愿湖南大學070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
	[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
	[考研] 化學308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
	[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5	jiajunX 2026-03-22	5/250	2026-03-23 18:07 by YMU施老師
	[考研] 308求調(diào)劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler
	[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
	[考研] 0703化學297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
	[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
	[考研] 材料學碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad