前提・実現したいこと
ディレクトリ内にあるファイル内の数値データ(カンマ区切り)を読み込み、
中心差分の計算を行った後、カンマ区切りでそれぞれの数値をファイルに出力していく
というプログラムを作成したいと考えています。
数値データを読み込み、計算するまでは行えたのですが、
ファイルへの出力が上手く出来ない状態です。
ご協力・コメント等いただけたら幸いです。
発生している問題・エラーメッセージ
特になし
しかし、np.savetxt()で新たにファイルを作成して書き込みを行っているためか
最後の計算結果のみしか書き込みがされていません。
改行して次の計算結果を同じファイル内に追記していくという方法で出力したい
と考えています。
該当のソースコード
import numpy as np
dt = 0.0003
f = 100
f2 = 200
f3 = 300
a = 'He_re1000_vor-0'
'''ファイル名の指定'''
for file in range(0,5):
```data2 = np.loadtxt(a+str(f2),delimiter=',',skiprows=1)
data3 = np.loadtxt(a+str(f3),delimiter=',',skiprows=1)
'''中心差分の計算とファイル出力'''
``````dvorx = [((data3[line,4]-data2[line,4])-(data2[line,4]-data1[line,4]))/(2*dt)]
```dvory =[((data3[line,5]-data2[line,5])-(data2[line,5]-data1[line,5]))/(2*dt)]
``````np.savetxt(dvorx)
```line += 1
f += 100
```f3 += 100
file += 1
試したこと
以前はstr()やnp.float64関係でのエラーが頻発していたので
文字型の問題かと思い、List化してから出力するようにしてから
書き込みが出来るようになりました。
補足情報(言語/FW/ツール等のバージョンなど)
プログラム環境 Python 3.6 Spyder
操作したいファイル名
He_re1000_vor-0100
He_re1000_vor-0200
He_re1000_vor-0300
・・・
He_re1000_vor-10900
※それぞれのファイル内には
点番号,座標データx,y,z,数値データ①,②,③
というデータが1,500,000行入っています。
-
気になる質問をクリップする
クリップした質問は、後からいつでもマイページで確認できます。
またクリップした質問に回答があった際、通知やメールを受け取ることができます。
クリップを取り消します
-
良い質問の評価を上げる
以下のような質問は評価を上げましょう
- 質問内容が明確
- 自分も答えを知りたい
- 質問者以外のユーザにも役立つ
評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。
質問の評価を上げたことを取り消します
-
評価を下げられる数の上限に達しました
評価を下げることができません
- 1日5回まで評価を下げられます
- 1日に1ユーザに対して2回まで評価を下げられます
質問の評価を下げる
teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。
- プログラミングに関係のない質問
- やってほしいことだけを記載した丸投げの質問
- 問題・課題が含まれていない質問
- 意図的に内容が抹消された質問
- 過去に投稿した質問と同じ内容の質問
- 広告と受け取られるような投稿
評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。
質問の評価を下げたことを取り消します
この機能は開放されていません
評価を下げる条件を満たしてません
質問の評価を下げる機能の利用条件
この機能を利用するためには、以下の事項を行う必要があります。
- 質問回答など一定の行動
-
メールアドレスの認証
メールアドレスの認証
-
質問評価に関するヘルプページの閲覧
質問評価に関するヘルプページの閲覧
checkベストアンサー
+3
まず、numpy.savetxt()
はファイルにデータを追記するのではなく、行列を一気にファイルに書き込みますので、ループで回すと最終結果だけが毎回上書きされてしまいます。
で、対応する方法ですが・・・、
順番に説明します。
現状のコードでは(ソースコードが崩れて記載されているので若干想定が入りますが)、data1
,data2
,data3
にそれぞれ (1,500,000 x 7) の行列が入っているかと思います。
それを
for line in range(3):
と行ごとにループを行い、更にループ内で
dvorx = [((data3[line,4]-data2[line,4])-(data2[line,4]-data1[line,4]))/(2*dt)]
dvory = [((data3[line,5]-data2[line,5])-(data2[line,5]-data1[line,5]))/(2*dt)]
dvorz = [((data3[line,6]-data2[line,6])-(data2[line,6]-data1[line,6]))/(2*dt)]
と5列目(index=4),6列目(index=5),7列目(index=6)それぞれで dvorx
,dvory
,dvorz
のデータを算出しているようです。
しかしながら、この部分は行列演算で使うことで、
dvor = ((data3[:,4:7]-data2[:,4:7])-(data2[:,4:7]-data1[:,4:7]))/(2*dt)
と一行で全データの計算を行うことができます。
data1[:,4:7]
は行列data1
より5列目・6列目・7列目のみを切り出した (1,500,000 x 3) の行列となります。
その結果dvor
には (1,500,000 x 3)の(結果が入った)行列を得ることができます。
よって、この行列を
np.savetxt("output_file.csv", dvor, delimiter=',')
とすることで、全結果が入ったCSVファイルを作成することができます。
【補足】
非常に細かい事ですが
for file in range(0,5):
は、file
に 0,1,2,3,4 の値を順番に入れてループがまわります。
よって
file += 1
は不要です。
line += 1
も同様です。
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
15分調べてもわからないことは、teratailで質問しよう!
- ただいまの回答率 88.32%
- 質問をまとめることで、思考を整理して素早く解決
- テンプレート機能で、簡単に質問をまとめられる
質問への追記・修正、ベストアンサー選択の依頼
pashango2
2017/07/07 18:31
中心差分の計算とファイル出力の部分を見ると、最後の結果を繰り返しファイルに書き込んでいます、コードのインデントがなくなっているので、コードを```(Shift+@)で囲んで下さい