43 call
crcrot1(2,t(j,j),c(j),s(j))
47 call
crcrot2(j+1,t(1,j),t(1,j+1),c(j),s(j))
64 call
crcrot2(n,u(1,j),u(1,j+1),c(j),s(j))
78 t = x(i)*c(i) - x(i+1)*
cmplx(0,s(i))
79 x(i+1) = x(i+1)*c(i) - x(i)*
cmplx(0,s(i))
92 t = x(i)*c + y(i)*
cmplx(0,s)
93 y(i) = y(i)*c + x(i)*
cmplx(0,s)
subroutine crsf2csf(n, t, u, c, s)
subroutine crcrot2(n, x, y, c, s)
subroutine crcrot1(n, x, c, s)
octave_value sqrt(void) const