25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/*
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER START
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis *
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * The contents of this file are subject to the terms of the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Common Development and Distribution License (the "License").
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * You may not use this file except in compliance with the License.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis *
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * or http://www.opensolaris.org/os/licensing.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * See the License for the specific language governing permissions
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * and limitations under the License.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis *
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * When distributing Covered Code, include this CDDL HEADER in each
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * If applicable, add the following below this CDDL HEADER, with the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * fields enclosed by brackets "[]" replaced with your own identifying
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * information: Portions Copyright [yyyy] [name of copyright owner]
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis *
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER END
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/*
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/*
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2005 Sun Microsystems, Inc. All rights reserved.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Use is subject to license terms.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis .file "remainder.s"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#include "libm.h"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr JasiukajtisLIBM_ANSI_PRAGMA_WEAK(remainder,function)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#include "libm_protos.h"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ENTRY(remainder)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis push %rbp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movq %rsp,%rbp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis subq $16,%rsp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movlpd %xmm1,-16(%rbp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movlpd %xmm0,-8(%rbp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ucomisd %xmm0,%xmm1 / if x or y is NaN, use fprem1
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jp 1f
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl -12(%rbp),%eax / eax <-- hi_32(y)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis andl $0x7fffffff,%eax / eax <-- hi_32(|y|)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis orl -16(%rbp),%eax / eax <-- lo_32(y)|hi_32(|y|)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .yzero_or_xinf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl -4(%rbp),%eax / eax <-- hi_32(x)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis andl $0x7fffffff,%eax / eax <-- hi_32(|x|)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0x7ff00000,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jne 1f
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0,-8(%rbp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .yzero_or_xinf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis1:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fldl -16(%rbp) / y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fldl -8(%rbp) / x
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.rem_loop:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fprem1 / partial remainder
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstsw %ax / store status word
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis andw $0x400,%ax / check for incomplete reduction
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jne .rem_loop / while incomplete, do fprem1 again
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstpl -8(%rbp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movsd -8(%rbp),%xmm0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstp %st(0)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis leave
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ret
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.yzero_or_xinf:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_SETUP(1)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl $28,%edi
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl $2,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis call PIC_F(_SVID_libm_err)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_WRAPUP
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis leave
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ret
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis .align 4
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis SET_SIZE(remainder)