shr.c 文件参考

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

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

浏览源代码.

函数
mp_limb_t	lmmp_shr1add_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_shr1add_nc_ (mp_ptr dst, mp_srcptr numa, mp_srcptr numb, mp_size_t n, mp_limb_t c)
	带进位加法后右移1位 [dst,n] = ([numa,n] + [numb,n] + c) >> 1
mp_limb_t	lmmp_shr1sub_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_shr1sub_nc_ (mp_ptr dst, mp_srcptr numa, mp_srcptr numb, mp_size_t n, mp_limb_t c)
	带借位减法后右移1位 [dst,n] = ([numa,n] - [numb,n] - c) >> 1
mp_limb_t	lmmp_shr_ (mp_ptr dst, mp_srcptr numa, mp_size_t na, mp_size_t shr)
	Copyright (C) 2026 HJimmyK(Jericho Knox)
mp_limb_t	lmmp_shr_c_ (mp_ptr dst, mp_srcptr numa, mp_size_t na, mp_size_t shr, mp_limb_t c)
	带进位的大数右移操作 [dst,na] = [numa,na]>>shr，dst的高shr位填充c的高shr位

函数说明

lmmp_shr1add_n_()

mp_limb_t lmmp_shr1add_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	第二个加数指针
n	limb长度

警告

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

返回

右移操作产生的进位值 [0|1]

在文件 shr.c 第 62 行定义.

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

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

被这些函数引用 lmmp_mul_toom32_(), lmmp_toom_interp5_(), lmmp_toom_interp6_() , 以及 lmmp_toom_interp7_().

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

lmmp_shr1add_nc_()

mp_limb_t lmmp_shr1add_nc_	(	mp_ptr	dst,
		mp_srcptr	numa,
		mp_srcptr	numb,
		mp_size_t	n,
		mp_limb_t	c
	)

带进位加法后右移1位 [dst,n] = ([numa,n] + [numb,n] + c) >> 1

参数

dst	结果输出指针
numa	第一个加数指针
numb	第二个加数指针
n	limb长度
c	初始进位值 [0|1]

警告

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

返回

右移操作产生的进位值 [0|1]

在文件 shr.c 第 89 行定义.

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

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

被这些函数引用 lmmp_mul_fft_(), lmmp_mul_fft_cache_() , 以及 lmmp_mullo_fft_().

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

lmmp_shr1sub_n_()

mp_limb_t lmmp_shr1sub_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	减数指针
n	操作数的位数（limb数量）

警告

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

返回

右移操作产生的进位值 (0或1)

在文件 shr.c 第 116 行定义.

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

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

被这些函数引用 lmmp_mul_toom32_(), lmmp_toom_interp5_(), lmmp_toom_interp6_() , 以及 lmmp_toom_interp7_().

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

lmmp_shr1sub_nc_()

mp_limb_t lmmp_shr1sub_nc_	(	mp_ptr	dst,
		mp_srcptr	numa,
		mp_srcptr	numb,
		mp_size_t	n,
		mp_limb_t	c
	)

带借位减法后右移1位 [dst,n] = ([numa,n] - [numb,n] - c) >> 1

参数

dst	结果输出指针
numa	被减数指针
numb	减数指针
n	limb长度
c	初始借位值 [0|1]

警告

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

返回

右移操作产生的进位值 [0|1]

在文件 shr.c 第 143 行定义.

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

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

lmmp_shr_()

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

Copyright (C) 2026 HJimmyK(Jericho Knox)

大数右移操作 [dst,na] = [numa,na]>>shr，dst的高shr位填充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/.

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

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

引用了 LIMB_BITS, lmmp_copy , 以及 n.

被这些函数引用 lmmp_bninv_(), lmmp_div_(), lmmp_fft_extract_coef_(), lmmp_invsqrt_newton_(), lmmp_sqrt_(), lmmp_sqrt_newton_(), lmmp_to_str_divide_(), lmmp_toom_interp6_() , 以及 lmmp_toom_interp7_().

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

lmmp_shr_c_()

mp_limb_t lmmp_shr_c_	(	mp_ptr	dst,
		mp_srcptr	numa,
		mp_size_t	na,
		mp_size_t	shr,
		mp_limb_t	c
	)

带进位的大数右移操作 [dst,na] = [numa,na]>>shr，dst的高shr位填充c的高shr位

参数

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

警告

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

返回

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

在文件 shr.c 第 40 行定义.

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

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

被这些函数引用 lmmp_ifft_bfy_(), lmmp_inv_() , 以及 lmmp_sqrt_divide_().

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