◆ kron()

template<typename T , std::size_t M, std::size_t N, std::size_t P, std::size_t Q>

constexpr matrix<T, M * P, N * Q> cotila::kron	(	const matrix< T, M, N > &	a,
		const matrix< T, P, Q > &	b
	)

Parameters

a	an $M \times N$ matrix
b	an $P \times Q$ matrix

Returns: An $MP \times NQ$ matrix $\textbf{a}\otimes\textbf{b}$ of type T such that $\left(\textbf{a}\otimes\textbf{b}\right)_{ij} = \textbf{a}_{\lfloor i/P \rfloor,\lfloor j/Q \rfloor}\textbf{b}_{i\textrm{%}P,j\textrm{%}Q}$ where $i \textrm{%} P$ is the remainder of

Computes the kronecker tensor product of two matrices.