binvert.c 文件参考

#include "../../../include/lammp/numth.h"
#include "../../../include/lammp/lmmpn.h"
#include "../../../include/lammp/impl/mparam.h"
#include "../../../include/lammp/impl/tmp_alloc.h"

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

浏览源代码.

宏定义
#define	ahi   (numa + halfn)
#define	ahi   (numa + halfn)
#define	alo   (numa)
#define	alo   (numa)
#define	k   (tp)
#define	k   (tp)
#define	scratch   (tp + 3 * halfn)
#define	scratch   (tp + 3 * halfn)
#define	xhi   (dst + halfn)
#define	xhi   (dst + halfn)
#define	xlo   (dst)
#define	xlo   (dst)
#define	xlo_sqr   (tp + halfn)
#define	xlo_sqr   (tp + halfn)
#define	xlo_sqr_mul_ahi   (tp + 2 * halfn)
#define	xlo_sqr_mul_ahi   (tp + 2 * halfn)

函数
static void	binvert_mulhi_ (mp_ptr dst, mp_srcptr xp, mp_srcptr ap, mp_size_t n, mp_ptr tp)
void	lmmp_binvert_n_dc_ (mp_ptr restrict dst, mp_srcptr restrict numa, mp_size_t n, mp_ptr restrict tp)
static void	lmmp_mullo_n_ (mp_ptr restrict dst, mp_srcptr restrict numa, mp_srcptr restrict numb, mp_size_t n, mp_ptr restrict tp)
static void	lmmp_sqrlo_n_ (mp_ptr restrict dst, mp_srcptr restrict numa, mp_size_t n, mp_ptr restrict tp)

宏定义说明

ahi [1/2]

#define ahi   (numa + halfn)

ahi [2/2]

#define ahi   (numa + halfn)

alo [1/2]

#define alo   (numa)

alo [2/2]

#define alo   (numa)

k [1/2]

#define k   (tp)

k [2/2]

#define k   (tp)

scratch [1/2]

#define scratch   (tp + 3 * halfn)

scratch [2/2]

#define scratch   (tp + 3 * halfn)

xhi [1/2]

#define xhi   (dst + halfn)

xhi [2/2]

#define xhi   (dst + halfn)

xlo [1/2]

#define xlo   (dst)

xlo [2/2]

#define xlo   (dst)

xlo_sqr [1/2]

#define xlo_sqr   (tp + halfn)

xlo_sqr [2/2]

#define xlo_sqr   (tp + halfn)

xlo_sqr_mul_ahi [1/2]

#define xlo_sqr_mul_ahi   (tp + 2 * halfn)

xlo_sqr_mul_ahi [2/2]

#define xlo_sqr_mul_ahi   (tp + 2 * halfn)

函数说明

binvert_mulhi_()

static void binvert_mulhi_ ( mp_ptr  dst, mp_srcptr  xp, mp_srcptr  ap, mp_size_t  n, mp_ptr  tp  )

inlinestatic

在文件 binvert.c 第 13 行定义.

                                                                                                  {
    if (n < BINVERT_MULHI_MERSENNE_THRESHOLD) {
        lmmp_mul_n_(tp, xp, ap, n);
        lmmp_copy(dst, tp + n, n);
    } else {
        mp_size_t m = lmmp_fft_next_size_((n * 2 + 1) >> 1);
        lmmp_debug_assert(n * 2 > m && m >= n);
        lmmp_mul_mersenne_(tp, m, xp, n, ap, n);
        lmmp_dec(tp);
        mp_size_t fn = m - n;   // 从 tp+n 开始的长度
        mp_size_t sn = n - fn;  // 从 tp 开始的长度
        lmmp_copy(dst, tp + n, fn);
        lmmp_copy(dst + fn, tp, sn);
    }
}

引用了 BINVERT_MULHI_MERSENNE_THRESHOLD, lmmp_copy, lmmp_debug_assert, lmmp_dec, lmmp_fft_next_size_(), lmmp_mul_mersenne_(), lmmp_mul_n_() , 以及 tp.

被这些函数引用 lmmp_binvert_n_dc_().

函数调用图:

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

lmmp_binvert_n_dc_()

void lmmp_binvert_n_dc_	(	mp_ptr restrict	dst,
		mp_srcptr restrict	numa,
		mp_size_t	n,
		mp_ptr restrict	tp
	)

