log2_exp2.h 文件参考

#include <stdint.h>

log2_exp2.h 的引用(Include)关系图:

此图展示该文件直接或间接的被哪些文件引用了:

浏览源代码.

函数
void	exp2_fixed_128 (uint64_t *dst, uint64_t high, uint64_t low)
	floor(exp2(x/B)*B-B), B=2^128
uint64_t	exp2_fixed_64 (uint64_t x)
	floor(exp2(x/B)*B-B), B=2^64
void	log2_fixed_128 (uint64_t *dst, uint64_t high, uint64_t low)
	floor(log2(1+x/B)*B), B=2^128
uint64_t	log2_fixed_64 (uint64_t x)
	Copyright (C) 2026 HJimmyK(Jericho Knox)

函数说明

exp2_fixed_128()

void exp2_fixed_128	(	uint64_t *	dst,
		uint64_t	high,
		uint64_t	low
	)

floor(exp2(x/B)*B-B), B=2^128

参数

high	输入的小数部分高64位
low	输入的小数部分低64位
dst	输出数组（2个元素）：dst[0] 为低64位，dst[1] 为高64位

注解

x = high * 2^64 + low

在文件 log2_exp2.c 第 409 行定义.

                                                                {
    uint64_t res[3];
    uint64_t x[2] = {low, high};
    uint64_t coeff[3];
    res[0] = exp2_coeffs[EXP2_COEFFS_SIZE - 1][0];
    res[1] = exp2_coeffs[EXP2_COEFFS_SIZE - 1][1];
    res[2] = exp2_coeffs[EXP2_COEFFS_SIZE - 1][2];
    for (int i = 1; i < EXP2_COEFFS_SIZE; i++) {
        coeff[0] = exp2_coeffs[EXP2_COEFFS_SIZE - 1 - i][0];
        coeff[1] = exp2_coeffs[EXP2_COEFFS_SIZE - 1 - i][1];
        coeff[2] = exp2_coeffs[EXP2_COEFFS_SIZE - 1 - i][2];
        umul192x128_tohi192(res, res, x);
        _u192add(res, coeff);
    }
    dst[0] = res[1];
    dst[1] = res[2];
}

引用了 _u192add, coeff, exp2_coeffs, EXP2_COEFFS_SIZE, n , 以及 umul192x128_tohi192().

被这些函数引用 lmmp_cbrtapprox_6_().

函数调用图:

这是这个函数的调用关系图:

exp2_fixed_64()

uint64_t exp2_fixed_64

(

uint64_t

x

)

floor(exp2(x/B)*B-B), B=2^64

参数

x

输入的小数部分

返回

floor(exp2(x/B)*B-B)

在文件 log2_exp2.c 第 485 行定义.

                                   {
    uint64_t res[2];
    uint64_t coeff[2];
    res[0] = exp2_coeffs_64bit[EXP2_COEFFS_SIZE_64BIT - 1] [0];
    res[1] = exp2_coeffs_64bit[EXP2_COEFFS_SIZE_64BIT - 1] [1];
    for (int i = 1; i < EXP2_COEFFS_SIZE_64BIT; i++) {
        coeff[0] = exp2_coeffs_64bit[EXP2_COEFFS_SIZE_64BIT - 1 - i][0];
        coeff[1] = exp2_coeffs_64bit[EXP2_COEFFS_SIZE_64BIT - 1 - i][1];
        umul128x64_tohi128(res, res, x);
        _u128add(res, res, coeff);
    }
    return res[1];
}

引用了 _u128add, coeff, exp2_coeffs_64bit, EXP2_COEFFS_SIZE_64BIT, n , 以及 umul128x64_tohi128().

被这些函数引用 lmmp_cbrtapprox_3_().

函数调用图:

这是这个函数的调用关系图:

log2_fixed_128()

void log2_fixed_128	(	uint64_t *	dst,
		uint64_t	high,
		uint64_t	low
	)

floor(log2(1+x/B)*B), B=2^128

参数

high	输入的小数部分高64位
low	输入的小数部分低64位
dst	输出数组（2个元素）：dst[0] 为低64位，dst[1] 为高64位

注解

x = high * 2^64 + low

