最快的Hadamard产品代码

Question

有两个维度为2的复杂数组我想计算一个点乘法（Hadamard乘积）：

complex(8) :: A(N,N), B(N,N), C(N,N)
...
do j = 1, N
  do i = 1, N
    C(i,j) = A(i,j)*B(i,j)
  enddo
enddo

是否有任何BLAS例程来优化它，或者这实际上已经是编写Hadamard产品的最有效方式了吗？或者编译器是否在这么简单的情况下为我做了工作？

我在Fortran中编码，因此第一个索引是快速索引。