1 Star 0 Fork 0

zhouxs1023/eigenmath_pratt

加入 Gitee
与超过 1200万 开发者一起发现、参与优秀开源项目,私有仓库也完全免费 :)
免费加入
文件
该仓库未声明开源许可证文件(LICENSE),使用请关注具体项目描述及其代码上游依赖。
克隆/下载
mcmp.cpp 1.31 KB
一键复制 编辑 原始数据 按行查看 历史
Calin Barbat 提交于 2018-02-24 10:52 . Initial commit.
// Bignum compare
//
// returns
//
// -1 a < b
//
// 0 a = b
//
// 1 a > b
#include "stdafx.h"
#include "defs.h"
int
mcmp(unsigned int *a, unsigned int *b)
{
int i;
if (MSIGN(a) == -1 && MSIGN(b) == 1)
return -1;
if (MSIGN(a) == 1 && MSIGN(b) == -1)
return 1;
// same sign
if (MLENGTH(a) < MLENGTH(b)) {
if (MSIGN(a) == 1)
return -1;
else
return 1;
}
if (MLENGTH(a) > MLENGTH(b)) {
if (MSIGN(a) == 1)
return 1;
else
return -1;
}
// same length
for (i = MLENGTH(a) - 1; i > 0; i--)
if (a[i] != b[i])
break;
if (a[i] < b[i]) {
if (MSIGN(a) == 1)
return -1;
else
return 1;
}
if (a[i] > b[i]) {
if (MSIGN(a) == 1)
return 1;
else
return -1;
}
return 0;
}
int
mcmpint(unsigned int *a, int n)
{
int t;
unsigned int *b;
b = mint(n);
t = mcmp(a, b);
mfree(b);
return t;
}
#if SELFTEST
void
test_mcmp(void)
{
int i, j, k;
unsigned int *x, *y;
logout("testing mcmp\n");
for (i = -1000; i < 1000; i++) {
x = mint(i);
for (j = -1000; j < 1000; j++) {
y = mint(j);
k = mcmp(x, y);
if (i == j && k != 0) {
logout("failed\n");
errout();
}
if (i < j && k != -1) {
logout("failed\n");
errout();
}
if (i > j && k != 1) {
logout("failed\n");
errout();
}
mfree(y);
}
mfree(x);
}
logout("ok\n");
}
#endif
马建仓 AI 助手
尝试更多
代码解读
代码找茬
代码优化
C
1
https://gitee.com/zhouxs1023/eigenmath_pratt.git
git@gitee.com:zhouxs1023/eigenmath_pratt.git
zhouxs1023
eigenmath_pratt
eigenmath_pratt
master

搜索帮助