Kronecker product: Difference between revisions

m
C++ - faster Kronecker product calculation
m (C++ - faster Kronecker product calculation)
Line 914:
size_t brows = b.rows();
size_t bcolumns = b.columns();
size_t rows =matrix<scalar_type> c(arows * brows, acolumns * bcolumns);
for (size_t columnsi = acolumns0; *i bcolumns< arows; ++i)
for (size_t j = 0; j < acolumns; ++j)
matrix<scalar_type> c(rows, columns);
for (size_t ik = 0; ik < rowsbrows; ++ik)
for (size_t jl = 0; jl < columnsbcolumns; ++jl)
c(i*brows + k, j*bcolumns + l) = a(i/brows, j/bcolumns) * b(i % browsk, j % bcolumnsl);
return c;
}
1,777

edits