shupdn.f

gsac程序包

	subroutine shupdn(s,saa)
	complex*16 s(2), saa(2,2)
	complex*16 tmp
		tmp  = saa(1,1)*s(1) + saa(1,2)*s(2)
		s(2) = saa(2,1)*s(1) + saa(2,2)*s(2)
		s(1) = tmp
	return
	end