Return to acos87.m CVS log

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

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

1.1       root        1: //////////
                      2: / libm 8087
                      3: / acos(d)
                      4: //////////
                      5: 
                      6: #include "larges.h"
                      7: #include "ifno8087.h"
                      8: 
                      9:        .globl  acos_
                     10:        .globl  asin
                     11:        .globl  cfcc87
                     12:        .globl  edomain
                     13: 
                     14: //////////
                     15: / double
                     16: / acos(d)
                     17: / double d;
                     18: //////////
                     19: 
                     20: d      =       RASIZE          / arg offset
                     21: 
                     22: acos_:
                     23:        ifno8087(_acos_)
                     24:        mov     bx, sp
                     25:        fdld    Pss d(bx)       / Load argument d.
                     26: 
                     27:        fld     st              / d, d
                     28:        fabs                    / abs(d), d
                     29:        fld1                    / 1, abs(d), d
                     30:        fcompp                  / Compare 1 to abs(d) and pop twice.
                     31:        Gcall   cfcc87
                     32:        jae     0f
                     33:        Gjmp    edomain         / d < -1 or 1 < d, issue EDOM error.
                     34: 0:
                     35:        Gcall   asin            / asin(d)
                     36:        fild    cs:minus1       / -1, asin(d)
                     37:        fldpi                   / pi, -1, asin(d)
                     38:        fscale                  / pi/2, -1, asin(d)
                     39:        fstp    st1             / pi/2, asin(d)
                     40:        fsubr                   / acos(d) = pi/2 - asin(d)
                     41:        Gret
                     42: 
                     43: / Data.
                     44: minus1:        .word   -1
                     45: 
                     46: / end of acos87.m

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