xref: /illumos-gate/usr/src/lib/libm/amd64/src/remquol.S (revision 55fea89d)
125c28e83SPiotr Jasiukajtis/*
225c28e83SPiotr Jasiukajtis * CDDL HEADER START
325c28e83SPiotr Jasiukajtis *
425c28e83SPiotr Jasiukajtis * The contents of this file are subject to the terms of the
525c28e83SPiotr Jasiukajtis * Common Development and Distribution License (the "License").
625c28e83SPiotr Jasiukajtis * You may not use this file except in compliance with the License.
725c28e83SPiotr Jasiukajtis *
825c28e83SPiotr Jasiukajtis * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
925c28e83SPiotr Jasiukajtis * or http://www.opensolaris.org/os/licensing.
1025c28e83SPiotr Jasiukajtis * See the License for the specific language governing permissions
1125c28e83SPiotr Jasiukajtis * and limitations under the License.
1225c28e83SPiotr Jasiukajtis *
1325c28e83SPiotr Jasiukajtis * When distributing Covered Code, include this CDDL HEADER in each
1425c28e83SPiotr Jasiukajtis * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
1525c28e83SPiotr Jasiukajtis * If applicable, add the following below this CDDL HEADER, with the
1625c28e83SPiotr Jasiukajtis * fields enclosed by brackets "[]" replaced with your own identifying
1725c28e83SPiotr Jasiukajtis * information: Portions Copyright [yyyy] [name of copyright owner]
1825c28e83SPiotr Jasiukajtis *
1925c28e83SPiotr Jasiukajtis * CDDL HEADER END
2025c28e83SPiotr Jasiukajtis */
2125c28e83SPiotr Jasiukajtis/*
2225c28e83SPiotr Jasiukajtis * Copyright 2011 Nexenta Systems, Inc.  All rights reserved.
2325c28e83SPiotr Jasiukajtis */
2425c28e83SPiotr Jasiukajtis/*
2525c28e83SPiotr Jasiukajtis * Copyright 2006 Sun Microsystems, Inc.  All rights reserved.
2625c28e83SPiotr Jasiukajtis * Use is subject to license terms.
2725c28e83SPiotr Jasiukajtis */
2825c28e83SPiotr Jasiukajtis
2925c28e83SPiotr Jasiukajtis	.file	"remquol.s"
3025c28e83SPiotr Jasiukajtis
3125c28e83SPiotr Jasiukajtis#include "libm.h"
3225c28e83SPiotr JasiukajtisLIBM_ANSI_PRAGMA_WEAK(remquol,function)
3325c28e83SPiotr Jasiukajtis#include "libm_protos.h"
3425c28e83SPiotr Jasiukajtis	ENTRY(remquol)
3525c28e83SPiotr Jasiukajtis	fldt	24(%rsp)		/ load arg y
3625c28e83SPiotr Jasiukajtis	fldt	8(%rsp)			/ load arg x
3725c28e83SPiotr Jasiukajtis.Lreml_loop:
3825c28e83SPiotr Jasiukajtis	fprem1				/ partial remainder
3925c28e83SPiotr Jasiukajtis	fstsw	%ax			/ store status word
4025c28e83SPiotr Jasiukajtis	andw	$0x400,%ax		/ check whether reduction complete
4125c28e83SPiotr Jasiukajtis	jne	.Lreml_loop		/ while reduction incomplete, do fprem1
4225c28e83SPiotr Jasiukajtis	fstsw	%ax
4325c28e83SPiotr Jasiukajtis	fwait
4425c28e83SPiotr Jasiukajtis	fstp	%st(1)
4525c28e83SPiotr Jasiukajtis	movw	%ax,%dx
4625c28e83SPiotr Jasiukajtis	andw	$0x4000,%dx		/ get C3
4725c28e83SPiotr Jasiukajtis	sarw	$13,%dx
4825c28e83SPiotr Jasiukajtis	movw	%ax,%cx
4925c28e83SPiotr Jasiukajtis	andw	$0x100,%cx		/ get C0
5025c28e83SPiotr Jasiukajtis	sarw	$6,%cx
5125c28e83SPiotr Jasiukajtis	addw	%cx,%dx
5225c28e83SPiotr Jasiukajtis	andw	$0x200,%ax		/ get C1
5325c28e83SPiotr Jasiukajtis	sarw	$9,%ax
5425c28e83SPiotr Jasiukajtis	addw	%dx,%ax
55*55fea89dSDan Cross	cwtl
5625c28e83SPiotr Jasiukajtis	movl	16(%rsp),%edx		/ sign and bexp of x
5725c28e83SPiotr Jasiukajtis	movl	32(%rsp),%ecx		/ sign and bexp of y
5825c28e83SPiotr Jasiukajtis	andl	$0x8000,%edx		/ edx <- sign(x)
5925c28e83SPiotr Jasiukajtis	andl	$0x8000,%ecx		/ ecx <- sign(y)
6025c28e83SPiotr Jasiukajtis	cmpl	%edx,%ecx
6125c28e83SPiotr Jasiukajtis	je	1f
6225c28e83SPiotr Jasiukajtis	negl	%eax			/ negative n
6325c28e83SPiotr Jasiukajtis1:
6425c28e83SPiotr Jasiukajtis	movl	%eax,(%rdi)		/ last 3 significant bits of quotient
6525c28e83SPiotr Jasiukajtis	ret
6625c28e83SPiotr Jasiukajtis	.align	16
6725c28e83SPiotr Jasiukajtis	SET_SIZE(remquol)
68