5261 libm should stop using synonyms.h 5298 fabs is 0-sized, confuses dis(1) and others Reviewed by: Josef 'Jeff' Sipek <jeffpc@josefsipek.net> Approved by: Gordon Ross <gwr@nexenta.com>
1 /* 2 * CDDL HEADER START 3 * 4 * The contents of this file are subject to the terms of the 5 * Common Development and Distribution License (the "License"). 6 * You may not use this file except in compliance with the License. 7 * 8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE 9 * or http://www.opensolaris.org/os/licensing. 10 * See the License for the specific language governing permissions 11 * and limitations under the License. 12 * 13 * When distributing Covered Code, include this CDDL HEADER in each 14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE. 15 * If applicable, add the following below this CDDL HEADER, with the 16 * fields enclosed by brackets "[]" replaced with your own identifying 17 * information: Portions Copyright [yyyy] [name of copyright owner] 18 * 19 * CDDL HEADER END 20 */ 21 /* 22 * Copyright 2011 Nexenta Systems, Inc. All rights reserved. 23 */ 24 /* 25 * Copyright 2005 Sun Microsystems, Inc. All rights reserved. 26 * Use is subject to license terms. 27 */ 28 29 .file "remainder.s" 30 31 #include "libm.h" 32 LIBM_ANSI_PRAGMA_WEAK(remainder,function) 33 #include "libm_protos.h" 34 35 ENTRY(remainder) 36 push %rbp 37 movq %rsp,%rbp 38 subq $16,%rsp 39 movlpd %xmm1,-16(%rbp) 40 movlpd %xmm0,-8(%rbp) 41 42 ucomisd %xmm0,%xmm1 / if x or y is NaN, use fprem1 43 jp 1f 44 45 movl -12(%rbp),%eax / eax <-- hi_32(y) 46 andl $0x7fffffff,%eax / eax <-- hi_32(|y|) 47 orl -16(%rbp),%eax / eax <-- lo_32(y)|hi_32(|y|) 48 je .yzero_or_xinf 49 50 movl -4(%rbp),%eax / eax <-- hi_32(x) 51 andl $0x7fffffff,%eax / eax <-- hi_32(|x|) 52 cmpl $0x7ff00000,%eax 53 jne 1f 54 cmpl $0,-8(%rbp) 55 je .yzero_or_xinf 56 1: 57 fldl -16(%rbp) / y 58 fldl -8(%rbp) / x 59 .rem_loop: 60 fprem1 / partial remainder 61 fstsw %ax / store status word 62 andw $0x400,%ax / check for incomplete reduction 63 jne .rem_loop / while incomplete, do fprem1 again 64 fstpl -8(%rbp) 65 movsd -8(%rbp),%xmm0 66 fstp %st(0) 67 leave 68 ret 69 70 .yzero_or_xinf: 71 PIC_SETUP(1) 72 movl $28,%edi 73 movl $2,%eax 74 call PIC_F(_SVID_libm_err) 75 PIC_WRAPUP 76 leave 77 ret 78 .align 4 79 SET_SIZE(remainder) --- EOF ---