s_cbrtl.s

Glibc 2.3.2源代码(解压后有100多M)

data4 0x9d8eeccb, 0x9ea3d011, 0x1f231fdf, 0x9f1dbdfa
data4 0x1e7507a3, 0x1ec42614, 0x9e8693cb, 0x9ec68398
data4 0x1d5b05fb, 0x1de32119, 0x9f003429, 0x9ec16d92
data4 0x9f095315, 0x9f119d2c, 0x9ed0c984, 0x9f090662
data4 0x9e59aa1f, 0x9ed4e64a, 0x9f2798a7, 0x9f23624d
data4 0x1e0467d9, 0x1f22e7e7, 0x1e915256, 0x9cb4df70
data4 0x9e6f687c, 0x9e3c35e5, 0x9e5757ab, 0x9f031fa1
data4 0x1f25bff7, 0x1f0e58c2, 0x1ef3ce04, 0x1f002ecb
data4 0x9ebdc836, 0x9ed657dd, 0x9f149441, 0x9e8544b2
data4 0x1cd8ff1e, 0x1e9bb463, 0x1eaa1c5c, 0x1f200c1a
data4 0x1edbfbaf, 0x1f18724d, 0x9ed63c22, 0x9f08e045
data4 0x1f13ad07, 0x9e949311, 0x9f0c50d4, 0x1e824516
data4 0x1d5e52ba, 0x1d583fbd, 0x1e3b60a9, 0x9effe6d3
data4 0x1f0d0508, 0x1f00be77, 0x9e404bfa, 0x9e1ca381
data4 0x9f084dd8, 0x9e6db85d, 0x1db698e4, 0x9ebd1871
data4 0x9ecc2679, 0x1ee68442, 0x1edb1050, 0x9dbc96a4
data4 0x9f27c1f4, 0x1c99b756, 0x1eb4400a, 0x9f24390a
data4 0x1d927875, 0x9f074faa, 0x1e9dc2c3, 0x1f13c0d2
data4 0x1e3c9685, 0x9e6b6f75, 0x9db9cb31, 0x1ea5f3aa
data4 0x9d992c61, 0x1f1015e4, 0x1f194f70, 0x9e19d2b3
data4 0x9d89116c, 0x1f23cd35, 0x1e33d3a2, 0x1ee331b8
data4 0x1d5ba7ec, 0x9f273788, 0x9e6907f4, 0x9ed5f912
data4 0x9edd458d, 0x1e2ca7b2, 0x1ef81fe4, 0x1dc7ade6
data4 0x1e876e51, 0x9f04ec89, 0x1f1da63a, 0x1ec02bd0
data4 0x9e71326f, 0x1e7847b4, 0x1f0de618, 0x9e036cb6
data4 0x1eec61e2, 0x1ef1758b, 0x9ee880a3, 0x1ed269d7
data4 0x1e27edd3, 0x9e8a81a1, 0x1eacb84d, 0x9e1aad37
data4 0x1f1aa8f7, 0x1e9bbd90, 0x1ea1b61f, 0x9ed41c2f
data4 0x1dbb5dd6, 0x1f0ec733, 0x9df06b1b, 0x1e06fef1
data4 0x9edede3a, 0x1edeb5e2, 0x1f0e63ee, 0x9db316bb
data4 0x9efc1ad3, 0x1f01fbb5, 0x9cc0d078, 0x1ea28b36
data4 0x9e9dd205, 0x9e791534, 0x1da1c8d5, 0x9e8195cc
data4 0x1f0681a4, 0x1eeaf1e2, 0x9ef83b37, 0x9f22a92b
data4 0x1eabc4ce, 0x1f10eefb, 0x1e06d9aa, 0x1e7cacd5
data4 0x1f1ea087, 0x1eb21983, 0x9f100c78, 0x1e840abe
data4 0x9efab66c, 0x1f183fa8, 0x9e84ee68, 0x9eea083d
data4 0x9ee23a74, 0x1f1351d7, 0x9ec5d42a, 0x9f071f57
data4 0x9ef578d9, 0x9f1aa7e7, 0x1eb02044, 0x1f151a2e
data4 0x9c0dc8b2, 0x9ef4087a, 0x1ec12b93, 0x1c1a946b
data4 0x1e89946f, 0x9dafe8c3, 0x1d295288, 0x9e8497ab
data4 0x1ec000c6, 0x1e102f29, 0x1e542256, 0x1e67d44d
data4 0x1ef688d8, 0x1f0e0f29, 0x1e67861f, 0x1e869748
data4 0x1ee6aa6e, 0x9e4d228b, 0x9e50be5b, 0x1e9fe225
data4 0x9ea34102, 0x9e628a3b, 0x9ed9fd83, 0x1ecd7109
data4 0x1f1864ff, 0x1ea19b76, 0x1db0d1c9, 0x9dff519b
data4 0x1e8fea71, 0x9ee82e9a, 0x9f08919b, 0x9ef5c8ae
data4 0x9ee446a4, 0x1ea59444, 0x1eb74230, 0x1ea13fbf
data4 0x9ea6a3ea, 0x1e5f2797, 0x9e0adb07, 0x9d3adadd
data4 0x1ebf2ee2, 0x1da19bfa, 0x1e8dea6d, 0x1ec4fea9
data4 0x1e669f22, 0x1dc5f919, 0x9ed25caa, 0x1ee475b1
data4 0x1ed0603e, 0x9eacb35c, 0x1dc00b27, 0x1e2f9991
data4 0x1e7b0406, 0x1eaa3387, 0x9d865bde, 0x1eb78a48
data4 0x1c40ae2e, 0x1ee9838b, 0x9f0f0d7f, 0x1e3e5d26
data4 0x1e99e7a6, 0x9e681ccf, 0x9e93ed65, 0x9eeb6a66
data4 0x1e29e9af, 0x9e96f923, 0x9e74f11d, 0x9f1474da
data4 0x1eec2ea7, 0x1ebf7aa3, 0x9c25dcca, 0x9f0553c2
data4 0x9e599efd, 0x1d2ab490, 0x1e95d7cd, 0x9ee4b20e
data4 0x9d988ce5, 0x9ef9787e, 0x9dbbba5b, 0x9f12c304
data4 0x1e3b9d70, 0x1e7bcae8, 0x9d98bb6e, 0x9e8e6b01
data4 0x9f07d03b, 0x9d67c822, 0x9f0ef69e, 0x1c7c0fe3
data4 0x9e9bfbb9, 0x9e83b84b, 0x1efbf15e, 0x9ecfa6a6
data4 0x9c91158e, 0x9ecf6770, 0x1ee1e3a8, 0x9dc95ec0
data4 0x1ef603f7, 0x1d5e52ba, 0x1c477d1b, 0x9e955cd8
data4 0x1ed665b0, 0x9e8376c4, 0x9c0ee88e, 0x1e8c989e
data4 0x1ea2df29, 0x9d961e5c, 0x1e101813, 0x1e7fffff
data4 0x9e5abff4, 0x1dbddd71, 0x1eb69100, 0x1e71f114
data4 0x1e9ca798, 0x1ef62c8d, 0x9db4e55a, 0x1dbe69ce
data4 0x9ef1c01f, 0x1f044a2a, 0x9eb9e0d7, 0x9ee59745
data4 0x9e874803, 0x1ea0b418, 0x9e13572a, 0x1ddbb3a2
data4 0x9ec0e391, 0x1e89fba1, 0x1ee8b261, 0x9e5d25f0
data4 0x9ef222cb, 0x9ef135ec, 0x1ea04b9a, 0x9f04291f
data4 0x9e969254, 0x9ee32f08, 0x9ed909d3, 0x9e362640
data4 0x9ec20735, 0x1e50131b, 0x9ed4e049, 0x1ee8e817
data4 0x1e1e09c0, 0x9ea643c5, 0x9e5a1ab6, 0x9e389059
data4 0x1e560947, 0x1d02b877, 0x1e4475ab, 0x9ea9aaf6
data4 0x1e95bc5e, 0x1eaf6afd, 0x1d43067d, 0x9d043821
data4 0x9e97baa9, 0x1de5c4f9, 0x9e9a0069, 0x9e1b9944
data4 0x1eb13686, 0x9eb907eb, 0x1e059589, 0x1cbd0f93
data4 0x9eb7e6ae, 0x1e9fa175, 0x1ee5bdf4, 0x1e8052f7
data4 0x9c80d1e3, 0x1bfbe28e, 0x9e672b3b, 0x9ecacf19
data4 0x9e3c04be, 0x1dfe8c5c, 0x1e1ba9cb, 0x1eb40b1e
data4 0x1ec7e7f6, 0x9d0d45b3, 0x1ef0113b, 0x9a155fa3
data4 0x1e28ec3b, 0x1e7ca8df, 0x9d2f91b4, 0x1eccd9ed
data4 0x9ed943bc, 0x9ccaab19, 0x9e8a5c58, 0x1ec3bca8
data4 0x1ed78dc7, 0x9ed391a8, 0x9e938f6e, 0x9ec4a030
data4 0x9e80346e, 0x1e7a4686, 0x9e284315, 0x9e39584c
data4 0x1ebdc9b4, 0x9e9cfce5, 0x9ef55c65, 0x1e2941e7
data4 0x9efbe59f, 0x1d87c41b, 0x1e40befc, 0x1e3d05b5
data4 0x1de9ea67, 0x1ec9a21c, 0x1decb69a, 0x1df6e75a
data4 0x9e8030ab, 0x9db20540, 0x9ef1e977, 0x1e3cdc43
data4 0x1e0492b0, 0x9e91d872, 0x1e775346, 0x9e939978
data4 0x1eb2714e, 0x1e49a203, 0x9e10195a, 0x1ef1ffc3
data4 0x9ea8b709, 0x9e832e27, 0x1ed5ac3b, 0x1edb20a6
data4 0x1e4dbd4e, 0x1efbb932, 0x1d8170ec, 0x1e6c4849
data4 0x1f008e17, 0x1e8000c4, 0x1d855ecf, 0x9e37cb85
data4 0x1ecffdf5, 0x1eba6519, 0x9edbe600, 0x1ea3e5e7
data4 0x1ed4fb39, 0x1f00be77, 0x1e6f4484, 0x9e9e7107
data4 0x9e30b29d, 0x9ee6e174, 0x1e3a2656, 0x9dd72f3f
data4 0x9ee12138, 0x1ed16fed, 0x9ece8a02, 0x9ca5b249
data4 0x9eafd508, 0x9ef0e9fc, 0x1d1307ac, 0x1eecee20
data4 0x1cf60c6f, 0x9d556216, 0x9eaed175, 0x9ec919f4
data4 0x1ec2c988, 0x1cd82772, 0x9dc99456, 0x1eab0467
data4 0x1e89b36f, 0x1c757944, 0x1eef9abd, 0x9e98664d
ASM_SIZE_DIRECTIVE(D_table)





