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 "log.s"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#include "libm.h"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr JasiukajtisLIBM_ANSI_PRAGMA_WEAK(log,function)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#include "libm_protos.h"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ENTRY(log)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fldln2 / loge(2)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 8(%esp),%eax / eax <-- hi_32(x)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis testl $0x80000000,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jnz .maybe_0_or_less
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis testl $0x7fffffff,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jz .maybe_0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fldl 4(%esp) / arg, loge(2)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fyl2x / loge(2)*log2(arg); ln(arg)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ret
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.maybe_0:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 4(%esp),%ecx / ecx <-- lo_32(x)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0,%ecx
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis je .zero / no branch if x is +denormal
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.neg_nan_reentry:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fldl 4(%esp) / arg, loge(2)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fyl2x / loge(2)*log2(arg); ln(arg)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ret
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.zero_or_less:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis / x =< 0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis testl $0x7fffffff,%eax
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jnz .less_than_0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 4(%esp),%ecx / ecx <-- lo_32(x)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0,%ecx
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jne .less_than_0 / branch if x is -denormal
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.zero:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis / x = +/-0
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl %ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl %esp,%ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_SETUP(1)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl $16
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jmp .merge
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.maybe_0_or_less:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0xfff00000,%eax / -INF below hi_32(x)?
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ja .neg_nan_reentry
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jb .zero_or_less
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl 4(%esp),%ecx / ecx <-- lo_32(x)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis cmpl $0,%ecx / is x NaN or -INF?
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis jne .neg_nan_reentry / branch if x is NaN with signbit = 1
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis / x = -INF
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.less_than_0:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl %ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis movl %esp,%ebp
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis PIC_SETUP(2)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl $17
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis.merge:
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis fstp %st(0) / stack empty
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl 12(%ebp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl 8(%ebp)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis pushl 12(%ebp)
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(log)