日本极品中文字幕不卡,亚洲男人天堂色噜噜av,免费一区二区三区日本

nicer

支持

11樓2010-04-07 07:13:03

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

tjyl

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3218.1
紅花: 2
帖子: 576
在線: 156.6小時
蟲號: 765184
注冊: 2009-05-07
專業(yè): 生物無機(jī)化學(xué)
管轄: 程序語言

★ ★ ★ ★ ★
qli0218(金幣+10):非常感謝��！ 2010-04-07 11:23
hometownlove(金幣+5):學(xué)習(xí)了 2010-04-07 20:17

計算軟件：abinit
應(yīng)用軟件：C,Gnuplot
作用：直接畫一個簡單能帶圖，方便看結(jié)果
方法：運(yùn)行時提取出xxx_EIG里的值，產(chǎn)生一個gnuplot用的腳本，然后把圖輸出到一個ps文件里。
使用：
以tutorial #3里3.5的輸出 t3xo_DS2_EIG為例子
gcc eig.c -lm
./a.out t3xo_DS2_EIG
輸入 kptopt 的絕對值
3
然后一側(cè)輸入ndivk的幾個值，每輸入一個按回車
10
12
17

程序（可選）：

CODE:

 #include

#include

#include

#include

double **matrix (int R, int L);        //為數(shù)組M[R][L]開辟空間。

void freem (double **m, int R, int L);

double emin (double **m, int R, int L);

double emax (double **m, int R, int L);

int

main (int argc, char **argv)

{

  int knum;

  int nband;

  int i, j;

  FILE *fp;

  FILE *pt, *pd;

  double **eig;

  double ymax, ymin;

  int kptopt = 3;

  int *ndivk;

  if (argc != 2)

    {

      printf ("this program only need 2 parameters.\n");

      return -1;

    }

  printf ("please input the kptopt:\n");

  scanf ("%d", &kptopt);

  ndivk = malloc (sizeof (int) * kptopt);

  printf ("please input the ndivk:\n");

  for (i = 0; i < kptopt; i++)

    scanf ("%d", ndivk + i);

  if ((fp = fopen (argv[1], "r")) == NULL)

    {

      printf ("can't open %s.\n", argv[1]);

      return -1;

    }

  fscanf (fp, "%*[^0-9]%d%*[^\n]\n", &knum);

  fscanf (fp, "%*[^=]=%d%*[^\n]\n", &nband);

  rewind (fp);

  eig = matrix (knum, nband);

  if (NULL == eig)

    {

      printf ("can't get memory for eig.\n");

      return -1;

    }

  fscanf (fp, "%*[^\n]\n");

  for (i = 0; i < knum; i++)

    {

      fscanf (fp, "%*[^\n]\n");

      for (j = 0; j < nband; j++)

        fscanf (fp, "%lf", &eig[i][j]);

      fscanf (fp, "\n");

    }

  pd = fopen ("plot.data", "w");

  if (NULL == pd)

    {

      printf ("fail to open plot.data.\n");

      return -1;

    }

  for (i = 0; i < knum; i++)

    {

      fprintf (pd, "%4d", i + 1);

      for (j = 0; j < nband; j++)

        fprintf (pd, " %.5f", eig[i][j]);

      fprintf (pd, "\n");

    }

  fclose (fp);

  fclose (pd);

  ymax = emax (eig, knum, nband);

  ymin = emin (eig, knum, nband);

  ymax = ymax + fabs (0.2 * ymax);        //設(shè)置畫圖所用的Y軸坐標(biāo)上下限

  ymin = ymin - fabs (0.2 * ymin);

  if ((pt = fopen ("plot", "w")) == NULL)

    {

      printf ("fail to open file 'plot'.\n");

      return -1;

    }

  fprintf (pt, "gnuplot<
        set term post landscape enhanced color dashed defaultplex \"Helvetica\" 14\n\

        set output \'pic.ps\'\n\

        set key outside\n\

        set title \"Energy Band\"\n\

        set ylabel \"Energy (Ev)\"\n");

  fprintf (pt, "set yrange[%f:%f]\n", ymin, ymax);

  fprintf (pt, "set xrange[%d:%d]\n", 1, knum);

  j = 0;

  for (i = 1; i < kptopt; i++)

    {

      j = j + ndivk[i - 1];

      fprintf (pt, "set arrow from %d,%f to %d,%f\n", j + 1, ymin, j + 1, ymax);        //在特殊點(diǎn)處畫箭頭

    }

  fprintf (pt,

           " plot 'plot.data' u 1:2  w l lw 2 smooth cspline,'plot.data' u 1:2  w p  ps 0.3 ");

  for (i = 1; i < nband; i++)

    fprintf (pt,

             " ,'plot.data' u 1:%d  w l lw 2 smooth cspline,'plot.data' u 1:%d w p  ps 0.3 ",

             i + 2, i + 2);

  fprintf (pt, "\nEOF");

  fclose (pt);

  chmod ("plot", S_IRUSR | S_IWUSR | S_IXUSR | S_IXGRP | S_IXOTH);

  system ("./plot");

  remove ("plot");

  remove ("plot.data");                //結(jié)束后刪除中間產(chǎn)生的數(shù)據(jù)和gnuplot腳本.

  freem (eig, knum, nband);

  return 0;

}

