昨晚深夜檢查了一下計(jì)算任務(wù)的結(jié)果胰舆，然后想做兩件事：
1骚露、畫出對(duì)數(shù)坐標(biāo)的彩圖蹬挤，如此頁中上圖所示http://matplotlib.org/examples/pylab_examples/pcolor_log.html?highlight=lognorm
2缚窿、換不同的K點(diǎn)設(shè)置檢查透射系數(shù)和電導(dǎo)。
結(jié)果快氣死我了：
1焰扳、matplotlib的一些模板我這邊的電腦上貌似沒有倦零，例如Times New Roman字體竟然沒有，這也算了吨悍，能忍扫茅。最不能忍的是我把LogNorm和ticker.LogLocator都試了，結(jié)果畫出來的圖的colorbar變成一格一格的育瓜，例如說原來的colormap是60個(gè)顏色葫隙，用了對(duì)數(shù)坐標(biāo)后只有6個(gè)顏色，畫出來的圖就變成一塊一塊的躏仇。像http://matplotlib.org/examples/pylab_examples/contourf_log.html?highlight=lognorm
而且colorbar上的標(biāo)簽沒了恋脚。折騰了很久。最后感覺應(yīng)該是在contourf的levels設(shè)計(jì)的問題焰手。但是試了半天色塊變得更少更大了糟描。似乎可以用numpy.logspace，但是看它的說明半天進(jìn)不了腦子书妻。不過最不爽的是人家用的也是默認(rèn)參數(shù)來畫圖船响，為什么色塊過渡就這么自然。
而且即使解決了分塊的問題，colorbar仍然不顯示標(biāo)簽见间。
花了幾小時(shí)在這種無聊斃了的問題上聊闯，搞得整個(gè)人都煩躁了。不用內(nèi)置的log屬性米诉，自己寫一個(gè)腳本來轉(zhuǎn)換畫圖的值馅袁，倒也可以，但是我對(duì)matplotlib的語法不熟荒辕，主要是不知怎么把colorbar上對(duì)數(shù)化后的值換回原值汗销，又得查半天。
只好用回線性坐標(biāo)了抵窒。浪費(fèi)了四小時(shí)左右在上面弛针，什么都沒做出來。
這時(shí)看到lx師兄N年前的文章里畫圖用的就是對(duì)數(shù)坐標(biāo)李皇。削茁。唉，要是師兄還在就好了掉房。不過就算他在我估計(jì)也不想or不敢問他茧跋。
（現(xiàn)在頭腦清醒時(shí)再去看了一下logspace卓囚，一下子就看懂了瘾杭。。真是得好好休息才行了哪亿。）
2粥烁、寫了個(gè)批量處理的腳本來運(yùn)行atk，結(jié)果一提交等一會(huì)就失敗蝇棉，說license版本對(duì)不上讨阻。我很奇怪，因?yàn)檫@個(gè)腳本我不大可能對(duì)之改動(dòng)篡殷，既然是沿用的钝吮，為什么以前行，現(xiàn)在不行呢板辽。因?yàn)閳?jiān)信這個(gè)腳本和能成功運(yùn)行的腳本是一樣的奇瘦，我竟然沒去看關(guān)鍵的版本號(hào)。檢查半天不知問題在哪戳气。當(dāng)時(shí)也實(shí)在太困了链患，坐在電腦前都快要睡著了，來回看著同幾行腳本一個(gè)字都沒看進(jìn)去瓶您。最后受不了麻捻，只好降低腳本的自動(dòng)化程度纲仍，手動(dòng)處理了一下，讓服務(wù)器先跑一部分贸毕，然后就去睡了郑叠。
今天睡了一覺腦袋清醒點(diǎn)再看，檢查了一下明棍，發(fā)現(xiàn)版本號(hào)竟然真的改動(dòng)過乡革。。摊腋。沸版。。兴蒸。视粮。我死活想不通是什么時(shí)候改的。我為什么會(huì)把一個(gè)不能用的腳本放到backup下橙凳。蕾殴。。岛啸。钓觉。。坚踩。荡灾。完全一頭霧水。
另外還有個(gè)bug很annoying堕虹。mendeley的merge citation失效率高達(dá)90%卧晓。選中想merge的citation后，這個(gè)merge citation按鈕基本不會(huì)出現(xiàn)赴捞，你只能手動(dòng)一篇一篇修改。如果那位置引用的文獻(xiàn)只有三四篇還好郁稍，如果是十篇以上赦政，就呵呵了。事實(shí)上耀怜，用了mendeley兩個(gè)月恢着，我只在一星期前的某天成功見到過merge按鈕的出現(xiàn)。去他們的論壇看财破，好多人反映merge功能的bug掰派，管理員說已經(jīng)修改好了，但實(shí)際上根本沒有左痢。我也留了個(gè)帖子靡羡，沒人鳥我系洛。
想想人家畢竟是免費(fèi)軟件，人家就是不鳥你你能如何略步。
這個(gè)bug真是非常讓人沮喪描扯，因?yàn)閙erge citation是一個(gè)非常基礎(chǔ)的功能趟薄，連這都做不到绽诚，就相當(dāng)于一臺(tái)手機(jī)的短信功能出了問題一樣。如果三個(gè)月后仍然是這樣杭煎，我應(yīng)該會(huì)轉(zhuǎn)去用別的軟件了恩够，雖然我還是挺喜歡它的。
順便吐槽一下羡铲，現(xiàn)在貌似還是比較流行用tex寫科技論文的吧∶的瘢現(xiàn)在被迫一直用word，感覺有點(diǎn)不大好犀勒。雖然我也不見得多喜歡tex屎飘，這貨也是老有bug。
總之就是基本上罕有“一切順利發(fā)展”的體驗(yàn)贾费。老是不是這里壞了就是那里壞了钦购，這種生活很討厭。
前幾天花了不少精力褂萧，終于畫出美好的圖了押桃，挑幾張展示一下：
誒？這才發(fā)現(xiàn)美麗的圖全是ud的majority导犹。
由于腳本改了又改唱凯，我都忘記關(guān)鍵改動(dòng)是什么了。目前最終的腳本的關(guān)鍵部分是這樣的：
levels = numpy.logspace(math.log10(min_dd), math.log10(max_dd), 40)
levels_bar = numpy.logspace(math.log10(min_dd), math.log10(max_dd), 5)
CF = pylab.contourf(K_A, K_B, T_dd,
interpolation = 'bilinear',
levels = levels,
cmap = pylab.cm.Spectral,
norm = LogNorm(clip=True))
pylab.colorbar(CF,
ticks = levels_bar,
format = '%5.1e')
不過有一點(diǎn)很關(guān)鍵的是谎痢，用atk算透射系數(shù)時(shí)最好不要用KrylovSelfEnergy算法磕昼，之前用它，得出一堆為負(fù)的透射系數(shù)节猿，害我浪費(fèi)不少時(shí)間去想它的物理意義票从，然后查了一下才知道只是算法的問題。然后換了RecursionSelfEnergy滨嘱，速度也不見得比前者慢峰鄙，但是數(shù)據(jù)就準(zhǔn)確很多。
另外太雨，Mendeley兩天后聯(lián)系我了吟榴，讓我發(fā)出問題的文檔過去。我做了個(gè)test文件發(fā)了過去囊扳，現(xiàn)在還沒回吩翻。