在文件 binvert.c 第 91 行定义.

                                                                                                       {
    lmmp_param_assert(dst != NULL && tp != NULL);
    lmmp_param_assert(numa != NULL && n > 0);
    lmmp_param_assert(numa[0] % 2 == 1);
    if (n == 1) {
        dst[0] = lmmp_binvert_ulong_(numa[0]);
    } else if (n == 2) {
        lmmp_binvert_2_(dst, numa);
    } else if (n == 3) {
        lmmp_binvert_3_(dst, numa);
    } else if (n == 4) {
        lmmp_binvert_4_(dst, numa);
    } else if (n % 2 == 0) {
        mp_size_t halfn = n / 2;
 
#define k               (tp)              // [tp,          halfn]
#define alo             (numa)            // [numa,        halfn]
#define ahi             (numa + halfn)    // [numa+halfn,  halfn]
#define xlo             (dst)             // [dst,         halfn]
#define xhi             (dst + halfn)     // [dst+halfn,   halfn]
#define xlo_sqr         (tp + halfn)      // [tp+halfn,    halfn]
#define xlo_sqr_mul_ahi (tp + 2 * halfn)  // [tp+2*halfn,  halfn]
#define scratch         (tp + 3 * halfn)  // [tp+3*halfn,2*halfn]
//      ________________________________________________________________
// tp : |_________________________5*(n+1)/2____________________________|
//      |   k   | xlo_sqr | xlo_sqr_mul_ahi |   scratch   | remaining  |
//      |_halfn_|__halfn__|______halfn______|___2*halfn___|            |
 
        lmmp_binvert_n_dc_(xlo, alo, halfn, tp);
        binvert_mulhi_(k, xlo, alo, halfn, tp + halfn);
        lmmp_sqrlo_n_(xlo_sqr, xlo, halfn, scratch);
        lmmp_mullo_n_(xlo_sqr_mul_ahi, ahi, xlo_sqr, halfn, scratch);
        lmmp_mullo_n_(xhi, xlo, k, halfn, scratch);
        lmmp_add_n_(xhi, xhi, xlo_sqr_mul_ahi, halfn);
        lmmp_not_(xhi, xhi, halfn);
        lmmp_inc(xhi);
    } else {
        mp_size_t halfn = n / 2 + 1;
        mp_size_t ahin = n - halfn;
 
#define k               (tp)              // [tp,          halfn]
#define alo             (numa)            // [numa,        halfn]
#define ahi             (numa + halfn)    // [numa+halfn,   ahin]
#define xlo             (dst)             // [dst,         halfn]
#define xhi             (dst + halfn)     // [dst+halfn,    ahin]
#define xlo_sqr         (tp + halfn)      // [tp+halfn,     ahin]
#define xlo_sqr_mul_ahi (tp + 2 * halfn)  // [tp+2*halfn,   ahin]
#define scratch         (tp + 3 * halfn)  // [tp+3*halfn, 2*ahin]
//      ________________________________________________________________
// tp : |_________________________5*(n+1)/2____________________________|
//      |    k    | xlo_sqr | xlo_sqr_mul_ahi |   scratch  | remaining |
//      |__halfn__|__halfn__|______halfn______|___2*ahin___|           |
 
        lmmp_binvert_n_dc_(xlo, alo, halfn, tp);
        binvert_mulhi_(k, xlo, alo, halfn, tp + halfn);
        lmmp_sqrlo_n_(xlo_sqr, xlo, ahin, scratch);
        lmmp_mullo_n_(xlo_sqr_mul_ahi, ahi, xlo_sqr, ahin, scratch);
        lmmp_mullo_n_(xhi, xlo, k, ahin, scratch);
        lmmp_add_n_(xhi, xhi, xlo_sqr_mul_ahi, ahin);
        lmmp_not_(xhi, xhi, ahin);
        lmmp_inc(xhi);
    }
}

引用了 ahi, alo, binvert_mulhi_(), k, lmmp_add_n_(), lmmp_binvert_2_(), lmmp_binvert_3_(), lmmp_binvert_4_(), lmmp_binvert_n_dc_(), lmmp_binvert_ulong_(), lmmp_inc, lmmp_mullo_n_(), lmmp_not_(), lmmp_param_assert, lmmp_sqrlo_n_(), scratch, tp, xhi, xlo, xlo_sqr , 以及 xlo_sqr_mul_ahi.

被这些函数引用 lmmp_binvert_n_dc_().

函数调用图:

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

lmmp_mullo_n_()

static void lmmp_mullo_n_ ( mp_ptr restrict  dst, mp_srcptr restrict  numa, mp_srcptr restrict  numb, mp_size_t  n, mp_ptr restrict  tp  )

inlinestatic

在文件 binvert.c 第 42 行定义.

  {
    if (n < MULLO_DC_THRESHOLD) {
        lmmp_mullo_dc_(dst, numa, numb, tp, n);
    } else {
        lmmp_mullo_fft_(dst, numa, numb, n, tp);
    }
}

引用了 lmmp_mullo_dc_(), lmmp_mullo_fft_(), MULLO_DC_THRESHOLD , 以及 tp.

被这些函数引用 lmmp_binvert_n_dc_().

函数调用图:

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

lmmp_sqrlo_n_()

static void lmmp_sqrlo_n_ ( mp_ptr restrict  dst, mp_srcptr restrict  numa, mp_size_t  n, mp_ptr restrict  tp  )

inlinestatic

在文件 binvert.c 第 29 行定义.

  {
    if (n < MULLO_DC_THRESHOLD) {
        lmmp_sqrlo_dc_(dst, numa, tp, n);
    } else {
        lmmp_mullo_fft_(dst, numa, numa, n, tp);
    }
}

引用了 lmmp_mullo_fft_(), lmmp_sqrlo_dc_(), MULLO_DC_THRESHOLD , 以及 tp.

被这些函数引用 lmmp_binvert_n_dc_().

函数调用图:

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