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

   1 /
   2 /   C E X P
   3 /   - - - -
   4 /
   5 /  COMPLEX EXPONENT ROUTINE
   6 /Z=X+IY
   7 /
   8 / VERSION 5A 4-25-77 MH
   9 /
  10 /
  11 /CEXP(Z)=EXP(X)*(COS(Y)+I*SIN(Y))
  12 /ENTER+EXIT IN COMPLEX
  13 /EXTERNAL EXP,SIN,COS
  14 /
  15         SECT    CEXP
  16         JA      #CEXP
  17         DPCHK
  18         TEXT    +CEXP  +
  19 CEXPXR, SETX    XR
  20         SETB    BP
  21         JA      .+3
  22 BP,     F 0.0
  23 XR,     F 0.0
  24 ARG,    F 0.0
  25         F 0.0
  26         ORG     10*3+BP
  27         FNOP
  28         JA      CEXPXR
  29         0
  30 RT,     JA      .
  31         BASE    0
  32 #CEXP,  STARTD
  33         FLDA    10*3
  34         FSTA    RT
  35         FLDA    0
  36         SETB    BP
  37         SETX    XR
  38         BASE    BP
  39         LDX     1,1
  40         FSTA    BP
  41         FLDA%   BP,1
  42         FSTA    BP
  43         STARTE
  44         FLDA%   BP
  45         FSTA    ARG
  46         STARTF
  47         EXTERN  EXP
  48         JSR     EXP             /EXP(X)
  49         JA      CEX1
  50         JA      ARG
  51 CEX1,   FSTA    ETEMP
  52         EXTERN  COS
  53         JSR     COS             /COS(Y)
  54         JA      CEX2
  55         JA      ARG+3
  56 CEX2,   FSTA    ARG
  57         EXTERN  SIN
  58         JSR     SIN             /SIN(Y)
  59         JA      CEX3
  60         JA      ARG+3
  61 CEX3,   FSTA    ARG+3
  62         FLDA    ETEMP
  63         FMULM   ARG
  64         FMULM   ARG+3
  65         STARTE
  66         FLDA    ARG
  67         FSTA    #CAC
  68         JA      RT
  69         EXTERN  #CAC
  70 ETEMP,  F 0.0
  71 \f