scipyの疎行列モジュール scipy.sparse を用いて行列A(n*n) とB(n*n)の要素同士の積を計算したいです。目標としては、C = A○Bとすると,cij = aij * bij のようにしたいのですが,なぜかできず,numpyと同様の操作をするとどうしても内積になってしまいます.
初歩だとは思うのですが,なかなか検索で探し当てることができず,教えていただけると助かります.
具体的には,
===
print("lil_matrix")
A = sp.lil_matrix([[1,1,1],[1,1,1],[1,1,1]])
B = sp.lil_matrix([[0,1,0],[0,1,0],[0,1,0]])
C = A*B
print(A.toarray())
print(B.toarray())
print(C.toarray())
===
として,本来であれば
C = [[0,1,0],[0,1,0],[0,1,0]] となってほしいところなのですが,
結果には
===
lil_matrix
[[1 1 1]
[1 1 1]
[1 1 1]]
[[0 1 0]
[0 1 0]
[0 1 0]]
[[0 3 0]
[0 3 0]
[0 3 0]]
===
と表示されてしまいます.
CSR格納形式の際も同様で,
===
print("csr_matrix")
A_csr = A.tocsr()
B_csr = B.tocsr()
C_csr = A_csr*B_csr
print(A_csr.toarray())
print(B_csr.toarray())
print(C_csr.toarray())
===
の結果が,
===
csr_matrix
[[1 1 1]
[1 1 1]
[1 1 1]]
[[0 1 0]
[0 1 0]
[0 1 0]]
[[0 3 0]
[0 3 0]
[0 3 0]]
===
ちなみに,np.dot(A,B)を用いても同様の結果となるので,要素同士の積をどのように計算したらいいのかわからず,基本的なところで詰まってしまっています.
どなたか教えていただけると助かります。
No.2ベストアンサー
- 回答日時:
A.multiply(B)
No.1
- 回答日時:
> numpyと同様の操作をする
numpy.matrixでも A*B は0 3 0;0 3 0 ;0 3 0 になります。
行列の積、ということでは、こちらの方が「普通」でしょう。
> 検索で探し当てることができず
http://docs.scipy.org/doc/scipy/reference/genera …
は読みましたか?
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
関連するカテゴリからQ&Aを探す
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
【至急!!!】python言語で本を見...
-
16進の10進変換について
-
VBAのNext forでエラーが出て困...
-
パイソン文法で(インデント)
-
数値かどうかの判定方法
-
変数名に変数を使用
-
シェルスクリプト+perl
-
pythonの*
-
ランダム文字列
-
Perl でワードファイルのヘッダ...
-
Perlにおける<?phpinfo();?>の...
-
シェルスクリプトで、空白(ス...
-
パイソンのプログラムで contin...
-
JSONで文字列が長い時
-
Application.ScreenUpdating = ...
-
FindFirst を複数条件で検索
-
テキスト入力後、エンターを押...
-
文字の横にプルダウンを表示さ...
-
formで特定のinputを送信しない...
-
sendmailで複数アドレスに送信...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
16進の10進変換について
-
pythonで演算子を変数に代入す...
-
【至急!!!】python言語で本を見...
-
pythonの*
-
c言語 16進数の2進数への変換
-
awkの正規表現での最左最短マッ...
-
数値かどうかの判定方法
-
関数「exists」と「defined」の...
-
シェルで、ランダム関数を変数...
-
I2C接続のLCDディスプレイを使う
-
シェルスクリプトで、空白(ス...
-
perでメールが送信出来ない場合...
-
UWSCのSELECT文の記述方法
-
int(input("○○"))の使い方
-
switch文のエラーについて
-
pythonで改行を含めてfindする方法
-
シェルスクリプト+perl
-
PythonでUTF8の文字列をリスト...
-
少数点以下の桁数について
-
qq!!の意味を教えてください。
おすすめ情報