シミュレーションによって得られた計算結果を自動的にgnuplotで
表示させるプログラミングを作成しました。シミュレーション自体は
fortranで行っていてfort.61というファイルが出力されます。
c言語でこのファイルを読みとって、それに少し細工をしfort.62
というファイルに新たに書き込んで、gnuplotで表示させたいのですが、
gnuplotの画面で no valid dataというエラーが表示されてしまいます。
ところが、cのプログラムを終了させてからgnuplotでこのfort.62を
通常の方法でplotすると、きちんとグラフが表示されます。
何が原因がわかるかたいらっしゃいましたら、教えて下さい。
ちなみにコンパイルはcygwin上でgccを使用しています。以下にソース
を添付します。gnuplotで読み込んでいるファイルplot.pltは
「plot "fort.62"」と記述してあるtextファイルです。
よろしくお願いします。
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int main(void) {
FILE *fp,*fp2;
int i,N;
double x,y;
double volt[100],curr[100];
double norm,para;
para = 0.2442E-03;
N=14;
fp=fopen("fort.61","r");
for(i=0;i<N;i++){
fscanf(fp,"%lf %lf",&x,&y);
volt[i] = x;
curr[i] = y;
if(volt[i] == -1.0){
norm = para/curr[i];
}
}
for(i=0;i<N;i++){
curr[i] = norm*curr[i];
}
close(fp);
fp2 = fopen("fort.62","w");
for(i=0;i<N;i++){
fprintf(fp2,"%lf %lf\n",volt[i],curr[i]);
}
close(fp2);
system("wgnuplot -persist plot.plt");
}
No.2ベストアンサー
- 回答日時:
私はgnuplotについて何も知らないので、
とんちんかんな解答だったらすみません。
> double norm,para;
> if(volt[i] == -1.0){
> norm = para/curr[i];
> }
volt[i]がたまたま-1.0になるとき「だけ」normを計算し、
volt[i]が-1.0になることがない場合はnormの値は不定(ゴミ)です。
それでよいのでしょうか?
> fp=fopen("fort.61","r");
> fp2 = fopen("fort.62","w");
ファイルをオープンできなかったときの処理が書いてありません。
それでよいのでしょうか?
> close(fp);
> close(fp2);
closeという関数は標準関数に見あたりません。
fcloseか何かの間違いでしょうか。
どうもありがとうございます。
いつもcloseという関数を使用していました、fcloseだったのですね。
コンパイル中にエラーが出ないので気づきませんでした、貴重な意見
どうもありがとうございました。
No.3
- 回答日時:
No.2の方が指摘されているfclose関数をclose関数に間違えているのが原因だと思います。
fprintf関数は出力をバッファリングするため、fprintf関数から復帰してきた時点で
ファイルへの出力が行われていない可能性があります。
正しくfclose関数を呼び出せばバッファ内のデータもファイルへ出力されます。
また、プログラムの終了時には開いたままのファイルに対する処理が行われ
その際にバッファ中のデータもファイルに書き込まれるのだと思います。
gnuloptのエラーメッセージが「no valid data」なので、ファイルは存在していて、
データが全く無いか正しい形式のデータが無いのだと思います。
これは、close関数を誤って使ったためにファイルに全データが書き込まれていないため
発生しているのではないでしょうか。
まったくご指摘の通りでした、僕もファイルの書き込みが正常に終了
していないのではと予想して、close関数の位置をいじったりしてみたの
ですが、まさか、関数自体が間違っているとは気づきませんでした。
大変勉強になりました。どうもありがとうございます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# このプログラミング誰か教えてくれませんか 1 2022/06/02 15:27
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# c言語でユーザ関数を利用して複素数のべき乗と絶対値の数列を計算するプログラムが作りたいです。 3 2023/01/29 22:13
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- Excel(エクセル) VBA フォルダ見える化のコードについて 2 2023/06/19 15:04
- C言語・C++・C# LU分解法のピボット選択機能実装について(C言語・gcc-9) 1 2022/07/22 15:20
- その他(プログラミング・Web制作) Fortranでの出力ファイル 2 2023/03/21 21:25
- C言語・C++・C# C言語(構造体) 3 2022/07/05 20:08
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
VBSで指定行に挿入
-
テキストファイルの行数を取得...
-
freadとfwrite
-
C/C++ ファイル入出力操作 (長...
-
改行までの一文字ずつのファイ...
-
バイナリファイルをテキストフ...
-
VS2010 MFC CStdioFileについて
-
C言語 バイナリファイルの読み...
-
fopenで開いたファイルのサイズ...
-
テキストファイルの先頭への文...
-
ハフマン符号のプログラム
-
複数テキストファイルを読み込...
-
winsock recvでの文字化け
-
fseekの使い方
-
0バイトファイルの作成
-
EOF判定されない
-
ファイル内のデータを1行削除...
-
【VB.Net】バイト型配列に読み...
-
C言語初心者の質問失礼します。
-
どんなプログラムを書いても指...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
ファイル内のデータを1行削除...
-
バイナリファイルをテキストフ...
-
VBSで指定行に挿入
-
fgetsで2行目から文字化け
-
テキストファイルの行数を取得...
-
ファイルサイズ指定し、ファイ...
-
巨大なテキストファイル(可変...
-
【VB.Net】バイト型配列に読み...
-
fopenで開いたファイルのサイズ...
-
改行までの一文字ずつのファイ...
-
freadとfwrite
-
c言語 2つのファイルを行ご...
-
0バイトファイルの作成
-
fopen(書き込みモード)でファイ...
-
テキストファイルの先頭への文...
-
テキストファイルの文字列の削除
-
C言語での採番について
-
VS2010 MFC CStdioFileについて
-
続・EOF判定されない
-
C言語での改行コードの扱いにつ...
おすすめ情報