shupdn.f

来自「gsac程序包」· F 代码 · 共 11 行

	subroutine shupdn(s,saa)
	implicit none
	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