sqrt_1.c 文件参考

#include <math.h>
#include "../../../include/lammp/impl/mparam.h"
#include "../../../include/lammp/numth.h"

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

浏览源代码.

函数
mp_limb_t	lmmp_sqrt_1_ (mp_ptr dsts, mp_limb_t x)
	计算算术平方根 floor(sqrt(x))
mp_limb_t	lmmp_sqrt_2_ (mp_ptr dsts, mp_ptr dstr, mp_srcptr numa)
	计算算术平方根 floor(sqrt([numa,2]))
ulong	lmmp_sqrt_ulong_ (ulong a)
	Copyright (C) 2026 HJimmyK(Jericho Knox)

函数说明

lmmp_sqrt_1_()

mp_limb_t lmmp_sqrt_1_	(	mp_ptr	dsts,
		mp_limb_t	x
	)

计算算术平方根 floor(sqrt(x))

参数

dsts	结果指针（1个limb）
x	被开方数

警告

x>=B/4, dsts!=NULL

注解

[dsts,1]=floor(sqrt(x)), return remainder

返回

余数

在文件 sqrt_1.c 第 33 行定义.

                                                 {
    lmmp_param_assert(x >= LIMB_B_4);
    mp_limb_t s = lmmp_sqrt_ulong_(x);
    *dsts = s;
    return x - s * s;
}

引用了 LIMB_B_4, lmmp_param_assert, lmmp_sqrt_ulong_(), n , 以及 s.

被这些函数引用 lmmp_sqrt_() , 以及 lmmp_sqrt_2_().

函数调用图:

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

lmmp_sqrt_2_()

mp_limb_t lmmp_sqrt_2_	(	mp_ptr	dsts,
		mp_ptr	dstr,
		mp_srcptr	numa
	)

计算算术平方根 floor(sqrt([numa,2]))

参数

dsts	结果指针（1个limb）
dstr	余数指针（1个limb）
numa	被开方数指针

警告

numa[1]>=B/4, dsts!=NULL, dstr!=NULL, numa!=NULL

注解

[dsts,1]=floor(sqrt([numa,2])), rh:[dstr,1]=remainder, return rh

返回

余数的高位

在文件 sqrt_1.c 第 40 行定义.

                                                                 {
    lmmp_param_assert(numa[1] >= LIMB_B_4);
    mp_limb_t rl, s, q, al, u;
    mp_slimb_t rh;
 
    rl = lmmp_sqrt_1_(&s, numa[1]);
    al = numa[0];
 
    //(r:alh)/2
    rl = rl << 31 | al >> 33;
    q = rl / s;
    q -= q >> 32;
 
    u = rl - s * q;
    s = s << 32 | q;
    rh = u >> 31;
    rl = (u << 33) | (al & (((mp_limb_t)1 << 33) - 1));
 
    q *= q;
    rh -= rl < q;
    rl -= q;
    if (rh < 0) {
        rl += s;
        rh += rl < s;
        --s;
        rl += s;
        rh += rl < s;
    }
 
    dsts[0] = s;
    dstr[0] = rl;
    return rh;
}

引用了 LIMB_B_4, lmmp_param_assert, lmmp_sqrt_1_(), n , 以及 s.

被这些函数引用 lmmp_sqrt_divide_().

函数调用图:

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

lmmp_sqrt_ulong_()

ulong lmmp_sqrt_ulong_

(

ulong

a

)

Copyright (C) 2026 HJimmyK(Jericho Knox)

计算算术平方根 floor(sqrt(a))

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

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

                                {
    ulong is;
 
    is = (ulong)sqrt((double)a);
 
    is -= (is * is > a);
    if (is == (1ULL << 32))
        is--;
    return is;
}

引用了 n.

被这些函数引用 lmmp_nthroot_ulong_() , 以及 lmmp_sqrt_1_().

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