Discussion :

[quiqui191]

Bonjour,
j'ai un problème : le produit matrice-vecteur en fortran.
Voila mon code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
program produit
implicit none
integer,parameter :: n=4
integer,parameter :: dp=kind(0.d00)
real(kind=dp), dimension (n,n):: M
real(kind=dp), dimension (n):: xvec=(/1,2,3,4/)
real(kind=dp), dimension (n):: y
real(kind=dp), 
integer :: i,j
call produit_matrice_vecteur(M,x,y,n)
end program produit 
 
 
subroutine produit_matrice_vecteur(mat,a,b,ndim)
implicit none
integer,parameter :: dp =kind(0.d00)
integer, intent(in) :: ndim
real(kind(0.0d00)), dimension(ndim,ndim):: mat 
real(kind(0.0d00)), dimension(ndim):: a
real(kind(0.0d00)), intent(out) :: b
integer :: i,prod
prod=0.0d00
prod=prod+b
do i=1,ndim
prod=dot_product(mat(i,:),a)
end do
write(*,*)"prod=b",b
 
return 
end subroutine produit_matrice_vecteur

[bobbyboy]

Bonjour,

Un petit effort s'il te plait, juste en relisant ton programme (ou au moins les erreurs du compilateur, tu peux trouver que la 8ème ligne de ton programme est clairement fausse et ne sert à rien ... a savoir :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
real(kind=dp),

Ensuite il y a multes incohérences dans ta subroutine :

1) je suppose que dans le main, la variable "y" est le vecteur contenant le résultat de la multiplication, mais l'argument correspondant dans ta subroutine, "b", n'est pas déclaré comme un tableau.

2) ensuite à quoi sert la variable "prod" ? Tu la déclares comme un entier mais tu lui affectes un réel double :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
integer :: i,prod
prod=0.0d00

3) Outre le problème de déclaration, regarde ce que tu fais de cette variable prod :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
prod=0.0d00
prod=prod+b

en gros tu fais : prod = 0 et apres : prod = 0 + b ... n'est-ce pas plus simple de faire simplement : prod = b ?

4) Surtout qu'après, dans la boucle, la valeur de prod est ndim fois écrasée ...

S'agissant apparement d'un travail scolaire, je t'invite à regarder un peu plus attentivement ton code et te laisse le soin de le débugger.