home *** CD-ROM | disk | FTP | other *** search
/ minnie.tuhs.org / unixen.tar / unixen / PDP-11 / Trees / V6 / usr / source / fort / rt1 / cabs.s < prev    next >
Encoding:
Text File  |  1975-07-17  |  545 b   |  58 lines
  1. /
  2. /
  3.  
  4. / cabs fortran function
  5.  
  6. .globl    cabs.
  7. .globl    dcabs.
  8.  
  9. .globl    sqrt
  10. .globl    retrn
  11. .globl    rerr
  12. .globl    temp
  13.  
  14. dcabs.:    temp
  15.     .+2
  16.     setd
  17.     br    1f
  18.  
  19. cabs.:    temp
  20.     .+2
  21.     setf
  22. 1:
  23.     mov    2(r3),r1
  24.     movf    (r1)+,r0
  25.     movf    (r1)+,r1
  26.     jsr    pc,hypot
  27.     bes    1f
  28.     movf    r0,temp
  29.     jmp    retrn
  30. 1:
  31.     jsr    r5,rerr; 4
  32. hypot:
  33.     movf    fr1,-(sp)
  34.     absf    fr0
  35.     absf    fr1
  36.     cmpf    fr0,fr1
  37.     cfcc
  38.     bne    1f
  39.     movf    $one,fr0
  40.     br    3f
  41. 1:
  42.     blt    2f
  43.     movf    fr0,-(sp)
  44.     movf    fr1,fr0
  45.     movf    (sp)+,fr1
  46. 2:
  47.     divf    fr1,fr0
  48.     mulf    fr0,fr0
  49. 3:
  50.     addf    $one,fr0
  51.     jsr    pc,sqrt
  52.     mulf    fr1,fr0
  53.     cfcc
  54.     movf    (sp)+,fr1
  55.     rts    pc
  56.  
  57. one = 40200
  58.