Return to exp.s CVS log

Up to [CSRG BSD Unix] / 43BSD / usr.lib / liboldnm

Annotation of 43BSD/usr.lib/liboldnm/exp.s, revision 1.1

1.1     ! root        1: #
        !             2: # Copyright (c) 1980 Regents of the University of California.
        !             3: # All rights reserved.  The Berkeley software License Agreement
        !             4: # specifies the terms and conditions for redistribution.
        !             5: #
        !             6: #      @(#)exp.s       5.1 (Berkeley) 5/8/85
        !             7: # 
        !             8: #
        !             9: # double exp(arg)
        !            10: # double arg; 0<= arg< 88
        !            11: # double exp10(arg)
        !            12: # method: range reduction to [-0.5,0.5], Hart&Chenet EXPB 1067
        !            13: # J F Jarvis, August 5, 1978
        !            14: .set   ERANGE,34
        !            15: .globl _exp
        !            16: .globl _exp10
        !            17: .globl _errno
        !            18: .text
        !            19: .align 1
        !            20: _exp10:
        !            21:        .word   0x07c0
        !            22:        bispsw  $0xe0
        !            23:        muld3   4(ap),loge10,r0
        !            24:        jbr     argtst
        !            25: .align 1
        !            26: _exp:
        !            27:        .word   0x07c0
        !            28:        bispsw  $0xe0
        !            29:        movd    4(ap),r0
        !            30: argtst:
        !            31:        jnequ   smlarg
        !            32:        movd    $0d1.0e+0,r0
        !            33:        ret
        !            34: smlarg:
        !            35:        cmpd    r0,minarg
        !            36:        jgeq    lrgarg
        !            37:        mnegd   huge,r0
        !            38:        ret
        !            39: lrgarg:
        !            40:        cmpd    r0,maxarg
        !            41:        jleq    range
        !            42:        movl    $ERANGE,_errno
        !            43:        movd    huge,r0
        !            44:        ret
        !            45: range:
        !            46:        emodd   log2e,log2ex,r0,r10,r6  # r10=int(arg), r6=frac(arg)
        !            47:        tstd    r0
        !            48:        jgtr    l1
        !            49:        addd2   $0d0.5e+0,r6
        !            50:        subw2   $1,r10
        !            51:        jbr     l2
        !            52: l1:    subd2   $0d0.5e+0,r6
        !            53: l2:                    # Hart&Cheney EXPB 1067
        !            54:        muld3   r6,r6,r8        # range [-.5,.5] D=18.1
        !            55:        polyd   r8,$2,pcoef
        !            56:        muld2   r0,r6
        !            57:        polyd   r8,$2,qcoef
        !            58:        subd3   r6,r0,r2
        !            59:        addd2   r6,r0
        !            60:        divd2   r2,r0
        !            61:        muld2   sqrt2,r0
        !            62:        extzv $7,$8,r0,r2
        !            63:        addl2   r2,r10
        !            64:        insv    r10,$7,$8,r0    # load correct exponent
        !            65:        ret
        !            66: .data
        !            67: .align 2
        !            68: pcoef:
        !            69:        .double 0d0.23093347753750233624e-1
        !            70:        .double 0d0.20202065651286927227886e+2
        !            71:        .double 0d0.15139067990543389159e+4
        !            72: qcoef:
        !            73:        .double 0d0.1e+1
        !            74:        .double 0d0.23318421142748162379e+3
        !            75:        .double 0d0.43682116627275584985e+4
        !            76: minarg:
        !            77:        .double 0d-88.028e+0
        !            78: maxarg:
        !            79:        .double 0d88.028e+0
        !            80: huge:  .double 0d1.7e+38
        !            81: loge10:
        !            82:        .double 0d2.30258509299404568401e+0
        !            83: sqrt2:
        !            84:        .double 0d1.41421356237309504880e+0
        !            85: log2e:
        !            86:        .double 0d1.44269504088896340735e+0
        !            87: log2ex:
        !            88:        .byte   0xbb

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