vo/votools/gasplib/varsvd.f


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

subroutine varsvd (v,ma,w,cvm,ncvm)
	parameter (mmax=20)
	real*8 v(ma,ma),w(ma),cvm(ncvm,ncvm),wti(mmax)
	real*8 sum

	do i=1,ma
		wti(i)=0.
		if(w(i).ne.0.0d0) wti(i)=1./(w(i)*w(i))
	enddo
	do i=1,ma
		do j=1,i
			sum=0.
			do k=1,ma
				sum=sum+v(i,k)*v(j,k)*wti(k)
			enddo
			cvm(i,j)=sum
			cvm(j,i)=sum
		enddo
	enddo
	return
	end