两个矩阵的直积

hongcha

[Fortran] 纯文本查看 复制代码

subroutine product_AB(H,A,B,m,n,r,s,index1,index2)
  implicit none
  integer(8)::m,n,r,s,ii,jj,kk,ll,row,column,index1,index2       
  complex(8)::A(m,n),B(r,s),H(index1,index2)
  row=0
  do ii=1,m
    do kk=1,r             
      row=row+1         
      column=0
      do jj=1,n
        do ll=1,s
          column=column+1
          h(row,column)=a(ii,jj)*b(kk,ll)
        enddo
      end do
    enddo
  enddo
end subroutine product_AB