在FORTRAN

Question

乘以兩個數組我試圖乘用FORTRAN兩個數組，下面是我的代碼

program multiplyarray 
    implicit none 
! Declare Variables 
    integer, dimension(2,3 :: first_array 
    integer, dimension(3,2) ::second_array 
    integer, dimension(2,2) :: resultant_array 
    integer:: i, j, a, b, k, product 


! Generating values for array 
    a = 2 
    do i = 1,2 
     do j = 1,3 

     first_array(i,j) = a 
     a = a + 2 
     end do 
    end do 

    print*,first_array 

    do i = 1,3 
     do j = 1,2 

      b = 0 
      second_array(i,j) = b 
      b = b + 3 
     end do 
    end do 

! Multiplying both arrays 
     do i = 1,2 
      do j = 1,2 
       do k = 1,3 
        product = product + first_array(i,k) * second_array(k,j) 
       end do 
        resultant_array (j,i) = product 
      end do 
     end do 

! Displaying the resultant matrix 
     do i = 1,2 
      do j = 1,2 
        print *,resultant_array(i,j) 
      end do 
     end do 

end program multiplyarray 

可能是什麼錯誤？我沒有得到個體矩陣，所以我可以乘他們得到結果矩陣。第一矩陣是0，而第二矩陣的結果也爲0

Answer 1

嘗試使用MATMUL功能

MATMUL（first_array，second_array）= resultant_array 打印*，resultant_array

Answer 2

第二個數組將始終爲零

b = 0 
    second_array(i,j) = b 
    b = b + 3 

你總是設置b爲零。行b = b + 3不起作用，因爲無論如何b始終設置爲零。

當second_array爲零時，產品也將爲零。

---

由於HighPerformanceMark建議，您可以輕鬆地乘用matmul或更不容易，但很快使用BLAS庫兩個矩陣。

result = matmul(first_array, second_array)