Return to atan.s CVS log

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

Annotation of 43BSD/usr.lib/liboldnm/atan.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: #      @(#)atan.s      5.1 (Berkeley) 5/8/85
        !             7: # 
        !             8: #
        !             9: # double atan(arg1); -pi/2 < atan < pi/2
        !            10: #double atan2(arg1,arg2); -pi< atan2<pi
        !            11: #double arg1,arg2
        !            12: #method: range reduction to [sqrt(2)-1,sqrt(2)+1]
        !            13: # followed by Hart&Cheney ARCTN 5076 D=17.55
        !            14: # J. F. Jarvis August 8, 1978
        !            15: .globl _atan
        !            16: .globl _atan2
        !            17: .text
        !            18: .align 1
        !            19: _atan:
        !            20:        .word   0x03c0
        !            21:        bispsw  $0xe0
        !            22:        movd    4(ap),r0
        !            23:        jgtr    a1
        !            24:        mnegd   r0,r0   # atan(arg1), arg1<0
        !            25:        bsbb    satan
        !            26:        mnegd   r0,r0
        !            27:        ret
        !            28: a1:    bsbb    satan   # atan(arg1), arg1>=0
        !            29:        ret
        !            30: #
        !            31: .align 1
        !            32: _atan2:
        !            33:        .word   0x03c0
        !            34:        bispsw  $0xe0
        !            35:        movd    4(ap),r0        # atan(arg1/arg2)
        !            36:        movd    12(ap),r2
        !            37:        addd3   r0,r2,r4
        !            38:        cmpd    r0,r4
        !            39:        jneq    b1
        !            40:        tstd    r0
        !            41:        jgeq    b2
        !            42:        mnegd   pio2,r0
        !            43:        ret
        !            44: b2:    movd    pio2,r0
        !            45:        ret
        !            46: #
        !            47: b1:    tstd    r2
        !            48:        jgeq    b3
        !            49:        divd2   r2,r0
        !            50:        jleq    b4
        !            51:        bsbb    satan   # arg1<0, arg2<0
        !            52:        subd2   pi,r0
        !            53:        ret
        !            54: b4:    mnegd   r0,r0   # arg1>0, arg2<0
        !            55:        bsbb    satan
        !            56:        subd3   r0,pi,r0
        !            57:        ret
        !            58: #
        !            59: b3:    divd2   r2,r0
        !            60:        jleq    b5
        !            61:        bsbb    satan   # arg1>0, arg2>0
        !            62:        ret
        !            63: b5:    mnegd   r0,r0   # arg1<0, arg2>0
        !            64:        bsbb    satan
        !            65:        mnegd   r0,r0
        !            66:        ret
        !            67: #
        !            68: .globl satan
        !            69: satan: # range reduction on positive arg(r0)
        !            70:        cmpd    r0,sq2m1
        !            71:        jgeq    c1
        !            72:        bsbb    xatan
        !            73:        rsb
        !            74: c1:    cmpd    r0,sq2p1
        !            75:        jleq    c2
        !            76:        divd3   r0,$0d1.0e+0,r0
        !            77:        bsbb    xatan
        !            78:        subd3   r0,pio2,r0
        !            79:        rsb
        !            80: c2:    addd3   $0d1.0e+0,r0,r2
        !            81:        subd2   $0d1.0e+0,r0
        !            82:        divd2   r2,r0
        !            83:        bsbb    xatan
        !            84:        addd2   pio4,r0
        !            85:        rsb
        !            86: #
        !            87: xatan: # compute arctan(r0) for:sqrt(2)-1<r0<sqrt(2)+1
        !            88: # Hart&Cheney ARCTN 5076 is evaluated
        !            89:        movd    r0,r8
        !            90:        muld3   r0,r0,r6
        !            91:        polyd   r6,$4,pcoef
        !            92:        muld2   r0,r8
        !            93:        polyd   r6,$4,qcoef
        !            94:        divd3   r0,r8,r0
        !            95:        rsb
        !            96: .data
        !            97: .align 2
        !            98: pcoef:
        !            99:        .double 0d0.1589740288482307048e+0
        !           100:        .double 0d0.66605790170092626575e+1
        !           101:        .double 0d0.40969264832102256374e+2
        !           102:        .double 0d0.77477687719204208616e+2
        !           103:        .double 0d0.44541340059290680319e+2
        !           104: qcoef:
        !           105:        .double 0d1.0e+0
        !           106:        .double 0d0.15503977551421987525e+2
        !           107:        .double 0d0.62835930511032376833e+2
        !           108:        .double 0d0.92324801072300974840e+2
        !           109:        .double 0d0.44541340059290680444e+2
        !           110: pio4: .double  0d0.78539816339744830961e+0
        !           111: pio2: .double  0d1.57079632679489661923e+0
        !           112: sq2p1: .double 0d2.41421356237309504880e+0
        !           113: sq2m1: .double 0d0.41421356237309504880e+0
        !           114: pi: .double 0d3.14159265358979323846e+0

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