| [d9be488] | 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 |
|
|---|
| [321a8e6] | 29 | #include <abi/asmtool.h>
|
|---|
| [d9be488] | 30 | #include <libarch/x87.h>
|
|---|
| 31 |
|
|---|
| 32 | .text
|
|---|
| 33 |
|
|---|
| [321a8e6] | 34 | FUNCTION_BEGIN(sin_f64)
|
|---|
| [d9be488] | 35 | # compute sine (no stack frame)
|
|---|
| 36 |
|
|---|
| 37 | fldl 4(%esp)
|
|---|
| 38 | fsin
|
|---|
| 39 |
|
|---|
| 40 | # detect if source operand is out of range
|
|---|
| 41 |
|
|---|
| 42 | fnstsw %ax
|
|---|
| 43 | andw $X87_STATUS_WORD_C2_MASK, %ax
|
|---|
| 44 | jnz fix_range
|
|---|
| 45 |
|
|---|
| 46 | ret
|
|---|
| 47 |
|
|---|
| 48 | # argument reduction
|
|---|
| 49 |
|
|---|
| 50 | fix_range:
|
|---|
| 51 | fldpi
|
|---|
| 52 | fadd %st(0)
|
|---|
| 53 | fxch %st(1)
|
|---|
| 54 |
|
|---|
| 55 | reduce:
|
|---|
| 56 | fprem1
|
|---|
| 57 | fnstsw %ax
|
|---|
| 58 | andw $X87_STATUS_WORD_C2_MASK, %ax
|
|---|
| 59 | jnz reduce
|
|---|
| 60 |
|
|---|
| [8f5e80be] | 61 | fstp %st(1)
|
|---|
| [d9be488] | 62 | fsin
|
|---|
| 63 |
|
|---|
| 64 | ret
|
|---|
| [321a8e6] | 65 | FUNCTION_END(sin_f64)
|
|---|
