1// Copyright (C) 1988-1994 Sun Microsystems, Inc. 2550 Garcia Avenue 2// Mountain View, California 94043 All rights reserved. 3// 4// Any person is hereby authorized to download, copy, use, create bug fixes, 5// and distribute, subject to the following conditions: 6// 7// 1. the software may not be redistributed for a fee except as 8// reasonable to cover media costs; 9// 2. any copy of the software must include this notice, as well as 10// any other embedded copyright notices; and 11// 3. any distribution of this software or derivative works thereof 12// must comply with all applicable U.S. export control laws. 13// 14// THE SOFTWARE IS MADE AVAILABLE "AS IS" AND WITHOUT EXPRESS OR IMPLIED 15// WARRANTY OF ANY KIND, INCLUDING BUT NOT LIMITED TO THE IMPLIED 16// WARRANTIES OF DESIGN, MERCHANTIBILITY, FITNESS FOR A PARTICULAR 17// PURPOSE, NON-INFRINGEMENT, PERFORMANCE OR CONFORMANCE TO 18// SPECIFICATIONS. 19// 20// BY DOWNLOADING AND/OR USING THIS SOFTWARE, THE USER WAIVES ALL CLAIMS 21// AGAINST SUN MICROSYSTEMS, INC. AND ITS AFFILIATED COMPANIES IN ANY 22// JURISDICTION, INCLUDING BUT NOT LIMITED TO CLAIMS FOR DAMAGES OR 23// EQUITABLE RELIEF BASED ON LOSS OF DATA, AND SPECIFICALLY WAIVES EVEN 24// UNKNOWN OR UNANTICIPATED CLAIMS OR LOSSES, PRESENT AND FUTURE. 25// 26// IN NO EVENT WILL SUN MICROSYSTEMS, INC. OR ANY OF ITS AFFILIATED 27// COMPANIES BE LIABLE FOR ANY LOST REVENUE OR PROFITS OR OTHER SPECIAL, 28// INDIRECT AND CONSEQUENTIAL DAMAGES, EVEN IF IT HAS BEEN ADVISED OF THE 29// POSSIBILITY OF SUCH DAMAGES. 30// 31// This file is provided with no support and without any obligation on the 32// part of Sun Microsystems, Inc. ("Sun") or any of its affiliated 33// companies to assist in its use, correction, modification or 34// enhancement. Nevertheless, and without creating any obligation on its 35// part, Sun welcomes your comments concerning the software and requests 36// that they be sent to fdlibm-comments@sunpro.sun.com. 37// atan2s(+-0,+anything but nan) is +-0 38T(RN, 0x0p+0, 0x0p+0, 0x0p+0, 0x0p+0, 0) 39T(RN, 0x0p+0, 0x1p-149, 0x0p+0, 0x0p+0, 0) 40T(RN, 0x0p+0, 0x1p-129, 0x0p+0, 0x0p+0, 0) 41T(RN, 0x0p+0, 0x1.fffffep+127, 0x0p+0, 0x0p+0, 0) 42T(RN, 0x0p+0, inf, 0x0p+0, 0x0p+0, 0) 43T(RN, -0x0p+0, 0x0p+0, -0x0p+0, 0x0p+0, 0) 44T(RN, -0x0p+0, 0x1p-149, -0x0p+0, 0x0p+0, 0) 45T(RN, -0x0p+0, 0x1p-129, -0x0p+0, 0x0p+0, 0) 46T(RN, -0x0p+0, 0x1.fffffep+127, -0x0p+0, 0x0p+0, 0) 47T(RN, -0x0p+0, inf, -0x0p+0, 0x0p+0, 0) 48// atan2s(+-0,-anything but nan) is +-pi 49T(RN, 0x0p+0, -inf, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 50T(RN, 0x0p+0,-0x1.fffffep+127, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 51T(RN, 0x0p+0, -0x1p-129, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 52T(RN, 0x0p+0, -0x1p-149, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 53T(RN, 0x0p+0, -0x0p+0, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 54T(RN, -0x0p+0, -inf, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 55T(RN, -0x0p+0,-0x1.fffffep+127, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 56T(RN, -0x0p+0, -0x1p-129, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 57T(RN, -0x0p+0, -0x1p-149, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 58T(RN, -0x0p+0, -0x0p+0, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 59// atan2s(+-anything but 0 and nan, 0) is +- pi/2 60T(RN, inf, 0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 61T(RN, inf, -0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 62T(RN, 0x1.fffffep+127, 0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 63T(RN, 0x1.fffffep+127, -0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 64T(RN, 0x1p-126, 0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 65T(RN, 0x1p-126, -0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 66T(RN, 0x1p-149, 0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 67T(RN, 0x1p-149, -0x0p+0, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 68T(RN, -inf, 0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 69T(RN, -inf, -0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 70T(RN,-0x1.fffffep+127, 0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 71T(RN,-0x1.fffffep+127, -0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 72T(RN, -0x1p-126, 0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 73T(RN, -0x1p-126, -0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 74T(RN, -0x1p-149, 0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 75T(RN, -0x1p-149, -0x0p+0, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 76// atan2s(big,small) :=: +-pi/2 77T(RN, 0x1.fffffep+127, 0x1p-126, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 78T(RN, 0x1.fffffep+127, -0x1p-126, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 79T(RN,-0x1.fffffep+127, 0x1p-126, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 80T(RN,-0x1.fffffep+127, -0x1p-126, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 81T(RZ, 0x1.fffffep+127, 0x1p-126, 0x1.921fb4p+0, -0x1.4442d2p-1, INEXACT) 82T(RZ, 0x1.fffffep+127, -0x1p-126, 0x1.921fb4p+0, -0x1.4442d2p-1, INEXACT) 83T(RZ,-0x1.fffffep+127, 0x1p-126, -0x1.921fb4p+0, 0x1.4442d2p-1, INEXACT) 84T(RZ,-0x1.fffffep+127, -0x1p-126, -0x1.921fb4p+0, 0x1.4442d2p-1, INEXACT) 85T(RU, 0x1.fffffep+127, 0x1p-126, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 86T(RU, 0x1.fffffep+127, -0x1p-126, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 87T(RU,-0x1.fffffep+127, 0x1p-126, -0x1.921fb4p+0, 0x1.4442d2p-1, INEXACT) 88T(RU,-0x1.fffffep+127, -0x1p-126, -0x1.921fb4p+0, 0x1.4442d2p-1, INEXACT) 89T(RD, 0x1.fffffep+127, 0x1p-126, 0x1.921fb4p+0, -0x1.4442d2p-1, INEXACT) 90T(RD, 0x1.fffffep+127, -0x1p-126, 0x1.921fb4p+0, -0x1.4442d2p-1, INEXACT) 91T(RD,-0x1.fffffep+127, 0x1p-126, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 92T(RD,-0x1.fffffep+127, -0x1p-126, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 93// atan2s(small,big) = small/big (big>0) 94T(RN, 0x1p-126, 0x1.fffffep+127, 0x0p+0,-0x1.000002p-105, INEXACT|UNDERFLOW) 95T(RN, -0x1p-126, 0x1.fffffep+127, -0x0p+0, 0x1.000002p-105, INEXACT|UNDERFLOW) 96T(RN, 0x1p-126,-0x1.fffffep+127, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 97T(RN, -0x1p-126,-0x1.fffffep+127, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 98T(RZ, 0x1p-126, 0x1.fffffep+127, 0x0p+0,-0x1.000002p-105, INEXACT|UNDERFLOW) 99T(RZ, -0x1p-126, 0x1.fffffep+127, -0x0p+0, 0x1.000002p-105, INEXACT|UNDERFLOW) 100T(RZ, 0x1p-126,-0x1.fffffep+127, 0x1.921fb4p+1, -0x1.4442d2p-1, INEXACT) 101T(RZ, -0x1p-126,-0x1.fffffep+127, -0x1.921fb4p+1, 0x1.4442d2p-1, INEXACT) 102T(RU, 0x1p-126, 0x1.fffffep+127, 0x1p-149, 0x1p+0, INEXACT|UNDERFLOW) 103T(RU, -0x1p-126, 0x1.fffffep+127, -0x0p+0, 0x1.000002p-105, INEXACT|UNDERFLOW) 104T(RU, 0x1p-126,-0x1.fffffep+127, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 105T(RU, -0x1p-126,-0x1.fffffep+127, -0x1.921fb4p+1, 0x1.4442d2p-1, INEXACT) 106T(RD, 0x1p-126, 0x1.fffffep+127, 0x0p+0,-0x1.000002p-105, INEXACT|UNDERFLOW) 107T(RD, -0x1p-126, 0x1.fffffep+127, -0x1p-149, -0x1p+0, INEXACT|UNDERFLOW) 108T(RD, 0x1p-126,-0x1.fffffep+127, 0x1.921fb4p+1, -0x1.4442d2p-1, INEXACT) 109T(RD, -0x1p-126,-0x1.fffffep+127, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 110// atan2s(+-x,+x) = +-pi/4 for normal x 111T(RN, 0x1p-126, 0x1p-126, 0x1.921fb6p-1, 0x1.777a5cp-2, INEXACT) 112T(RN, -0x1p-126, 0x1p-126, -0x1.921fb6p-1, -0x1.777a5cp-2, INEXACT) 113T(RN, 0x1.fffffep+127, 0x1.fffffep+127, 0x1.921fb6p-1, 0x1.777a5cp-2, INEXACT) 114T(RN,-0x1.fffffep+127, 0x1.fffffep+127, -0x1.921fb6p-1, -0x1.777a5cp-2, INEXACT) 115// atan2s(+-x,-x) = +-3pi/4 for normal x 116T(RN, 0x1p-126, -0x1p-126, 0x1.2d97c8p+1, 0x1.99bc5cp-6, INEXACT) 117T(RN, -0x1p-126, -0x1p-126, -0x1.2d97c8p+1, -0x1.99bc5cp-6, INEXACT) 118T(RN, 0x1p+0, -0x1p+0, 0x1.2d97c8p+1, 0x1.99bc5cp-6, INEXACT) 119T(RN, -0x1p+0, -0x1p+0, -0x1.2d97c8p+1, -0x1.99bc5cp-6, INEXACT) 120T(RN, 0x1.fffffep+127,-0x1.fffffep+127, 0x1.2d97c8p+1, 0x1.99bc5cp-6, INEXACT) 121T(RN,-0x1.fffffep+127,-0x1.fffffep+127, -0x1.2d97c8p+1, -0x1.99bc5cp-6, INEXACT) 122// random arguments between -2.0 and 2.0 123T(RN, -0x1.13284cp-1, 0x1.6ca8ep+0, -0x1.716d2p-2, 0x1.90111p-3, INEXACT) 124T(RN, 0x1.c2ca6p+0, -0x1.55f12p+0, 0x1.1c206ep+1, 0x1.1ac042p-2, INEXACT) 125T(RN, -0x1.15679ep-1, -0x1.41e132p-4, -0x1.b6ff44p+0, -0x1.2ecaf4p-4, INEXACT) 126T(RN, 0x1.281b0ep+0, 0x1.b5ce34p+0, 0x1.30789cp-1, 0x1.17e61ap-3, INEXACT) 127T(RN, -0x1.583482p-2, -0x1.ea8224p+0, -0x1.7be508p+1, 0x1.e4c608p-3, INEXACT) 128T(RN, -0x1.aae18p-1, 0x1.7a9da2p-4, -0x1.75db7cp+0, 0x1.483b7cp-6, INEXACT) 129T(RN, -0x1.845p+0, 0x1.6ca322p+0, -0x1.a237b2p-1, -0x1.e3262p-3, INEXACT) 130T(RN, -0x1.fe09bep+0, -0x1.ff6c7ep+0, -0x1.2dc43cp+1, 0x1.348c12p-4, INEXACT) 131T(RN, -0x1.d24c82p-1, 0x1.29b682p+0, -0x1.5428a2p-1, 0x1.e59fe4p-3, INEXACT) 132T(RN, 0x1.25ea7ep-1, 0x1.f99598p-1, 0x1.0d9c96p-1, -0x1.3ef8ecp-7, INEXACT) 133// atan2s involve nan 134T(RN, 0x0p+0, nan, nan, 0x0p+0, 0) 135T(RN, 0x1p+0, nan, nan, 0x0p+0, 0) 136T(RN, nan, 0x1p-149, nan, 0x0p+0, 0) 137T(RN, nan,-0x1.fffffep+127, nan, 0x0p+0, 0) 138T(RN, nan, nan, nan, 0x0p+0, 0) 139T(RN, nan, nan, nan, 0x0p+0, 0) 140T(RZ, 0x0p+0, nan, nan, 0x0p+0, 0) 141T(RZ, 0x1p+0, nan, nan, 0x0p+0, 0) 142T(RZ, nan, 0x1p-149, nan, 0x0p+0, 0) 143T(RZ, nan,-0x1.fffffep+127, nan, 0x0p+0, 0) 144T(RZ, nan, nan, nan, 0x0p+0, 0) 145T(RZ, nan, nan, nan, 0x0p+0, 0) 146T(RU, 0x0p+0, nan, nan, 0x0p+0, 0) 147T(RU, 0x1p+0, nan, nan, 0x0p+0, 0) 148T(RU, nan, 0x1p-149, nan, 0x0p+0, 0) 149T(RU, nan,-0x1.fffffep+127, nan, 0x0p+0, 0) 150T(RU, nan, nan, nan, 0x0p+0, 0) 151T(RU, nan, nan, nan, 0x0p+0, 0) 152T(RD, 0x0p+0, nan, nan, 0x0p+0, 0) 153T(RD, 0x1p+0, nan, nan, 0x0p+0, 0) 154T(RD, nan, 0x1p-149, nan, 0x0p+0, 0) 155T(RD, nan,-0x1.fffffep+127, nan, 0x0p+0, 0) 156T(RD, nan, nan, nan, 0x0p+0, 0) 157T(RD, nan, nan, nan, 0x0p+0, 0) 158// atan2s(+-(anything but inf and nan), +inf) is +-0 159T(RN, 0x1p-149, inf, 0x0p+0, 0x0p+0, 0) 160T(RN, 0x1.fffffep+127, inf, 0x0p+0, 0x0p+0, 0) 161T(RN, -0x1p-149, inf, -0x0p+0, 0x0p+0, 0) 162T(RN,-0x1.fffffep+127, inf, -0x0p+0, 0x0p+0, 0) 163// atan2s(+-(anything but inf and nan), -inf) is +-pi 164T(RN, 0x1p-149, -inf, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 165T(RN, 0x1.fffffep+127, -inf, 0x1.921fb6p+1, 0x1.777a5cp-2, INEXACT) 166T(RN, -0x1p-149, -inf, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 167T(RN,-0x1.fffffep+127, -inf, -0x1.921fb6p+1, -0x1.777a5cp-2, INEXACT) 168// atan2s(+-inf,+inf ) is +-pi/4 169T(RN, inf, inf, 0x1.921fb6p-1, 0x1.777a5cp-2, INEXACT) 170T(RN, -inf, inf, -0x1.921fb6p-1, -0x1.777a5cp-2, INEXACT) 171// atan2s(+-inf,-inf ) is +-3pi/4 172T(RN, inf, -inf, 0x1.2d97c8p+1, 0x1.99bc5cp-6, INEXACT) 173T(RN, -inf, -inf, -0x1.2d97c8p+1, -0x1.99bc5cp-6, INEXACT) 174// atan2s(+-inf, (anything but,0,nan, and inf)) is +-pi/2 175T(RN, inf, 0x1p-149, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 176T(RN, inf, -0x1p-149, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 177T(RN, inf, 0x1.fffffep+127, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 178T(RN, inf,-0x1.fffffep+127, 0x1.921fb6p+0, 0x1.777a5cp-2, INEXACT) 179T(RN, -inf, 0x1p-149, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 180T(RN, -inf, -0x1p-149, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 181T(RN, -inf, 0x1.fffffep+127, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 182T(RN, -inf,-0x1.fffffep+127, -0x1.921fb6p+0, -0x1.777a5cp-2, INEXACT) 183