在文件 log2_exp2.c 第 292 行定义.

                                                                {
    uint64_t res[3];
    uint64_t x[3] = {0, low, high};
    uint64_t coeff[3];
    if (high < 0x4000000000000000ULL) {
        bool sign = log2_coeffs_bool_1[LOG2_COEFFS_SIZE_1 - 1];
        bool coeff_sign;
        res[0] = log2_coeffs_1[LOG2_COEFFS_SIZE_1 - 1][0];
        res[1] = log2_coeffs_1[LOG2_COEFFS_SIZE_1 - 1][1];
        res[2] = log2_coeffs_1[LOG2_COEFFS_SIZE_1 - 1][2];
        for (int i = 1; i < LOG2_COEFFS_SIZE_1; i++) {
            coeff[0] = log2_coeffs_1[LOG2_COEFFS_SIZE_1 - 1 - i][0];
            coeff[1] = log2_coeffs_1[LOG2_COEFFS_SIZE_1 - 1 - i][1];
            coeff[2] = log2_coeffs_1[LOG2_COEFFS_SIZE_1 - 1 - i][2];
            coeff_sign = log2_coeffs_bool_1[LOG2_COEFFS_SIZE_1 - 1 - i];
 
            umul192x128_tohi192(res, res, x + 1);
            if (sign == coeff_sign) {
                _u192add(res, coeff);
            } else {
                if (leq_192(res, coeff)) {
                    _u192sub(coeff, res);
                    res[0] = coeff[0];
                    res[1] = coeff[1];
                    res[2] = coeff[2];
                    sign = coeff_sign;
                } else {
                    _u192sub(res, coeff);
                }
            }
        }
    } else if (high < 0x8000000000000000ULL) {
        bool sign = log2_coeffs_bool_2[LOG2_COEFFS_SIZE_2 - 1];
        bool coeff_sign;
        res[0] = log2_coeffs_2[LOG2_COEFFS_SIZE_2 - 1][0];
        res[1] = log2_coeffs_2[LOG2_COEFFS_SIZE_2 - 1][1];
        res[2] = log2_coeffs_2[LOG2_COEFFS_SIZE_2 - 1][2];
        for (int i = 1; i < LOG2_COEFFS_SIZE_2; i++) {
            coeff[0] = log2_coeffs_2[LOG2_COEFFS_SIZE_2 - 1 - i][0];
            coeff[1] = log2_coeffs_2[LOG2_COEFFS_SIZE_2 - 1 - i][1];
            coeff[2] = log2_coeffs_2[LOG2_COEFFS_SIZE_2 - 1 - i][2];
            coeff_sign = log2_coeffs_bool_2[LOG2_COEFFS_SIZE_2 - 1 - i];
 
            umul192x128_tohi192(res, res, x + 1);
            if (sign == coeff_sign) {
                _u192add(res, coeff);
            } else {
                if (leq_192(res, coeff)) {
                    _u192sub(coeff, res);
                    res[0] = coeff[0];
                    res[1] = coeff[1];
                    res[2] = coeff[2];
                    sign = coeff_sign;
                } else {
                    _u192sub(res, coeff);
                }
            }
        }
    } else if (high < 0xC000000000000000ULL) {
        bool sign = log2_coeffs_bool_3[LOG2_COEFFS_SIZE_3 - 1];
        bool coeff_sign;
        res[0] = log2_coeffs_3[LOG2_COEFFS_SIZE_3 - 1][0];
        res[1] = log2_coeffs_3[LOG2_COEFFS_SIZE_3 - 1][1];
        res[2] = log2_coeffs_3[LOG2_COEFFS_SIZE_3 - 1][2];
        for (int i = 1; i < LOG2_COEFFS_SIZE_3; i++) {
            coeff[0] = log2_coeffs_3[LOG2_COEFFS_SIZE_3 - 1 - i][0];
            coeff[1] = log2_coeffs_3[LOG2_COEFFS_SIZE_3 - 1 - i][1];
            coeff[2] = log2_coeffs_3[LOG2_COEFFS_SIZE_3 - 1 - i][2];
            coeff_sign = log2_coeffs_bool_3[LOG2_COEFFS_SIZE_3 - 1 - i];
 
            umul192x128_tohi192(res, res, x + 1);
            if (sign == coeff_sign) {
                _u192add(res, coeff);
            } else {
                if (leq_192(res, coeff)) {
                    _u192sub(coeff, res);
                    res[0] = coeff[0];
                    res[1] = coeff[1];
                    res[2] = coeff[2];
                    sign = coeff_sign;
                } else {
                    _u192sub(res, coeff);
                }
            }
        }
    } else {
        bool sign = log2_coeffs_bool_4[LOG2_COEFFS_SIZE_4 - 1];
        bool coeff_sign;
        res[0] = log2_coeffs_4[LOG2_COEFFS_SIZE_4 - 1][0];
        res[1] = log2_coeffs_4[LOG2_COEFFS_SIZE_4 - 1][1];
        res[2] = log2_coeffs_4[LOG2_COEFFS_SIZE_4 - 1][2];
        for (int i = 1; i < LOG2_COEFFS_SIZE_4; i++) {
            coeff[0] = log2_coeffs_4[LOG2_COEFFS_SIZE_4 - 1 - i][0];
            coeff[1] = log2_coeffs_4[LOG2_COEFFS_SIZE_4 - 1 - i][1];
            coeff[2] = log2_coeffs_4[LOG2_COEFFS_SIZE_4 - 1 - i][2];
            coeff_sign = log2_coeffs_bool_4[LOG2_COEFFS_SIZE_4 - 1 - i];
 
            umul192x128_tohi192(res, res, x + 1);
            if (sign == coeff_sign) {
                _u192add(res, coeff);
            } else {
                if (leq_192(res, coeff)) {
                    _u192sub(coeff, res);
                    res[0] = coeff[0];
                    res[1] = coeff[1];
                    res[2] = coeff[2];
                    sign = coeff_sign;
                } else {
                    _u192sub(res, coeff);
                }
            }
        }
    }
    dst[0] = (res[1] << 1) | (res[0] >> 63);
    dst[1] = (res[2] << 1) | (res[1] >> 63);
}

