mul_1.c 文件参考

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

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

浏览源代码.

函数
mp_limb_t	lmmp_addmul_1_ (mp_ptr restrict numa, mp_srcptr restrict numb, mp_size_t n, mp_limb_t b)
mp_limb_t	lmmp_mul_1_ (mp_ptr restrict dst, mp_srcptr restrict numa, mp_size_t na, mp_limb_t x)
	Copyright (C) 2026 HJimmyK(Jericho Knox)
void	lmmp_mullo_basecase_ (mp_ptr restrict dst, mp_srcptr numa, mp_srcptr numb, mp_size_t n)
mp_limb_t	lmmp_submul_1_ (mp_ptr restrict numa, mp_srcptr restrict numb, mp_size_t n, mp_limb_t b)

函数说明

lmmp_addmul_1_()

mp_limb_t lmmp_addmul_1_	(	mp_ptr restrict	numa,
		mp_srcptr restrict	numb,
		mp_size_t	n,
		mp_limb_t	b
	)

在文件 mul_1.c 第 106 行定义.

                                                                                                  {
    mp_limb_t cl = 0;
    mp_size_t i = 0;
 
    mp_limb_t ul0, ul1, ul2, ul3;
    mp_limb_t lpl0, hpl0, rl0;
    mp_limb_t lpl1, hpl1, rl1;
    mp_limb_t lpl2, hpl2, rl2;
    mp_limb_t lpl3, hpl3, rl3;
 
    if (numa == numb) {
        for (; i + 4 <= n; i += 4) {
            ul0 = numa[i + 0];
            ul1 = numa[i + 1];
            ul2 = numa[i + 2];
            ul3 = numa[i + 3];
 
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = ul0 + lpl0;
            cl += (lpl0 < ul0);
 
            _umul64to128_(ul1, b, &lpl1, &hpl1);
            lpl1 += cl;
            cl = (lpl1 < cl) + hpl1;
            lpl1 = ul1 + lpl1;
            cl += (lpl1 < ul1);
 
            _umul64to128_(ul2, b, &lpl2, &hpl2);
            lpl2 += cl;
            cl = (lpl2 < cl) + hpl2;
            lpl2 = ul2 + lpl2;
            cl += (lpl2 < ul2);
 
            _umul64to128_(ul3, b, &lpl3, &hpl3);
            lpl3 += cl;
            cl = (lpl3 < cl) + hpl3;
            lpl3 = ul3 + lpl3;
            cl += (lpl3 < ul3);
 
            numa[i + 0] = lpl0;
            numa[i + 1] = lpl1;
            numa[i + 2] = lpl2;
            numa[i + 3] = lpl3;
        }
        for (; i < n; i++) {
            ul0 = numa[i];
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = ul0 + lpl0;
            cl += (lpl0 < ul0);
            numa[i] = lpl0;
        }
    } else {
        /* seq(numa,numb) */
        for (; i + 4 <= n; i += 4) {
            ul0 = numb[i + 0];
            ul1 = numb[i + 1];
            ul2 = numb[i + 2];
            ul3 = numb[i + 3];
            rl0 = numa[i + 0];
            rl1 = numa[i + 1];
            rl2 = numa[i + 2];
            rl3 = numa[i + 3];
 
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = rl0 + lpl0;
            cl += (lpl0 < rl0);
 
            _umul64to128_(ul1, b, &lpl1, &hpl1);
            lpl1 += cl;
            cl = (lpl1 < cl) + hpl1;
            lpl1 = rl1 + lpl1;
            cl += (lpl1 < rl1);
 
            _umul64to128_(ul2, b, &lpl2, &hpl2);
            lpl2 += cl;
            cl = (lpl2 < cl) + hpl2;
            lpl2 = rl2 + lpl2;
            cl += (lpl2 < rl2);
 
            _umul64to128_(ul3, b, &lpl3, &hpl3);
            lpl3 += cl;
            cl = (lpl3 < cl) + hpl3;
            lpl3 = rl3 + lpl3;
            cl += (lpl3 < rl3);
 
            numa[i + 0] = lpl0;
            numa[i + 1] = lpl1;
            numa[i + 2] = lpl2;
            numa[i + 3] = lpl3;
        }
        for (; i < n; i++) {
            ul0 = numb[i];
            rl0 = numa[i];
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = rl0 + lpl0;
            cl += (lpl0 < rl0);
            numa[i] = lpl0;
        }
    }
    return cl;
}

