-1

PythonでBigFloatライブラリを使用してexpを計算しています。しかし、BigFloats の逆行列を計算する必要があります。関数を使用していますnumpy.linalg.invが、次のエラーが発生します。

ufunc inv の指定された署名とキャストに一致するループが見つかりませんでした

したがって、 BigFloat を他の型に変換する必要があると思います。どうすればいいですか?

4

1 に答える 1

0

Bigfloat 値の 2D リスト (exp 付き) を考えると、次のようになります。

from bigfloat import *
import numpy as np

row = 4
col = 4
mat = []
with precision(1000):
    for i in range(0,row):
        row = []
        for j in range(0,col):
            row.append(exp(1./(i+j+1)))
        mat.append(row)
print "A single element in bigfloat: \n",mat[0][0]
mat_numpivized = np.array([[np.float(col) for col in row] for row in mat])
print "\nNumpived matrix of bigfloat: \n",mat_numpivized
inversed_mat = np.linalg.inv(mat_numpivized)
print "\nInverse of numpivized bigfloat matrix: \n", inversed_mat

結果:

A single element in bigfloat: 
    2.71828182845904523536028747135266249775724709369995957496696762772407663035354759457138217852516642742746639193200305992181741359662904357290033429526059563073813232862794349076323382988075319525101901157383418793070215408914993488416750924476146066808226480016847741185374234544243710753907774499206950

Numpived matrix of bigfloat: 
[[ 2.71828183  1.64872127  1.39561243  1.28402542]
 [ 1.64872127  1.39561243  1.28402542  1.22140276]
 [ 1.39561243  1.28402542  1.22140276  1.18136041]
 [ 1.28402542  1.22140276  1.18136041  1.15356499]]

Inverse of numpivized bigfloat matrix: 
[[    3.82106765   -30.68951855    61.37113204   -34.60880961]
 [  -30.68951855   385.05574805  -890.19682475   538.10687764]
 [   61.37113204  -890.19682475  2211.11911647 -1390.16166038]
 [  -34.60880961   538.10687764 -1390.16166038   893.29628089]]
于 2016-10-23T16:05:20.447 に答える