home *** CD-ROM | disk | FTP | other *** search
/ Simtel MSDOS 1992 June / SIMTEL_0692.cdr / msdos / ddjmag / ddj8601.arc / BARR.JAN next >
Encoding:
Text File  |  1986-01-31  |  1.8 KB  |  63 lines
  1. cseg    segment    para public 'CODE'
  2.  
  3.     org    100h
  4.  
  5.     assume    cs:cseg,ds:cseg,es:cseg,ss:cseg
  6.  
  7. main    proc    near
  8.     nop
  9.     call sqrt 
  10.     jmp    main
  11. main    endp
  12.  
  13.  
  14. ; Michael Barr's 32-bit Square Root Routine
  15. ; 
  16. ; Call with CX:BX = argument
  17. ; Returns   BX    = result
  18. ;
  19.  
  20. sqrt      proc      near           ;CX:BX = argument
  21.           push      ax             ;save other registers
  22.           push      dx
  23.           push      di
  24.           jcxz      sqrt3
  25.           mov       dx,cx          ;prepare first try
  26.           mov       di,-1
  27. sqrt1:    shl       dx,1           ;estimating size of arg.
  28.           jc        sqrt2          ;to guess initial try
  29.           shl       dx,1
  30.           jc        sqrt2
  31.           shr       di,1
  32.           jmp       sqrt1
  33. sqrt2:    mov       dx,cx          ;restore argument
  34.           mov       ax,bx
  35.           cmp       dx,di          ;prevent overflow
  36.           jae       sqrt4
  37.           div       di
  38.           cmp       ax,di          ;comp quotient and divisor
  39.           jae       sqrt4
  40.           add       di,ax          ;average them
  41.           rcr       di,1
  42.           jmp       sqrt2          ;and do it again
  43. sqrt3:    mov       dx,bx          ;prepare first try
  44.           mov       di,0ffh
  45.           or        bx,bx          ;lower half zero?
  46.           jnz       sqrt1          ;no, jump
  47.           mov       di,bx
  48. sqrt4:    mov       bx,di          ;return result in BX
  49.           pop       di
  50.           pop       dx
  51.           pop       ax
  52.           ret
  53. sqrt      endp
  54.  
  55.  
  56. cseg    ends
  57.  
  58.     end    main
  59. 
  60. cseg    ends
  61.  
  62.     end    main
  63.