引用了 _umul64to128_(), n , 以及 numb.

被这些函数引用 lmmp_mullo_basecase_().

函数调用图:

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

lmmp_mul_1_()

mp_limb_t lmmp_mul_1_	(	mp_ptr restrict	dst,
		mp_srcptr restrict	numa,
		mp_size_t	na,
		mp_limb_t	x
	)

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

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

                                                                                               {
    mp_limb_t cl = 0;
    mp_size_t i = 0;
 
    mp_limb_t ul0, ul1, ul2, ul3;
    mp_limb_t lpl0, hpl0, lpl1, hpl1;
    mp_limb_t lpl2, hpl2, lpl3, hpl3;
 
    if (dst == numa) {
        for (; i + 4 <= na; i += 4) {
            ul0 = dst[i + 0];
            ul1 = dst[i + 1];
            ul2 = dst[i + 2];
            ul3 = dst[i + 3];
 
            _umul64to128_(ul0, x, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
 
            _umul64to128_(ul1, x, &lpl1, &hpl1);
            lpl1 += cl;
            cl = (lpl1 < cl) + hpl1;
 
            _umul64to128_(ul2, x, &lpl2, &hpl2);
            lpl2 += cl;
            cl = (lpl2 < cl) + hpl2;
 
            _umul64to128_(ul3, x, &lpl3, &hpl3);
            lpl3 += cl;
            cl = (lpl3 < cl) + hpl3;
 
            dst[i + 0] = lpl0;
            dst[i + 1] = lpl1;
            dst[i + 2] = lpl2;
            dst[i + 3] = lpl3;
        }
 
        for (; i < na; i++) {
            ul0 = dst[i];
            _umul64to128_(ul0, x, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            dst[i] = lpl0;
        }
    } else {
        /* seq(dst,numa) */
        for (; i + 4 <= na; i += 4) {
            ul0 = numa[i + 0];
            ul1 = numa[i + 1];
            ul2 = numa[i + 2];
            ul3 = numa[i + 3];
 
            _umul64to128_(ul0, x, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
 
            _umul64to128_(ul1, x, &lpl1, &hpl1);
            lpl1 += cl;
            cl = (lpl1 < cl) + hpl1;
 
            _umul64to128_(ul2, x, &lpl2, &hpl2);
            lpl2 += cl;
            cl = (lpl2 < cl) + hpl2;
 
            _umul64to128_(ul3, x, &lpl3, &hpl3);
            lpl3 += cl;
            cl = (lpl3 < cl) + hpl3;
 
            dst[i + 0] = lpl0;
            dst[i + 1] = lpl1;
            dst[i + 2] = lpl2;
            dst[i + 3] = lpl3;
        }
 
        for (; i < na; i++) {
            ul0 = numa[i];
            _umul64to128_(ul0, x, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            dst[i] = lpl0;
        }
    }
 
    return cl;
}

引用了 _umul64to128_() , 以及 n.

被这些函数引用 lmmp_mullo_basecase_().

函数调用图:

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

lmmp_mullo_basecase_()

void lmmp_mullo_basecase_	(	mp_ptr restrict	dst,
		mp_srcptr	numa,
		mp_srcptr	numb,
		mp_size_t	n
	)

在文件 mul_1.c 第 326 行定义.

                                                                                            {
    mp_limb_t h;
    h = numa[0] * numb[n - 1];
    if (n != 1) {
        mp_size_t i;
        mp_limb_t b0;
 
        b0 = *numb++;
        h += numa[n - 1] * b0 + lmmp_mul_1_(dst, numa, n - 1, b0);
        dst++;
 
        for (i = n - 2; i > 0; i--) {
            b0 = *numb++;
            h += numa[i] * b0 + lmmp_addmul_1_(dst, numa, i, b0);
            dst++;
        }
    }
 
    dst[0] = h;
}

引用了 b0, lmmp_addmul_1_(), lmmp_mul_1_(), n , 以及 numb.

函数调用图:

lmmp_submul_1_()

mp_limb_t lmmp_submul_1_	(	mp_ptr restrict	numa,
		mp_srcptr restrict	numb,
		mp_size_t	n,
		mp_limb_t	b
	)

在文件 mul_1.c 第 216 行定义.

                                                                                                  {
    mp_limb_t cl = 0;
    mp_size_t i = 0;
 
    mp_limb_t ul0, ul1, ul2, ul3;
    mp_limb_t lpl0, hpl0, rl0;
    mp_limb_t lpl1, hpl1, rl1;
    mp_limb_t lpl2, hpl2, rl2;
    mp_limb_t lpl3, hpl3, rl3;
 
    if (numa == numb) {
        for (; i + 4 <= n; i += 4) {
            ul0 = numa[i + 0];
            ul1 = numa[i + 1];
            ul2 = numa[i + 2];
            ul3 = numa[i + 3];
 
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = ul0 - lpl0;
            cl += (lpl0 > ul0);
 
            _umul64to128_(ul1, b, &lpl1, &hpl1);
            lpl1 += cl;
            cl = (lpl1 < cl) + hpl1;
            lpl1 = ul1 - lpl1;
            cl += (lpl1 > ul1);
 
            _umul64to128_(ul2, b, &lpl2, &hpl2);
            lpl2 += cl;
            cl = (lpl2 < cl) + hpl2;
            lpl2 = ul2 - lpl2;
            cl += (lpl2 > ul2);
 
            _umul64to128_(ul3, b, &lpl3, &hpl3);
            lpl3 += cl;
            cl = (lpl3 < cl) + hpl3;
            lpl3 = ul3 - lpl3;
            cl += (lpl3 > ul3);
 
            numa[i + 0] = lpl0;
            numa[i + 1] = lpl1;
            numa[i + 2] = lpl2;
            numa[i + 3] = lpl3;
        }
        for (; i < n; i++) {
            ul0 = numa[i];
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = ul0 - lpl0;
            cl += (lpl0 > ul0);
            numa[i] = lpl0;
        }
    } else {
        /* seq(numa,numb) */
        for (; i + 4 <= n; i += 4) {
            ul0 = numb[i + 0];
            ul1 = numb[i + 1];
            ul2 = numb[i + 2];
            ul3 = numb[i + 3];
            rl0 = numa[i + 0];
            rl1 = numa[i + 1];
            rl2 = numa[i + 2];
            rl3 = numa[i + 3];
 
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = rl0 - lpl0;
            cl += (lpl0 > rl0);
 
            _umul64to128_(ul1, b, &lpl1, &hpl1);
            lpl1 += cl;
            cl = (lpl1 < cl) + hpl1;
            lpl1 = rl1 - lpl1;
            cl += (lpl1 > rl1);
 
            _umul64to128_(ul2, b, &lpl2, &hpl2);
            lpl2 += cl;
            cl = (lpl2 < cl) + hpl2;
            lpl2 = rl2 - lpl2;
            cl += (lpl2 > rl2);
 
            _umul64to128_(ul3, b, &lpl3, &hpl3);
            lpl3 += cl;
            cl = (lpl3 < cl) + hpl3;
            lpl3 = rl3 - lpl3;
            cl += (lpl3 > rl3);
 
            numa[i + 0] = lpl0;
            numa[i + 1] = lpl1;
            numa[i + 2] = lpl2;
            numa[i + 3] = lpl3;
        }
        for (; i < n; i++) {
            ul0 = numb[i];
            rl0 = numa[i];
            _umul64to128_(ul0, b, &lpl0, &hpl0);
            lpl0 += cl;
            cl = (lpl0 < cl) + hpl0;
            lpl0 = rl0 - lpl0;
            cl += (lpl0 > rl0);
            numa[i] = lpl0;
        }
    }
    return cl;
}

引用了 _umul64to128_(), n , 以及 numb.

函数调用图: