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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
subroutine lspdsp(ne,ind,m,n,maxc,mode,ll,lunit,cw)
c!but
c dspdsp visualise une matrice booleenne creuse
c!liste d'appel
c
c subroutine lspdsp(ne,ind,m,n,maxc,mode,ll,lunit,cw)
c
c double precision x(*)
c integer ind(*)
c integer nx,m,n,maxc,mode,ll,lunit
c character cw*(*)
c
c ind : indices specifiant la position des elements vrais
c m : nombre de ligne de la matrice
c n : nombre de colonnes de la matrice
c maxc : nombre de caracteres maximum autorise pour
c representer un nombre
c mode : inutilise
c ll : longueur de ligne maximum admissible
c lunit : etiquette logique du support d'edition
c cw : chaine de caracteres de travail de longueur au moins ll
c!
c Copyright INRIA
integer maxc
integer ind(*)
character cw*(*)
c
if(ne.eq.0) then
write(cw,'(''('',i5,'','',i5,'') False sparse matrix'')') m,n
call basout(io,lunit,cw(1:33))
call basout(io,lunit,' ')
goto 99
else
write(cw,'(''('',i5,'','',i5,'') sparse matrix'')') m,n
call basout(io,lunit,cw(1:27))
call basout(io,lunit,' ')
if(io.eq.-1) goto 99
endif
ilr=1
ilc=m+1
cw=' '
i0=0
i1=i0
l=1
do 20 k=1,ne
cw=' '
11 i0=i0+1
if(i0-i1.gt.ind(l)) then
i1=i0
l=l+1
goto 11
endif
i=l
j=ind(ilc-1+k)
write(cw,'(''('',i5,'','',i5,'') T'')') i,j
l1=19
call basout(io,lunit,cw(1:l1) )
if (io.eq.-1) goto 99
20 continue
99 continue
c
end
|
