Return to atan87.m CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libm / i8087

Annotation of coherent/b/lib/libm/i8087/atan87.m, revision 1.1.1.1

1.1       root        1: //////////
                      2: / libm 8087
                      3: / atan(d)
                      4: //////////
                      5: 
                      6: #include "larges.h"
                      7: #include "ifno8087.h"
                      8: 
                      9:        .globl  atan_
                     10:        .globl  atan
                     11:        .globl  cfcc87
                     12:        .globl  tstcc
                     13: 
                     14: //////////
                     15: / double
                     16: / atan(d)
                     17: / double d;
                     18: /
                     19: / Uses fpatan, which requires 0 < st1 < st.
                     20: //////////
                     21: 
                     22: d      =       RASIZE          / arg offset
                     23: 
                     24: atan_:
                     25:        ifno8087(_atan_)
                     26:        mov     bx, sp
                     27:        fdld    Pss d(bx)       / Load argument d.
                     28: 
                     29: atan:                          / d
                     30:        Gcall   tstcc
                     31:        je      3f              / d = 0, just return it.
                     32:        pushf                   / Save the flags.
                     33:        ja      0f              / d is positive.
                     34:        fchs                    / d is negative, negate it.
                     35: 
                     36: 0:     fld1                    / Load a 1.
                     37:        fcom
                     38:        Gcall   cfcc87          / Compare 1 to d.
                     39:        jae     1f              / 1 >= d, ok.
                     40:        fdivr   st1, st         / and set d = 1/d.
                     41: 
                     42: 1:     fpatan                  / Result = atan(d/1).
                     43:        jae     2f              / 1 >= d, OK as is.
                     44:        fild    cs:minus1       / -1, result
                     45:        fldpi                   / pi, -1, result
                     46:        fscale                  / pi/2, -1, result
                     47:        fstp    st1             / pi/2, result
                     48:        fsubr                   / Result = pi/2 - result.
                     49: 
                     50: 2:     popf                    / Restore the sign flag.
                     51:        jae     3f              / Positive, done.
                     52:        fchs                    / Result = -result.
                     53: 
                     54: 3:     Gret
                     55: 
                     56: / Data.
                     57: minus1:        .word   -1
                     58: 
                     59: / end of atan87.m

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