DataMuseum.dk

top - metrics - download

    Length: 6144 (0x1800)
    Types: TextFile
    Names: »gonimtridql«

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦1248b0c55⟧ »gobib« 
            └─⟦this⟧ 

TextFile

;gosav time.90 lines.5000
nimtridql=set 80
permanent nimtridql.15
dqltxt=set 20
dqltxt=edit gonimfreq0
l./zone ev/,d./nev:=/,l./in,khi)/,d,l./<:khi/,d1,i/
      <:deltaM:>,<<ddddd>,dM,nl,1);
/,l2,d./u1(/,i/
read(in,J,gamma,m0);
if J<0 then goto stop;
/,l./freq:= 0;/,d./N:= ((J/,i?
Jeven:= J//2*2=J; JJ:=(J+1)*J;
Mmax:= m0+dM; pm:= Mmax//3;
if pm*3=Mmax then begin
  plus:= gamma=0;
  write(res,nl,2,<:J,  gamma, Mmax:>,<<dddd>,J,gamma,m0);
  Mmin:= 0; meven:= true; M1:= J; M2:= 1-gamma;
end else begin
  write(res,nl,2,<:J, Mmax:>,<<dddd>,J,m0);
  Mmin:= Mmax-pm*3-dM;
  meven:= Mmin//2*2=Mmin;
  if meven then begin
    M1:= J; M2:= (J+2)//2; 
  end else begin
    M1:= J+1; M2:= M1//2
end end;
N:= ((J*2+1)*(Mmax-Mmin)//3+M1)//2+M2;
?,l./m:= if m0=0/,s,i?
tridql(N,u,h);
m0:= if Mmin=0 then 0 else Mmin+dM;
for m:= m0 step 3 until Mmax-dM do begin
  Kmax:= if Jeven==m//2*2=m then J else J-1;
  k:= if m=0 then 1 else index(m,-Kmax);
  l:= index(m,Kmax); j:= 0;
  k01:= if m=0 then gamma*2 else -Kmax;
  for k0:= k01 step 2 until Kmax do begin
    p:= 0;
    for j:=j+1 while p<=0.5 and j<=N do begin
      p:= 0;
      for i:=k step 1 until l do p:= p+h(j,i)**2;
    end;
    if p>0.5 then begin
      j:= j-1; go:= u(j);
      write(res,nl,1,<<-ddd>,m,k0,<<  -dddddddd.ddd>,go+XJ);
      if k0>k01 then write(res,<<  -dddddddd.ddd>,go-gu);
      gu:= go;
end end end;
end;  goto rep end;
stop: write(res,<:<25>:>); close(res,closeres) end;
?,f
nimtridql=algol dqltxt index.no list.yes
▶EOF◀