sysdeps/ia64/fpu/e_log2f.S

   1 .file "log2f.s"
   2
   3
   4 // Copyright (c) 2000 - 2003, Intel Corporation
   5 // All rights reserved.
   6 //
   7 // Contributed 2000 by the Intel Numerics Group, Intel Corporation
   8 //
   9 // Redistribution and use in source and binary forms, with or without
  10 // modification, are permitted provided that the following conditions are
  11 // met:
  12 //
  13 // * Redistributions of source code must retain the above copyright
  14 // notice, this list of conditions and the following disclaimer.
  15 //
  16 // * Redistributions in binary form must reproduce the above copyright
  17 // notice, this list of conditions and the following disclaimer in the
  18 // documentation and/or other materials provided with the distribution.
  19 //
  20 // * The name of Intel Corporation may not be used to endorse or promote
  21 // products derived from this software without specific prior written
  22 // permission.
  23
  24 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
  25 // "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
  26 // LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
  27 // A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL INTEL OR ITS
  28 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
  29 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
  30 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
  31 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
  32 // OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY OR TORT (INCLUDING
  33 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
  34 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  35 //
  36 // Intel Corporation is the author of this code, and requests that all
  37 // problem reports or change requests be submitted to it directly at
  38 // http://www.intel.com/software/products/opensource/libraries/num.htm.
  39 //
  40 // History
  41 //==============================================================
  42 // 09/11/00 Initial version
  43 // 05/20/02 Cleaned up namespace and sf0 syntax
  44 // 02/10/03 Reordered header: .section, .global, .proc, .align
  45 //
  46 // API
  47 //==============================================================
  48 // float log2f(float)
  49 //
  50 // Overview of operation
  51 //==============================================================
  52 // Background
  53 //
  54 // Implementation
  55 //
  56 // Let x = 2^l * m, where     m=1.b1 b2 ... b8 b9 ... b52
  57 //     y=frcpa(m),   r=m*y-1, f=b1 b2 .. b8 (table index)
  58 // j=0 if f<128; j=1 if f>=128
  59 // T is a table that stores log2(1/y) (in entries 1..255) rounded to
  60 //   double extended precision; f is used as an index; T[255]=0
  61 //
  62 // If f=0 and b9=0, r is set to 2^{-8}* 0.b9 b10 ... b52 = m-1 (fractional part of m),
  63 //                  and 0 is used instead of T[0]
  64 //                  (polynomial evaluation only, for m=1+r, 0<=r<2^{-9})
  65 // If f=255, r is set to (m-2)/2  (T[255]=0, and only polynomial evaluation is used
  66 //                                 for m=2(1-r'), 0<=r'<2^{-9})
  67 //
  68 // log2f(x) is approximated as
  69 //     (l-j) + T[f] + (c1*r+c2*r^2+...+c6*r^6), if f>0
  70 //
  71
  72
  73 // Special values
  74 //==============================================================
  75 //  log2f(0)=-inf, raises Divide by Zero
  76 //  log2f(+inf)=inf
  77 //  log2f(x)=NaN,  raises Invalid if x<0
  78 //
  79
  80
  81 // Registers used
  82 //==============================================================
  83 //   f6-f14
  84 //   r2-r3, r23-r30
  85 //   p6,p7,p8,p12
  86 //
  87
  88
  89 GR_SAVE_B0                    = r33
  90 GR_SAVE_PFS                   = r34
  91 GR_SAVE_GP                    = r35 // This reg. can safely be used
  92 GR_SAVE_SP                    = r36
  93
  94 GR_Parameter_X                = r37
  95 GR_Parameter_Y                = r38
  96 GR_Parameter_RESULT           = r39
  97 GR_Parameter_TAG              = r40
  98
  99 FR_X             = f10
 100 FR_Y             = f1
 101 FR_RESULT        = f8
 102
 103
 104
 105
 106 // Data tables
 107 //==============================================================
 108
 109 RODATA
 110
 111 .align 16
 112
 113 LOCAL_OBJECT_START(poly_coeffs)
 114
 115 data8 0x3fdec709dc3a03fd, 0xbfd71547652b82fe //C_3 and C_4
 116 data8 0xb8aa3b295c17f0bc, 0x00003fff  // C_1
 117 data8 0xb8aa3b295c17f0bc, 0x0000bffe  // C_2
 118 LOCAL_OBJECT_END(poly_coeffs)
 119
 120
 121 LOCAL_OBJECT_START(T_table)
 122
 123 data8 0x3f671b0ea42e5fda, 0x3f815cfe8eaec830
 124 data8 0x3f8cfee70c5ce5dc, 0x3f94564a62192834
 125 data8 0x3f997723ace35766, 0x3f9f5923c69b54a1
 126 data8 0x3fa2a094a085d693, 0x3fa538941776b01e
 127 data8 0x3fa8324c9b914bc7, 0x3faacf54ce07d7e9
 128 data8 0x3fadced958dadc12, 0x3fb0387efbca869e
 129 data8 0x3fb18ac6067479c0, 0x3fb30edd3e13530d
 130 data8 0x3fb463c15936464e, 0x3fb5b9e13c3fa21d
 131 data8 0x3fb7113f3259e07a, 0x3fb869dd8d1b2035
 132 data8 0x3fb9c3bea49d3214, 0x3fbb1ee4d7961701
 133 data8 0x3fbc7b528b70f1c5, 0x3fbdd90a2c676ed4
 134 data8 0x3fbf05d4976c2028, 0x3fc032fbbaee6d65
 135 data8 0x3fc0e3b5a9f3284a, 0x3fc195195c7d125b
 136 data8 0x3fc22dadc2ab3497, 0x3fc2e050231df57d
 137 data8 0x3fc379f79c2b255b, 0x3fc42ddd2ba1b4a9
 138 data8 0x3fc4c89b9e6807f5, 0x3fc563dc29ffacb2
 139 data8 0x3fc619a25f5d798d, 0x3fc6b5ffbf367644
 140 data8 0x3fc752e1f660f8d6, 0x3fc7f049e753e7cf
 141 data8 0x3fc8a8980abfbd32, 0x3fc94724cca657be
 142 data8 0x3fc9e63a24971f46, 0x3fca85d8feb202f7
 143 data8 0x3fcb2602497d5346, 0x3fcbc6b6f5ee1c9b
 144 data8 0x3fcc67f7f770a67e, 0x3fcceec4b2234fba
 145 data8 0x3fcd91097ad13982, 0x3fce33dd57f3d335
 146 data8 0x3fced74146bc7b10, 0x3fcf7b3646fef683
 147 data8 0x3fd00223a943dc19, 0x3fd054a474bf0eb7
 148 data8 0x3fd0999d9b9259a1, 0x3fd0eca66d3b2581
 149 data8 0x3fd13ffa2e85b475, 0x3fd185a444fa0a7b
 150 data8 0x3fd1cb8312f27eff, 0x3fd21fa1441ce5e8
 151 data8 0x3fd265f526e603cb, 0x3fd2baa0c34be1ec
 152 data8 0x3fd3016b45de21ce, 0x3fd3486c38aa29a8
 153 data8 0x3fd38fa3efaa8262, 0x3fd3e562c0816a02
 154 data8 0x3fd42d141f53b646, 0x3fd474fd543f222c
 155 data8 0x3fd4bd1eb680e548, 0x3fd505789e234bd1
 156 data8 0x3fd54e0b64003b70, 0x3fd596d761c3c1f0
 157 data8 0x3fd5dfdcf1eeae0e, 0x3fd6291c6fd9329c
 158 data8 0x3fd6729637b59418, 0x3fd6bc4aa692e0fd
 159 data8 0x3fd7063a1a5fb4f2, 0x3fd75064f1ed0715
 160 data8 0x3fd79acb8cf10390, 0x3fd7d67c1e43ae5c
 161 data8 0x3fd8214f4068afa7, 0x3fd86c5f36dea3dc
 162 data8 0x3fd8b7ac64dd7f9d, 0x3fd8f4167a0c6f92
 163 data8 0x3fd93fd2d5e1bf1d, 0x3fd98bcd84296946
 164 data8 0x3fd9c8c333e6e9a5, 0x3fda152f142981b4
 165 data8 0x3fda527fd95fd8ff, 0x3fda9f5e3edeb9e6
 166 data8 0x3fdadd0b2b5755a7, 0x3fdb2a5d6f51ff83
 167 data8 0x3fdb686799b00be3, 0x3fdbb62f1b887cd8
 168 data8 0x3fdbf4979f666668, 0x3fdc332a6e8399d4
 169 data8 0x3fdc819dc2d45fe4, 0x3fdcc0908e19b7bd
 170 data8 0x3fdcffae611ad12b, 0x3fdd3ef776d43ff4
 171 data8 0x3fdd8e5002710128, 0x3fddcdfb486cb9a1
 172 data8 0x3fde0dd294245fe4, 0x3fde4dd622a28840
 173 data8 0x3fde8e06317114f0, 0x3fdece62fe9a9915
 174 data8 0x3fdf1f164a15389a, 0x3fdf5fd8a9063e35
 175 data8 0x3fdfa0c8937e7d5d, 0x3fdfe1e649bb6335
 176 data8 0x3fe011990641535a, 0x3fe032560e91e59e
 177 data8 0x3fe0532a5ebcd44a, 0x3fe0741617f5fc28
 178 data8 0x3fe08cd653f38839, 0x3fe0adeb55c1103b
 179 data8 0x3fe0cf181d5d1dd0, 0x3fe0f05ccd0aced7
 180 data8 0x3fe111b9875788ab, 0x3fe1332e6f1bcf73
 181 data8 0x3fe154bba77c2088, 0x3fe16df59bfa06c1
 182 data8 0x3fe18fadb6e2d3c2, 0x3fe1b17e849adc26
 183 data8 0x3fe1caeb6a0de814, 0x3fe1ece7c830eec9
 184 data8 0x3fe20efd3dae01df, 0x3fe2289de375d901
 185 data8 0x3fe24adf9b6a6fe0, 0x3fe26d3ad1aebcfc
 186 data8 0x3fe287100c2771f4, 0x3fe2a9983b3c1b28
 187 data8 0xbfda78e146f7bef4, 0xbfda33760a7f6051
 188 data8 0xbfd9ff43476fb5f7, 0xbfd9b97c3c4eec8f
 189 data8 0xbfd98504431717fc, 0xbfd93ee07535f967
 190 data8 0xbfd90a228d5712b2, 0xbfd8c3a104cb24f5
 191 data8 0xbfd88e9c72e0b226, 0xbfd847bc33d8618e
 192 data8 0xbfd812703988bb69, 0xbfd7dd0569c04bff
 193 data8 0xbfd7959c202292f1, 0xbfd75fe8d2c5d48f
 194 data8 0xbfd72a1637cbc183, 0xbfd6e221cd9d0cde
 195 data8 0xbfd6ac059985503b, 0xbfd675c99ce81f92
 196 data8 0xbfd63f6db2590482, 0xbfd5f6c138136489
 197 data8 0xbfd5c01a39fbd688, 0xbfd58952cf519193
 198 data8 0xbfd5526ad18493ce, 0xbfd51b6219bfe6ea
 199 data8 0xbfd4d1cdf8b4846f, 0xbfd49a784bcd1b8b
 200 data8 0xbfd4630161832547, 0xbfd42b6911cf5465
 201 data8 0xbfd3f3af3461e1c4, 0xbfd3bbd3a0a1dcfb
 202 data8 0xbfd383d62dac7ae7, 0xbfd34bb6b2546218
 203 data8 0xbfd313750520f520, 0xbfd2db10fc4d9aaf
 204 data8 0xbfd2a28a6dc90387, 0xbfd269e12f346e2c
 205 data8 0xbfd2311515e2e855, 0xbfd1f825f6d88e13
 206 data8 0xbfd1bf13a6c9c69f, 0xbfd185ddfa1a7ed0
 207 data8 0xbfd14c84c4dd6128, 0xbfd11307dad30b76
 208 data8 0xbfd0d9670f6941fe, 0xbfd09fa235ba2020
 209 data8 0xbfd0790adbb03009, 0xbfd03f09858c55fb
 210 data8 0xbfd004e3a7c97cbd, 0xbfcf9532288fcf69
 211 data8 0xbfcf205339208f27, 0xbfceab2a23a5b83e
 212 data8 0xbfce5ce55fdd37a5, 0xbfcde73fe3b1480f
 213 data8 0xbfcd714f44623927, 0xbfccfb1321b8c400
 214 data8 0xbfccac163c770dc9, 0xbfcc355b67195dd0
 215 data8 0xbfcbbe540a3f036f, 0xbfcb6ecf175f95e9
 216 data8 0xbfcaf74751e1be33, 0xbfca7f71fb7bab9d
 217 data8 0xbfca2f632320b86b, 0xbfc9b70ba539dfae
 218 data8 0xbfc93e6587910444, 0xbfc8edcae8352b6c
 219 data8 0xbfc874a0db01a719, 0xbfc7fb27199df16d
 220 data8 0xbfc7a9fec7d05ddf, 0xbfc72fff456ac70d
 221 data8 0xbfc6de7d66023dbc, 0xbfc663f6fac91316
 222 data8 0xbfc6121ac74813cf, 0xbfc5970c478fff4a
 223 data8 0xbfc51bab907a5c8a, 0xbfc4c93d33151b24
 224 data8 0xbfc44d527fdadf55, 0xbfc3fa87be0f3a1b
 225 data8 0xbfc3a797cd35d959, 0xbfc32ae9e278ae1a
 226 data8 0xbfc2d79c6937efdd, 0xbfc25a619370d9dc
 227 data8 0xbfc206b5bde2f8b8, 0xbfc188ecbd1d16be
 228 data8 0xbfc134e1b489062e, 0xbfc0b6894488e95f
 229 data8 0xbfc0621e2f556b5c, 0xbfc00d8c711a12cc
 230 data8 0xbfbf1cd21257e18c, 0xbfbe72ec117fa5b2
 231 data8 0xbfbdc8b7c49a1ddb, 0xbfbcc8d5e467b710
 232 data8 0xbfbc1ddc9c39c7a1, 0xbfbb7294093cdd0f
 233 data8 0xbfba7111df348494, 0xbfb9c501cdf75872
 234 data8 0xbfb918a16e46335b, 0xbfb81579a73e83c6
 235 data8 0xbfb7684f39f4ff2d, 0xbfb6bad3758efd87
 236 data8 0xbfb60d060d7e41ac, 0xbfb507b836033bb7
 237 data8 0xbfb4591d6310d85a, 0xbfb3aa2fdd27f1c3
 238 data8 0xbfb2faef55ccb372, 0xbfb1f3723b4ae6db
 239 data8 0xbfb14360d6136ffa, 0xbfb092fb594145c1
 240 data8 0xbfafc482e8b48a7e, 0xbfae6265ace11ae4
 241 data8 0xbfacff9e5c4341d0, 0xbfaaea3316095f72
 242 data8 0xbfa985bfc3495194, 0xbfa820a01ac754cb
 243 data8 0xbfa6bad3758efd87, 0xbfa554592bb8cd58
 244 data8 0xbfa3ed3094685a26, 0xbfa2855905ca70f6
 245 data8 0xbfa11cd1d5133413, 0xbf9dfd78881399f1
 246 data8 0xbf9b28f618cc85df, 0xbf98530faa3c087b
 247 data8 0xbf957bc3dddcd7fa, 0xbf92a3115322f9e6
 248 data8 0xbf8f91ed4eef8370, 0xbf89dae4ec6b8b2e
 249 data8 0xbf842106b1499209, 0xbf7cc89f97d67594
 250 data8 0xbf71497accf7e11d, 0x0000000000000000
 251 LOCAL_OBJECT_END(T_table)
 252
 253
 254 .section .text
 255 GLOBAL_LIBM_ENTRY(__log2f)
 256
 257 { .mfi
 258   alloc r32=ar.pfs,1,4,4,0
 259   // y=frcpa(x)
 260   frcpa.s1 f6,p0=f1,f8
 261   // will form significand of 1.5 (to test whether the index is 128 or above)
 262   mov r24=0xc
 263 }
 264 {.mfi
 265   nop.m 0
 266   // normalize x
 267   fma.s1 f7=f8,f1,f0
 268   // r2 = pointer to C_1...C_6 followed by T_table
 269   addl r2 = @ltoff(poly_coeffs), gp;;
 270 }
 271 {.mfi
 272   // get significand
 273   getf.sig r25=f8
 274   // f8 denormal ?
 275   fclass.m p8,p10=f8,0x9
 276   // will form significand of 1.5 (to test whether the index is 128 or above)
 277   shl r24=r24,60
 278 }
 279 {.mfi
 280   mov r26=0x804
 281   nop.f 0
 282   // r23=bias-1
 283   mov r23=0xfffe;;
 284 }
 285
 286 {.mmf
 287   getf.exp r29=f8
 288   // load start address for C_1...C_6 followed by T_table
 289   ld8 r2=[r2]
 290   // will continue only for positive normal/denormal numbers
 291   fclass.nm.unc p12,p7 = f8, 0x19 ;;
 292 }
 293
 294 .pred.rel "mutex",p8,p10
 295 {.mfi
 296   // denormal input, repeat get significand (after normalization)
 297   (p8) getf.sig r25=f7
 298   // x=1 ?
 299   fcmp.eq.s0 p6,p0=f8,f1
 300   // get T_index
 301   (p10) shr.u r28=r25,63-8
 302 }
 303 {.mfi
 304   // f12=0.5
 305   setf.exp f12=r23
 306   nop.f 0
 307   // r27=bias
 308   mov r27=0xffff;;
 309 }
 310
 311 {.mfb
 312   // denormal input, repeat get exponent (after normalization)
 313   (p8) getf.exp r29=f7
 314   nop.f 0
 315   (p12) br.cond.spnt SPECIAL_log2f
 316 }
 317 {.mfi
 318   cmp.geu p12,p0=r25,r24
 319   nop.f 0
 320   mov r23=0xff;;
 321 }
 322
 323 {.mfi
 324   add r3=32,r2
 325   // r=1-x*y
 326   fms.s1 f6=f6,f8,f1
 327   // r26=0x80400...0 (threshold for using polynomial approximation)
 328   shl r26=r26,64-12
 329 }
 330 {.mfi
 331   // load C_3, C_4
 332   ldfpd f10,f11=[r2],16
 333   nop.f 0
 334   // r27=bias-1 (if index >=128, will add exponent+1)
 335   (p12) mov r27=0xfffe;;
 336 }
 337
 338 {.mfi
 339   // load C_1
 340   ldfe f14=[r2],32
 341   // x=1, return 0
 342   (p6) fma.s.s0 f8=f0,f0,f0
 343   (p8) shr.u r28=r25,63-8
 344 }
 345 {.mib
 346   // load C_2
 347   ldfe f13=[r3]
 348   // r29=exponent-bias
 349   sub r29=r29,r27
 350   // x=1, return
 351   (p6) br.ret.spnt b0;;
 352 }
 353
 354
 355 {.mfi
 356   // get T_index
 357   and r28=r28,r23
 358   fmerge.se f7=f1,f7
 359   // if first 9 bits after leading 1 are all zero, then p8=1
 360   cmp.ltu p8,p12=r25,r26;;
 361 }
 362 {.mfi
 363   // f8=expon - bias
 364   setf.sig f8=r29
 365   nop.f 0
 366   // get T address
 367   shladd r2=r28,3,r2
 368 }
 369 {.mfi
 370   // first 8 bits after leading 1 are all ones ?
 371   cmp.eq p10,p0=r23,r28
 372   // if first 8 bits after leading bit are 0, use polynomial approx. only
 373   (p8) fms.s1 f6=f7,f1,f1
 374   nop.i 0;;
 375 }
 376 {.mfi
 377   //r26=1
 378   mov r26=1
 379   // if first 8 bits after leading 1 are all ones, use polynomial approx. only
 380   (p10) fms.s1 f6=f7,f12,f1
 381   nop.i 0;;
 382 }
 383
 384 .pred.rel "mutex",p8,p12
 385 {.mmf
 386   // load T (unless first 9 bits after leading 1 are 0)
 387   (p12) ldfd f12=[r2]
 388   nop.m 0
 389   // set T=0 (if first 9 bits after leading 1 are 0)
 390   (p8) fma.s1 f12=f0,f0,f0;;
 391 }
 392
 393 {.mfi
 394   nop.m 0
 395   // P34=C_3+C_4*r
 396   fma.s1 f10=f11,f6,f10
 397   // r26=2^{63}
 398   shl r26=r26,63
 399 }
 400 {.mfi
 401   nop.m 0
 402   // r2=r*r
 403   fma.s1 f11=f6,f6,f0
 404   nop.i 0;;
 405 }
 406 {.mfi
 407   // significand of x is 1 ?
 408   cmp.eq p0,p6=r25,r26
 409   // P12=C_1+C_2*r
 410   fma.s1 f14=f13,f6,f14
 411   nop.i 0;;
 412 }
 413 {.mfi
 414   nop.m 0
 415   // normalize additive term (l=exponent of x)
 416   fcvt.xf f8=f8
 417   // if significand(x)=1, return exponent (l)
 418   nop.i 0;;
 419 }
 420 {.mfi
 421   nop.m 0
 422   // add T+l
 423   (p6) fma.s1 f8=f8,f1,f12
 424   nop.i 0
 425 }
 426 {.mfi
 427   nop.m 0
 428   // P14=P12+r2*P34
 429   (p6) fma.s1 f13=f10,f11,f14
 430   nop.i 0;;
 431 }
 432
 433 {.mfb
 434   nop.m 0
 435   // result=T+l+r*P14
 436   (p6) fma.s.s0 f8=f13,f6,f8
 437   // return
 438   br.ret.sptk b0;;
 439 }
 440
 441
 442 SPECIAL_log2f:
 443 {.mfi
 444   nop.m 0
 445   // x=+Infinity ?
 446   fclass.m p7,p0=f8,0x21
 447   nop.i 0;;
 448 }
 449 {.mfi
 450   nop.m 0
 451   // x=+/-Zero ?
 452   fclass.m p8,p0=f8,0x7
 453   nop.i 0;;
 454 }
 455 {.mfi
 456   nop.m 0
 457   // x=-Infinity, -normal, -denormal ?
 458   fclass.m p6,p0=f8,0x3a
 459   nop.i 0;;
 460 }
 461 {.mfb
 462   nop.m 0
 463   // log2f(+Infinity)=+Infinity
 464   nop.f 0
 465   (p7) br.ret.spnt b0;;
 466 }
 467 {.mfi
 468   (p8) mov GR_Parameter_TAG = 172
 469   // log2f(+/-0)=-infinity, raises Divide by Zero
 470   // set f8=-0
 471   (p8) fmerge.ns f8=f0,f8
 472   nop.i 0;;
 473 }
 474 {.mfb
 475   nop.m 0
 476   (p8) frcpa.s0 f8,p0=f1,f8
 477   (p8) br.cond.sptk __libm_error_region;;
 478 }
 479 {.mfb
 480   (p6) mov GR_Parameter_TAG = 173
 481   // x<0: return NaN, raise Invalid
 482   (p6) frcpa.s0 f8,p0=f0,f0
 483   (p6) br.cond.sptk __libm_error_region;;
 484 }
 485
 486
 487 {.mfb
 488   nop.m 0
 489   // Remaining cases: NaNs
 490   fma.s.s0 f8=f8,f1,f0
 491   br.ret.sptk b0;;
 492 }
 493
 494 GLOBAL_LIBM_END(__log2f)
 495 #ifdef SHARED
 496 .symver __log2f,log2f@@GLIBC_2.27
 497 .weak __log2f_compat
 498 .set __log2f_compat,__log2f
 499 .symver __log2f_compat,log2f@GLIBC_2.2
 500 #endif
 501
 502
 503 LOCAL_LIBM_ENTRY(__libm_error_region)
 504 .prologue
 505 { .mfi
 506         add   GR_Parameter_Y=-32,sp             // Parameter 2 value
 507         nop.f 0
 508 .save   ar.pfs,GR_SAVE_PFS
 509         mov  GR_SAVE_PFS=ar.pfs                 // Save ar.pfs
 510 }
 511 { .mfi
 512 .fframe 64
 513         add sp=-64,sp                           // Create new stack
 514         nop.f 0
 515         mov GR_SAVE_GP=gp                       // Save gp
 516 };;
 517 { .mmi
 518         stfs [GR_Parameter_Y] = FR_Y,16         // STORE Parameter 2 on stack
 519         add GR_Parameter_X = 16,sp              // Parameter 1 address
 520 .save   b0, GR_SAVE_B0
 521         mov GR_SAVE_B0=b0                       // Save b0
 522 };;
 523 .body
 524 { .mib
 525         stfs [GR_Parameter_X] = FR_X                  // STORE Parameter 1 on stack
 526         add   GR_Parameter_RESULT = 0,GR_Parameter_Y  // Parameter 3 address
 527         nop.b 0
 528 }
 529 { .mib
 530         stfs [GR_Parameter_Y] = FR_RESULT             // STORE Parameter 3 on stack
 531         add   GR_Parameter_Y = -16,GR_Parameter_Y
 532         br.call.sptk b0=__libm_error_support#         // Call error handling function
 533 };;
 534 { .mmi
 535         nop.m 0
 536         nop.m 0
 537         add   GR_Parameter_RESULT = 48,sp
 538 };;
 539 { .mmi
 540         ldfs  f8 = [GR_Parameter_RESULT]       // Get return result off stack
 541 .restore sp
 542         add   sp = 64,sp                       // Restore stack pointer
 543         mov   b0 = GR_SAVE_B0                  // Restore return address
 544 };;
 545 { .mib
 546         mov   gp = GR_SAVE_GP                  // Restore gp
 547         mov   ar.pfs = GR_SAVE_PFS             // Restore ar.pfs
 548         br.ret.sptk     b0                     // Return
 549 };;
 550
 551 LOCAL_LIBM_END(__libm_error_region)
 552 .type   __libm_error_support#,@function
 553 .global __libm_error_support#