extract_bits.c 文件参考

#include "../../../include/lammp/lmmpn.h"
#include "../../../include/lammp/impl/inlines.h"

extract_bits.c 的引用(Include)关系图:

浏览源代码.

函数
mp_bitcnt_t	lmmp_extract_bits_ (mp_srcptr restrict num, mp_size_t n, mp_limb_t *restrict ext, int bits)
	Copyright (C) 2026 HJimmyK(Jericho Knox)

函数说明

lmmp_extract_bits_()

mp_bitcnt_t lmmp_extract_bits_	(	mp_srcptr restrict	num,
		mp_size_t	n,
		mp_limb_t *restrict	ext,
		int	bits
	)

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/.

在文件 extract_bits.c 第 20 行定义.

                                                                                                       {
    lmmp_param_assert(bits <= LIMB_BITS && bits > 0);
    lmmp_param_assert(n > 0);
    if (n == 1) {
        int lb = lmmp_limb_bits_(num[0]);
        if (lb <= bits) {
            *ext = num[0];
            return 0;
        } else {
            *ext = num[0] >> (lb - bits);
            return lb - bits;
        }
    } else {
        int lb = lmmp_limb_bits_(num[n - 1]);
        if (lb < bits) {
            *ext = num[n - 1] << (bits - lb);
            *ext |= num[n - 2] >> (LIMB_BITS - bits + lb);
            return LIMB_BITS * (n - 1) - (bits - lb);
        } else if (lb == bits) {
            *ext = num[n - 1];
            return LIMB_BITS * (n - 1);
        } else {
            *ext = num[n - 1] >> (lb - bits);
            return LIMB_BITS * (n - 1) + (lb - bits);
        }
    }
}

引用了 LIMB_BITS, lmmp_limb_bits_, lmmp_param_assert , 以及 n.