trialdiv.c 文件参考

#include "../../../include/lammp/numth.h"
#include "../../../include/lammp/lmmpn.h"
#include "../../../include/lammp/impl/prime_table.h"
#include "../../../include/lammp/impl/tmp_alloc.h"
#include "../../../include/lammp/impl/ele_mul.h"

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

浏览源代码.

宏定义
#define	MAX_T   0xffffffffffffull
	Copyright (C) 2026 HJimmyK(Jericho Knox)

函数
ushortp	lmmp_trialdiv_ (mp_srcptr restrict np, mp_size_t nn, ushort N, ushort *rn)
static ushortp	lmmp_trialdiv_short_ (mp_srcptr restrict np, mp_size_t nn, ushort N, ushort *rn)

宏定义说明

MAX_T

#define MAX_T   0xffffffffffffull

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

在文件 trialdiv.c 第 22 行定义.

函数说明

lmmp_trialdiv_()

ushortp lmmp_trialdiv_	(	mp_srcptr restrict	np,
		mp_size_t	nn,
		ushort	N,
		ushort *	rn
	)

在文件 trialdiv.c 第 80 行定义.

                                                                                  {
    ushort primen = lmmp_prime_cnt16_(N);
    if (primen > 4 * nn) {
        return lmmp_trialdiv_short_(np, nn, N, rn);
    } else {
    /*
    if [np,nn] is too large, we calculate the product of primes up to N,
    and divide [np,nn] by the product to get remainder.
    */
        TEMP_S_DECL;
        ulongp restrict pp = SALLOC_TYPE(primen / 4 + 1, ulong);
        ulong t = 1;
        mp_size_t pn = 0;
        for (ushort i = 0; i < primen; i++) {
            t *= prime_short_table[i];
            if (t > MAX_T) {
                pp[pn++] = t;
                t = 1;
            }
        }
        if (t > 1) {
            pp[pn++] = t;
        }
        mp_ptr restrict prod = SALLOC_TYPE(pn * 2, mp_limb_t);
        pn = lmmp_elem_mul_ulong_(prod, pp, pn, prod + pn);
 
        lmmp_div_(NULL, prod, np, nn, prod, pn);
        while (pn > 0 && prod[pn - 1] == 0) --pn;
        ushortp restrict ret;
        if (pn == 0) {
            ret = ALLOC_TYPE(primen, ushort);
            for (ushort i = 0; i < primen; i++) {
                ret[i] = prime_short_table[i];
            }
            *rn = primen;
        } else {
            ret = lmmp_trialdiv_short_(prod, pn, N, rn);
        } 
        TEMP_S_FREE;
        return ret;
    }
}

引用了 ALLOC_TYPE, lmmp_div_(), lmmp_elem_mul_ulong_(), lmmp_prime_cnt16_(), lmmp_trialdiv_short_(), MAX_T, n, prime_short_table, SALLOC_TYPE, t, TEMP_S_DECL , 以及 TEMP_S_FREE.

函数调用图:

lmmp_trialdiv_short_()

static ushortp lmmp_trialdiv_short_ ( mp_srcptr restrict  np, mp_size_t  nn, ushort  N, ushort *  rn  )

static

在文件 trialdiv.c 第 24 行定义.

                                                                                               {
    if (np == NULL || nn == 0) {
        *rn = 0;
        return NULL;
    }
    ulong t = 1;
    ushort idx[20];
    uint idx_cnt = 0;
    ushort retn_max = 10;
    ushortp ret = ALLOC_TYPE(retn_max, ushort);
    ushort retn = 0;
    if (!(np[0] & 1)) {
        ret[retn++] = 2;
    }
    for (mp_size_t i = 1;; i++) {
        ushort p = prime_short_table[i];
        if (p > N || i >= PRIME_SHORT_TABLE_SIZE - 1) break;
        t *= p;
        idx[idx_cnt++] = p;
        if (t > MAX_T || idx_cnt == 20) {
            mp_limb_t r = lmmp_mod_1_(np, nn, t);
            for (uint j = 0; j < idx_cnt; j++) {
                if (r == 0 || r % idx[j] == 0) {
                    ret[retn++] = idx[j];
                    if (retn == retn_max) {
                        retn_max = retn_max * 12 / 10;
                        ret = REALLOC_TYPE(ret, retn_max, ushort);
                    }
                }
            }
            idx_cnt = 0;
            t = 1;
        }
    }
    if (idx_cnt > 0) {
        mp_limb_t r = lmmp_mod_1_(np, nn, t);
        for (uint j = 0; j < idx_cnt; j++) {
            if (r == 0 || r % idx[j] == 0) {
                ret[retn++] = idx[j];
                if (retn == retn_max) {
                    retn_max += 10;
                    ret = REALLOC_TYPE(ret, retn_max, ushort);
                }
            }
        }
    }
    if (retn == 0) {
        lmmp_free(ret);
        *rn = 0;
        return NULL;
    } else {
        *rn = retn;
        return ret;
    }
}

引用了 ALLOC_TYPE, lmmp_free(), lmmp_mod_1_(), MAX_T, n, prime_short_table, PRIME_SHORT_TABLE_SIZE, REALLOC_TYPE , 以及 t.

被这些函数引用 lmmp_trialdiv_().

函数调用图:

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