sw/os8/v3d/sources/fortran/dectapes/dectape2/exp3.ra

   1 /
   2 /
   3 /
   4 /       E  X  P  3
   5 /       -  -  -  -
   6 /
   7 /SUBROUTINE     EXP3(B,E)
   8 /
   9 / VERSION 5A 4-26-77 MH
  10 /
  11         SECT    EXP3
  12         JA      #EXP3
  13         EXTERN  #ARGER
  14 EXP3ER, TRAP4   #ARGER
  15         TEXT    +EXP3  +
  16 EXP3XR, SETX    XREXP3
  17         SETB    BPEXP3
  18 BPEXP3, FNOP
  19         0
  20         0
  21 XREXP3, F 0.0
  22 EXP31,  F 0.0
  23 EXP32,  F 0.0
  24 FP1XP3, F 1.
  25         ORG     10*3+BPEXP3
  26         FNOP
  27         JA      EXP3XR
  28         0
  29 XP3RTN, JA      .
  30         BASE            0
  31 #EXP3,  STARTD
  32         FLDA    10*3
  33         FSTA    XP3RTN
  34         FLDA    0
  35         SETX    XREXP3
  36         SETB    BPEXP3
  37         BASE    BPEXP3
  38         LDX     1,1
  39         FSTA    BPEXP3
  40         FLDA%   BPEXP3,1        /ADDR OF B
  41         FSTA    EXP31
  42         FLDA%   BPEXP3,1+       /ADDR OF E
  43         FSTA    EXP32
  44         STARTF
  45         FLDA%   EXP31   /GET B
  46         JEQ     XP3RTN  /0 ^ X = 0
  47         FSTA    EXP31   /SAVE BASE
  48         FLDA%   EXP32   /GET E
  49         JEQ     EXP3ON  /X ^ 0 = 1
  50         FSTA    EXP32   /SAVE EXPONENT
  51         FLDA    EXP31
  52         JLT     EXP3ER  /ALL IS NOT WELL
  53         EXTERN  ALOG
  54         JSR     ALOG    /CALL LOG
  55         JA      .+4     /TAKE LOG (B)
  56         JA      EXP31
  57         FMUL    EXP32   /MULT BY EXPONENT-E
  58         FSTA    EXP31
  59         EXTERN  EXP
  60         JSR     EXP     /CALL EXP.
  61         JA      XP3RTN
  62         JA      EXP31
  63 EXP3ON, FLDA    FP1XP3
  64         JA      XP3RTN
  65 \f