Hi!
I would like to multiply two arrays in Fortran using DGEMM (BLAS procedure).
I have written a simple program:
And I get this result:
Could you point me where I have a mistake?
I will be very grateful for your help.
I would like to multiply two arrays in Fortran using DGEMM (BLAS procedure).
I have written a simple program:
Code:
program matrix
implicit none
double precision mac(m,n),mac2(n,k),mac3(m,k)
integer R,L
integer m,n,k
m=4
n=2
k=3
call FILLMATRIX(m,n,mac)
call FILLMATRIX(n,k,mac2)
call DGEMM("N","N",m,k,n,1.d0,mac,m,mac2,n,0.d0,mac3,m)
WRITE(*,*)'Matrix C:'
DO R=1,m
write(*,*) (mac3(R,L),L=1,k)
end do
END
SUBROUTINE FILLMATRIX (M,N,MATRIX)
INTEGER M,N
DOUBLE PRECISION MATRIX(M,N)
do I=1,M
do J=1,N
MATRIX(I,J) = I+2*J
end do
end do
WRITE(*,*) 'Matrix:'
do I=1,M
write(*,*) (MATRIX(I,J),J=1,N)
end do
END
And I get this result:
Code:
Matrix C:
29.00000000000000 0.0000000000000000 0.0000000000000000
36.00000000000000 0.0000000000000000 0.0000000000000000
27.00000000000000 7.2315834086755357E-308 0.0000000000000000
34.00000000000000 1.0037966311151816E-317 0.0000000000000000
Segmentation fault
I will be very grateful for your help.