fold.f

来自「这是一个希尔伯特变换的程序包」· F 代码 · 共 11 行

      SUBROUTINE FOLD(LA,A,LB,B,LC,C)
      DIMENSION A(*),B(*),C(*)
      LC=LA+LB-1
      CALL MOVE(0,C,0,LC*4)
      DO 10 I=1,LA
      DO 10 J=1,LB
      K=I+J-1
   10 C(K)=A(I)*B(J)+C(K)
      RETURN
      END