svupdn.f

gsac程序包

	subroutine svupdn(s,haa)
	implicit none
	complex*16 s(4), haa(4,4)
	complex*16 t(4), tmp
	integer i, j
	do 100 i=1,4
		tmp = dcmplx(0.0d+00,0.0d+00)
		do 101 j=1,4
			tmp = tmp + haa(i,j)*s(j)
  101		continue
		t(i) = tmp
  100	continue
	do 200 i=1,4
		s(i) = t(i)
  200	continue
	return
	end