问题标签 [pari-gp]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
3 回答
568 浏览

computer-algebra-systems - 有没有办法使 Pari/GP 中的小数合理化?

我正在寻找一种在 PARI/GP 计算器 (gp) 中将十进制值转换为分数的自动化方法。比如我想把0.759765625变成389/512。我知道我可以通过输入“759765625/10^9”手动执行此操作,但我想要一种更简单的方法,例如“合理化(0.759765625);” 将在千里马工作。

如果它不能直接做到这一点,也许有一个计算小数位数的功能?然后我可以将该函数的结果提高 10。(该函数必须计算 0 的前导小数位才有用。)

0 投票
3 回答
510 浏览

pari-gp - Pari/GP 将整数与实数进行比较

我想测试一个整数是否是 Pari-gp 中的完美幂。该测试sqrt(n)==floor(sqrt(n))适用于测试正方形,但对于其他所有幂都失败:sqrtn(n,k)==floor(sqrtn(n,k))使用k >=3.

我想可能是因为一个数字是实数,另一个是整数。该测试仍然适用于正方形。我究竟做错了什么?

0 投票
1 回答
505 浏览

pari-gp - 在 PARI/GP 中舍入一个数字

我知道“\p x”命令设置所有计算的精度,但我正在寻找一些不同的东西。给定一个以高精度计算的给定数字,我想仅针对我的代码的一部分将其舍入到较低的精度。这样做的原因是取复数的绝对值往往会产生一些舍入错误,从而导致以后出现问题。我希望这些绝对值是相同的,并且没有他们现在的微小错误。

换句话说,我怎样才能让 1.xxxxxxxxxxxxxx -> 1.xxxxxxx 只为一个特定的数字?

0 投票
1 回答
1040 浏览

vector - 如何在 PARI/GP 中获取行向量的总长度

从一定数量的除数中,我试图随机选择一个除数。

但是我无法随机获得一个数字,因为我找不到任何可以让我获得行向量长度的函数。

例子:

希望这里有人可以帮助我。谢谢!:)

0 投票
1 回答
2323 浏览

pari-gp - Pari/GP:用户定义函数

我已经定义了一些 arity 1 的函数,比如func1(-)func2(-)。我已经对它们进行了测试,发现它们实际上做了他们应该做的事情。

我希望定义第三个函数,比如func3(-),输出 和 的func1(-)差异func2(-)。这就是我所做的

然而,它并没有返回它应该返回的东西。让我们假设func1(5)输出 10 和func2(5)输出 2。那么,func3(5)应该输出 8,对吗?相反,它返回func1(5)一行中的输出,func2(2)另一行中的输出,然后返回零(即使相应输出的差异不是 0)。

你知道 的定义有什么问题func3(-)吗?

0 投票
3 回答
1944 浏览

python - 从 Python 调用 PARI/GP

我想从 Python 调用PARI/GP只是为了计算我定义nextprime(n)的不同s 的函数。n不幸的是,我无法安装pari-pythonos.system ,所以我想我会在 Python 中使用命令行调用它。但是,我在手册页中看不到如何让 PARI/GP 在非交互模式下运行。有没有办法做到这一点?

0 投票
1 回答
3328 浏览

pari-gp - 如何在 PARI/GP 中制作多行嵌套 for 循环?

如何在 PARI/GP 中创建跨越每个级别的多行的嵌套循环?我经常需要在 for 循环中做很多事情,而且为了可读性,我不喜欢将循环写在一行上。对于一个变量的循环,我一直在这样做:

但是,对于嵌套循环,我只设法将换行符放在一个级别。这有效:

这也有效:

但是,这是我真正想做的:

最后一个例子不起作用;它给出了一个错误:

我在 OS X 10.8.3 上使用 PARI/GP 2.5.3。我将脚本写入文件并使用Bashnested_for.gp运行它们。gp ./nested_for.gp

0 投票
0 回答
441 浏览

exception - Pari/GP 异常/中断循环

ellpow(E, P, m) 总是会抛出异常:

其中 x 和 y 是整数。

我想捕获值 x,而不是完成程序以便以后使用它。

代码片段是:

整个程序是:

0 投票
1 回答
1235 浏览

pari - 如何使用 PARI C 库

我已经搜索了一些关于如何在C 程序中使用 PARI 库的教程,但不幸的是,我只找到了在命令行界面中使用 GP 包的通用教程。

有人可以帮忙吗?例如,我想初始化一个 3x3 矩阵并找到它的 1000 次方。在 gp(CLI)中,这可以通过键入以下内容来实现:

但我发现在C 源文件中做这个安静简单的事情没有任何帮助。是否有提供此类示例的综合教程或文档?任何评论将不胜感激!

0 投票
1 回答
162 浏览

pari-gp - 是否有一种已知的电力塔算法以管理所有情况的数字为模?

我想在 PARI/GP 中实现计算

a_1 ^ a_2 ^ ... ^ a_n (mod m)

它管理所有情况,尤其是在 phi 链中出现高权力的情况。

有谁知道这样的实现?