1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
program intrinsics_109
integer, parameter :: N = 2
real :: A(N, N), B(N, N), C(N, N), D(N,N), E(N,N)
! Initialize matrices A and B as 2x2 identity matrices
A = reshape([1.0, 0.0, 0.0, 1.0], shape(A))
B = reshape([1.0, 0.0, 0.0, 1.0], shape(B))
! Perform matrix multiplication: C = A * B
C = matmul( matrix_a = A, matrix_b = B )
D = matmul( A, matrix_b = B )
E = matmul( A, B )
! Print the matrices
print *, "Matrix A (2x2 identity matrix):"
print *, A
print *, "Matrix B (2x2 identity matrix):"
print *, B
print *, "Matrix C (result of A * B):"
print *, C
end
|
