DataMuseum.dk

Presents historical artifacts from the history of:

RC4000/8000/9000

This is an automatic "excavation" of a thematic subset of
artifacts from Datamuseum.dk's BitArchive.

See our Wiki for more about RC4000/8000/9000

Excavated with: AutoArchaeologist - Free & Open Source Software.


top - metrics - download

⟦2004c2817⟧ TextFile

    Length: 2304 (0x900)
    Types: TextFile
    Names: »discreteinp«

Derivation

└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
    └─⟦b2ec5d50f⟧ 
        └─⟦this⟧ »discreteinp« 

TextFile


<*generates input to discrete
1982-01-08
Anders Lindgård
*>
if discrete then
begin
  <*generate input for discrete*>
  f:=0;
  readinfp(inp.f,1);
  connectlso;
  write(out,"nl",1,<:mode list.yes
lkj=lookup discrete
if ok.no
i taldiscr
mode list.yes
(
o allist
head cpu
discrete
head cpu
o c
convert allist
mode list.no):>);
  <*card 1*>
  write(out,"nl",1,<: beam foil set :>,inp,"sp",1);
  writeatsym(out,S,atno,Z);
  write(out,"sp",2);
  writestate(out,L,n,l,J);
  write(out,<: MeV :>,<< d.dd>,MeV,<: slit :>,slitw,<: tau :>,tau);
  if false then write(out,"sp",10);
  <*card 2*>
  write(out,"nl",2,
   if continue then <:F:> else <:T:>,
   if -,regint   then <:F:> else <:T:>,
   if true     then <:F:> else <:T:>, <*NOBASEline*>
   if true     then <:F:> else <:T:>, <*NONNEGative exponents*>
   if false    then <:F:> else <:T:>, <*PRint Y values*>
   <:F:>,                             <*PRintREL*>
   if false    then <:F:> else <:T:>, <*PRintFinal results*>
   if false    then <:F:> else <:T:>, <*PLOTReSults*>
   if true     then <:F:> else <:T:>, <*REPEAT summary*>
   "nl",1,<< dd>,maxexp,
   "nl",1,2<*weighted as sqrt(y)*>,
   "nl",1,maxtries);
  if -,regint then
  begin
    <*card 3a*>
    write(out,"nl",1,<< ddd>,no-max+1);
    <*card set 4a*>
    for i:=max step 1 until no do
      write(out,"nl",if (i-max+1) mod 5=1 then 1 else 0,
       << dddd.dd>,datax(i)*tau);
  end 3a and 4a;
  <*card 3b*>
  if regint then write(out,"nl",1,<<ddd>,1, <*one interval*>
      "nl",1,<< dddd.dd>,dt*tau,(no-max+1)*dt*tau,<< dddd>,(no-max+1));
  <*card set 5*>
  for i:=max step 1 until no do
    write(out,"nl",if (i-max+1) mod 5=1 then 1 else 0,
       << dddd.dd>,datay(i));
  <*card set 6*>
  <*only if IWT=4*>
  write(out,"nl",2);
  if fpout then closeout;
end discrete;
▶EOF◀