618 DOUBLE PRECISION sfac
621 COMPLEX*16 ccomp(len), csize(len), ctrue(len)
625 DOUBLE PRECISION scomp(20), ssize(20), strue(20)
629 INTRINSIC dimag, dble
632 scomp(2*i-1) = dble(ccomp(i))
633 scomp(2*i) = dimag(ccomp(i))
634 strue(2*i-1) = dble(ctrue(i))
635 strue(2*i) = dimag(ctrue(i))
636 ssize(2*i-1) = dble(csize(i))
637 ssize(2*i) = dimag(csize(i))
640 CALL stest(2*len,scomp,strue,ssize,sfac)
subroutine stest(LEN, SCOMP, STRUE, SSIZE, SFAC)