double **

matrix (int R, int L)

{

  double **M;

  int i;

  if ((M = malloc (sizeof (double *) * R)) == NULL)

    {

      printf ("can't get memory1.\n");

      return NULL;

    }

  for (i = 0; i < R; i++)

    {

      M[i] = calloc (sizeof (double), L);

      if (NULL == M[i])

        {

          printf ("can't get memory2.\n");

          return NULL;

        }

    }

  return M;

}

void

freem (double **m, int R, int L)

{

  int i;

  for (i = 0; i < R; i++)

    free ( m[i] );

  free (m);

}

double

emax (double **m, int R, int L)

{

  int i, j;

  double max;

  max = m[0][0];

  for (i = 0; i < R; i++)

    for (j = 0; j < L; j++)

      max = (max > m[i][j]) ? max : m[i][j];

  return max;

}

double

emin (double **m, int R, int L)

{

  int i, j;

  double min;

  min = m[0][0];

  for (i = 0; i < R; i++)

    for (j = 0; j < L; j++)

      min = (min < m[i][j]) ? min : m[i][j];

  return min;

}

引用回帖:

Originally posted by qli0218 at 2010-04-06 14:11:26:
[box=white]現(xiàn)在做第一性原理計算越來越多，但是如何發(fā)表高層次的文章卻與處理、挖掘所得數(shù)據(jù)的深度有著密切的聯(lián)系。所以希望大家互幫互助，分享處理數(shù)據(jù)的經(jīng)驗。

大家所用軟件都不一樣， ...

[ Last edited by tjyl on 2010-4-7 at 10:23 ]

贊一下(1人)

回復(fù)此樓

12樓2010-04-07 10:14:18

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

tjyl

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3218.1
紅花: 2
帖子: 576
在線: 156.6小時
蟲號: 765184
注冊: 2009-05-07
專業(yè): 生物無機(jī)化學(xué)
管轄: 程序語言

★ ★ ★ ★ ★ ★ ★
aylayl08(金幣+5):非常感謝，第一個提供程序的，鼓勵一下 2010-04-07 13:07
hometownlove(金幣+2):謝謝交流 2010-04-07 20:18

前面的代碼里freem函數(shù)貼出來后顯示的有問題。

CODE:

void freem (double **m, int R, int L)

{

  int i;

  for (i = 0; i < R; i++)

    free ( m[i]);

  free (m);

}

如果能帶數(shù)量多，就加
unset key\n\

CODE:

