//// Polynomial coefficients for the tgammal(x), 1.75 <= |x| < 2.0 //A5 = 2.6698206874501426502654943818e-04 + 3.4033756836921062797887300000e-20data8 0x3F317F3740FE2A68, 0x3BE417093234B06E//A4 = 7.4249010753513894345090307070e-02 + 3.9810018444482764697014200000e-18data8 0x3FB301FBB0F25A92, 0x3C525BEFFABB622F//A3 = -8.1576919247086265851720554565e-02+ -5.2716624487804746360745000000e-19data8 0xBFB4E239984650AC, 0xBC2372F1C4F276FFdata8 0xFEF3AEE71038E9A3, 0x00003FEB //A20 = 1.8995395865421509009969188571e-06data8 0xA11CFA2672BF876A, 0x0000BFEB //A19 =-1.2003868221414015771269244270e-06data8 0xF8E107215DAE2164, 0x00003FEC //A18 = 3.7085863210303833432006027217e-06data8 0xBCDDD3FC011EF7D6, 0x00003FEC //A17 = 2.8143303971756051015245433043e-06data8 0x8683C4687FA22E68, 0x00003FEE //A16 = 8.0177018464360416764308252462e-06data8 0xFDA09E5D33E32968, 0x00003FEE //A15 = 1.5117372062443781157389064848e-05data8 0xFFB00D0CFF4089B4, 0x00003FEF //A14 = 3.0480348961227424242198174995e-05data8 0xFEF6C39566785085, 0x00003FF0 //A13 = 6.0788135974125244644334004947e-05//A2 = 4.1184033042643969357854416558e-01 + 1.2103396182129232634761000000e-18data8 0x3FDA5B978B96BEBF, 0x3C3653AAD0A139E4//A1 = -4.2278433509846713445057275749e-01+ -4.9429151528135657430413000000e-18data8 0xBFDB0EE6072093CE, 0xBC56CB907027554F//A0 = 1.0000000000000000000000000000e+00 + 1.0969171200000000000000000000e-31data8 0x3FF0000000000000, 0x3981CC6A5B20B4D5data8 0xFF2B7BA9A8D68C37, 0x00003FF1 //A12 = 1.2167446884801403650547161615e-04data8 0xFCA53468E3692EF1, 0x00003FF2 //A11 = 2.4094136329542400976250900707e-04data8 0x808D698A9C993615, 0x00003FF4 //A10 = 4.9038845704938303659791698883e-04data8 0xF10F8E3FB8BB4AFB, 0x00003FF4 //A9 = 9.1957383840999861214472423976e-04data8 0x89E224E42F93F005, 0x00003FF6 //A8 = 2.1039333407187324139473634747e-03data8 0xBAF374824937A323, 0x00003FF6 //A7 = 2.8526458211545152218493600470e-03data8 0xB6BF7564F52140C6, 0x00003FF8 //A6 = 1.1154045718131014476684982178e-02//// Polynomial coefficients for the tgammal(x), 0.0 <= |x| < 0.125 //A5 = -9.8199506890314514073736518185e-01+ -5.9363811993837985890950900000e-17data8 0xBFEF6C80EC38B67A, 0xBC911C46B447C81F//A4 = 9.8172808683440015986576554496e-01 + 2.7457414262802803699834200000e-17data8 0x3FEF6A51055096B5, 0x3C7FA7FF90ACAD1F//A3 = -9.0747907608088618225394839101e-01 + -1.0676255850934306734701780000e-16data8 0xBFED0A118F324B62, 0xBC9EC5AFB633438Ddata8 0x9217E83FA207CB80, 0x00003FFD //A20 = 2.8533864762086088781083621561e-01data8 0xA8DABFA52FDF03EC, 0x0000BFFE //A19= -6.5958783896337186303285832783e-01data8 0xE331ED293AF39F9B, 0x00003FFE //A18 = 8.8748056656454687449654731184e-01data8 0xF9163C5DDB52419D, 0x0000BFFE //A17= -9.7299554149078295602977718525e-01data8 0xFEC0A1C672CB9265, 0x00003FFE //A16 = 9.9512683005268190987854104489e-01data8 0xFFD2D65B8EA7B5F4, 0x0000BFFE //A15= -9.9931087241443958201592847861e-01data8 0xFFF93AA39EE53445, 0x00003FFE //A14 = 9.9989668364186884793382816496e-01data8 0xFFFB99A9A3F5F480, 0x0000BFFE //A13= -9.9993286506283835663204999212e-01//A2 = 9.8905599532797250361682017683e-01 + 5.1778575360788420716540100000e-17data8 0x3FEFA658C23B1578, 0x3C8DD92B45408D07//A1 = -5.7721566490153275452712478000e-01+ -1.0607938730998824663273110000e-16data8 0xBFE2788CFC6FB618, 0xBC9E9346F8FDE55B//A0 = 9.9999999999999988897769753748e-01 + 1.1102230246251564036631420000e-16data8 0x3FEFFFFFFFFFFFFF, 0x3C9FFFFFFFFFFFFFdata8 0xFFF7FEBB545812C1, 0x00003FFE //A12 = 9.9987785409425126648628395084e-01data8 0xFFF00C02E943A3F2, 0x0000BFFE //A11= -9.9975657530855116454438747397e-01data8 0xFFE0420AADC53820, 0x00003FFE //A10 = 9.9951565514290485919027183699e-01data8 0xFFC01EB42EF27EEB, 0x0000BFFE //A9 = -9.9902526759155739377365522320e-01data8 0xFF83DAD0BF23FF12, 0x00003FFE //A8 = 9.9810569378236378800364235948e-01data8 0xFEF9F8ABDBCDB2F3, 0x0000BFFE //A7 = -9.9600176044241699109053158187e-01data8 0xFE3F05375988491D, 0x00003FFE //A6 = 9.9314911462127599008937257662e-01LOCAL_OBJECT_END(Constants_Tgammal_poly_splitted).align 64LOCAL_OBJECT_START(Constants_Tgammal_common)// Positive overflow valuedata8 0x3FE0000000000000 // 0.5data8 0x3FF8000000000000 // 1.5data8 0x3FD0000000000000 // 0.25data8 0x0000000000000000 // 0data8 0xDB718C066B352E21, 0x00004009 // Positive overflow valueLOCAL_OBJECT_END(Constants_Tgammal_common)//=======================================================// Lgamma registers// General Purpose RegistersGR_l_Log_Table        = r33GR_l_Log_Table1       = r34GR_l_BIAS             = r34GR_l_Index1           = r35GR_l_Index2           = r36GR_l_signif_Z         = r37 GR_l_X_0              = r38GR_l_X_1              = r39GR_l_X_2              = r40GR_l_Z_1              = r41GR_l_Z_2              = r42GR_l_N                = r43GR_l_Index3           = r44GR_l_Stirling_Table   = r45GR_l_N_Unbiased       = r46// Floating Point RegistersFR_l_logl_X            = f8FR_l_h_3              = f10FR_l_poly_hi          = f10FR_l_W                = f11FR_l_S                = f12FR_l_GS_hi            = f13FR_l_Y_lo             = f13FR_l_r_cor            = f14FR_l_G_1              = f15FR_l_G                = f15FR_l_H_1              = f32FR_l_H                = f32FR_l_h                = f33FR_l_h_1              = f33FR_l_N                = f33FR_l_G_2              = f34FR_l_H_2              = f35FR_l_h_2              = f36FR_l_G_3              = f37FR_l_log2_hi          = f38FR_l_GS_lo            = f39FR_l_H_3              = f40FR_l_float_N          = f41FR_l_Q_4              = f42FR_l_Q_3              = f43FR_l_Q_2              = f44FR_l_Q_1              = f45FR_l_Q_5              = f46FR_l_Q_6              = f47FR_l_log2_lo          = f48FR_l_r                = f49FR_l_poly_lo          = f50FR_l_poly             = f51FR_l_rsq              = f52FR_l_Y_lo_res         = f53FR_l_Y0               = f55  FR_l_Q0               = f56  FR_l_E0               = f57  FR_l_E2               = f58  FR_l_E1               = f59  FR_l_Y1               = f60  FR_l_E3               = f61  FR_l_Y2               = f62  FR_l_Z                = f63    FR_l_Z2               = f64    FR_l_Z4               = f65    FR_l_Z8               = f66                                   FR_l_CH               = f67    FR_l_CL               = f68    FR_l_B2               = f69  FR_l_B4               = f70  FR_l_B6               = f71  FR_l_B8               = f72  FR_l_B10              = f73  FR_l_B12              = f74  FR_l_B14              = f75  FR_l_B16              = f76  FR_l_B18              = f77  FR_l_Half             = f78  FR_l_SS               = f79FR_l_AbsX_m_Half      = f80FR_l_CXH              = f81FR_l_CXL              = f82FR_l_SSCXH            = f83FR_l_SSCXL            = f84FR_l_XYH              = f85FR_l_XYL              = f86FR_l_Temp             = f87FR_l_logl_YHi         = f88FR_l_logl_YLo         = f89FR_l_SignedXYH        = f123FR_l_AbsX             = f127//=======================================================// Negative part registers// General Purpose RegistersGR_n_sin_Table        = r47GR_n_XN               = r48// Float point registersFR_n_IXNS             = f125FR_n_IXN              = f126FR_n_XNS              = f90FR_n_XS               = f91FR_n_XS2              = f92FR_n_XS2L             = f93FR_n_XS4              = f94FR_n_XS7              = f95FR_n_XS8              = f96FR_n_TT               = f97FR_n_TH               = f98FR_n_TL               = f99FR_n_A2H              = f100     FR_n_A2L              = f101     FR_n_A1H              = f102     FR_n_A1L              = f103     FR_n_A9               = f104     FR_n_A8               = f105     FR_n_A7               = f106     FR_n_A6               = f107     FR_n_A5               = f108     FR_n_A4               = f109     FR_n_A3               = f110     FR_n_PolyH            = f111FR_n_PolyL            = f112FR_n_Poly1H           = f113FR_n_SinxH            = f113 // the same as FR_n_Poly1HFR_n_Poly1L           = f114FR_n_SinxL            = f114 // the same as FR_n_Poly1LFR_n_Tail             = f115FR_n_NegOne           = f116FR_n_Y0               = f117FR_n_Q0               = f118   FR_n_E0               = f119                                  FR_n_E2               = f120   FR_n_E1               = f121                                  FR_n_Y1               = f55    FR_n_E3               = f56                                   FR_n_Y2               = f57    FR_n_R0               = f58    FR_n_E4               = f59FR_n_RcpResH          = f60FR_n_Y3               = f61FR_n_R1               = f62FR_n_Temp             = f63FR_n_RcpResL          = f64FR_n_ResH             = f65FR_n_ResL             = f66//=======================================================// Exp registers// General Purpose RegistersGR_e_ad_Arg           = r33GR_e_ad_A             = r34GR_e_signexp_x        = r35GR_e_exp_x            = r35GR_e_exp_mask         = r36GR_e_ad_W1            = r37GR_e_ad_W2            = r38GR_e_M2               = r39GR_e_M1               = r40GR_e_K                = r41GR_e_exp_2_mk         = r42GR_e_exp_2_k          = r43GR_e_ad_T1            = r44GR_e_ad_T2            = r45GR_e_N_fix            = r46GR_e_one              = r47GR_e_exp_bias         = r48GR_e_sig_inv_ln2      = r49GR_e_rshf_2to51       = r50GR_e_exp_2tom51       = r51GR_e_rshf             = r52// Floating Point RegistersFR_e_RSHF_2TO51       = f10FR_e_INV_LN2_2TO63    = f11FR_e_W_2TO51_RSH      = f12FR_e_2TOM51           = f13FR_e_RSHF             = f14FR_e_Y_hi             = f15FR_e_Y_lo             = f32FR_e_scale            = f33FR_e_float_N          = f34FR_e_N_signif         = f35FR_e_L_hi             = f36FR_e_L_lo             = f37FR_e_r                = f38FR_e_W1               = f39FR_e_T1               = f40FR_e_W2               = f41FR_e_T2               = f42FR_e_W1_p1            = f43FR_e_rsq              = f44FR_e_A2               = f45FR_e_r4               = f46FR_e_A3               = f47FR_e_poly             = f48FR_e_T                = f49FR_e_W                = f50FR_e_Wp1              = f51FR_e_r6               = f52FR_e_2_mk             = f53FR_e_A1               = f54FR_e_T_scale          = f55FR_e_result_lo        = f56FR_e_W_T_scale        = f57FR_e_Wp1_T_scale      = f58FR_e_expl_Input_X     = f123FR_e_expl_Input_Y     = f124FR_e_expl_Output_X    = f123FR_e_expl_Output_Y    = f124FR_e_expl_Input_AbsX  = f122//=======================================================// Common registers// General Purpose RegistersGR_c_Table            = r53GR_c_NegUnderflow     = r54GR_c_NegSingularity   = r55GR_c_X                = r56GR_c_SignBit          = r57