add_n.c 文件参考

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

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

浏览源代码.

函数
mp_limb_t	lmmp_add_n_ (mp_ptr dst, mp_srcptr numa, mp_srcptr numb, mp_size_t n)
	无进位的n位加法 [dst,n] = [numa,n] + [numb,n]
mp_limb_t	lmmp_add_nc_ (mp_ptr dst, mp_srcptr numa, mp_srcptr numb, mp_size_t n, mp_limb_t c)
	Copyright (C) 2026 HJimmyK(Jericho Knox)

函数说明

lmmp_add_n_()

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

无进位的n位加法 [dst,n] = [numa,n] + [numb,n]

参数

dst	结果输出指针
numa	第一个加数指针
numb	第二个加数指针
n	limb长度

警告

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

返回

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

在文件 add_n.c 第 81 行定义.

                                                                               {
    mp_size_t i = 0;
    mp_limb_t cy = 0;
 
    for (; i + 4 <= n; i += 4) {
        mp_limb_t a0, b0, r0;
        mp_limb_t a1, b1, r1;
        mp_limb_t a2, b2, r2;
        mp_limb_t a3, b3, r3;
 
        a0 = numa[i + 0];
        b0 = numb[i + 0];
 
        a1 = numa[i + 1];
        b1 = numb[i + 1];
 
        a2 = numa[i + 2];
        b2 = numb[i + 2];
 
        a3 = numa[i + 3];
        b3 = numb[i + 3];
 
        r0 = a0 + cy;
        cy = (r0 < cy);
        r0 += b0;
        cy += (r0 < b0);
 
        r1 = a1 + cy;
        cy = (r1 < cy);
        r1 += b1;
        cy += (r1 < b1);
 
        r2 = a2 + cy;
        cy = (r2 < cy);
        r2 += b2;
        cy += (r2 < b2);
 
        r3 = a3 + cy;
        cy = (r3 < cy);
        r3 += b3;
        cy += (r3 < b3);
 
        dst[i + 0] = r0;
        dst[i + 1] = r1;
        dst[i + 2] = r2;
        dst[i + 3] = r3;
    }
 
    for (; i < n; i++) {
        mp_limb_t a, b, r;
        a = numa[i];
        b = numb[i];
        r = a + cy;
        cy = (r < cy);
        r += b;
        cy += (r < b);
        dst[i] = r;
    }
 
    return cy;
}

引用了 a0, a1, a2, a3, b0, b1, b2, b3, n, numb, r0, r1, r2 , 以及 r3.

被这些函数引用 lmmp_add_(), lmmp_binvert_n_dc_(), lmmp_binvert_unbalanced_(), lmmp_div_(), lmmp_div_basecase_(), lmmp_div_divide_n_(), lmmp_div_mulinv_(), lmmp_div_s_(), lmmp_from_str_divide_(), lmmp_invappr_newton_(), lmmp_mul_(), lmmp_mul_fft_unbalance_(), lmmp_mul_toom22_(), lmmp_mul_toom32_(), lmmp_mul_toom33_(), lmmp_mul_toom42_(), lmmp_mul_toom42_cache_(), lmmp_mul_toom42_cache_init_(), lmmp_mul_toom42_unbalance_(), lmmp_mul_toom43_(), lmmp_mul_toom62_(), lmmp_mul_toom62_cache_(), lmmp_mul_toom62_cache_init_(), lmmp_mul_toom62_unbalance_(), lmmp_mullo_dc_(), lmmp_sqr_toom2_(), lmmp_sqrlo_dc_(), lmmp_sqrt_divide_(), lmmp_toom_eval_dgr3_pm1_(), lmmp_toom_eval_dgr3_pm2_(), lmmp_toom_eval_pm1_(), lmmp_toom_interp5_(), lmmp_toom_interp6_() , 以及 lmmp_toom_interp7_().

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

lmmp_add_nc_()

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

Copyright (C) 2026 HJimmyK(Jericho Knox)

带进位的n位加法 [dst,n] = [numa,n] + [numb,n] + c

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

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

                                                                                             {
    mp_size_t i = 0;
    mp_limb_t cy = c;
 
    for (; i + 4 <= n; i += 4) {
        mp_limb_t a0, b0, r0;
        mp_limb_t a1, b1, r1;
        mp_limb_t a2, b2, r2;
        mp_limb_t a3, b3, r3;
 
        a0 = numa[i + 0];
        b0 = numb[i + 0];
 
        a1 = numa[i + 1];
        b1 = numb[i + 1];
 
        a2 = numa[i + 2];
        b2 = numb[i + 2];
 
        a3 = numa[i + 3];
        b3 = numb[i + 3];
 
        r0 = a0 + cy;
        cy = (r0 < cy);
        r0 += b0;
        cy += (r0 < b0);
 
        r1 = a1 + cy;
        cy = (r1 < cy);
        r1 += b1;
        cy += (r1 < b1);
 
        r2 = a2 + cy;
        cy = (r2 < cy);
        r2 += b2;
        cy += (r2 < b2);
 
        r3 = a3 + cy;
        cy = (r3 < cy);
        r3 += b3;
        cy += (r3 < b3);
 
        dst[i + 0] = r0;
        dst[i + 1] = r1;
        dst[i + 2] = r2;
        dst[i + 3] = r3;
    }
 
    for (; i < n; i++) {
        mp_limb_t a, b, r;
        a = numa[i];
        b = numb[i];
        r = a + cy;
        cy = (r < cy);
        r += b;
        cy += (r < b);
        dst[i] = r;
    }
 
    return cy;
}

引用了 a0, a1, a2, a3, b0, b1, b2, b3, c, n, numb, r0, r1, r2 , 以及 r3.

被这些函数引用 lmmp_add_n_sub_n_(), lmmp_fft_bfy_(), lmmp_ifft_bfy_() , 以及 lmmp_invappr_newton_().

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