配列の中身が変わってしまう。

Question

前提・実現したいこと

gnuplot で配列の中身を描画したい。だが、なぜか配列の中身が変わって、困ってしまった。

/* プロット開始 */
    ///*
gp = popen("gnuplot -persist", "w"); // パイプを開き、gnuplotの立ち上げ

の行の前後で配列の中身が変わってしまっている。
「ここからおかしい」の配列の1個目と2個目が
100, 75
から
-0.0, 0.0

に変わってしまっている。

どうすれば正しい値が出力されるのでしょうか？

発生している問題・エラーメッセージ

エラーは起きていない

該当のソースコード

#include<stdio.h>
#include<math.h> // 累乗計算で使う
#define N 5 //#define マクロ名 定数や式

/* PIM throughput                            */
/* Find the PIM throughput when N=2 to N=50. */

/*  Programming                                             */
/* Write a program to obtain the result when N=2 to N=50.   */
/* Plot the graph.                                          */


/*
再起関数で組み合わせnCrを作る
forで2~50回回す
毎回出力する
*/



float combination(float n, float r){
    /* nCrの計算 */

    if(n == r) return(1.0);        // (ex) 3C3 = 1
    else if(n == 0.0) return(1.0); // (ex) 3C0 = 1
    else if(r == 1.0) return(n);   // (ex) 3C1 = 3

    /* ------------------- */
    // 組み合わせは公式
    //   nCr = n-1Cr + n-1Cr-1
    // より、再起関数で表せる。
    /* ------------------- */
    return combination(n-1, r) + combination(n-1, r-1);
}


float caluculation(int n){
    /* --------------------------------------------- */
    /* 今回計算したいやつ                               */
    /* ∑{n}{k = 1} nCk * (1 / n)^k * (1 - 1/n)^(n-k) */
    /* --------------------------------------------- */

    float count_box = 0.0; // 一時的に保管する変数

    for(int r=1; r<=n; r++){
        // (ex) 0.75 = 2C1 * (1/2)^1 * (1 - 1/2)^(2-1) + 2C2 * (1/2)^2 * (1 - 1/2)^(0)
        count_box += combination(n, r) * pow(1.0/n, r) * pow(1.0 - 1.0/n, n-r);
    }
    return count_box;
}

int x[]         = {};

