应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 矩阵求特征值和特征向量用arpack和lapack哪个好些？
121/2返回列表12下一页
查看: 2485  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

田山东

捐助贵宾 (著名写手)

[求助] 矩阵求特征值和特征向量用arpack和lapack哪个好些？

要解一个很大的矩阵，10000*10000，求特征值和特征向量，大家有没有什么好办法可以算的速度快点？arpack好些还是lapack好些？还有没有其他办法？有没有可能把大矩阵化成小矩阵啊？
arpack虽然可以求解，但是一直没找到它的子程序。有同学说netlib上有，但是也没找到相关子程序。很希望得到进一步的帮助。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
everythinghasitsseason.enjoyyourlife.
1楼 2012-10-09 10:55:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lurencyj

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
田山东: 金币+5, ★★★★★最佳答案, 能否再进一步？ 2012-10-09 17:27:17
xzhdty: 金币+2, 谢谢参与 2012-10-11 09:09:31
arpack和lapack的算法是不一样的，选哪个需要看具体需求。
对于大型的矩阵（太大以致没法在内存中存下所有元素），那得用arpack。
对于不需要全部本征值的，arpack和lapack都能适用。
如果需要全部本征值，那得用lapack。
arpack的子程序实例非常实用，楼主好好找找。Examples. Good luck!

矩阵大小：一个10000*10000双精度矩阵全部存下来需要的内存是10000×10000×8/1024/1024 ～ 763MB空间，用lapack没有问题吧。

[ 发自手机版 http://muchong.com/3g ]
赞一下(1人) 回复此楼
很女子很弓虽大
2楼2012-10-09 16:35:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

田山东

捐助贵宾 (著名写手)
引用回帖:
2楼: Originally posted by lurencyj at 2012-10-09 16:35:18
arpack和lapack的算法是不一样的，选哪个需要看具体需求。
对于大型的矩阵（太大以致没法在内存中存下所有元素），那得用arpack。
对于不需要全部本征值的，arpack和lapack都能适用。
如果需要全部本征值，那得用 ...

老兄，哪个网址有介绍arpack子程序的？我找好久一直没找到。可以帮下忙啊？
赞一下 回复此楼
everythinghasitsseason.enjoyyourlife.
3楼2012-10-09 17:28:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)
引用回帖:
2楼: Originally posted by lurencyj at 2012-10-09 16:35:18
arpack和lapack的算法是不一样的，选哪个需要看具体需求。
对于大型的矩阵（太大以致没法在内存中存下所有元素），那得用arpack。
对于不需要全部本征值的，arpack和lapack都能适用。
如果需要全部本征值，那得用 ...

我的是双精度复数，要成衣4的吧？
赞一下 回复此楼
everythinghasitsseason.enjoyyourlife.
4楼2012-10-09 17:47:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lurencyj

木虫 (著名写手)
进步一的帮助啊。。。

只能跟你说，arpack看似简陋，功能暴强，matlab里面的eigs函数都是用的arpack子程求解。
赞一下 回复此楼
很女子很弓虽大
5楼2012-10-09 19:33:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lurencyj

木虫 (著名写手)
引用回帖:
4楼: Originally posted by 田山东 at 2012-10-09 17:47:33
我的是双精度复数，要成衣4的吧？...

双精度复数，要乘以16. 现在的机器。
赞一下 回复此楼
很女子很弓虽大
6楼2012-10-09 19:34:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)
引用回帖:
5楼: Originally posted by lurencyj at 2012-10-09 19:33:07
进步一的帮助啊。。。

只能跟你说，arpack看似简陋，功能暴强，matlab里面的eigs函数都是用的arpack子程求解。

你用arpack解矩阵特征值和特征向量没？在Lunix下解的吗？有没有Lunix下安装arpack的方法？
赞一下 回复此楼
everythinghasitsseason.enjoyyourlife.
7楼2012-10-10 09:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lurencyj

木虫 (著名写手)
用过，apt-get install libarpack...
赞一下 回复此楼
很女子很弓虽大
8楼2012-10-10 13:07:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)
引用回帖:
8楼: Originally posted by lurencyj at 2012-10-10 13:07:03
用过，apt-get install libarpack...

今天找同学装上了arpack，parpack96.tar.gz和patch.tar.z，结果在编译arpack自带的一个程序的时候，提示： Can't open included file 'debug.h'，是怎么回事呢？
赞一下 回复此楼
everythinghasitsseason.enjoyyourlife.
9楼2012-10-10 14:07:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)
你还记得当时怎么装的不？我用的是lunix  gfortran编辑器。能否给提供点安装步骤啊？
赞一下 回复此楼
everythinghasitsseason.enjoyyourlife.
10楼2012-10-10 15:14:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 田山东 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了8-15 +11 archvillain 2026-08-15 20/1000 2026-08-15 17:47 by 医学老男孩
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 6/300 2026-08-15 17:34 by 7huP8NPKZCp9
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 15:10 by QXtWNz7PI7wZ
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +7 majunge000 2026-08-11 9/450 2026-08-15 13:34 by kiz2011
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
信息提示
关闭
请填处理意见
关闭