mirror of
https://sourceware.org/git/glibc.git
synced 2024-11-22 21:10:07 +00:00
548 lines
15 KiB
ArmAsm
548 lines
15 KiB
ArmAsm
.file "log2f.s"
|
|
|
|
|
|
// Copyright (c) 2000 - 2003, Intel Corporation
|
|
// All rights reserved.
|
|
//
|
|
// Contributed 2000 by the Intel Numerics Group, Intel Corporation
|
|
//
|
|
// Redistribution and use in source and binary forms, with or without
|
|
// modification, are permitted provided that the following conditions are
|
|
// met:
|
|
//
|
|
// * Redistributions of source code must retain the above copyright
|
|
// notice, this list of conditions and the following disclaimer.
|
|
//
|
|
// * Redistributions in binary form must reproduce the above copyright
|
|
// notice, this list of conditions and the following disclaimer in the
|
|
// documentation and/or other materials provided with the distribution.
|
|
//
|
|
// * The name of Intel Corporation may not be used to endorse or promote
|
|
// products derived from this software without specific prior written
|
|
// permission.
|
|
|
|
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
|
|
// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
|
|
// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
|
|
// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL INTEL OR ITS
|
|
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
|
|
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
|
|
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
|
|
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
|
|
// OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY OR TORT (INCLUDING
|
|
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
|
|
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
//
|
|
// Intel Corporation is the author of this code, and requests that all
|
|
// problem reports or change requests be submitted to it directly at
|
|
// http://www.intel.com/software/products/opensource/libraries/num.htm.
|
|
//
|
|
// History
|
|
//==============================================================
|
|
// 09/11/00 Initial version
|
|
// 05/20/02 Cleaned up namespace and sf0 syntax
|
|
// 02/10/03 Reordered header: .section, .global, .proc, .align
|
|
//
|
|
// API
|
|
//==============================================================
|
|
// float log2f(float)
|
|
//
|
|
// Overview of operation
|
|
//==============================================================
|
|
// Background
|
|
//
|
|
// Implementation
|
|
//
|
|
// Let x = 2^l * m, where m=1.b1 b2 ... b8 b9 ... b52
|
|
// y=frcpa(m), r=m*y-1, f=b1 b2 .. b8 (table index)
|
|
// j=0 if f<128; j=1 if f>=128
|
|
// T is a table that stores log2(1/y) (in entries 1..255) rounded to
|
|
// double extended precision; f is used as an index; T[255]=0
|
|
//
|
|
// If f=0 and b9=0, r is set to 2^{-8}* 0.b9 b10 ... b52 = m-1 (fractional part of m),
|
|
// and 0 is used instead of T[0]
|
|
// (polynomial evaluation only, for m=1+r, 0<=r<2^{-9})
|
|
// If f=255, r is set to (m-2)/2 (T[255]=0, and only polynomial evaluation is used
|
|
// for m=2(1-r'), 0<=r'<2^{-9})
|
|
//
|
|
// log2f(x) is approximated as
|
|
// (l-j) + T[f] + (c1*r+c2*r^2+...+c6*r^6), if f>0
|
|
//
|
|
|
|
|
|
// Special values
|
|
//==============================================================
|
|
// log2f(0)=-inf, raises Divide by Zero
|
|
// log2f(+inf)=inf
|
|
// log2f(x)=NaN, raises Invalid if x<0
|
|
//
|
|
|
|
|
|
// Registers used
|
|
//==============================================================
|
|
// f6-f14
|
|
// r2-r3, r23-r30
|
|
// p6,p7,p8,p12
|
|
//
|
|
|
|
|
|
GR_SAVE_B0 = r33
|
|
GR_SAVE_PFS = r34
|
|
GR_SAVE_GP = r35 // This reg. can safely be used
|
|
GR_SAVE_SP = r36
|
|
|
|
GR_Parameter_X = r37
|
|
GR_Parameter_Y = r38
|
|
GR_Parameter_RESULT = r39
|
|
GR_Parameter_TAG = r40
|
|
|
|
FR_X = f10
|
|
FR_Y = f1
|
|
FR_RESULT = f8
|
|
|
|
|
|
|
|
|
|
// Data tables
|
|
//==============================================================
|
|
|
|
RODATA
|
|
|
|
.align 16
|
|
|
|
LOCAL_OBJECT_START(poly_coeffs)
|
|
|
|
data8 0x3fdec709dc3a03fd, 0xbfd71547652b82fe //C_3 and C_4
|
|
data8 0xb8aa3b295c17f0bc, 0x00003fff // C_1
|
|
data8 0xb8aa3b295c17f0bc, 0x0000bffe // C_2
|
|
LOCAL_OBJECT_END(poly_coeffs)
|
|
|
|
|
|
LOCAL_OBJECT_START(T_table)
|
|
|
|
data8 0x3f671b0ea42e5fda, 0x3f815cfe8eaec830
|
|
data8 0x3f8cfee70c5ce5dc, 0x3f94564a62192834
|
|
data8 0x3f997723ace35766, 0x3f9f5923c69b54a1
|
|
data8 0x3fa2a094a085d693, 0x3fa538941776b01e
|
|
data8 0x3fa8324c9b914bc7, 0x3faacf54ce07d7e9
|
|
data8 0x3fadced958dadc12, 0x3fb0387efbca869e
|
|
data8 0x3fb18ac6067479c0, 0x3fb30edd3e13530d
|
|
data8 0x3fb463c15936464e, 0x3fb5b9e13c3fa21d
|
|
data8 0x3fb7113f3259e07a, 0x3fb869dd8d1b2035
|
|
data8 0x3fb9c3bea49d3214, 0x3fbb1ee4d7961701
|
|
data8 0x3fbc7b528b70f1c5, 0x3fbdd90a2c676ed4
|
|
data8 0x3fbf05d4976c2028, 0x3fc032fbbaee6d65
|
|
data8 0x3fc0e3b5a9f3284a, 0x3fc195195c7d125b
|
|
data8 0x3fc22dadc2ab3497, 0x3fc2e050231df57d
|
|
data8 0x3fc379f79c2b255b, 0x3fc42ddd2ba1b4a9
|
|
data8 0x3fc4c89b9e6807f5, 0x3fc563dc29ffacb2
|
|
data8 0x3fc619a25f5d798d, 0x3fc6b5ffbf367644
|
|
data8 0x3fc752e1f660f8d6, 0x3fc7f049e753e7cf
|
|
data8 0x3fc8a8980abfbd32, 0x3fc94724cca657be
|
|
data8 0x3fc9e63a24971f46, 0x3fca85d8feb202f7
|
|
data8 0x3fcb2602497d5346, 0x3fcbc6b6f5ee1c9b
|
|
data8 0x3fcc67f7f770a67e, 0x3fcceec4b2234fba
|
|
data8 0x3fcd91097ad13982, 0x3fce33dd57f3d335
|
|
data8 0x3fced74146bc7b10, 0x3fcf7b3646fef683
|
|
data8 0x3fd00223a943dc19, 0x3fd054a474bf0eb7
|
|
data8 0x3fd0999d9b9259a1, 0x3fd0eca66d3b2581
|
|
data8 0x3fd13ffa2e85b475, 0x3fd185a444fa0a7b
|
|
data8 0x3fd1cb8312f27eff, 0x3fd21fa1441ce5e8
|
|
data8 0x3fd265f526e603cb, 0x3fd2baa0c34be1ec
|
|
data8 0x3fd3016b45de21ce, 0x3fd3486c38aa29a8
|
|
data8 0x3fd38fa3efaa8262, 0x3fd3e562c0816a02
|
|
data8 0x3fd42d141f53b646, 0x3fd474fd543f222c
|
|
data8 0x3fd4bd1eb680e548, 0x3fd505789e234bd1
|
|
data8 0x3fd54e0b64003b70, 0x3fd596d761c3c1f0
|
|
data8 0x3fd5dfdcf1eeae0e, 0x3fd6291c6fd9329c
|
|
data8 0x3fd6729637b59418, 0x3fd6bc4aa692e0fd
|
|
data8 0x3fd7063a1a5fb4f2, 0x3fd75064f1ed0715
|
|
data8 0x3fd79acb8cf10390, 0x3fd7d67c1e43ae5c
|
|
data8 0x3fd8214f4068afa7, 0x3fd86c5f36dea3dc
|
|
data8 0x3fd8b7ac64dd7f9d, 0x3fd8f4167a0c6f92
|
|
data8 0x3fd93fd2d5e1bf1d, 0x3fd98bcd84296946
|
|
data8 0x3fd9c8c333e6e9a5, 0x3fda152f142981b4
|
|
data8 0x3fda527fd95fd8ff, 0x3fda9f5e3edeb9e6
|
|
data8 0x3fdadd0b2b5755a7, 0x3fdb2a5d6f51ff83
|
|
data8 0x3fdb686799b00be3, 0x3fdbb62f1b887cd8
|
|
data8 0x3fdbf4979f666668, 0x3fdc332a6e8399d4
|
|
data8 0x3fdc819dc2d45fe4, 0x3fdcc0908e19b7bd
|
|
data8 0x3fdcffae611ad12b, 0x3fdd3ef776d43ff4
|
|
data8 0x3fdd8e5002710128, 0x3fddcdfb486cb9a1
|
|
data8 0x3fde0dd294245fe4, 0x3fde4dd622a28840
|
|
data8 0x3fde8e06317114f0, 0x3fdece62fe9a9915
|
|
data8 0x3fdf1f164a15389a, 0x3fdf5fd8a9063e35
|
|
data8 0x3fdfa0c8937e7d5d, 0x3fdfe1e649bb6335
|
|
data8 0x3fe011990641535a, 0x3fe032560e91e59e
|
|
data8 0x3fe0532a5ebcd44a, 0x3fe0741617f5fc28
|
|
data8 0x3fe08cd653f38839, 0x3fe0adeb55c1103b
|
|
data8 0x3fe0cf181d5d1dd0, 0x3fe0f05ccd0aced7
|
|
data8 0x3fe111b9875788ab, 0x3fe1332e6f1bcf73
|
|
data8 0x3fe154bba77c2088, 0x3fe16df59bfa06c1
|
|
data8 0x3fe18fadb6e2d3c2, 0x3fe1b17e849adc26
|
|
data8 0x3fe1caeb6a0de814, 0x3fe1ece7c830eec9
|
|
data8 0x3fe20efd3dae01df, 0x3fe2289de375d901
|
|
data8 0x3fe24adf9b6a6fe0, 0x3fe26d3ad1aebcfc
|
|
data8 0x3fe287100c2771f4, 0x3fe2a9983b3c1b28
|
|
data8 0xbfda78e146f7bef4, 0xbfda33760a7f6051
|
|
data8 0xbfd9ff43476fb5f7, 0xbfd9b97c3c4eec8f
|
|
data8 0xbfd98504431717fc, 0xbfd93ee07535f967
|
|
data8 0xbfd90a228d5712b2, 0xbfd8c3a104cb24f5
|
|
data8 0xbfd88e9c72e0b226, 0xbfd847bc33d8618e
|
|
data8 0xbfd812703988bb69, 0xbfd7dd0569c04bff
|
|
data8 0xbfd7959c202292f1, 0xbfd75fe8d2c5d48f
|
|
data8 0xbfd72a1637cbc183, 0xbfd6e221cd9d0cde
|
|
data8 0xbfd6ac059985503b, 0xbfd675c99ce81f92
|
|
data8 0xbfd63f6db2590482, 0xbfd5f6c138136489
|
|
data8 0xbfd5c01a39fbd688, 0xbfd58952cf519193
|
|
data8 0xbfd5526ad18493ce, 0xbfd51b6219bfe6ea
|
|
data8 0xbfd4d1cdf8b4846f, 0xbfd49a784bcd1b8b
|
|
data8 0xbfd4630161832547, 0xbfd42b6911cf5465
|
|
data8 0xbfd3f3af3461e1c4, 0xbfd3bbd3a0a1dcfb
|
|
data8 0xbfd383d62dac7ae7, 0xbfd34bb6b2546218
|
|
data8 0xbfd313750520f520, 0xbfd2db10fc4d9aaf
|
|
data8 0xbfd2a28a6dc90387, 0xbfd269e12f346e2c
|
|
data8 0xbfd2311515e2e855, 0xbfd1f825f6d88e13
|
|
data8 0xbfd1bf13a6c9c69f, 0xbfd185ddfa1a7ed0
|
|
data8 0xbfd14c84c4dd6128, 0xbfd11307dad30b76
|
|
data8 0xbfd0d9670f6941fe, 0xbfd09fa235ba2020
|
|
data8 0xbfd0790adbb03009, 0xbfd03f09858c55fb
|
|
data8 0xbfd004e3a7c97cbd, 0xbfcf9532288fcf69
|
|
data8 0xbfcf205339208f27, 0xbfceab2a23a5b83e
|
|
data8 0xbfce5ce55fdd37a5, 0xbfcde73fe3b1480f
|
|
data8 0xbfcd714f44623927, 0xbfccfb1321b8c400
|
|
data8 0xbfccac163c770dc9, 0xbfcc355b67195dd0
|
|
data8 0xbfcbbe540a3f036f, 0xbfcb6ecf175f95e9
|
|
data8 0xbfcaf74751e1be33, 0xbfca7f71fb7bab9d
|
|
data8 0xbfca2f632320b86b, 0xbfc9b70ba539dfae
|
|
data8 0xbfc93e6587910444, 0xbfc8edcae8352b6c
|
|
data8 0xbfc874a0db01a719, 0xbfc7fb27199df16d
|
|
data8 0xbfc7a9fec7d05ddf, 0xbfc72fff456ac70d
|
|
data8 0xbfc6de7d66023dbc, 0xbfc663f6fac91316
|
|
data8 0xbfc6121ac74813cf, 0xbfc5970c478fff4a
|
|
data8 0xbfc51bab907a5c8a, 0xbfc4c93d33151b24
|
|
data8 0xbfc44d527fdadf55, 0xbfc3fa87be0f3a1b
|
|
data8 0xbfc3a797cd35d959, 0xbfc32ae9e278ae1a
|
|
data8 0xbfc2d79c6937efdd, 0xbfc25a619370d9dc
|
|
data8 0xbfc206b5bde2f8b8, 0xbfc188ecbd1d16be
|
|
data8 0xbfc134e1b489062e, 0xbfc0b6894488e95f
|
|
data8 0xbfc0621e2f556b5c, 0xbfc00d8c711a12cc
|
|
data8 0xbfbf1cd21257e18c, 0xbfbe72ec117fa5b2
|
|
data8 0xbfbdc8b7c49a1ddb, 0xbfbcc8d5e467b710
|
|
data8 0xbfbc1ddc9c39c7a1, 0xbfbb7294093cdd0f
|
|
data8 0xbfba7111df348494, 0xbfb9c501cdf75872
|
|
data8 0xbfb918a16e46335b, 0xbfb81579a73e83c6
|
|
data8 0xbfb7684f39f4ff2d, 0xbfb6bad3758efd87
|
|
data8 0xbfb60d060d7e41ac, 0xbfb507b836033bb7
|
|
data8 0xbfb4591d6310d85a, 0xbfb3aa2fdd27f1c3
|
|
data8 0xbfb2faef55ccb372, 0xbfb1f3723b4ae6db
|
|
data8 0xbfb14360d6136ffa, 0xbfb092fb594145c1
|
|
data8 0xbfafc482e8b48a7e, 0xbfae6265ace11ae4
|
|
data8 0xbfacff9e5c4341d0, 0xbfaaea3316095f72
|
|
data8 0xbfa985bfc3495194, 0xbfa820a01ac754cb
|
|
data8 0xbfa6bad3758efd87, 0xbfa554592bb8cd58
|
|
data8 0xbfa3ed3094685a26, 0xbfa2855905ca70f6
|
|
data8 0xbfa11cd1d5133413, 0xbf9dfd78881399f1
|
|
data8 0xbf9b28f618cc85df, 0xbf98530faa3c087b
|
|
data8 0xbf957bc3dddcd7fa, 0xbf92a3115322f9e6
|
|
data8 0xbf8f91ed4eef8370, 0xbf89dae4ec6b8b2e
|
|
data8 0xbf842106b1499209, 0xbf7cc89f97d67594
|
|
data8 0xbf71497accf7e11d, 0x0000000000000000
|
|
LOCAL_OBJECT_END(T_table)
|
|
|
|
|
|
.section .text
|
|
GLOBAL_LIBM_ENTRY(log2f)
|
|
|
|
{ .mfi
|
|
alloc r32=ar.pfs,1,4,4,0
|
|
// y=frcpa(x)
|
|
frcpa.s1 f6,p0=f1,f8
|
|
// will form significand of 1.5 (to test whether the index is 128 or above)
|
|
mov r24=0xc
|
|
}
|
|
{.mfi
|
|
nop.m 0
|
|
// normalize x
|
|
fma.s1 f7=f8,f1,f0
|
|
// r2 = pointer to C_1...C_6 followed by T_table
|
|
addl r2 = @ltoff(poly_coeffs), gp;;
|
|
}
|
|
{.mfi
|
|
// get significand
|
|
getf.sig r25=f8
|
|
// f8 denormal ?
|
|
fclass.m p8,p10=f8,0x9
|
|
// will form significand of 1.5 (to test whether the index is 128 or above)
|
|
shl r24=r24,60
|
|
}
|
|
{.mfi
|
|
mov r26=0x804
|
|
nop.f 0
|
|
// r23=bias-1
|
|
mov r23=0xfffe;;
|
|
}
|
|
|
|
{.mmf
|
|
getf.exp r29=f8
|
|
// load start address for C_1...C_6 followed by T_table
|
|
ld8 r2=[r2]
|
|
// will continue only for positive normal/denormal numbers
|
|
fclass.nm.unc p12,p7 = f8, 0x19 ;;
|
|
}
|
|
|
|
.pred.rel "mutex",p8,p10
|
|
{.mfi
|
|
// denormal input, repeat get significand (after normalization)
|
|
(p8) getf.sig r25=f7
|
|
// x=1 ?
|
|
fcmp.eq.s0 p6,p0=f8,f1
|
|
// get T_index
|
|
(p10) shr.u r28=r25,63-8
|
|
}
|
|
{.mfi
|
|
// f12=0.5
|
|
setf.exp f12=r23
|
|
nop.f 0
|
|
// r27=bias
|
|
mov r27=0xffff;;
|
|
}
|
|
|
|
{.mfb
|
|
// denormal input, repeat get exponent (after normalization)
|
|
(p8) getf.exp r29=f7
|
|
nop.f 0
|
|
(p12) br.cond.spnt SPECIAL_log2f
|
|
}
|
|
{.mfi
|
|
cmp.geu p12,p0=r25,r24
|
|
nop.f 0
|
|
mov r23=0xff;;
|
|
}
|
|
|
|
{.mfi
|
|
add r3=32,r2
|
|
// r=1-x*y
|
|
fms.s1 f6=f6,f8,f1
|
|
// r26=0x80400...0 (threshold for using polynomial approximation)
|
|
shl r26=r26,64-12
|
|
}
|
|
{.mfi
|
|
// load C_3, C_4
|
|
ldfpd f10,f11=[r2],16
|
|
nop.f 0
|
|
// r27=bias-1 (if index >=128, will add exponent+1)
|
|
(p12) mov r27=0xfffe;;
|
|
}
|
|
|
|
{.mfi
|
|
// load C_1
|
|
ldfe f14=[r2],32
|
|
// x=1, return 0
|
|
(p6) fma.s.s0 f8=f0,f0,f0
|
|
(p8) shr.u r28=r25,63-8
|
|
}
|
|
{.mib
|
|
// load C_2
|
|
ldfe f13=[r3]
|
|
// r29=exponent-bias
|
|
sub r29=r29,r27
|
|
// x=1, return
|
|
(p6) br.ret.spnt b0;;
|
|
}
|
|
|
|
|
|
{.mfi
|
|
// get T_index
|
|
and r28=r28,r23
|
|
fmerge.se f7=f1,f7
|
|
// if first 9 bits after leading 1 are all zero, then p8=1
|
|
cmp.ltu p8,p12=r25,r26;;
|
|
}
|
|
{.mfi
|
|
// f8=expon - bias
|
|
setf.sig f8=r29
|
|
nop.f 0
|
|
// get T address
|
|
shladd r2=r28,3,r2
|
|
}
|
|
{.mfi
|
|
// first 8 bits after leading 1 are all ones ?
|
|
cmp.eq p10,p0=r23,r28
|
|
// if first 8 bits after leading bit are 0, use polynomial approx. only
|
|
(p8) fms.s1 f6=f7,f1,f1
|
|
nop.i 0;;
|
|
}
|
|
{.mfi
|
|
//r26=1
|
|
mov r26=1
|
|
// if first 8 bits after leading 1 are all ones, use polynomial approx. only
|
|
(p10) fms.s1 f6=f7,f12,f1
|
|
nop.i 0;;
|
|
}
|
|
|
|
.pred.rel "mutex",p8,p12
|
|
{.mmf
|
|
// load T (unless first 9 bits after leading 1 are 0)
|
|
(p12) ldfd f12=[r2]
|
|
nop.m 0
|
|
// set T=0 (if first 9 bits after leading 1 are 0)
|
|
(p8) fma.s1 f12=f0,f0,f0;;
|
|
}
|
|
|
|
{.mfi
|
|
nop.m 0
|
|
// P34=C_3+C_4*r
|
|
fma.s1 f10=f11,f6,f10
|
|
// r26=2^{63}
|
|
shl r26=r26,63
|
|
}
|
|
{.mfi
|
|
nop.m 0
|
|
// r2=r*r
|
|
fma.s1 f11=f6,f6,f0
|
|
nop.i 0;;
|
|
}
|
|
{.mfi
|
|
// significand of x is 1 ?
|
|
cmp.eq p0,p6=r25,r26
|
|
// P12=C_1+C_2*r
|
|
fma.s1 f14=f13,f6,f14
|
|
nop.i 0;;
|
|
}
|
|
{.mfi
|
|
nop.m 0
|
|
// normalize additive term (l=exponent of x)
|
|
fcvt.xf f8=f8
|
|
// if significand(x)=1, return exponent (l)
|
|
nop.i 0;;
|
|
}
|
|
{.mfi
|
|
nop.m 0
|
|
// add T+l
|
|
(p6) fma.s1 f8=f8,f1,f12
|
|
nop.i 0
|
|
}
|
|
{.mfi
|
|
nop.m 0
|
|
// P14=P12+r2*P34
|
|
(p6) fma.s1 f13=f10,f11,f14
|
|
nop.i 0;;
|
|
}
|
|
|
|
{.mfb
|
|
nop.m 0
|
|
// result=T+l+r*P14
|
|
(p6) fma.s.s0 f8=f13,f6,f8
|
|
// return
|
|
br.ret.sptk b0;;
|
|
}
|
|
|
|
|
|
SPECIAL_log2f:
|
|
{.mfi
|
|
nop.m 0
|
|
// x=+Infinity ?
|
|
fclass.m p7,p0=f8,0x21
|
|
nop.i 0;;
|
|
}
|
|
{.mfi
|
|
nop.m 0
|
|
// x=+/-Zero ?
|
|
fclass.m p8,p0=f8,0x7
|
|
nop.i 0;;
|
|
}
|
|
{.mfi
|
|
nop.m 0
|
|
// x=-Infinity, -normal, -denormal ?
|
|
fclass.m p6,p0=f8,0x3a
|
|
nop.i 0;;
|
|
}
|
|
{.mfb
|
|
nop.m 0
|
|
// log2f(+Infinity)=+Infinity
|
|
nop.f 0
|
|
(p7) br.ret.spnt b0;;
|
|
}
|
|
{.mfi
|
|
(p8) mov GR_Parameter_TAG = 172
|
|
// log2f(+/-0)=-infinity, raises Divide by Zero
|
|
// set f8=-0
|
|
(p8) fmerge.ns f8=f0,f8
|
|
nop.i 0;;
|
|
}
|
|
{.mfb
|
|
nop.m 0
|
|
(p8) frcpa.s0 f8,p0=f1,f8
|
|
(p8) br.cond.sptk __libm_error_region;;
|
|
}
|
|
{.mfb
|
|
(p6) mov GR_Parameter_TAG = 173
|
|
// x<0: return NaN, raise Invalid
|
|
(p6) frcpa.s0 f8,p0=f0,f0
|
|
(p6) br.cond.sptk __libm_error_region;;
|
|
}
|
|
|
|
|
|
{.mfb
|
|
nop.m 0
|
|
// Remaining cases: NaNs
|
|
fma.s.s0 f8=f8,f1,f0
|
|
br.ret.sptk b0;;
|
|
}
|
|
|
|
GLOBAL_LIBM_END(log2f)
|
|
|
|
|
|
LOCAL_LIBM_ENTRY(__libm_error_region)
|
|
.prologue
|
|
{ .mfi
|
|
add GR_Parameter_Y=-32,sp // Parameter 2 value
|
|
nop.f 0
|
|
.save ar.pfs,GR_SAVE_PFS
|
|
mov GR_SAVE_PFS=ar.pfs // Save ar.pfs
|
|
}
|
|
{ .mfi
|
|
.fframe 64
|
|
add sp=-64,sp // Create new stack
|
|
nop.f 0
|
|
mov GR_SAVE_GP=gp // Save gp
|
|
};;
|
|
{ .mmi
|
|
stfs [GR_Parameter_Y] = FR_Y,16 // STORE Parameter 2 on stack
|
|
add GR_Parameter_X = 16,sp // Parameter 1 address
|
|
.save b0, GR_SAVE_B0
|
|
mov GR_SAVE_B0=b0 // Save b0
|
|
};;
|
|
.body
|
|
{ .mib
|
|
stfs [GR_Parameter_X] = FR_X // STORE Parameter 1 on stack
|
|
add GR_Parameter_RESULT = 0,GR_Parameter_Y // Parameter 3 address
|
|
nop.b 0
|
|
}
|
|
{ .mib
|
|
stfs [GR_Parameter_Y] = FR_RESULT // STORE Parameter 3 on stack
|
|
add GR_Parameter_Y = -16,GR_Parameter_Y
|
|
br.call.sptk b0=__libm_error_support# // Call error handling function
|
|
};;
|
|
{ .mmi
|
|
nop.m 0
|
|
nop.m 0
|
|
add GR_Parameter_RESULT = 48,sp
|
|
};;
|
|
{ .mmi
|
|
ldfs f8 = [GR_Parameter_RESULT] // Get return result off stack
|
|
.restore sp
|
|
add sp = 64,sp // Restore stack pointer
|
|
mov b0 = GR_SAVE_B0 // Restore return address
|
|
};;
|
|
{ .mib
|
|
mov gp = GR_SAVE_GP // Restore gp
|
|
mov ar.pfs = GR_SAVE_PFS // Restore ar.pfs
|
|
br.ret.sptk b0 // Return
|
|
};;
|
|
|
|
LOCAL_LIBM_END(__libm_error_region)
|
|
.type __libm_error_support#,@function
|
|
.global __libm_error_support#
|