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