私は、2つの正方行列のクロネッカー積を計算する手書きのメソッドを確認するための次のコードを持っています。最初の部分は実際に私の繰り返しとタイリングの方法がそれぞれ同じ出力をもたらすことを確認しています。aとb import pyopencl as cl
import numpy
from time import time
N = 3
num_iter = 1
a = numpy.ra
私は行列回ベクトルプログラムを作成しようとしましたが、私はプログラムを実行すると、私は出力が 少し奇妙であることがわかった。ここに私のコードは #include <stdio.h>
main()
{
int i, j, n, a[3][3], b[3], c[3];
n = 3;
printf("Matrix A\n");
for (i = 0;
密行列Aとベクトルe、fのA * diag（e）* A^T * fの高速乗算の提案はありませんか？ これは私が今持っているものです。 v[:] = 0
for i in range(N):
for j in range(N):
v[i] = v[i]+A[i,j]*e[j]*np.dot(A[:,j],f)
おかげで、@rubenvbによって作ら