当前时区为 UTC + 8 小时



发表新帖 回复这个主题  [ 6 篇帖子 ] 
作者 内容
1 楼 
 文章标题 : [问题]新得利软件包里面的libarpack如何调用呢?(已解决)
帖子发表于 : 2008-01-11 10:32 

注册: 2006-12-15 18:17
帖子: 33
送出感谢: 0 次
接收感谢: 0 次
因为要计算矩阵特征值,要调用arpack,安装了源里面的libarpack,却不知道怎么调用,请教达人

参考中科大bbs上一个帖子把问题解决了
1、在新得利里面搜索arpack,下载后安装
2、在终端输入whereis libarpack
我的机器上得的结果:
libarpack: /usr/lib/libarpack.a /usr/lib/libarpack.so
用第一个路径
3、去http://www.caam.rice.edu/software/ARPACK/ 下载arpack96.tar.gz
4、终端中输入tar xzvf arpack96.tar.gz
5、cd ARPACK/EXAMPLES/NONSYM
6、输入gfortran dndrv1.f -L/usr/lib -larpack -o ss
举例:如果程序aaa.f要调用/dir/lib/lib***.a (我的机器上路径是/usr/lib/libarpack.a)
编译时写作 gfortran aaa.f -L/dir/lib -l*** -o -ss
7、执行./ss
8、得到结果
Ritz values (Real,Imag) and relative residuals
----------------------------------------------
Col 1 Col 2 Col 3
Row 1: 1.96054E+01 0.00000E+00 1.94458E-14
Row 2: 7.68333E+01 0.00000E+00 4.23099E-15
Row 3: 4.82193E+01 0.00000E+00 7.80421E-15
Row 4: 4.82193E+01 0.00000E+00 6.98407E-15


_NDRV1
======

Size of the matrix is 100
The number of Ritz values requested is 4
The number of Arnoldi vectors generated (NCV) is 20
What portion of the spectrum: SM
The number of converged Ritz values is 4
The number of Implicit Arnoldi update iterations taken is 9
The number of OP*x is 139
*********************************************************************************
说明调用成功了
附注:1、网上有关于安装arpack的帖子,不过比较复杂,不推荐使用,源里面有现成的可以用
2、其他软件包没有尝试过,希望拍砖


页首
 用户资料  
 
2 楼 
 文章标题 :
帖子发表于 : 2008-04-06 18:25 

注册: 2006-10-18 15:48
帖子: 700
送出感谢: 0 次
接收感谢: 0 次
连arpack你都不自己编译,
估计是因为实验室没有师兄师姐指点你吧

:P


页首
 用户资料  
 
3 楼 
 文章标题 :
帖子发表于 : 2008-05-30 9:02 

注册: 2006-12-15 18:17
帖子: 33
送出感谢: 0 次
接收感谢: 0 次
之前实验室里面没有人用这个东西,都是现成的有限元软件计算,这回自己写程序又要算流固耦合,所以想用arpack


页首
 用户资料  
 
4 楼 
 文章标题 :
帖子发表于 : 2008-05-30 9:20 

注册: 2006-10-18 15:48
帖子: 700
送出感谢: 0 次
接收感谢: 0 次
自己编译arpack吧
arpack的效率依赖于blas和lapack
debian估计也是atlus_blas这个库
效率不高
你安装intel mkl库(包含高效的blas和lapack)
再在编译arpack时连接mkl
计算效率提高N倍.

编译arpack不难
关键是修改arpack96.tar.gz --> ARPACK/ARMAKES/下的配置文件。


页首
 用户资料  
 
5 楼 
 文章标题 :
帖子发表于 : 2008-05-30 10:22 

注册: 2006-12-15 18:17
帖子: 33
送出感谢: 0 次
接收感谢: 0 次
尝试过自己编译,配置文件修改无从下手
目前例子调通了,可是改驱动还是头疼的


页首
 用户资料  
 
6 楼 
 文章标题 :
帖子发表于 : 2008-05-30 16:22 

注册: 2006-10-18 15:48
帖子: 700
送出感谢: 0 次
接收感谢: 0 次
以前有没有在linux下编程的经验?
如果没有的话
周围又没人教你
那估计是很不爽
因为得附带着学习很多别的东西

不过如果这个arpack在你的程序中占有极大的重要性
那还是静下心来好好学习啊!


页首
 用户资料  
 
显示帖子 :  排序  
发表新帖 回复这个主题  [ 6 篇帖子 ] 

当前时区为 UTC + 8 小时


在线用户

正在浏览此版面的用户:没有注册用户 和 3 位游客


不能 在这个版面发表主题
不能 在这个版面回复主题
不能 在这个版面编辑帖子
不能 在这个版面删除帖子
不能 在这个版面提交附件

前往 :  
本站点为公益性站点,用于推广开源自由软件,由 DiaHosting VPSBudgetVM VPS 提供服务。
我们认为:软件应可免费取得,软件工具在各种语言环境下皆可使用,且不会有任何功能上的差异;
人们应有定制和修改软件的自由,且方式不受限制,只要他们自认为合适。

Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group
简体中文语系由 王笑宇 翻译