質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.48%
NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

1回答

14029閲覧

numpyでnp.linalg.solveをつかい、計算したいのですが、エラーが発生してしまいます。

date-07

総合スコア13

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2018/01/30 05:45

編集2018/01/30 05:51

python3.6.3でnumpyを使っています。

Decimal(0.295535291)このような数値をfloat64にcastし、np.linalg.solveを使い、結果を出したいのですが、
TypeError: No loop matching the specified signature and casting
was found for ufunc solve1
とエラーが出てしまいます。
cast時に数値が、-0.70446470900000002このようになってしまっています。これが原因でしょうか。

質問は、このエラーはなぜ発生してしまっているのか、Decimalからfloat64へcast時、正しい値へcastするにはどうしたらいいか。の二つです。

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

TypeError: No loop matching the specified signature and casting was found for ufunc solve1

該当のソースコード

val1 = np.array([ [Decimal(-0.704464709), Decimal(0.295535291), Decimal(0.295535291), Decimal(0.295535291), Decimal(0.295535291), Decimal(0.295535291), Decimal(0.295535291), Decimal(1)], [Decimal(0.0000822818198462602), Decimal(-1), Decimal(0), Decimal(0), Decimal(0), Decimal(0), Decimal(0), Decimal(1)], [Decimal(0.063900458), Decimal(0.06398274), Decimal(-0.93601726), Decimal(0.06398274), Decimal(0.06398274), Decimal(0.06398274), Decimal(0.06398274), Decimal(1)], [Decimal(0.185757364), Decimal(0.185757364), Decimal(0.185757364), Decimal(-0.814242636), Decimal(0.185757364), Decimal(0.185757364), Decimal(0.185757364), Decimal(1)], [Decimal(0.030625048), Decimal(0), Decimal(0), Decimal(0), Decimal(-1), Decimal(0), Decimal(0), Decimal(1)], [Decimal(0.420372803), Decimal(0.450997851), Decimal(0.450997851), Decimal(0.450997851), Decimal(0.450997851), Decimal(-0.549002149), Decimal(0.450997851), Decimal(1)], [Decimal(0.003726755), Decimal(0.003726755), Decimal(0.003726755), Decimal(0.003726755), Decimal(0.003726755), Decimal(0.003726755), Decimal(-0.996273245), Decimal(1)], [Decimal(1), Decimal(1), Decimal(1), Decimal(1), Decimal(1), Decimal(1), Decimal(1), Decimal(0)] ]) val2 = np.array( [0, 0, 0, 0, 0, 0, 0, 1] ) for i, data in enumerate(val1): val1[i] = [np.float64(x) for x in data] ans = np.linalg.solve(val1, val2) print(ans)

試したこと

Decimalをfloat64にcastしたとき、桁数が多くなっていることが原因かと思い、roundを使い四捨五入してみましたが同じエラーでした。
型が違うのかと思い、四捨五入後にfloat32にcast、float64にcastなども試しましたが同じエラーでした。

val1 = np.array([ [-0.704464709, 0.295535291, 0.295535291, 0.295535291, 0.295535291, 0.295535291, 0.295535291, 1], [8.23E-05, -1, 0, 0, 0, 0, 0, 1], [0.063900458, 0.06398274, -0.93601726, 0.06398274, 0.06398274, 0.06398274, 0.06398274, 1], [0.185757364, 0.185757364, 0.185757364, -0.814242636, 0.185757364, 0.185757364, 0.185757364, 1], [0.030625048, 0, 0, 0, -1, 0, 0, 1], [0.420372803, 0.450997851, 0.450997851, 0.450997851, 0.450997851, -0.549002149, 0.450997851, 1], [0.003726755, 0.003726755, 0.003726755, 0.003726755, 0.003726755, 0.003726755, -0.996273245, 1], [1, 1, 1, 1, 1, 1, 1, 0] ]) val2 = np.array( [0, 0, 0, 0, 0, 0, 0, 1] ) ans = np.linalg.solve(val1, val2)

最初から、np.arrayに入力した場合はエラーにはなりませんでした。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

コード詳細は確認できていませんがval1 = val1.astype(np.float64)で変換でよいです。

投稿2018/01/30 06:25

can110

総合スコア38266

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

date-07

2018/01/30 06:29

ありがとうございます! できました!
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問