26 lines
595 B
Fortran
26 lines
595 B
Fortran
SUBROUTINE dsvbksb(u,w,v,m,n,mp,np,b,x)
|
|
INTEGER m,mp,n,np,NMAX
|
|
DOUBLE PRECISION b(mp),u(mp,np),v(np,np),w(np),x(np)
|
|
PARAMETER (NMAX=500)
|
|
INTEGER i,j,jj
|
|
DOUBLE PRECISION s,tmp(NMAX)
|
|
do 12 j=1,n
|
|
s=0.0d0
|
|
if(w(j).ne.0.0d0)then
|
|
do 11 i=1,m
|
|
s=s+u(i,j)*b(i)
|
|
11 continue
|
|
s=s/w(j)
|
|
endif
|
|
tmp(j)=s
|
|
12 continue
|
|
do 14 j=1,n
|
|
s=0.0d0
|
|
do 13 jj=1,n
|
|
s=s+v(j,jj)*tmp(jj)
|
|
13 continue
|
|
x(j)=s
|
|
14 continue
|
|
return
|
|
END
|