613
  614
  615
  616
  617
  618      DOUBLE PRECISION SFAC
  619      INTEGER          LEN
  620
  621      COMPLEX*16       CCOMP(LEN), CSIZE(LEN), CTRUE(LEN)
  622
  623      INTEGER          I
  624
  625      DOUBLE PRECISION SCOMP(20), SSIZE(20), STRUE(20)
  626
  628
  629      INTRINSIC        dimag, dble
  630
  631      DO 20 i = 1, len
  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))
  638   20 CONTINUE
  639
  640      CALL stest(2*len,scomp,strue,ssize,sfac)
 
  641      RETURN
subroutine stest(len, scomp, strue, ssize, sfac)