int main(void){
    /* --------------------------------------------- */
    /* nを2~50まで以下の式を回す。                       */
    /* ∑{n}{k = 1} nCk * (1 / n)^k * (1 - 1/n)^(n-k) */
    /* --------------------------------------------- */

    FILE *gp;

    float throughput[] = {};

    printf("\n");

    for(int i=1; i<=N; i++){
        x[i] = i;
        printf("x[%d] = %d\n", i, x[i]);

        throughput[i] = caluculation(i) * 100;

        printf("N= %d : Throughput=%.2f%%\n", i, throughput[i]);

        if(i%10 == 0) printf("\n");
    }
    printf("\n");


    for (int j=1; j <= N; j++) {
        // 確認
        printf("j=%d\tx[%d]=%d\tthroughput[%d]=%.2f\n", j, j, x[j], j, throughput[j]);
        /*
        j=1     x[1]=1  throughput[1]=100.00
        j=2     x[2]=2  throughput[2]=75.00
        j=3     x[3]=3  throughput[3]=70.37
        j=4     x[4]=4  throughput[4]=68.36
        j=5     x[5]=5  throughput[5]=67.23
        */
    }
    printf("\n");


    for (int m=1; m <= N; m++) {
        // 確認
        printf("m3=%d\tx[%d]=%d\tthroughput[%d]=%.2f\n", m, m, x[m], m, throughput[m]);
        /*
       m3=1    x[1]=1  throughput[1]=100.00
       m3=2    x[2]=2  throughput[2]=75.00
       m3=3    x[3]=3  throughput[3]=70.37
       m3=4    x[4]=4  throughput[4]=68.36
       m3=5    x[5]=5  throughput[5]=67.23
        */
    }
    printf("\n");

    /* プロット開始 */
    ///*
    gp = popen("gnuplot -persist", "w"); // パイプを開き、gnuplotの立ち上げ


    for (int m=1; m <= N; m++) {
        // 確認
        printf("ここからおかしいm2=%d\tx[%d]=%d\tthroughput[%d]=%.2f\n", m, m, x[m], m, throughput[m]);
        /*
        こからおかしいm2=1    x[1]=1  throughput[1]=-0.00
        ここからおかしいm2=2    x[2]=2  throughput[2]=0.00
        ここからおかしいm2=3    x[3]=3  throughput[3]=70.37
        ここからおかしいm2=4    x[4]=4  throughput[4]=68.36
        ここからおかしいm2=5    x[5]=5  throughput[5]=67.23
        */

    }
    printf("\n");



    fprintf(gp, "set multiplot\n"); // マルチプロットモード
    fprintf(gp, "set xrange [0:%d]\n", x[N]); // 範囲の指定
    fprintf(gp, "set yrange [60:100]\n");

    fprintf(gp, "set xlabel \"size n\"\n"); // ラベル表示
    fprintf(gp, "set ylabel \"Throughput\"\n");

    for (int m=1; m <= N; m++) {
        // 確認
        printf("m=%d\tx[%d]=%d\tthroughput[%d]=%.2f\n", m, m, x[m], m, throughput[m]);
        /*
        m=1     x[1]=1  throughput[1]=-0.00
        m=2     x[2]=2  throughput[2]=0.00
        m=3     x[3]=3  throughput[3]=70.37
        m=4     x[4]=4  throughput[4]=68.36
        m=5     x[5]=5  throughput[5]=67.23
        */

    }
    printf("\n");

    // 曲線のプロット
    fprintf(gp, "plot '-' with lines linetype 1\n");
    for (int k = 1; k <= N; k++) {
        printf("k=%d\tx[%d]=%d\tthroughput[%d]=%.2f\n", k, k, x[k], k, throughput[k]);
        /*
        k=1     x[1]=1  throughput[1]=-0.00
        k=2     x[2]=2  throughput[2]=0.00
        k=3     x[3]=3  throughput[3]=70.37
        k=4     x[4]=4  throughput[4]=68.36
        k=5     x[5]=5  throughput[5]=67.23
        */
        fprintf(gp, "%d\t%f\n", x[k], throughput[k]);
    }
    printf("\n");
    fprintf(gp, "e\n");




    // 点のプロット
    fprintf(gp, "plot '-' with points pointtype 1\n");
    for (int l = 1; l <= N; l++) {
        printf("l=%d\tx[%d]=%d\tthroughput[%d]=%.2f\n", l, l, x[l], l, throughput[l]);
        /*
        l=1     x[1]=1  throughput[1]=-0.00
        l=2     x[2]=2  throughput[2]=0.00
        l=3     x[3]=3  throughput[3]=70.37
        l=4     x[4]=4  throughput[4]=68.36
        l=5     x[5]=5  throughput[5]=67.23
        */
        fprintf(gp, "%d\t%f\n", x[l], throughput[l]);
    }
    fprintf(gp, "e\n");




    fprintf(gp, "set nomultiplot\n"); // マルチプロットモード終了
    fprintf(gp, "exit\n"); // gnuplotの終了
    fflush(gp);
    pclose(gp); // パイプを閉じる

    //*/
    return 0;
}

試したこと

ここに問題に対して試したことを記載してください。

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Answer 1

int x[]         = {};
float throughput[] = {};

要素数が0です。
範囲外アクセスしているため何が起こっても不思議はないです。
必要なサイズ(インデックスがすべて1開始になっているのでN+1)の領域を確保しましょう。

int x[N+1]         = {};
float throughput[N+1] = {};

Comments

[off]

回答ありがとうございます！修正したら、正しい値を配列に格納することができました！

配列を初期化する時点で、何個分必要なのかという枠を定義してあげる必要があるということですか。
大変参考になりました。ありがとうございました。