DataMuseum.dk

top - metrics - download

    Length: 768 (0x300)
    Types: TextFile
    Names: »tfft1«

└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
    └─⟦b2ec5d50f⟧ 
        └─⟦09b4e9619⟧ »thcømat« 
            └─⟦this⟧ 

TextFile

;klab3 1
lookup allist
if ok.yes
mode 10.yes
slet fourint.3
beskyt proc.fourint.3.2
if 10.yes
fourint=algol list.yes
fourint=algol
external real procedure fourint(RE,IM,n,si,dx);
value n,si,dx;
integer n,si; real dx;
array RE,IM;
begin
integer i,N;
real f0,ft,af,g0,gt,ag,dy,y,T;
N:=2**n-1;
T:=N*dx;
f0:=RE(0); ft:=RE(N); af:=(ft-f0)/N;
g0:=IM(0); gt:=IM(N); ag:=(gt-g0)/N;
for i:=0 step 1 until N do begin
  RE(i):=RE(i)-af*i-f0;
  IM(i):=IM(i)-ag*i-g0;
  end;
fft(RE,IM,n,si);
RE(0):=RE(0)*dx;
IM(0):=IM(0)*dx;
fourint:=dy:=2*pi/dx/N;
for i:=1 step 1 until N do begin
  y:=i*dy;
  RE(i):=RE(i)*dx;
  IM(i):=IM(i)*dx+si*(ft*(0/y-T*1)-f0*(-T+0/y))/y/T;
  end;
end;
end;
▶EOF◀