1
scipy.sparse.csr.csr_matrix
を指定すると、対角要素を素早く返すことができますか?私はこれをやりたいscipyスパース行列の対角要素を返す
理由はD
は対角エントリ(A
が対角線上の非ゼロを持つことが保証私のスパース行列、である)A
に同意する対角行列であるinv(D) A
を計算することです。
scipy.sparse.csr.csr_matrix
を指定すると、対角要素を素早く返すことができますか?私はこれをやりたいscipyスパース行列の対角要素を返す
理由はD
は対角エントリ(A
が対角線上の非ゼロを持つことが保証私のスパース行列、である)A
に同意する対角行列であるinv(D) A
を計算することです。
は、行列の主対角返し
例:
>>> import numpy as np
>>> from scipy.sparse import csr_matrix
>>> mymat = csr_matrix((4, 4), dtype=np.int8)
>>> mymat.diagonal()
array([0, 0, 0, 0], dtype=int8)
*ため息*。あなたがグーグルで完全に失敗したとき:-)答えてくれてありがとう! +1(10分のロックアウト後に受け入れます) – parsiad
あなたの例では、 'mymat'は密な配列です(その' toarray() ':))。まばらな対角線は密なものとはまったく一般的ではありません。 – hpaulj
@hpauljの提案に感謝します。私はウォーレンの提案でポストを更新しました。両方の人が私の答えを今好きになることを願っています! :)ところで、私はrefによって動かされましたが、あなたが言うことは理にかなっています。 – gsamaras