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 2006 Sun Microsystems, Inc. All rights reserved.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Use is subject to license terms.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis .file "nextafterf.s"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#include "libm.h"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr JasiukajtisLIBM_ANSI_PRAGMA_WEAK(nextafterf,function)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis .data
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis .align 4
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr JasiukajtisFmaxf: .long 0x7f7fffff
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr JasiukajtisFminf: .long 0x1
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisftmpf: .long 0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ENTRY(nextafterf)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl %ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl %esp,%ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl $0,%eax /// upper half of %eax must be initialized
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis flds 12(%ebp) / y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis subl $4,%esp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis flds 8(%ebp) / x, y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fucom / x : y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstsw %ax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis sahf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jp .NaN
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .equal
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstp %st(1) / x
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ja .bigger
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis / x < y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ftst / x : 0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl $0x1,-4(%ebp) / -4(%ebp) contains Fminf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fnstsw %ax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis sahf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .final
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ja .addulp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jb .subulp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.bigger:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis / x > y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ftst / x : 0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl $0x80000001,-4(%ebp) / -4(%ebp) contains -Fminf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fnstsw %ax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis sahf
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .final
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jb .addulp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.subulp:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 8(%ebp),%eax / x
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis subl $1,%eax / x - ulp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl %eax,-4(%ebp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jmp .final
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.addulp:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 8(%ebp),%eax / x
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis addl $1,%eax / x + ulp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl %eax,-4(%ebp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.final:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstp %st(0) / empty
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis flds -4(%ebp) / z
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis andl $0x7f800000,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jz .underflow
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0x7f800000,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .overflow
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jmp .return
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.overflow:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_SETUP(1)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis flds PIC_L(Fmaxf) / Fmaxf, z
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fmul %st(0),%st / overflow-to-Inf, z
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstps PIC_L(ftmpf) / z & create overflow signal
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_WRAPUP
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jmp .return
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.underflow:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_SETUP(2)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis flds PIC_L(Fminf) / Fminf, z
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fmul %st(0),%st / underflow-to-0, z
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstps PIC_L(ftmpf) / z & create underflow signal
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_WRAPUP
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jmp .return
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.equal:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstp %st(0) / C99 says to return y when x == y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jmp .return
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.NaN:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis faddp %st,%st(1) / x+y
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.return:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fwait
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis leave
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ret
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis .align 4
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis SET_SIZE(nextafterf)