LAMMP 4.2.0
Lamina High-Precision Arithmetic Library
载入中...
搜索中...
未找到
宏定义 | 函数
signed.h 文件参考
#include "../lmmpn.h"
+ signed.h 的引用(Include)关系图:
+ 此图展示该文件直接或间接的被哪些文件引用了:

浏览源代码.

宏定义

#define INLINE_   static inline
 Copyright (C) 2026 HJimmyK(Jericho Knox)
 

函数

static mp_ssize_t lmmp_add_signed_ (mp_ptr dst, mp_srcptr numa, mp_ssize_t na, mp_srcptr numb, mp_ssize_t nb)
 计算带符号数的加法
 
static mp_ssize_t lmmp_mul_signed_ (mp_ptr dst, mp_srcptr numa, mp_ssize_t na, mp_srcptr numb, mp_ssize_t nb)
 计算带符号数的乘法
 
static mp_ssize_t lmmp_sqr_signed_ (mp_ptr dst, mp_srcptr numa, mp_ssize_t na)
 计算带符号数的平方
 

宏定义说明

◆ INLINE_

#define INLINE_   static inline

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

在文件 signed.h 第 22 行定义.

函数说明

◆ lmmp_add_signed_()

static mp_ssize_t lmmp_add_signed_ ( mp_ptr  dst,
mp_srcptr  numa,
mp_ssize_t  na,
mp_srcptr  numb,
mp_ssize_t  nb 
)
inlinestatic

计算带符号数的加法

参数
dst结果指针，自行保证有足够的空间，需要的最多的空间为max(abs(na),abs(nb)) + 1
numa第一个数的指针
na第一个数的长度，为负数表示此数为负数，绝对值表示实际长度
numb第二个数的指针
nb第二个数的长度，为负数表示此数为负数，绝对值表示实际长度
返回
结果的长度，为负数表示此数为负数，绝对值表示实际长度
警告
dst!=NULL, eqsep(dst,[numa|numb])

在文件 signed.h 第 35 行定义.

35 {
36 if (na == 0 && nb != 0) {
37 lmmp_copy(dst, numb, LMMP_ABS(nb));
38 return nb;
39 } else if (na != 0 && nb == 0) {
40 lmmp_copy(dst, numa, LMMP_ABS(na));
41 return na;
42 } else if (na == 0 && nb == 0) {
43 dst[0] = 0;
44 return 0;
45 } else if (na < 0 && nb < 0) {
46 na = -na;
47 nb = -nb;
48 if (na > nb) {
49 dst[na] = lmmp_add_(dst, numa, na, numb, nb);
50 return dst[na] == 0 ? -na : -(na + 1);
51 } else {
52 dst[nb] = lmmp_add_(dst, numb, nb, numa, na);
53 return dst[nb] == 0 ? -nb : -(nb + 1);
54 }
55 } else if (na > 0 && nb > 0) {
56 if (na > nb) {
57 dst[na] = lmmp_add_(dst, numa, na, numb, nb);
58 return dst[na] == 0 ? na : na + 1;
59 } else {
60 dst[nb] = lmmp_add_(dst, numb, nb, numa, na);
61 return dst[nb] == 0 ? nb : nb + 1;
62 }
63 } else if (na < 0 && nb > 0) {
64 na = -na;
65 if (na < nb) {
66 lmmp_sub_(dst, numb, nb, numa, na);
67 while (dst[nb - 1] == 0 && nb > 0) {
68 --nb;
69 }
70 return nb;
71 } else if (na > nb) {
72 lmmp_sub_(dst, numa, na, numb, nb);
73 while (dst[na - 1] == 0 && na > 0) {
74 --na;
75 }
76 return -na;
77 } else {
78 int cmp = lmmp_cmp_(numa, numb, na);
79 if (cmp < 0) {
80 lmmp_sub_(dst, numb, nb, numa, na);
81 while (dst[nb - 1] == 0 && nb > 0) {
82 --nb;
83 }
84 return nb;
85 } else if (cmp > 0) {
86 lmmp_sub_(dst, numa, na, numb, nb);
87 while (dst[na - 1] == 0 && na > 0) {
88 --na;
89 }
90 return -na;
91 } else {
92 dst[0] = 0;
93 return 0;
94 }
95 }
96 } else {
97 /* na > 0 && nb < 0 */
98 nb = -nb;
99 if (na < nb) {
100 lmmp_sub_(dst, numb, nb, numa, na);
101 while (dst[nb - 1] == 0 && nb > 0) {
102 --nb;
103 }
104 return -nb;
105 } else if (na > nb) {
106 lmmp_sub_(dst, numa, na, numb, nb);
107 while (dst[na - 1] == 0 && na > 0) {
108 --na;
109 }
110 return na;
111 } else {
112 int cmp = lmmp_cmp_(numa, numb, na);
113 if (cmp < 0) {
114 lmmp_sub_(dst, numb, nb, numa, na);
115 while (dst[nb - 1] == 0 && nb > 0) {
116 --nb;
117 }
118 return -nb;
119 } else if (cmp > 0) {
120 lmmp_sub_(dst, numa, na, numb, nb);
121 while (dst[na - 1] == 0 && na > 0) {
122 --na;
123 }
124 return na;
125 } else {
126 dst[0] = 0;
127 return 0;
128 }
129 }
130 }
131}
lmmp_copy
#define lmmp_copy(dst, src, n)
Definition lmmp.h:367
LMMP_ABS
#define LMMP_ABS(x)
Definition lmmp.h:349
lmmp_add_
static mp_limb_t lmmp_add_(mp_ptr dst, mp_srcptr numa, mp_size_t na, mp_srcptr numb, mp_size_t nb)
大数加法静态内联函数 [dst,na]=[numa,na]+[numb,nb]
Definition lmmpn.h:1050
lmmp_cmp_
static int lmmp_cmp_(mp_srcptr numa, mp_srcptr numb, mp_size_t n)
大数比较函数（内联）
Definition lmmpn.h:996
lmmp_sub_
static mp_limb_t lmmp_sub_(mp_ptr dst, mp_srcptr numa, mp_size_t na, mp_srcptr numb, mp_size_t nb)
大数减法静态内联函数 [dst,na]=[numa,na]-[numb,nb]
Definition lmmpn.h:1064
numb
#define numb
n
#define n