.align 32
.global cbrtl#

.section .text
.proc  cbrtl#
.align 32
cbrtl: 


{ .mfi
  getf.sig r3=f8
     // will continue only for normal/denormal numbers          
(p0)  fclass.nm.unc p12,p7 = f8, 0x1b           
  // r2 = pointer to C_1...C_6 followed by T_table
  addl r2 = @ltoff(poly_coeffs), gp;;
}
{.mfi
  // r29=2/3*bias -63=0xaaaa-0x3f=0xaa6b
  mov r29=0xaa6b    
  // normalize a
  fma.s1 f14=f8,f1,f0
  // r27 = pointer to D table
  addl r27 = @ltoff(D_table), gp;;
}
{.mib
  nop.m 0
  (p7) cmp.eq p12,p0=r3,r0
  nop.b 0;;
}
{.mfb
  // load start address for C_1...C_6 followed by T_table
  ld8 r2=[r2]
  (p12) fma.s0 f8=f8,f1,f0
  (p12) br.ret.spnt b0;;
}
{.mmf 
  // load C_1
  ldfe f7=[r2],16
  // load start address of D table
  ld8 r27=[r27]
  // y=frcpa(a)
  frcpa.s0 f8,p6=f1,f8;;
}
{.mmi
  // load C_2
  ldfe f9=[r2],16;;
   // load C_3, C_4
  ldfpd f10,f11=[r2],16
  nop.i 0;;
}
{.mmi
  // get normalized significand
  getf.sig r23=f14
  // get exponent
  getf.exp r24=f14
  mov r25=0x20000;;
}
{.mii
  // get r26=sign
  and r26=r24,r25
  // eliminate leading 1 from r23=2nd table index
  shl r23=r23,1
  // eliminate sign from exponent (r25)
  andcm r25=r24,r25;;
}
{.mfi
  // load C_5,C_6
  (p6) ldfpd f12,f13=[r2],16
  // r=1-a*y
  (p6) fnma.s1 f6=f8,f14,f1
  // 1: exponent*=5;  // (2^{16}-1)/3=0x5555
  shladd r24=r25,2,r25;;
}
{.mib
  // r30=(5*expon)*16
  shladd r30=r24,4,r0
  // r28=3*exponent
  shladd r28=r25,1,r25
  nop.b 0;;
}
{.mmi
  // r28=6*exponent
  shladd r28=r28,1,r0
  // r24=17*expon
  add r24=r24,r30
  // r23=2nd table index (8 bits)
  shr.u r23=r23,56;;
}
{.mmi
  // adjust T_table pointer by 2nd index
  shladd r2=r23,3,r2
  // adjust D_table pointer by 2nd index
  shladd r27=r23,2,r27
  // r30=(17*expon)*16^2
  shl r30=r24,8;;
}
{.mmi
  // r24=expon*(2^16-1)/3
  add r24=r24,r30;;
  // r24=expon*(2^20+2)/3=expon*0x55556
  shladd r24=r24,4,r28
  nop.i 0;;
}
{.mii
  nop.m 0
  // r24=floor(expon/3)
  shr.u r24=r24,20
  nop.i 0;;
}
{.mmi
  nop.m 0
  // r28=3*exponent
  shladd r28=r24,1,r24
  // bias exponent
  add r24=r29,r24;;
}
{.mmi
  // get remainder of exponent/3
  sub r25=r25,r28;;
  // add sign to exponent
  or r24=r24,r26
  // remainder <<=8
  shl r25=r25,8;;
}  
{.mfi
  // adjust D_table pointer by 1st index
  shladd r27=r25,2,r27
  // P_1=C_1+C_2*r
  (p6) fma.s1 f7=f9,f6,f7
  // adjust T_table pointer by 1st index
  shladd r2=r25,3,r2
}
{.mfi
  // f14=sign*2^{exponent/3}
  (p6) setf.exp f14=r24
  // r2=r*r
  (p6) fma.s1 f9=f6,f6,f0
  nop.i 0;;
}
{.mfi
  // load D
  (p6) ldfs f15=[r27]
  // P_2=C_3+C_4*r
  (p6) fma.s1 f10=f11,f6,f10
  nop.i 0
}
{.mfi
  // load T
  (p6) ldf8 f8=[r2]
  // P_3=C_5+C_6*r
  (p6) fma.s1 f12=f13,f6,f12
  nop.i 0;;
}
{.mfi
  nop.m 0
  // P_4=D-r*P_1
  (p6) fnma.s1 f15=f6,f7,f15
  nop.i 0
}
{.mfi
  nop.m 0
  // r3=r*r2
  (p6) fma.s1 f6=f6,f9,f0
  nop.i 0;;
}
{.mfi
  nop.m 0
  // P_5=P_2+r2*P_3
  (p6) fma.s1 f10=f9,f12,f10
  nop.i 0;;
}
{.mfi
  nop.m 0
  // T=T*(sign*2^{exponent/3})
  (p6) fma.s1 f8=f8,f14,f0
  nop.i 0
}
{.mfi
  nop.m 0
  // P=P_4-r3*P_5
  (p6) fnma.s1 f6=f6,f10,f15
  nop.i 0;;
}
{.mfb
  nop.m 0
  // result=T+T*p
  (p6) fma.s0 f8=f8,f6,f8
  br.ret.sptk b0;;
}
.endp cbrtl
ASM_SIZE_DIRECTIVE(cbrtl)