vo/votools/gasplib/ksbsvd.f


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

subroutine ksbsvd (u,w,v,m,n,b,x)
	parameter (nmax=1000)
	real*8 u(m,n),w(n),v(n,n),b(m),x(n),tmp(nmax)
	real*8 s

	do j=1,n
		s=0.
		if(w(j).ne.0.)then
			do i=1,m
				s=s+u(i,j)*b(i)
			enddo
			s=s/w(j)
		endif
		tmp(j)=s
	enddo
	do j=1,n
		s=0.
		do jj=1,n
			s=s+v(j,jj)*tmp(jj)
		enddo
		x(j)=s
	enddo
	return
	end