shl.c 文件参考

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

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

浏览源代码.

函数
mp_limb_t	lmmp_addshl1_n_ (mp_ptr dst, mp_srcptr numa, mp_srcptr numb, mp_size_t n)
	加法结合左移1位操作 [dst,n] = [numa,n] + ([numb,n] << 1)
mp_limb_t	lmmp_shl_ (mp_ptr dst, mp_srcptr numa, mp_size_t na, mp_size_t shr)
	Copyright (C) 2026 HJimmyK(Jericho Knox)
mp_limb_t	lmmp_shl_c_ (mp_ptr dst, mp_srcptr numa, mp_size_t na, mp_size_t shr, mp_limb_t c)
	带进位的大数左移操作 [dst,na] = [numa,na]<<shl，dst的低shl位填充c的低shl位
mp_limb_t	lmmp_subshl1_n_ (mp_ptr dst, mp_srcptr numa, mp_srcptr numb, mp_size_t n)
	减法结合左移1位操作 [dst,n] = [numa,n] - ([numb,n] << 1)

函数说明

lmmp_addshl1_n_()

mp_limb_t lmmp_addshl1_n_	(	mp_ptr	dst,
		mp_srcptr	numa,
		mp_srcptr	numb,
		mp_size_t	n
	)

加法结合左移1位操作 [dst,n] = [numa,n] + ([numb,n] << 1)

参数

dst	结果输出指针
numa	被加数指针
numb	加数指针（先左移1位）
n	limb长度

警告

n>0, eqsep(dst,[numa|numb])

返回

运算后的进位值 [0|1|2]

在文件 shl.c 第 66 行定义.

                                                                                   {
    mp_size_t i, c = 0, mb = 0;
 
    for (i = 0; i < n; i++) {
        mp_limb_t a, b, r;
        a = numa[i];
        b = (numb[i] << 1) + mb;
        mb = numb[i] >> (LIMB_BITS - 1);
        r = a + c;
        c = (r < c);
        r += b;
        c += (r < b);
        dst[i] = r;
    }
    return c + mb;
}

引用了 c, LIMB_BITS, n , 以及 numb.

被这些函数引用 lmmp_mul_toom33_(), lmmp_mul_toom42_(), lmmp_mul_toom42_cache_(), lmmp_mul_toom42_cache_init_(), lmmp_mul_toom44_(), lmmp_mul_toom53_(), lmmp_mul_toom62_(), lmmp_mul_toom62_cache_(), lmmp_mul_toom62_cache_init_(), lmmp_sqr_toom3_(), lmmp_sqr_toom4_(), lmmp_sqrlo_dc_() , 以及 lmmp_sqrt_divide_().

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

lmmp_shl_()

mp_limb_t lmmp_shl_	(	mp_ptr	dst,
		mp_srcptr	numa,
		mp_size_t	na,
		mp_size_t	shr
	)

Copyright (C) 2026 HJimmyK(Jericho Knox)

大数左移操作 [dst,na] = [numa,na]<<shl，dst的低shl位填充0

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

在文件 shl.c 第 19 行定义.

                                                                             {
    if (shr == 0) {
        lmmp_copy(dst, numa, na);
        return 0;
    } else {
        const mp_limb_t rshr = LIMB_BITS - shr;
        mp_limb_t high_limb, low_limb;
        mp_limb_t retval;
        numa += na;
        dst += na;
        low_limb = *--numa;
        retval = low_limb >> rshr;
        high_limb = (low_limb << shr);
        while (--na != 0) {
            low_limb = *--numa;
            *--dst = high_limb | (low_limb >> rshr);
            high_limb = (low_limb << shr);
        }
        *--dst = high_limb;
        return retval;
    }
}

引用了 LIMB_BITS, lmmp_copy , 以及 n.

