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
31
32
33
34
|
subroutine dmmul(a,na,b,nb,c,nc,l,m,n)
c!but
c ce sous programme effectue le produit matriciel:
c c=a*b .
c!liste d'appel
c
c subroutine dmmul(a,na,b,nb,c,nc,l,m,n)
c double precision a(na,m),b(nb,n),c(nc,n)
c integer na,nb,nc,l,m,n
c
c a tableau de taille na*m contenant la matrice a
c na nombre de lignes du tableau a dans le programme appel
c b,nb,c,nc definitions similaires a celles de a,na
c l nombre de ligne des matrices a et c
c m nombre de colonnes de a et de lignes de b
c n nombre de colonnes de b et c
c!sous programmes utilises
c neant
c!
double precision a(*),b(*),c(*)
double precision ddot
integer na,nb,nc,l,m,n
integer i,j,ib,ic
c
ib=1
ic=0
do 30 j=1,n
do 20 i=1,l
20 c(ic+i)=ddot(m,a(i),na,b(ib),1)
ic=ic+nc
ib=ib+nb
30 continue
return
end
|
