From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- vo/votools/gasplib/varsvd.f | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) create mode 100644 vo/votools/gasplib/varsvd.f (limited to 'vo/votools/gasplib/varsvd.f') diff --git a/vo/votools/gasplib/varsvd.f b/vo/votools/gasplib/varsvd.f new file mode 100644 index 00000000..81e800f3 --- /dev/null +++ b/vo/votools/gasplib/varsvd.f @@ -0,0 +1,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 -- cgit