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
35
36
37
38
39
40
41
42
43
44
45
46
|
C/MEMBR ADD NAME=WMMUL,SSI=0
subroutine wmmul(ar,ai,na,br,bi,nb,cr,ci,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 wmmul(ar,ai,na,br,bi,nb,cr,ci,nc,l,m,n)
c double precision ar(*),ai(*),br(*),bi(*),cr(*),ci(*)
c integer i,j,k,ia,ib,ic
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 ar(*),ai(*),br(*),bi(*),cr(*),ci(*)
double precision tr,ti
integer na,nb,nc,l,m,n
integer i,j,k,ia,ib,ic
c
ib=0
ic=0
do 30 j=1,n
do 20 i=1,l
ia=i
tr=0.0d+0
ti=0.0d+0
do 10 k=1,m
ib=ib+1
tr=tr+ar(ia)*br(ib)-ai(ia)*bi(ib)
ti=ti+ar(ia)*bi(ib)+ai(ia)*br(ib)
ia=ia+na
10 continue
ib=ib-m
cr(ic+i)=tr
20 ci(ic+i)=ti
ic=ic+nc
ib=ib+nb
30 continue
return
end
|
