divexact.h 文件参考

#include "../lmmpn.h"

divexact.h 的引用(Include)关系图:

此图展示该文件直接或间接的被哪些文件引用了:

浏览源代码.

宏定义
#define	MODLIMB_INVERSE_15   ((mp_limb_t)0xEEEEEEEEEEEEEEEF)
#define	MODLIMB_INVERSE_3   ((mp_limb_t)0xAAAAAAAAAAAAAAAB)
	Copyright (C) 2026 HJimmyK(Jericho Knox)
#define	MODLIMB_INVERSE_9   ((mp_limb_t)0x8E38E38E38E38E39)

函数
static void	lmmp_divexact_by15_ (mp_ptr dst, mp_srcptr numa, mp_size_t na)
	精确除以15（[dst,na] = [numa,na] / 15）
static void	lmmp_divexact_by3_ (mp_ptr dst, mp_srcptr numa, mp_size_t na)
	精确除以3（[dst,na] = [numa,na] / 3）
static void	lmmp_divexact_by9_ (mp_ptr dst, mp_srcptr numa, mp_size_t na)
	精确除以9（[dst,na] = [numa,na] / 9）

宏定义说明

MODLIMB_INVERSE_15

#define MODLIMB_INVERSE_15   ((mp_limb_t)0xEEEEEEEEEEEEEEEF)

在文件 divexact.h 第 73 行定义.

MODLIMB_INVERSE_3

#define MODLIMB_INVERSE_3   ((mp_limb_t)0xAAAAAAAAAAAAAAAB)

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

在文件 divexact.h 第 21 行定义.

MODLIMB_INVERSE_9

#define MODLIMB_INVERSE_9   ((mp_limb_t)0x8E38E38E38E38E39)

在文件 divexact.h 第 47 行定义.

函数说明

lmmp_divexact_by15_()

static void lmmp_divexact_by15_ ( mp_ptr  dst, mp_srcptr  numa, mp_size_t  na  )

inlinestatic

精确除以15（[dst,na] = [numa,na] / 15）

参数

dst	结果存储位置
numa	被除数
na	被除数长度

警告

eqsep(dst,numa), na>0

在文件 divexact.h 第 82 行定义.

                                                                                 {
    mp_limb_t c = 0;
    mp_limb_t l, q, s, t, carry;
    mp_size_t i = 0;
    do {
        s = numa[i];
        l = s - c;
        c = l > s;
        q = l * MODLIMB_INVERSE_15;
        dst[i] = q;
        // carry from 15*q = (q<<4) - q
        t = q << 4;
        carry = (q >> (LIMB_BITS - 4)) - (t < q);
        c += carry;
    } while (++i < na);
}

引用了 c, l, LIMB_BITS, MODLIMB_INVERSE_15, s , 以及 t.

被这些函数引用 lmmp_toom_interp7_().

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

lmmp_divexact_by3_()

static void lmmp_divexact_by3_ ( mp_ptr  dst, mp_srcptr  numa, mp_size_t  na  )

inlinestatic

精确除以3（[dst,na] = [numa,na] / 3）

参数

dst	结果存储位置
numa	被除数
na	被除数长度

警告

eqsep(dst,numa), na>0

在文件 divexact.h 第 30 行定义.

                                                                                {
    mp_limb_t c = 0;
    mp_limb_t l, q, s;
    mp_size_t i = 0;
    do {
        s = numa[i];
        l = s - c;
        c = l > s;
        q = l * MODLIMB_INVERSE_3;
        dst[i] = q;
        l = q + q;
        c += l < q;
        l += q;
        c += l < q;
    } while (++i < na);
}

引用了 c, l, MODLIMB_INVERSE_3 , 以及 s.

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

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

lmmp_divexact_by9_()

static void lmmp_divexact_by9_ ( mp_ptr  dst, mp_srcptr  numa, mp_size_t  na  )

inlinestatic

精确除以9（[dst,na] = [numa,na] / 9）

参数

dst	结果存储位置
numa	被除数
na	被除数长度

警告

eqsep(dst,numa), na>0

在文件 divexact.h 第 56 行定义.

                                                                                {
    mp_limb_t c = 0;
    mp_limb_t l, q, s, t, carry;
    mp_size_t i = 0;
    do {
        s = numa[i];
        l = s - c;
        c = l > s;
        q = l * MODLIMB_INVERSE_9;
        dst[i] = q;
        // carry from 9*q = (q<<3) + q
        t = q << 3;
        carry = (q >> (LIMB_BITS - 3)) + ((t + q) < t);
        c += carry;
    } while (++i < na);
}

引用了 c, l, LIMB_BITS, MODLIMB_INVERSE_9, s , 以及 t.

被这些函数引用 lmmp_toom_interp7_().

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