diff options
Diffstat (limited to 'vo/votools/gasplib/ksbsvd.f')
| -rw-r--r-- | vo/votools/gasplib/ksbsvd.f | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/vo/votools/gasplib/ksbsvd.f b/vo/votools/gasplib/ksbsvd.f new file mode 100644 index 00000000..d93609e9 --- /dev/null +++ b/vo/votools/gasplib/ksbsvd.f @@ -0,0 +1,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 |