diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /noao/astcat/src/awcs/ksbsvd.f | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'noao/astcat/src/awcs/ksbsvd.f')
| -rw-r--r-- | noao/astcat/src/awcs/ksbsvd.f | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/noao/astcat/src/awcs/ksbsvd.f b/noao/astcat/src/awcs/ksbsvd.f new file mode 100644 index 00000000..3c78ec23 --- /dev/null +++ b/noao/astcat/src/awcs/ksbsvd.f @@ -0,0 +1,27 @@ +C This routines was copied from the stsdas$pkg/analysis/gasp/gasplib/ +C directory. See the file stsdas$copyright.stsdas for copyright +C restrictions. + 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 |