Return to ldexp.s CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libc / crt / i386

Annotation of coherent/b/lib/libc/crt/i386/ldexp.s, revision 1.1

1.1     ! root        1: //////////
        !             2: / libc/crt/i386/ldexp.s
        !             3: / i386 C runtime library.
        !             4: / IEEE software floating point support.
        !             5: / i8086 IEEE software floating point library.
        !             6: / ldexp()
        !             7: //////////
        !             8: 
        !             9: //////////
        !            10: / double
        !            11: / ldexp(dbl, exp) double dbl; int exp;
        !            12: /
        !            13: / Returns dbl * 2^^exp.
        !            14: //////////
        !            15: 
        !            16:        .globl  ldexp
        !            17: 
        !            18: dbl    .equ    4
        !            19: exp    .equ    dbl+8
        !            20: EXPMASK        .equ    0x7FF00000
        !            21: SGNMASK        .equ    0x80000000
        !            22: MAXEXP .equ    2047
        !            23: 
        !            24: ldexp:
        !            25:        movl    %ecx, dbl+4(%esp)               / high word of dbl to ECX
        !            26:        movl    %edx, $EXPMASK
        !            27:        andl    %edx, %ecx                      / shifted exponent to EDX
        !            28:        xorl    %ecx, %edx                      / mask off exponent in ECX
        !            29:        shrl    %edx, $20                       / exponent in EDX
        !            30:        addl    %edx, exp(%esp)                 / bump by exp
        !            31:        jle     ?underflow
        !            32:        cmpl    %edx, $MAXEXP
        !            33:        jge     ?overflow
        !            34:        shll    %edx, $20
        !            35:        orl     %edx, %ecx                      / hi mantissa bits to EDX
        !            36:        movl    %eax, dbl(%esp)                 / lo manitssa dword to EAX
        !            37:        ret
        !            38: 
        !            39: ?underflow:
        !            40:        andl    %ecx, $SGNMASK                  / preserve sign
        !            41:        jmp     ?done
        !            42: 
        !            43: ?overflow:
        !            44:        andl    %ecx, $SGNMASK                  / preserve sign
        !            45:        orl     %ecx, $EXPMASK                  / max exponent
        !            46: ?done:
        !            47:        movl    %edx, %ecx                      / sign/exp to EDX
        !            48:        subl    %eax, %eax                      / 0 lo mantissa to EAX
        !            49:        ret
        !            50: 
        !            51: / end of libc/crt/i386/ldexp.s

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