引用了 _u192add, _u192sub, coeff, leq_192(), log2_coeffs_1, log2_coeffs_2, log2_coeffs_3, log2_coeffs_4, log2_coeffs_bool_1, log2_coeffs_bool_2, log2_coeffs_bool_3, log2_coeffs_bool_4, LOG2_COEFFS_SIZE_1, LOG2_COEFFS_SIZE_2, LOG2_COEFFS_SIZE_3, LOG2_COEFFS_SIZE_4, n , 以及 umul192x128_tohi192().

被这些函数引用 lmmp_cbrtapprox_6_().

函数调用图:

这是这个函数的调用关系图:

log2_fixed_64()

uint64_t log2_fixed_64

(

uint64_t

x

)

Copyright (C) 2026 HJimmyK(Jericho Knox)

This file is part of LAMMP.

LAMMP is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License (LGPL) as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version.

This program is distributed WITHOUT ANY WARRANTY.

See https://www.gnu.org/licenses/.

floor(log2(1+x/B)*B), B=2^64

参数

x

输入的小数部分

返回

floor(log2(1+x/B)*B)

在文件 log2_exp2.c 第 458 行定义.

                                   {
    uint64_t res[2];
    uint64_t coeff[2];
    bool sign = log2_coeffs_bool_64bit[LOG2_COEFFS_SIZE_64BIT - 1];
    bool coeff_sign;
    res[0] = log2_coeffs_64bit[LOG2_COEFFS_SIZE_64BIT - 1][0];
    res[1] = log2_coeffs_64bit[LOG2_COEFFS_SIZE_64BIT - 1][1];
    for (int i = 1; i < LOG2_COEFFS_SIZE_64BIT; i++) {
        coeff[0] = log2_coeffs_64bit[LOG2_COEFFS_SIZE_64BIT - 1 - i][0];
        coeff[1] = log2_coeffs_64bit[LOG2_COEFFS_SIZE_64BIT - 1 - i][1];
        coeff_sign = log2_coeffs_bool_64bit[LOG2_COEFFS_SIZE_64BIT - 1 - i];
        umul128x64_tohi128(res, res, x);
 
        if (sign == coeff_sign) {
            _u128add(res, res, coeff);
        } else {
            if (_u128cmp(res, coeff)) {
                _u128sub(res, coeff, res);
                sign = coeff_sign;
            } else {
                _u128sub(res, res, coeff);
            }
        }
    }
    return (res[1] << 1) | (res[0] >> 63);
}

引用了 _u128add, _u128cmp, _u128sub, coeff, log2_coeffs_64bit, log2_coeffs_bool_64bit, LOG2_COEFFS_SIZE_64BIT, n , 以及 umul128x64_tohi128().

被这些函数引用 lmmp_cbrtapprox_3_().

函数调用图:

这是这个函数的调用关系图: