Return to fnorm.s CVS log

Up to [CSRG BSD Unix] / 43BSDTahoe / lib / libc / tahoe / fpe

Annotation of 43BSDTahoe/lib/libc/tahoe/fpe/fnorm.s, revision 1.1.1.1

1.1       root        1: #ifdef LIBC_SCCS
                      2:        .asciz  "@(#)fnorm.s    1.1 (Berkeley/CCI) 7/2/86"
                      3: #endif LIBC_SCCS
                      4: 
                      5: #include <tahoemath/fp.h>
                      6: #include "DEFS.h"
                      7: 
                      8: ENTRY(fnorm, R2|R3|R4|R5|R6)
                      9:        movl    r0,r4           # copy to temporary.
                     10:        jneq    inr0
                     11:        movl    r1,r5
                     12:        clrl    r3              # r3 - pos of m.s.b
                     13: inr1:  ffs     r5,r6
                     14:        incl    r6
                     15:        addl2   r6,r3
                     16:        shrl    r6,r5,r5
                     17:        jneq    inr1
                     18:        cmpl    $0,r3
                     19:        jeql    retzero
                     20:        jmp     cmpshift
                     21: inr0:  movl    $32,r3
                     22: inr00: ffs     r4,r6
                     23:        incl    r6
                     24:        addl2   r6,r3
                     25:        shrl    r6,r4,r4
                     26:        jneq    inr00
                     27: 
                     28: cmpshift:
                     29:                                # compute the shift (r4).
                     30:        subl3   r3,$HID_R0R1,r4
                     31:        jlss    shiftr          # if less then zero we shift right.
                     32:        shlq    r4,r0,r0        # else we shift left.
                     33:        subl2   r4,r2           # uodate exponent.
                     34:        jleq    underflow       # if less then 0 (biased) it is underflow.
                     35:        jmp     combine         # go to combine exponent and fraction.
                     36: shiftr:
                     37:        mnegl   r4,r4
                     38:        shrq    r4,r0,r0        # shift right.
                     39:        addl2   r4,r2           # update exponent
                     40:        cmpl    r2,$256
                     41:        jgeq    overflow        # check for overflow.
                     42: combine:
                     43:        andl2   $CLEARHID,r0    # clear the hidden bit.
                     44:        shal    $EXPSHIFT,r2,r2 # shift the exponent to its proper place.
                     45:        orl2    r2,r0
                     46:        ret
                     47: 
                     48: underflow:
                     49:        callf   $4,fpunder
                     50:        ret
                     51: 
                     52: overflow:
                     53:        callf   $4,fpover
                     54:        ret
                     55: retzero:
                     56:        clrl    r0
                     57:        clrl    r1
                     58:        ret

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.