fprintf(pt, "gnuplot<
        set term post landscape enhanced color dashed defaultplex \"Helvetica\" 14\n\

        set output \'pic.ps\'\n\

        set key outside\n\

        set title \"Energy Band\"\n\

        unset key\n\

        set ylabel \"Energy (Ev)\"\n");

引用回帖:

Originally posted by tjyl at 2010-04-07 10:14:18:
計算軟件：abinit
應(yīng)用軟件：C,Gnuplot
作用：直接畫一個簡單能帶圖，方便看結(jié)果
方法：運(yùn)行時提取出xxx_EIG里的值，產(chǎn)生一個gnuplot用的腳本，然后把圖輸出到一個ps文件里。
使用：
以tutorial #3里3.5的輸 ...

[ Last edited by tjyl on 2010-4-7 at 10:51 ]

贊一下(2人)

回復(fù)此樓

13樓2010-04-07 10:18:13

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

銳利的碎片

木蟲 (正式寫手)

star watcher

1ST強(qiáng)帖: 1
應(yīng)助: 136 (高中生)
金幣: 3637.1
散金: 252
紅花: 22
帖子: 988
在線: 1224.9小時
蟲號: 961933
注冊: 2010-03-05
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

其實scilab處理數(shù)據(jù)很方便，主要是可以進(jìn)行很多分析工作。語法和matlab十分接近。直接下載可執(zhí)行版加上環(huán)境變量就可以用了。當(dāng)然freemat也很好。

[ Last edited by 銳利的碎片 on 2010-4-7 at 10:54 ]

贊一下

回復(fù)此樓

14樓2010-04-07 10:53:06

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

tjyl

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3218.1
紅花: 2
帖子: 576
在線: 156.6小時
蟲號: 765184
注冊: 2009-05-07
專業(yè): 生物無機(jī)化學(xué)
管轄: 程序語言

octave用法和matlab基本一樣的。

引用回帖:

Originally posted by 銳利的碎片 at 2010-04-07 10:53:06:
其實scilab處理數(shù)據(jù)很方便，主要是可以進(jìn)行很多分析工作。語法和matlab十分接近。直接下載可執(zhí)行版加上環(huán)境變量就可以用了。當(dāng)然freemat也很好。

[ Last edited by 銳利的碎片 on 2010-4-7 at 10:54 ]

贊一下

回復(fù)此樓

15樓2010-04-07 11:36:36

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

楊婧

16樓2010-04-07 11:40:51

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

銳利的碎片

木蟲 (正式寫手)

star watcher

1ST強(qiáng)帖: 1
應(yīng)助: 136 (高中生)
金幣: 3637.1
散金: 252
紅花: 22
帖子: 988
在線: 1224.9小時
蟲號: 961933
注冊: 2010-03-05
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by tjyl at 2010-04-07 11:36:36:
octave用法和matlab基本一樣的。

單純功能而言，scilab要比octave好。要長得像freemat的幫助直接抄的matlab。

贊一下

回復(fù)此樓

17樓2010-04-07 11:49:31

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

RHOW

18樓2010-04-07 12:52:10

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

manin8769

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 908.8
帖子: 290
在線: 38小時
蟲號: 581569
注冊: 2008-07-18
性別: MM
專業(yè): 無機(jī)合成和制備化學(xué)

MATLAB

回復(fù)此樓

19樓2010-04-07 14:30:30

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

faith2007

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 42.9
帖子: 76
在線: 97.6小時
蟲號: 620425
注冊: 2008-10-08
專業(yè): 核技術(shù)及其應(yīng)用

占個沙發(fā)

回復(fù)此樓

20樓2010-04-07 20:07:33

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

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

[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 266求調(diào)劑 +4	陽陽哇塞 2026-03-14	8/400	2026-03-18 12:43 by Linda Hu
[考研] 化學(xué)工程321分求調(diào)劑 +14	大米飯！ 2026-03-15	17/850	2026-03-18 12:40 by Linda Hu
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 考研調(diào)劑 +4	芬達(dá)46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao