base_table.h

浏览该文件的文档.

/*
 * [LAMMP]
 * Copyright (C) [2025-2026] [HJimmyK(Jericho Knox)]
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public License
 * along with this program.  If not, see <https://www.gnu.org/licenses/>.
 */
 
#ifndef __LAMMP_BASE_TABLE_H__
#define __LAMMP_BASE_TABLE_H__
 
#include "../lmmpn.h"
 
typedef struct mp_base_t {
    // 单个limb能容纳的基数的最大幂次
    // 二的幂次存储 log2(base)
    // large_base = base ^ digits_in_limb
    mp_limb_t large_base;
    // ceiling(2^64*log2(base)/log2(2^64))
    // N 位 base 数最多需要 N * lg_base / 2^64 + 1 个 limb
    mp_limb_t lg_base;
    // N 位二进制数最多需要 N * inv_lg_base / 2^64 + 1 个 base 进制位
    // ceiling(2^64/log2(base))
    mp_limb_t inv_lg_base;
    // 单个limb可容纳的最大基数位数
    // floor(64/log2(base))
    int digits_in_limb;
    // 基数（2~256）
    int base;
} mp_base_t;
 
typedef struct mp_basepow_t {
    // 基数幂值(base^digits)
    mp_ptr p;
    // p的 limb 长度
    mp_size_t np;
    // 归一化p的逆元
    mp_ptr invp;
    // invp的有效长度
    mp_size_t ni;
    // 去除的末尾零 limb 长度
    mp_size_t zeros;
    // 基数幂的指数（log_base(p)）
    mp_size_t digits;
    // p归一化时的移位位数
    int norm_cnt;
    // 基数
    int base;
} mp_basepow_t;
 
extern const mp_base_t lmmp_bases_table[255];
 
#endif // __LAMMP_BASE_TABLE_H__