gnuplotをc言語で使用しようとしているのですが、エラーが出ます。

Question

gnuplotでline 0: warning: Skipping data file with no valid points とエラーが表示されます。

数値計算をしたデータをcsvファイルで出力し、そのファイルをgnuplotでグラフにしようとしているのですが、

line 0: warning: Skipping data file with no valid points
line 0: warning: Skipping data file with no valid points 
line 0: warning: Skipping data file with no valid points  

という表示が出てグラフが出力されません。
解決方法がわかる方いらっしゃいましたら、回答よろしくお願いします。

エラーが出るソースコードを以下に貼ります。

#include <stdio.h>
#include <stdlib.h>
#define GNUPLOT_PATH "C:/gnuplot/bin/gnuplot.exe"
#define FILE_PATH "./plot.csv"
#define J_MAX 100

int main()
{
    int n_max;
    double dx = 0.01, dt = 0.01, c = 0.5, time = 0.0;
    double c_fl, t_max;
    double u[J_MAX+1], u_nn[J_MAX+1], u_ex[J_MAX+1], px[J_MAX+1];
    printf("Input n_max: ");
    scanf("%d", &n_max);
    c_fl = c * dt / dx;
    t_max = (float)n_max * dt;

    for (int j = 0; j <= J_MAX; j++)
    {
        px[j] = dx * j;
        if (px[j] >= 0.2 && px[j] <= 0.4)
            u[j] = 1.0;
        else
            u[j] = 0.0;
    }

    for (int n = 0; n < n_max; n++)
    {
        for (int j = 0; j < J_MAX; j++)
        {
            u_nn[j] = u[j] + c_fl * (u[j+1] - u[j]);
            u[j] = u_nn[j];
        }
        time += dt;
    }

    printf("t = %.2f\n", time);

    for (int j = 0; j <= J_MAX; j++)
    {
        px[j] = j * dx;
        if (px[j] >= 0.2+c*t_max && px[j] <= 0.4+c*t_max)
            u_ex[j] = 1.0;
        else
            u_ex[j] = 0.0;
    }

    
    FILE *fp;
    fp = fopen(FILE_PATH, "w");
    if (fp == NULL)
    {
        fprintf(stderr, "Oops, i can't find %s.", FILE_PATH);
        exit(EXIT_FAILURE);
    }

    for (int j = 0; j < J_MAX; j++)
    {
        fprintf(fp, "%f, %f, %f\n", px[j], u[j], u_ex[j]);
    }


    FILE *gp;
    gp = _popen(GNUPLOT_PATH, "w");
    if (gp == NULL) 
    {
        fprintf(stderr, "Oops, I can't find %s.", GNUPLOT_PATH);
        exit(EXIT_FAILURE);
    }
    // gnuplotにコマンドを送る
    fprintf(gp, "set xrange [-2.0:2.0]\n");
    fprintf(gp, "set yrange [-2.0:2.0]\n");
    fprintf(gp, "set xlabel 'x'\n");
    fprintf(gp, "set ylabel 'u(x,t)'\n");
    fprintf(gp, "plot 'plot.csv' using 1:2 with lines\n");
    fprintf(gp, "replot 'plot.csv' using 1:3 with lines\n");
    fflush(gp); // バッファに格納されているデータを吐き出す（必須）
    getchar();	// 入力待ち
    _pclose(gp);
    fclose(fp);
    exit(EXIT_SUCCESS);
}

Answer 1

gnuplot では標準では CSV を受け付けません。

set datafile separator comma

するか、出力内容にスペースも入れたままなら

set datafile separator " ,"

とする。

separator を設定しないなら出力部分を

fprintf(fp, "%f, %f, %f\n", px[j], u[j], u_ex[j]);

ではなくて

fprintf(fp, "%f %f %f\n", px[j], u[j], u_ex[j]);

とすることでも解決しますね。

---

プログラムを動かすと fclose の位置が悪いためにファイルが出力される前の段階で gnuplot に読み込ませていますね。データ出力が完了したら fclose してしまいましょう。

	for (int j = 0; j < J_MAX; j++)
	{
		fprintf(fp, "%f %f %f\n", px[j], u[j], u_ex[j]);
	}
	fclose(fp);

それでデータを読み込むことはできるようになります。あとは最後のキー入力待ちが正しく動作しません。入力バッファをクリアした後で getchar で待ちましょう。

	rewind(stdin);
	getchar();    // 入力待ち

Comments

[chiku_soh]

datファイルに変更し、fprintfの部分を変更したのですが、まだ同じエラーが出ます。どうすればいいですか。

[atata0319]

回答を追記しました。

[chiku_soh]

エラーは消えましたが、グラフが一瞬で閉じてしまうのですが、sleepのような関数って使用可能ですか？

[atata0319]

回答の最後に記載してある通り、getchar の前に rewind してみましょう。または、fgets で止めるのもありかと思います。それなら Enter キー押すまで進まないので。

[atata0319]

fgets は引数の指定が面倒なので、再度 scanf("%d", &n_max); で止めても良いかと思います。

[chiku_soh]

ありがとうございました