[Fortran 95] Prodotto righe per colonne
December 29, 2008 – 1:44 amDescrizione
Calcola la matrice risultante dal prodotto righe per colonne di due matrici a scelta dell’utente.
Codice
program prodmatr
real :: a(100,100), b(100,100), g(100,100)
integer :: m,n,p,i,j
read*, n,p,m
do i=1,n
do j=1,m
read*, A(i,j)
enddo
enddo
do i=1,m
do j=1,p
read*, B(i,j)
enddo
enddo
call prod(a,b,g,m,n,p)
do i=1,n
do j=1, p
print*, g(i,j)
enddo
enddo
endsubroutine prod(a,b,g,m,n,p)
real :: a(100, m), b(100, p), g(100, p)
integer :: i,j,l,m,n,p
do i = 1, n
do j = 1, p
g(i,j)=0
do l = 1, m
g(i,j)=g(i,j)+a(i,l)*b(l,j)
enddo
enddo
enddo
end
![[Valid RSS]](http://v3rgil.altervista.org/Images/valid-rss.png "Validate my RSS feed")
2 Responses to “[Fortran 95] Prodotto righe per colonne”
io ho lo stesso algoritmo, ma tanto per cambiare non ho capito come funziona! è possibile che solo con questa cosa g(i,j)=g(i,j)+a(i,l)*b(l,j) si faccia il prodotto? non riesco a capacitarmene.
By Mary on Jan 28, 2009
Beh… il “prodotto” è dato da questa cosa “g(i,j)=g(i,j)+a(i,l)*b(l,j)” innestata nei tre cicli
.
By V3rgil on Jan 28, 2009