diff options
author | Igor Pashev <pashev.igor@gmail.com> | 2012-09-11 19:12:10 +0400 |
---|---|---|
committer | Igor Pashev <pashev.igor@gmail.com> | 2012-09-11 19:12:10 +0400 |
commit | 19700b860d9ec70d01e885d92c3d4f62fd052873 (patch) | |
tree | b6b265cd81fb97b4675f4a285f397479609fdf18 /usr/src/libm/src/i386/amd64/remainder.s | |
download | libm-19700b860d9ec70d01e885d92c3d4f62fd052873.tar.gz |
Imported Upstream version 20060131HEADupstream/20060131upstreammaster
Diffstat (limited to 'usr/src/libm/src/i386/amd64/remainder.s')
-rw-r--r-- | usr/src/libm/src/i386/amd64/remainder.s | 79 |
1 files changed, 79 insertions, 0 deletions
diff --git a/usr/src/libm/src/i386/amd64/remainder.s b/usr/src/libm/src/i386/amd64/remainder.s new file mode 100644 index 0000000..984b7ec --- /dev/null +++ b/usr/src/libm/src/i386/amd64/remainder.s @@ -0,0 +1,79 @@ +/* + * CDDL HEADER START + * + * The contents of this file are subject to the terms of the + * Common Development and Distribution License (the "License"). + * You may not use this file except in compliance with the License. + * + * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE + * or http://www.opensolaris.org/os/licensing. + * See the License for the specific language governing permissions + * and limitations under the License. + * + * When distributing Covered Code, include this CDDL HEADER in each + * file and include the License file at usr/src/OPENSOLARIS.LICENSE. + * If applicable, add the following below this CDDL HEADER, with the + * fields enclosed by brackets "[]" replaced with your own identifying + * information: Portions Copyright [yyyy] [name of copyright owner] + * + * CDDL HEADER END + */ +/* + * Copyright 2005 Sun Microsystems, Inc. All rights reserved. + * Use is subject to license terms. + */ + + .ident "@(#)remainder.s 1.2 06/01/23 SMI" + + .file "remainder.s" + +#include "libm.h" +LIBM_ANSI_PRAGMA_WEAK(remainder,function) +#include "libm_synonyms.h" +#include "libm_protos.h" + + ENTRY(remainder) + push %rbp + movq %rsp,%rbp + subq $16,%rsp + movlpd %xmm1,-16(%rbp) + movlpd %xmm0,-8(%rbp) + + ucomisd %xmm0,%xmm1 / if x or y is NaN, use fprem1 + jp 1f + + movl -12(%rbp),%eax / eax <-- hi_32(y) + andl $0x7fffffff,%eax / eax <-- hi_32(|y|) + orl -16(%rbp),%eax / eax <-- lo_32(y)|hi_32(|y|) + je .yzero_or_xinf + + movl -4(%rbp),%eax / eax <-- hi_32(x) + andl $0x7fffffff,%eax / eax <-- hi_32(|x|) + cmpl $0x7ff00000,%eax + jne 1f + cmpl $0,-8(%rbp) + je .yzero_or_xinf +1: + fldl -16(%rbp) / y + fldl -8(%rbp) / x +.rem_loop: + fprem1 / partial remainder + fstsw %ax / store status word + andw $0x400,%ax / check for incomplete reduction + jne .rem_loop / while incomplete, do fprem1 again + fstpl -8(%rbp) + movsd -8(%rbp),%xmm0 + fstp %st(0) + leave + ret + +.yzero_or_xinf: + PIC_SETUP(1) + movl $28,%edi + movl $2,%eax + call PIC_F(_SVID_libm_err) + PIC_WRAPUP + leave + ret + .align 4 + SET_SIZE(remainder) |