714
715
716
717
718
719 DOUBLE PRECISION SFAC
720 INTEGER LEN
721
722 COMPLEX*16 CCOMP(LEN), CSIZE(LEN), CTRUE(LEN)
723
724 INTEGER I
725
726 DOUBLE PRECISION SCOMP(20), SSIZE(20), STRUE(20)
727
729
730 INTRINSIC dimag, dble
731
732 DO 20 i = 1, len
733 scomp(2*i-1) = dble(ccomp(i))
734 scomp(2*i) = dimag(ccomp(i))
735 strue(2*i-1) = dble(ctrue(i))
736 strue(2*i) = dimag(ctrue(i))
737 ssize(2*i-1) = dble(csize(i))
738 ssize(2*i) = dimag(csize(i))
739 20 CONTINUE
740
741 CALL stest(2*len,scomp,strue,ssize,sfac)
742 RETURN
743
744
745
subroutine stest(len, scomp, strue, ssize, sfac)