matirix market データ読み込み方法は？

質問

今、現在連立一次方程式の解法に関する研究を行っています。そのために、MATRIX MARKETで行列データを取得して、それをC言語で読み込んで、配列に格納しようとしているのですが、うまく読み込めません。
格納しようとしているのは、
%%MatrixMarket matrix coordinate real general
% RBTitle: TUBULAR REACTOR MODEL
% RBKey: TUB100
100 100 396
1 1 -1.06376100E+03
2 1 1.80726500E-01
3 1 5.04700000E+02
1 2 -8.04861900E+00
2 2 -1.06287600E+03
4 2 5.04700000E+02
1 3 9.60400000E+02
3 3 -9.60783900E+02
4 3 1.91944900E-01
5 3 5.04700000E+02
2 4 9.60400000E+02
3 4 -8.37828500E+00
4 4 -9.59710900E+02
　　　・・・・・・・・・・・
　　　・・・・・・・・・・・
こういった非ゼロ要素のみが記述されている.mtxファイルです。
ご存知でしたら、よろしくお願いします。

回答No.2

echoes_x86 · ベストアンサー

こんばんは．

「何をしたいのか」について補足をお願いします．
>コマンドプロンプトに読み込んだデータを出力することはできるのですが、配列に格納…
この仰りようが既に矛盾しているということにお気づきでしょうか？

私が触ったのは結構前ですが，そのときからコードの内容が変わったのかと思い，
現在のコードを入手し，手元で再び実行してみました．
与えたデータファイルは
http://math.nist.gov/MatrixMarket/data/NEP/tubular/tub100.html
です．
結果，正しく標準出力にデータが出力されました．
これはexample_read.cの実行結果です．
つまり，この時点で「配列I，J，valにデータは格納されています」．

フォーマットを何か別の型式に変更したいということでしょうか？
それとも「example_read.cはコンパイルできたのだが，自分のコードはうまく行かない」ということでしょうか？
後者であればmatrix market云々ではなく，単純に質問者様の実装が不適切なだけ，
C言語そのものについて知識が不足しているということになります．
とりあえず，
fscanf(fp,"%d %d %d ",&a[i],&col_ind[i],&val[i]);
の部分は書式指定が不適切なのでこのままではデータを正しく読めません．
勝手に変更したりせずにサンプルコードをそのまま使えば良いのでは？

お礼

再度のご回答ありがとうございます。
結局、example_read.cにおいて既に配列に格納されていることに恥ずかしながら気付いていませんでした。。
回答者様のご指摘のとおり、Ｃ言語についてはまだまだ知識足らずで、自分で上記のコードで読み込もうとしていました。
丁寧に教えていただきありがとうございました。
無事、したいことができました。

回答No.1

echoes_x86

こんばんは．

ちょっと触ったことがある程度ですが…
matrix marketに読込の為のソースコード(参考URL)がありますよ．
ヘッダの解析等一通りできるはずです．
まぁ，私は「非零要素が何の順番で並んでいるか」
が分からなかったので読み込んだ後ソートして使ってましたが．
こんなことをしていたのは「行列の乗算」を「転置行列の乗算」きちんと行いたかったためです．

参考URL：

http://math.nist.gov/MatrixMarket/mmio-c.html

お礼

ご回答ありがとうございます。
おっしゃるとおり、そのソースコードを使ってみたのですが、うまくいかないのです。コマンドプロンプトに読み込んだデータを出力することはできるのですが、配列に格納、もしくは、別ファイルに出力することができません。
読み込みには、
fp=fopen("tub100.mtx","r")
while( 1)
{
　　fscanf(fp,"%d %d %d ",&a[i],&col_ind[i],&val[i]);
　　if( feof(fp)) break;
　　i++;
}
という感じで配列に格納しようとしています。