src/lexical/cached_float80.rs - master - germanubis/json

src/lexical/cached_float80.rs
206 строк6 KB
David Tolnay
Add brief note to top of lexical sources
08 июн 2020, 01:20
08 июн 2020, 01:20ddf5f17
// Adapted from https://github.com/Alexhuszagh/rust-lexical.

//! Cached exponents for basen values with 80-bit extended floats.
//!
//! Exact versions of base**n as an extended-precision float, with both
//! large and small powers. Use the large powers to minimize the amount
//! of compounded error.
//!
//! These values were calculated using Python, using the arbitrary-precision
//! integer to calculate exact extended-representation of each value.
//! These values are all normalized.

use super::cached::{ExtendedFloatArray, ModeratePathPowers};

// LOW-LEVEL
// ---------

// BASE10

const BASE10_SMALL_MANTISSA: [u64; 10] = [
    9223372036854775808,  // 10^0
    11529215046068469760, // 10^1
    14411518807585587200, // 10^2
    18014398509481984000, // 10^3
    11258999068426240000, // 10^4
    14073748835532800000, // 10^5
    17592186044416000000, // 10^6
    10995116277760000000, // 10^7
    13743895347200000000, // 10^8
    17179869184000000000, // 10^9
];
const BASE10_SMALL_EXPONENT: [i32; 10] = [
    -63, // 10^0
    -60, // 10^1
    -57, // 10^2
    -54, // 10^3
    -50, // 10^4
    -47, // 10^5
    -44, // 10^6
    -40, // 10^7
    -37, // 10^8
    -34, // 10^9
];
const BASE10_LARGE_MANTISSA: [u64; 66] = [
    11555125961253852697, // 10^-350
    13451937075301367670, // 10^-340
    15660115838168849784, // 10^-330
    18230774251475056848, // 10^-320
    10611707258198326947, // 10^-310
    12353653155963782858, // 10^-300
    14381545078898527261, // 10^-290
    16742321987285426889, // 10^-280
    9745314011399999080,  // 10^-270
    11345038669416679861, // 10^-260
    13207363278391631158, // 10^-250
    15375394465392026070, // 10^-240
    17899314949046850752, // 10^-230
    10418772551374772303, // 10^-220
    12129047596099288555, // 10^-210
    14120069793541087484, // 10^-200
    16437924692338667210, // 10^-190
    9568131466127621947,  // 10^-180
    11138771039116687545, // 10^-170
    12967236152753102995, // 10^-160
    15095849699286165408, // 10^-150
    17573882009934360870, // 10^-140
    10229345649675443343, // 10^-130
    11908525658859223294, // 10^-120
    13863348470604074297, // 10^-110
    16139061738043178685, // 10^-100
    9394170331095332911,  // 10^-90
    10936253623915059621, // 10^-80
    12731474852090538039, // 10^-70
    14821387422376473014, // 10^-60
    17254365866976409468, // 10^-50
    10043362776618689222, // 10^-40
    11692013098647223345, // 10^-30
    13611294676837538538, // 10^-20
    15845632502852867518, // 10^-10
    9223372036854775808,  // 10^0
    10737418240000000000, // 10^10
    12500000000000000000, // 10^20
    14551915228366851806, // 10^30
    16940658945086006781, // 10^40
    9860761315262647567,  // 10^50
    11479437019748901445, // 10^60
    13363823550460978230, // 10^70
    15557538194652854267, // 10^80
    18111358157653424735, // 10^90
    10542197943230523224, // 10^100
    12272733663244316382, // 10^110
    14287342391028437277, // 10^120
    16632655625031838749, // 10^130
    9681479787123295682,  // 10^140
    11270725851789228247, // 10^150
    13120851772591970218, // 10^160
    15274681817498023410, // 10^170
    17782069995880619867, // 10^180
    10350527006597618960, // 10^190
    12049599325514420588, // 10^200
    14027579833653779454, // 10^210
    16330252207878254650, // 10^220
    9505457831475799117,  // 10^230
    11065809325636130661, // 10^240
    12882297539194266616, // 10^250
    14996968138956309548, // 10^260
    17458768723248864463, // 10^270
    10162340898095201970, // 10^280
    11830521861667747109, // 10^290
    13772540099066387756, // 10^300
];
const BASE10_LARGE_EXPONENT: [i32; 66] = [
    -1226, // 10^-350
    -1193, // 10^-340
    -1160, // 10^-330
    -1127, // 10^-320
    -1093, // 10^-310
    -1060, // 10^-300
    -1027, // 10^-290
    -994,  // 10^-280
    -960,  // 10^-270
    -927,  // 10^-260
    -894,  // 10^-250
    -861,  // 10^-240
    -828,  // 10^-230
    -794,  // 10^-220
    -761,  // 10^-210
    -728,  // 10^-200
    -695,  // 10^-190
    -661,  // 10^-180
    -628,  // 10^-170
    -595,  // 10^-160
    -562,  // 10^-150
    -529,  // 10^-140
    -495,  // 10^-130
    -462,  // 10^-120
    -429,  // 10^-110
    -396,  // 10^-100
    -362,  // 10^-90
    -329,  // 10^-80
    -296,  // 10^-70
    -263,  // 10^-60
    -230,  // 10^-50
    -196,  // 10^-40
    -163,  // 10^-30
    -130,  // 10^-20
    -97,   // 10^-10
    -63,   // 10^0
    -30,   // 10^10
    3,     // 10^20
    36,    // 10^30
    69,    // 10^40
    103,   // 10^50
    136,   // 10^60
    169,   // 10^70
    202,   // 10^80
    235,   // 10^90
    269,   // 10^100
    302,   // 10^110
    335,   // 10^120
    368,   // 10^130
    402,   // 10^140
    435,   // 10^150
    468,   // 10^160
    501,   // 10^170
    534,   // 10^180
    568,   // 10^190
    601,   // 10^200
    634,   // 10^210
    667,   // 10^220
    701,   // 10^230
    734,   // 10^240
    767,   // 10^250
    800,   // 10^260
    833,   // 10^270
    867,   // 10^280
    900,   // 10^290
    933,   // 10^300
];
const BASE10_SMALL_INT_POWERS: [u64; 10] = [
    1, 10, 100, 1000, 10000, 100000, 1000000, 10000000, 100000000, 1000000000,
];
const BASE10_STEP: i32 = 10;
const BASE10_BIAS: i32 = 350;

// HIGH LEVEL
// ----------

const BASE10_POWERS: ModeratePathPowers = ModeratePathPowers {
    small: ExtendedFloatArray {
        mant: &BASE10_SMALL_MANTISSA,
        exp: &BASE10_SMALL_EXPONENT,
    },
    large: ExtendedFloatArray {
        mant: &BASE10_LARGE_MANTISSA,
        exp: &BASE10_LARGE_EXPONENT,
    },
    small_int: &BASE10_SMALL_INT_POWERS,
    step: BASE10_STEP,
    bias: BASE10_BIAS,
};

/// Get powers from base.
pub(crate) fn get_powers() -> &'static ModeratePathPowers {
    &BASE10_POWERS
}