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 pushl %ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl %esp,%ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fldl 16(%esp) / load arg y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fldl 8(%esp) / load arg x
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fucom
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fnstsw %ax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis sahf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jp .rem_loop / if x or y is NaN, use fprem1
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 20(%esp),%eax / eax <-- hi_32(y)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis andl $0x7fffffff,%eax / eax <-- hi_32(|y|)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis orl 16(%esp),%eax / eax <-- lo_32(y)|hi_32(|y|)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .yzero_or_xinf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 12(%esp),%eax / eax <-- hi_32(x)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis andl $0x7fffffff,%eax / eax <-- hi_32(|x|)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0x7ff00000,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jne .rem_loop
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0,8(%esp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .yzero_or_xinf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
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 fstp %st(1)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis leave
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ret
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.yzero_or_xinf:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_SETUP(1)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstp %st(0) / x
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstp %st(0) / empty NPX stack
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl $28 / case 28 in _SVID_libm_err
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl 20(%ebp) / pass y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl 16(%ebp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl 12(%ebp) / pass x
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl 8(%ebp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis call PIC_F(_SVID_libm_err)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis addl $20,%esp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_WRAPUP
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis leave
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ret
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis .align 4
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis SET_SIZE(remainder)