mul_toom_interp5.c 文件参考

#include "../../include/lammp/impl/divexact.h"
#include "../../include/lammp/impl/toom_interp.h"

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

浏览源代码.

宏定义
#define	r0   (dst)
#define	r1   (dst + n)
#define	r2   (dst + 2 * n)
#define	r3   (dst + 3 * n)
#define	r4   (dst + 4 * n)
#define	v0   r0
#define	v1   r2
#define	vinf   r4

函数
void	lmmp_toom_interp5_ (mp_ptr dst, mp_ptr v2, mp_ptr vm1, mp_size_t n, mp_size_t spt, int vm1_neg, mp_limb_t vinf0)
	Toom插值计算（5点插值），用于Toom-33和Toom-42乘法算法

宏定义说明

r0

#define r0   (dst)

r1

#define r1   (dst + n)

r2

#define r2   (dst + 2 * n)

r3

#define r3   (dst + 3 * n)

r4

#define r4   (dst + 4 * n)

v0

#define v0   r0

v1

#define v1   r2

vinf

#define vinf   r4

函数说明

lmmp_toom_interp5_()

void lmmp_toom_interp5_	(	mp_ptr	dst,
		mp_ptr	v2,
		mp_ptr	vm1,
		mp_size_t	n,
		mp_size_t	spt,
		int	vm1_neg,
		mp_limb_t	vinf0
	)

Toom插值计算（5点插值），用于Toom-33和Toom-42乘法算法

参数

dst	输出结果缓冲区，存储插值计算结果 v(0)储存在[dst,2n]，v(1)储存在[dst+2n,2n]
v2	v(2)插值点值，长度为 2n+1
vm1	v(-1)插值点值，长度为 2n+1
n	操作数的 limb 长度
spt	系数r4的长度
vm1_neg	符号标记：v(-1)是否为负数（1表示负，0表示正）
vinf0	无穷远点插值的低64位值

在文件 mul_toom_interp5.c 第 10 行定义.

                                                                                                                     {
    mp_limb_t cy, saved;
    mp_size_t dnp = 2 * n + 1;
 
#define r0 (dst)
#define r1 (dst + n)
#define r2 (dst + 2 * n)
#define r3 (dst + 3 * n)
#define r4 (dst + 4 * n)
#define v0 r0
#define v1 r2
#define vinf r4
 
    // v2=(v2-vm1)/3
    if (vm1_neg)
        lmmp_add_n_(v2, v2, vm1, dnp);
    else
        lmmp_sub_n_(v2, v2, vm1, dnp);
    lmmp_divexact_by3_(v2, v2, dnp);
 
    // vm1=(v1-vm1)/2
    if (vm1_neg)
        lmmp_shr1add_n_(vm1, v1, vm1, dnp);
    else
        lmmp_shr1sub_n_(vm1, v1, vm1, dnp);
 
    // v1=v1-v0
    v1[2 * n] -= lmmp_sub_n_(v1, v1, v0, 2 * n);
 
    // v2=(v2-v1)/2
    lmmp_shr1sub_n_(v2, v2, v1, dnp);
 
    // v1=v1-vm1
    lmmp_sub_n_(v1, v1, vm1, dnp);
 
    // add vm1 at correct place.
    cy = lmmp_add_n_(r1, r1, vm1, dnp);
    lmmp_inc_1(r3 + 1, cy);  // at most propagate to v1[2*n]
 
    saved = v1[2 * n];  // it is vinf[0]
    vinf[0] = vinf0;    // correct value of vinf
 
    // v2=v2-vinf*2
    cy = lmmp_shl_(vm1, vinf, spt, 1);
    cy += lmmp_sub_n_(v2, v2, vm1, spt);
    lmmp_dec_1(v2 + spt, cy);
 
    // vinf+=v2h, no overflow
    cy = lmmp_add_n_(vinf, vinf, v2 + n, n + 1);
    lmmp_inc_1(r3 + dnp, cy);
 
    // v1-=vinf, (same time vmh-=v2h)
    cy = lmmp_sub_n_(v1, v1, vinf, spt);
    vinf0 = vinf[0];
    v1[2 * n] = saved;  // correct value of v1
    lmmp_dec_1(v1 + spt, cy);
 
    // vml-=v2l
    cy = lmmp_sub_n_(r1, r1, v2, n);
    lmmp_dec_1(v1, cy);
 
    // last v2l
    cy = lmmp_add_n_(r3, r3, v2, n);
    v1[2 * n] += cy;  // no carry
    lmmp_inc_1(vinf, vinf0);
}

引用了 lmmp_add_n_(), lmmp_dec_1, lmmp_divexact_by3_(), lmmp_inc_1, lmmp_shl_(), lmmp_shr1add_n_(), lmmp_shr1sub_n_(), lmmp_sub_n_(), r1, r3, v0, v1, v2, vinf , 以及 vm1.

被这些函数引用 lmmp_mul_toom33_(), lmmp_mul_toom42_(), lmmp_mul_toom42_cache_(), lmmp_mul_toom42_cache_init_() , 以及 lmmp_sqr_toom3_().

函数调用图:

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