找回密码
 欢迎注册
查看: 10985|回复: 7

[求助] 关于pari/gp的效率,以及想问一下该如何手动编译……编译几乎总是不过

[复制链接]
发表于 2018-6-22 11:42:01 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
本帖最后由 .·.·. 于 2018-6-22 11:47 编辑

折腾了好久……终于发现,凭自己大概率是不可能成功的……
首先,想问一下有没有给pari/gp的循环加速的方法
一个简单的循环
for(a=1,100000000,c++)
需要跑10秒……
作为对比,R(一个超慢的统计软件)只需要1/3的时间就可以跑完b=0;for(a in 1:100000000)b=b+1
pari/gp究竟丧心病狂地做了什么?

以及,看到有一个gp2c很不错……windows强行编译……强行浪费了3个小时
cygwin版本一次成功……然而Cygwin多线程比windows单线程还慢,同时,gp2c的生成代码似乎也是减速而非加速……

换MSYS2,每次编译都费尽千辛万苦(比如在编GMP的时候替换掉gcc以过掉一些出了问题的汇编)
当然,最后编译是成功的……然而执行的时候不知为什么总提示某一个地址不是32位程序能用的……
我明明编的64位……
./configure --with-gmp=/d/pen/pari-2.9.5/msys2-ver/gmp-6.1.2 --with-gmp-include=/d/pen/pari-2.9.5/msys2-ver/gmp-6.1.2 --with-gmp-lib=/d/pen/pari-2.9.5/msys2-ver/gmp-6.1.2/.libs --with-readline --with-runtime-perl=/usr/bin --enable-tls --mt=pthread --host=msys2-x86_64

pa2c,由于最后链接的是人家的pari,最后的编译直接失败,于是跑不动……
./configure --with-perl=/usr/bin --with-paricfg=/d/pen/pari-2.9.5/Pari64-2-9-5/lib/pari/pari.cfg

有没有大牛分享一下编译pari/gp的经验呢?
难道只有用linux才能解决一切问题吗?

最后还有一个问题,如何判断某一个元素是否在集合中
比如我们有一个程序,每次会生成一个整数i
我们想判断这个整数i是否跟以往产生的整数重复
于是我们可以构造一个集合a=Set([])
之后每次执行
len=length(a)
a=Set([i,a])
if(length(a)-len,没重复,有重复)
有没有更快更方便的方法呢?比如python可以直接用i in a和a.add(i)来完成一切操作……然而萌新表示,看看这个程序的文档实在不知道该从哪里找更好的表述方法……
还望dalao不吝赐教
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2018-6-22 16:49:12 | 显示全部楼层
pari/gp的确效率不是很好,能够使用的内存也是预先分配的,不能使用太多内存。Linux比Windows效率明显要高,但是还不是特别好.

点评

BTW我也是醉了……明明知道Win10有linux子系统,非要作死折腾Cygwin跟MSYS2…… 如果一切顺利,大约明天我会写一下在windows自带的ubuntu上面运行pari/gp的感受:)  发表于 2018-6-22 19:52
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2018-6-22 18:36:07 | 显示全部楼层
mathe 发表于 2018-6-22 16:49
pari/gp的确效率不是很好,能够使用的内存也是预先分配的,不能使用太多内存。Linux比Windows效率明显要高 ...


弄了集合操作
差不多是用a=Set(compat(l,a))跟length(a)执行的
刚开始每次循环一秒就够
算了五个小时
现在一次循环十几秒
简直吓人

话说算数论的话有什么更快的程序吗
难道真的要上网下载盗版mathematica吗……

补充内容 (2018-6-23 02:44):
找到程序为什么从1.7s涨到17s的原因了……
原因是记错了(……)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2018-6-22 22:38:06 | 显示全部楼层
len=length(a)
a=Set([i,a])
if(length(a)-len,没重复,有重复)

将插入换成集合的搜索操作比较合适,即换成setsearch函数。


setsearch(S, x, {flag = 0})
Determines whether x belongs to the set S (see setisset).

We first describe the default behaviour, when flag is zero or omitted. If x belongs to the set S, returns the index j such that S[j] = x, otherwise returns 0.

  ? T = [7,2,3,5]; S = Set(T);
  ? setsearch(S, 2)
  %2 = 1
  ? setsearch(S, 4)      \\ not found
  %3 = 0
  ? setsearch(T, 7)      \\ search in a randomly sorted vector
  %4 = 0 \\ WRONG !
If S is not a set, we also allow sorted lists with respect to the cmp sorting function, without repeated entries, as per listsort(L,1); otherwise the result is undefined.



根据PARI/Gp的Set的setsearch的解释来看,是有序的集合,所以其实现应该是平衡树,插入和查找都是O(lgn)复杂度。
而根据你的需求来看,有两个高频操作,插入和查找。所以最理想的情况应该是字典来实现。插入和查找都是O(1)复杂度。如c++11 unordered_map结构。
======

点评

才发现大神你用的是2.10?为了跑gp2c特地装了2.9……然而2.9没这些函数……桑心  发表于 2018-6-24 18:18
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2018-6-23 01:41:54 | 显示全部楼层
本帖最后由 .·.·. 于 2018-6-23 02:46 编辑
wayne 发表于 2018-6-22 22:38
将插入换成集合的搜索操作比较合适,即换成setsearch函数。

……
找到程序为什么从1.7s涨到17s的原因了……
原因是记错了(……)

很抱歉虽然这一步是唯一可能会随时间变化变得复杂的步骤,但单凭一步的确不足以达到我描述的效果……
现在用linux,准备强行开多线程:)
我的系统差点就重装成ubuntu了……

悬崖勒马……

点评

嗯,说得对。 我这个也只是提出改进的空间,具体改进多少,得看具体的值  发表于 2018-6-23 09:53
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-3-29 14:56 , Processed in 0.044038 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表