DataMuseum.dk

top - metrics - download

    Length: 3072 (0xc00)
    Types: TextFile
    Names: »algrotcor«

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

TextFile

;gosav
algcor=set 30
algcor=edit algrotfrekv
d./GOS/,i/
æ12æ


rotcor

CALCULATION OF Pg-MATRIX ELEMENTS.
GOS.  1-2-1980.
/,l2,d1,m e
i?
comment
    The program calculates the square of J-normed Pg-matrix-
    elements defined by
          pg2 = (2*norm(<R1! Pg !R2>)/(J*(J+1)))**2.
    The results are written in binary form on the disc file
    <starkmat>.
    As input is given:
?,l./m= 100/,d./appropriate/,l./6)/
d./8)/,r/9/6/,l1,d1,i/
    7)  The smallest value of pg2 to be considered.
/,l1,r/11/ 8/,l3,d1,i/
    9)  Three integers = 0 or 1, indicating, by 1 for yes, if
        elements of Pa2, Pb2 or Pc2 must be calculated.
/,l1,d./;/,i/
   10)  An integer n <= 200, and n set of quantum numbers
        in order of increasing J. Only matrix elements con-
        nected with the levels specified in this way
        are transferred to the disc file <starkmat>;
/,l./integer array/,r/;/, qtn(1:200);/
l./procedure overskrift/,d./end overskrift/
l./centrifugal:=/,r/m mod 1000/100/
l./read(in,Jmin/,r/Jmin, JQmax, JPRmax,//,l1,i/
read(in,n); for i:=1 step 1 until n do begin
   read(in,J,k,m); if i=1 then Jmin:= J;
   qtn(i):= J shift 12 add (J+k-m) end;
JQmax:= J;  JPRmax:= qtn(n+1):= 0;
vægt(1):= vægt(2):= 1;
/,l./read(in,ki/,d1
l./Limits for F/,l2,d./end else begin/
l./Small/,l1,d,l1,d./Lscale/
l./overskrift/,r/ overskrift;//
l./J:= Jmin/,r/;/; Jmin:= 0;/,l./if sorter/,d2
l./FASE2:/,l9,d4
l./N1-q1-N2+q2/,i/
hb:= false; i:= Jmin;
for i:=i+1 while qtn(i) shift (-12) extract 12 = J do begin
   n:= (qtn(i) extract 12) - J; hb:= hb or
       n = tau1-4*q1 or n = tau2-4*q2 end;
/,l2,r/if/if hb and/
l./if intensitet/,d./end bereg/,d,l6,r/42*3/32*4/
l1,r/f;/W1(q1+84);/,r/A;/W2(q2+84);
   L(p+4):= A;/
l1,r/p+3=126/p+4=128/
l4,d./end udskrift/
l./SLUT:/,l2,d,l./-,Jlige;/,d,i/
i:= Jmin;
for i:=i+1 while qtn(i) shift (-12) extract 12 = J do Jmin:= i;
J1:= qtn(i) shift (-12) extract 12; J:= if J1=0 then J+1 else J1;
J_lige:= J mod 2 = 0;
/,l./end FASE2/,l1,i/
q:= JQmax;
/,l2,r/42*3/32*4/
l7,d./end sorter;/,l1,r/;/; l:= l+4;/
l./stop:/,i/
begin real a; integer J,k1,k2;
open(L,4,<:starkmat:>,0);
nylinie; nylinie; p:= 0; inrec(L,128); a:= L(1);
repeat:
nylinie;
for j:=-24,0 do begin
J:= a shift (j-16) extract 8;
k1:= a shift (j-8) extract 8;
k2:= a shift   j   extract 8;
write(res,<<ddd>,J,k1,k2)  end;
write(res,<<-dd ddd ddd.dd0>,L(p+2),L(p+3),
      <<    d.ddd ddd'-d>,L(p+4));
p:= (p+4) mod 128;  if p=0 then inrec(L,128);
a:= L(p+1);  if a shift (-16) extract 8 <=q then goto repeat;
close(L,true);
end;
/,f
rotcor=set 70
rotcor=algol algcor index.no
edit algcor
l./comment/,l1,p./;/,s,f
clear algcor
▶EOF◀