Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | A(:,:) | |||
integer, | intent(in) | :: | B(:,:) |
function i_kronecker_product(A,B) result(AxB) integer,intent(in) :: A(:,:), B(:,:) integer :: i,j integer :: rowA,colA integer :: rowB,colB integer :: AxB(size(A,1)*size(B,1),size(A,2)*size(B,2)) AxB = 0 rowA=size(A,1) ; colA=size(A,2) rowB=size(B,1) ; colB=size(B,2) forall(i=1:rowA,j=1:colA) AxB(1+rowB*(i-1):rowB*i,1+colB*(j-1):colB*j) = A(i,j)*B(:,:) end forall end function i_kronecker_product