|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 768 (0x300)
Types: TextFile
Names: »extsyminver«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt.
└─⟦0364f57e3⟧
└─⟦58ca399f1⟧ »extbib«
└─⟦this⟧
syminverse=algol list.yes index.no
external
procedure syminverse(A,L,n);
value n; integer n; array A,L;
begin
integer i,j,k,m;
real h,H;
array f(1:n);
for i:=1 step 1 until n do
for j:=i step 1 until n do L(i,j):= A(i,j);
comment A is inverted (by the method INVRS from CACM 66-P1-0);
m:= n-1;
for k:=0 step 1 until m do begin
H:= 1/L(1,1);
for i:=2 step 1 until n do f(i-1):= L(1,i);
for i:=1 step 1 until m do begin
L(i,n):= h:= -f(i)*H;
for j:=i step 1 until m do L(i,j):= L(i+1,j+1) + f(j)*h
end i;
L(n,n):= -H
end k;
for i:=1 step 1 until n do
for j:=i step 1 until n do L(i,j):= L(j,i):= -L(i,j);
end;
end
▶EOF◀