CSVファイルの読み込み

質問

以下のようなCSVファイルを読み込みたいと思っています。
0.575092,0.030525
0.565324,0.018315
0.555556,0.013431
0.553114,0.013431
0.54823,0.015873
0.538462,0.013431
0.531136,0.006105
0.52381,-0.003663
0.516484,-0.010989
0.501832,-0.015873
0.489622,-0.018315
0.477412,-0.020757
0.46276,-0.020757
・
・
・
50000＊2のデータで、以下のようなプログラムを使用しました。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

/* 確保するデータ保存領域の大きさ(N行×M列) */
#define N 50000
#define M 2
/* データの区切り文字 */
#define SEP_DATA ','

int csv_read(char filename[], double csv[N][M]) {
/* ファイルオープン */
FILE *fp;
if( (fp = fopen(filename, "r")) == NULL ) {
printf(" file open error!!\n");
return -1;
}

/* １行毎に読み出し */
char line[256], *ptr;
int i, j, k;
i=0;
while (fgets(line, 256, fp) != NULL) {
printf("*%s", line);
ptr = line; j=0;
do{
/* line[j]から次のタブ文字までを数値に変換 */
csv[i][j] = atof(ptr);
/* 次のタブ文字の位置を探す */
ptr = strchr(ptr, SEP_DATA);
/* タブ文字の次の文字を示す */
if (ptr!=NULL) { ptr++; }
j++;
}while(ptr!=NULL && j<M);
i++;
}

/* ファイルクローズ */
fclose(fp);
return 0;
}

int main(int argv, char *argc[]) {
char filename[256];
if( argv > 1){
strcpy(filename, argc[1]);
} else {
printf("Please Input Filename:");
scanf("%s", filename);
}

/* データ保存用の領域を確保 */
double (*csvdata)[M];
csvdata = (double(*)[M])malloc(sizeof(csvdata) * N);
if ( csvdata == NULL ){
return -1;
}
int i,j;
/* 配列の初期化 */
for( i=0; i<N; i++) {
for( j=0; j<M; j++) {
csvdata[i][j] = 0.0;
}
}

/* CSVデータの読み込み */
if( csv_read(filename, csvdata) < 0 ) {
return -1;
}

/* 配列の出力 */
for( i=0; i<N; i++) {
printf("%lf", csvdata[i][0]);
for( j=1; j<M; j++) {
printf("\t%lf", csvdata[i][j]);
}
printf("\n");
}
free(csvdata);
csvdata = NULL;
return 0;
}
しかし、コンパイルは通るのですが実行時にエラーが出てしまいます。
データの大きさに問題があると思うのですが、原因がわかりません。おかしな所のご指摘を頂きたいです。よろしくお願いします。

回答No.5

php504 · ベストアンサー

原因を書いてませんでしたね
csvdata = (double(*)[M])malloc(sizeof(csvdata) * N);
でsizeof(csvdata)はポインタのサイズなので4とか8の環境が多いでしょう
double2個だと普通16なので足りませんよね
csvdata = (double(*)[M])malloc(sizeof(double[M]) * N);
でいいと思いますが私はtypedefしたほうが好きです

補足

ありがとうございます。なるほど、原因がわかりました。
おっしゃる方法で無事に読み込みも成功しました。

回答No.4

Tacosan

あなたの想像通り, malloc で確保するサイズを間違えています.
よ～く考えればわかりますが, * を 1個忘れてます.

補足

ありがとうございます。＃３のかたのご指摘で確保するサイズが間違っているのは理解できたのですが、＊をいれるというのがどの部分かわかりません。よろしければご指南いただけないでしょうか。

回答No.3

php504

自動変数で10万個の配列は無謀ですが静的変数なら大丈夫でしょう
2次元配列（本当は配列ではなくてポインタですが）の動的確保にはいろいろありますが
typedef double DARRAY[M];
としておいて
DARRAY* csvdata;
csvdata = (DARRAY*)malloc(sizeof(DARRAY) * N);
というのはどうでしょう

回答No.2

kot_o

通りかかったものです。

ソースをざっと見て思ったことを、書いていきます。

/* データ保存用の領域を確保 */
>double (*csvdata)[M];

ここは何をやっているのですか？csvdataっという変数は本当に実体はありますか？
二次元配列の宣言ならなら
double *csvdata[N];
() はいらないかと思います。

>csvdata = (double(*)[M])malloc(sizeof(csvdata) * N);
これは、ちゃんと動作していますか？
forループを使い、配列ごとに動的配列を入れていってはどうでしょう。

for( int i = 0 ; i < N ; i++ )
{
csvdata[i] = (double*)malloc(sizeof(double) * M);
}

ここまで見た感じだと、MとMの使い方が逆ではありませんか？
そこを、デバッグすることをお勧めします。

お礼

ありがとうございます。
ご指摘の通り、メモリの確保の部分で引っかかっておりました。
自分で作ったプログラムではないので、自分もここが理解できなかったのですが、（）をとるとキャストのエラーが出てしまうためつけておりました。
おっしゃる方法で行ってみましたが、やはり関数呼び出しの段階で、
『 error C2664: 'csv_read' : 2 番目の引数を 'double *[2]' から 'double [][2]' に変換できません。(新しい機能 ; ヘルプを参照)
指示された型は関連がありません。変換には reinterpret_cast、C スタイル キャストまたは関数スタイルのキャストが必要です。』
のようなエラーが出てしまいます。

回答No.1

php504

確保する個数が前もってわかっているなら動的に確保する必要はないのでは
static double csvdata[N][M];
ではだめでしょうか

お礼

ありがとうございます。
上記は実験のサンプリングデータなのですが、サンプリング周波数や計測時間を変えた実験も行うので、データ長によっては、スタック領域ではいっぱいになってしまうと考え、このようにしました。