引用了 LMMP_ABS, lmmp_add_(), lmmp_cmp_(), lmmp_copy, lmmp_sub_() , 以及 numb.

被这些函数引用 lmmp_mat22_mul_basecase_(), lmmp_mat22_mul_strassen_(), lmmp_mat22_sqr_basecase_() , 以及 lmmp_mat22_sqr_strassen_().

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

◆ lmmp_mul_signed_()

static mp_ssize_t lmmp_mul_signed_ ( mp_ptr  dst,
mp_srcptr  numa,
mp_ssize_t  na,
mp_srcptr  numb,
mp_ssize_t  nb 
)
inlinestatic

计算带符号数的乘法

参数
dst结果指针，自行保证有足够的空间，需要空间为abs(na)+abs(nb)
numa第一个数的指针
na第一个数的长度，为负数表示此数为负数，绝对值表示实际长度
numb第二个数的指针
nb第二个数的长度，为负数表示此数为负数，绝对值表示实际长度
返回
结果的长度，为负数表示此数为负数，绝对值表示实际长度
警告
dst!=NULL, sep(dst,[numa|numb])

在文件 signed.h 第 143 行定义.

143 {
144 mp_ssize_t sign = ((na > 0) ^ (nb > 0)) ? -1 : 1;
145 na = LMMP_ABS(na);
146 nb = LMMP_ABS(nb);
147 if (na == 0 || nb == 0) {
148 dst[0] = 0;
149 return 0;
150 }
151 if (na < nb)
152 lmmp_mul_(dst, numb, nb, numa, na);
153 else
154 lmmp_mul_(dst, numa, na, numb, nb);
155 na += nb;
156 na -= (dst[na - 1] == 0);
157 return sign * na;
158}
mp_ssize_t
int64_t mp_ssize_t
Definition lmmp.h:79
lmmp_mul_
void lmmp_mul_(mp_ptr dst, mp_srcptr numa, mp_size_t na, mp_srcptr numb, mp_size_t nb)
不等长大数乘法操作 [dst,na+nb] = [numa,na] * [numb,nb]

引用了 LMMP_ABS, lmmp_mul_() , 以及 numb.

被这些函数引用 lmmp_mat22_mul_basecase_(), lmmp_mat22_mul_strassen_(), lmmp_mat22_sqr_basecase_() , 以及 lmmp_mat22_sqr_strassen_().

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

◆ lmmp_sqr_signed_()

static mp_ssize_t lmmp_sqr_signed_ ( mp_ptr  dst,
mp_srcptr  numa,
mp_ssize_t  na 
)
inlinestatic

计算带符号数的平方

参数
dst结果指针，自行保证有足够的空间，需要空间为abs(na)*2
numa第一个数的指针
na第一个数的长度，为负数表示此数为负数，绝对值表示实际长度
返回
结果的长度，平方一定为非负数
警告
dst!=NULL, sep(dst,numa)

在文件 signed.h 第 168 行定义.

168 {
169 if (na < 0) {
170 na = -na;
171 } else if (na == 0) {
172 dst[0] = 0;
173 return 0;
174 }
175 lmmp_sqr_(dst, numa, na);
176 na <<= 1;
177 na -= (dst[na - 1] == 0) ? 1 : 0;
178 return na;
179}
lmmp_sqr_
#define lmmp_sqr_
Definition inlines.h:166

引用了 lmmp_sqr_ , 以及 n.

被这些函数引用 lmmp_mat22_sqr_basecase_() , 以及 lmmp_mat22_sqr_strassen_().

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