660 DOUBLE PRECISION sfac
663 COMPLEX*16 ccomp(len), csize(len), ctrue(len)
667 DOUBLE PRECISION scomp(20), ssize(20), strue(20)
671 INTRINSIC dimag, dble
674 scomp(2*i-1) = dble(ccomp(i))
675 scomp(2*i) = dimag(ccomp(i))
676 strue(2*i-1) = dble(ctrue(i))
677 strue(2*i) = dimag(ctrue(i))
678 ssize(2*i-1) = dble(csize(i))
679 ssize(2*i) = dimag(csize(i))
682 CALL stest(2*len,scomp,strue,ssize,sfac)
subroutine stest(LEN, SCOMP, STRUE, SSIZE, SFAC)