被这些函数引用 lmmp_10pow_1_(), lmmp_12pow_1_(), lmmp_14pow_1_(), lmmp_2factorial_(), lmmp_6pow_1_(), lmmp_arith_seqprod_(), lmmp_bninv_(), lmmp_bninv_appr_newton_(), lmmp_div_(), lmmp_factorial_(), lmmp_fft_shl_coef_(), lmmp_hyperfac_(), lmmp_inv_(), lmmp_mul_fermat_recombine_(), lmmp_mul_mersenne_(), lmmp_mul_mersenne_single_(), lmmp_mul_toom43_(), lmmp_mul_toom44_(), lmmp_mul_toom53_(), lmmp_mul_toom62_(), lmmp_mul_toom62_cache_(), lmmp_mul_toom62_cache_init_(), lmmp_multinomial_(), lmmp_nCr_(), lmmp_nPr_(), lmmp_pow_1_(), lmmp_sqr_toom4_(), lmmp_sqrt_(), lmmp_sqrt_newton_(), lmmp_superfac_(), lmmp_to_str_(), lmmp_to_str_divide_(), lmmp_toom_eval_dgr3_pm2_(), lmmp_toom_eval_pm2_(), lmmp_toom_interp5_(), lmmp_toom_interp6_(), lmmp_toom_interp7_() , 以及 pow_nPr_().

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

lmmp_shl_c_()

mp_limb_t lmmp_shl_c_	(	mp_ptr	dst,
		mp_srcptr	numa,
		mp_size_t	na,
		mp_size_t	shl,
		mp_limb_t	c
	)

带进位的大数左移操作 [dst,na] = [numa,na]<<shl，dst的低shl位填充c的低shl位

参数

dst	结果输出指针
numa	源操作数指针
na	limb长度
shl	左移的位数 (0~63)
c	进位值（其高(64-shl)位必须为0）

警告

na>0, 0<=shl<64, eqsep(dst,numa) c的高(64-shl)位必须为0 允许dst指针地址大于numa（即支持原地长移位操作）

返回

其最低shl个比特位填充[numa,na]被移出的shl个最高位，其余比特位为0

在文件 shl.c 第 42 行定义.

                                                                                            {
    if (shr == 0) {
        lmmp_copy(dst, numa, na);
        return 0;
    } else {
        const mp_limb_t rshr = LIMB_BITS - shr;
        mp_limb_t high_limb, low_limb;
        mp_limb_t retval;
        numa += na;
        dst += na;
        low_limb = *--numa;
        retval = low_limb >> rshr;
        high_limb = (low_limb << shr);
        while (--na != 0) {
            low_limb = *--numa;
            *--dst = high_limb | (low_limb >> rshr);
            high_limb = (low_limb << shr);
        }
        c &= ((mp_limb_t)1 << shr) - 1;
        *--dst = high_limb | c;
        return retval;
    }
}

引用了 c, LIMB_BITS, lmmp_copy , 以及 n.

被这些函数引用 lmmp_fft_bfy_().

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

lmmp_subshl1_n_()

mp_limb_t lmmp_subshl1_n_	(	mp_ptr	dst,
		mp_srcptr	numa,
		mp_srcptr	numb,
		mp_size_t	n
	)

减法结合左移1位操作 [dst,n] = [numa,n] - ([numb,n] << 1)

参数

dst	结果输出指针
numa	被减数指针
numb	减数指针（先左移1位）
n	limb长度

警告

n>0, eqsep(dst,[numa|numb])

返回

运算后的借位值 [0|1|2]

在文件 shl.c 第 83 行定义.

                                                                                   {
    mp_size_t i, c = 0, mb = 0;
 
    for (i = 0; i < n; i++) {
        mp_limb_t a, b;
        a = numa[i];
        b = (numb[i] << 1) + mb;
        mb = numb[i] >> (LIMB_BITS - 1);
        b += c;
        c = (b < c);
        c += (a < b);
        dst[i] = a - b;
    }
    return c + mb;
}

引用了 c, LIMB_BITS, n , 以及 numb.