Return to exp.s CVS log

Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.lib / libm

Annotation of cci/usr/src/usr.lib/libm/exp.s, revision 1.1

1.1     ! root        1: /*     @(#)exp.s       4.1/4.2         10/31/84  CCI-CPG       */
        !             2: 
        !             3: 
        !             4: /*
        !             5:        Double-precision exponential function recoded with
        !             6:        new algorithm polynomial at CCI-CPG.
        !             7: 
        !             8:        The coefficients are #1067 from Hart and etc. (18.08D)
        !             9: */
        !            10:        .data
        !            11:        .comm   _errno,4
        !            12:        .align  2
        !            13: _p0:    .long  0x45bd3d04, 0x7f734dbd # .double 1513.9067990543389
        !            14:        .align  2
        !            15: _p1:   .long   0x42a19dd4, 0x989f60db # .double 20.202065651286927
        !            16:        .align  2
        !            17: _p2:   .long   0x3dbd2e42, 0xab70bdab # .double 0.023093347753750234
        !            18:        .align  2
        !            19: _q0:   .long   0x468881b1, 0x7c3a652b # .double 4368.2116627275585
        !            20:        .align  2
        !            21: _q1:   .long   0x44692f28, 0x7ae89543 # .double 233.18421142748162
        !            22:        .align  2
        !            23: _q2:   .long   0x40800000, 0x00000000 # .double 1
        !            24:        .align  2
        !            25: _log2e:        .long   0x40b8aa3b, 0x295c17f1 # .double 1.4426950408889634
        !            26:        .align  2
        !            27: _sqrt2:        .long   0x40b504f3, 0x33f9de65 # .double 1.4142135623730951
        !            28:        .align  2
        !            29: _maxf: .long   0x471c4000, 0x00000000 # .double 10000
        !            30:        .text
        !            31: LL0:   .align  1
        !            32:        .globl  _exp
        !            33:        .data
        !            34:        .align  2
        !            35: L56:   .long   0x40800000, 0x00000000 # .double 1
        !            36:        .text
        !            37:        .data
        !            38:        .align  2
        !            39: L59:   .long   0x7ffffffe, 0xfffffffe # .double 1.7014117331926443e+38
        !            40:        .text
        !            41:        .data
        !            42:        .align  2
        !            43: L60:   .long   0x40000000, 0x00000000 # .double 0.5
        !            44:        .text
        !            45:        .set    L51,0xc
        !            46:        .data
        !            47:        .text
        !            48: _exp:  .word   L51
        !            49:        subl3   $88,fp,sp
        !            50:        tstl    4(fp)
        !            51:        jneq    L55
        !            52:        movl    L56+4,r1
        !            53:        movl    L56,r0
        !            54:        ret
        !            55: L55:   lnd     _maxf
        !            56:        cmpd    4(fp)
        !            57:        jleq    L57
        !            58:        clrl    r1
        !            59:        clrl    r0
        !            60:        ret
        !            61: L57:   cmpd2   4(fp),_maxf
        !            62:        jleq    L58
        !            63:        movl    $34,_errno
        !            64:        movl    L59+4,r1
        !            65:        movl    L59,r0
        !            66:        ret
        !            67: L58:   ldd     4(fp)
        !            68:        muld    _log2e
        !            69:        std     4(fp)
        !            70:        pushl   8(fp)
        !            71:        pushl   4(fp)
        !            72:        callf   $12,_floor
        !            73:        ldd     r0
        !            74:        cvdl    -88(fp)
        !            75:        cvld    -88(fp)
        !            76:        std     r0
        !            77:        ldd     4(fp)
        !            78:        subd    r0
        !            79:        subd    L60
        !            80:        std     -60(fp)
        !            81:        muld    -60(fp)
        !            82:        std     -84(fp)
        !            83:        ldd     _p2
        !            84:        muld    -84(fp)
        !            85:        addd    _p1
        !            86:        muld    -84(fp)
        !            87:        addd    _p0
        !            88:        muld    -60(fp)
        !            89:        std     -68(fp)
        !            90:        ldd     _q2
        !            91:        muld    -84(fp)
        !            92:        addd    _q1
        !            93:        muld    -84(fp)
        !            94:        addd    _q0
        !            95:        std     -76(fp)
        !            96:        pushl   -88(fp)
        !            97:        addd    -68(fp)
        !            98:        muld    _sqrt2
        !            99:        std     r0
        !           100:        ldd     -76(fp)
        !           101:        subd    -68(fp)
        !           102:        std     r2
        !           103:        ldd     r0
        !           104:        divd    r2
        !           105:        pushd
        !           106:        callf   $16,_ldexp
        !           107:        ret
        !           108: 

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