| 1 | #
|
|---|
| 2 | # Copyright (c) 2014 Martin Decky
|
|---|
| 3 | # All rights reserved.
|
|---|
| 4 | #
|
|---|
| 5 | # Redistribution and use in source and binary forms, with or without
|
|---|
| 6 | # modification, are permitted provided that the following conditions
|
|---|
| 7 | # are met:
|
|---|
| 8 | #
|
|---|
| 9 | # - Redistributions of source code must retain the above copyright
|
|---|
| 10 | # notice, this list of conditions and the following disclaimer.
|
|---|
| 11 | # - Redistributions in binary form must reproduce the above copyright
|
|---|
| 12 | # notice, this list of conditions and the following disclaimer in the
|
|---|
| 13 | # documentation and/or other materials provided with the distribution.
|
|---|
| 14 | # - The name of the author may not be used to endorse or promote products
|
|---|
| 15 | # derived from this software without specific prior written permission.
|
|---|
| 16 | #
|
|---|
| 17 | # THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
|
|---|
| 18 | # IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
|
|---|
| 19 | # OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
|
|---|
| 20 | # IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
|
|---|
| 21 | # INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
|
|---|
| 22 | # NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
|
|---|
| 23 | # DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
|
|---|
| 24 | # THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
|
|---|
| 25 | # (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
|
|---|
| 26 | # THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|---|
| 27 | #
|
|---|
| 28 |
|
|---|
| 29 | #include <abi/asmtool.h>
|
|---|
| 30 | #include <libarch/x87.h>
|
|---|
| 31 |
|
|---|
| 32 | .text
|
|---|
| 33 |
|
|---|
| 34 | FUNCTION_BEGIN(cos_f64)
|
|---|
| 35 | pushq %rbp
|
|---|
| 36 | movq %rsp, %rbp
|
|---|
| 37 |
|
|---|
| 38 | # compute cosine (using red zone)
|
|---|
| 39 |
|
|---|
| 40 | movsd %xmm0, -8(%rbp)
|
|---|
| 41 | fldl -8(%rbp)
|
|---|
| 42 |
|
|---|
| 43 | fcos
|
|---|
| 44 |
|
|---|
| 45 | # detect if source operand is out of range
|
|---|
| 46 |
|
|---|
| 47 | fnstsw %ax
|
|---|
| 48 | andw $X87_STATUS_WORD_C2_MASK, %ax
|
|---|
| 49 | jnz fix_range
|
|---|
| 50 |
|
|---|
| 51 | fstpl -8(%rbp)
|
|---|
| 52 | movsd -8(%rbp), %xmm0
|
|---|
| 53 |
|
|---|
| 54 | leave
|
|---|
| 55 | retq
|
|---|
| 56 |
|
|---|
| 57 | # argument reduction
|
|---|
| 58 |
|
|---|
| 59 | fix_range:
|
|---|
| 60 | fldpi
|
|---|
| 61 | fadd %st(0)
|
|---|
| 62 | fxch %st(1)
|
|---|
| 63 |
|
|---|
| 64 | reduce:
|
|---|
| 65 | fprem1
|
|---|
| 66 | fnstsw %ax
|
|---|
| 67 | andw $X87_STATUS_WORD_C2_MASK, %ax
|
|---|
| 68 | jnz reduce
|
|---|
| 69 |
|
|---|
| 70 | fstp %st(1)
|
|---|
| 71 | fcos
|
|---|
| 72 |
|
|---|
| 73 | fstpl -8(%rbp)
|
|---|
| 74 | movsd -8(%rbp), %xmm0
|
|---|
| 75 |
|
|---|
| 76 | leave
|
|---|
| 77 | retq
|
|---|
| 78 | FUNCTION_END(cos_f64)
|
|---|
