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)