c言語を使ってgnuplotを起動してグラフを書かせる場合、次の2通りの方法が考えられる。

　(a) c言語でtxtファイルに出力してから、gnuplotを起動して書かせる

　(b) c言語でgnuplotにデータをそのまま渡して書かせる

参考ページ
C言語からgnuplotを操作する

Q : それではどちらの方が実行速度が早いのだろうか？

今回はこの疑問を解決する

---

プログラムで一番実行に時間がかかるのはデータIO(Input-Output、データの書き込みや読み込み)だと言われている
とすると、(a)はデータの書き込みをしてファイルを作るので(a)の方が時間がかかりそうな先入観が・・・

実際にやってみた(トリビア風)

---

(a)と(b)の実行時間を比較してみました。



赤が(a)、青が(b)
横軸はデータの長さ
縦軸は実行時間 (/usr/bin/time を使って5回測定してその平均した値)

これを見る限りでは、データを一度txtファイルとかに出力した方が早い・・・・
たぶんgnuplotに直接データを渡すとき、一度全部溜め込んでから書くとか微妙なことをしてるのかなー・・・

何にせよデータ数が300行を超えたあたりからは、一度出力した方が早そうですねー
(どのみち再現するためにもtxtには出力したいところなんだけどね)


ソースは右下の続きを読むに付けておきます。

(a) のソース

int main(int argc, char *argv[]){

FILE *gp, *fp;
int i;
double n, a, b;

if(argc > 1){
n = atof(argv[1]);
}
else{
n = 100;
}
fprintf(stdout, "n = %.0le\n",n);

if( (fp = fopen("test_data.txt","w")) == NULL){
perror("data fopen error");
exit(1);
fclose(fp);
}

for(i=1; i<=n; i++){
a = 1/n*i;
b = sin(2*M_PI*10*a);
fprintf(fp, "%5.5lf %5.5lf\n", a, b);
}

fclose(fp);

gp = popen("gnuplot -persist","w");
fprintf(gp, "set terminal pdf\n");
fprintf(gp, "set output '/dev/null'\n");

fprintf(gp, "plot \"test_data.txt\" w l lw 2\n");
fprintf(gp, "set output \"hoge.pdf\"\n");
fprintf(gp, "replot\n");

pclose(gp);
return 0;
}

(b) のソース

int main(int argc, char *argv[]){

FILE *gp, *fp;
int i;
double n, a, b;

if(argc > 1){
n = atof(argv[1]);
}
else{
n = 100;
}
fprintf(stdout, "n = %.0le\n",n);

gp = popen("gnuplot -persist","w");

fprintf(gp, "set terminal pdf\n");
fprintf(gp, "set output \"hoge.pdf\"\n");
fprintf(gp, "plot '-' w l\n");
for(i=0; i a = 1/n*i;
b = sin(2*M_PI*10*a);
fprintf(gp, "%5.5lf %5.5lf\n", a, b);
}
fprintf(gp, "e\n");

pclose(gp);
return 0;
}