/* * Copyright (c) 2002 Apple Computer, Inc. All rights reserved. * * @APPLE_LICENSE_HEADER_START@ * * The contents of this file constitute Original Code as defined in and * are subject to the Apple Public Source License Version 1.1 (the * "License"). You may not use this file except in compliance with the * License. Please obtain a copy of the License at * http://www.apple.com/publicsource and read it before using this file. * * This Original Code and all software distributed under the License are * distributed on an "AS IS" basis, WITHOUT WARRANTY OF ANY KIND, EITHER * EXPRESS OR IMPLIED, AND APPLE HEREBY DISCLAIMS ALL SUCH WARRANTIES, * INCLUDING WITHOUT LIMITATION, ANY WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE OR NON-INFRINGEMENT. Please see the * License for the specific language governing rights and limitations * under the License. * * @APPLE_LICENSE_HEADER_END@ */ /* * Written by J.T. Conklin <jtc@netbsd.org>. * Public domain. * Adapted for exp2 by Stephen C. Peters */ #include <machine/asm.h> #include "abi.h" #warning scp: edge case +-Inf #warning obsolete ENTRY(exp2l) XMM_ONE_ARG_LONG_DOUBLE_PROLOGUE fldt ARG_LONG_DOUBLE_ONE fld %st(0) frndint /* int(x) */ fxch %st(1) fsub %st(1),%st /* fract(x) */ f2xm1 /* 2^(fract(x)) - 1 */ fld1 faddp /* 2^(fract(x)) */ fscale /* 2^x */ fstp %st(1) XMM_LONG_DOUBLE_EPILOGUE ret