24小时热门版块排行榜    

查看: 4498  |  回复: 7

uuv2010

荣誉版主 (职业作家)

优秀版主

[求助] fortran程序怎样调用lapack库进行矩阵对角化?

在linux下想应用lapack库来进行矩阵的对角化,在fortran的程序中怎样实现lapack库的调用?非常感谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jjdg(金币+1): 感谢参与 2011-11-12 00:41:36
ben_ladeng: 专家考核存档 2011-11-13 11:12:45
余泽成(金币+10): 替楼主发放奖励! 2011-11-14 11:17:32
引用回帖:
1楼: Originally posted by uuv2010 at 2011-11-11 22:29:52:
在linux下想应用lapack库来进行矩阵的对角化,在fortran的程序中怎样实现lapack库的调用?非常感谢!

首先你需要确定你安装lapack函数库了吧。假设你的lapck安装在目录home/dicp/lib中,这个目录下会有静态库lapack.a,或者liblapack.a libblas.a
假设你的代码有三个a.f90 b.f90 c.f90
调用方式有两种,
写个简单的makefile 做个例子:
LIB = -L home/dicp/lib
a.out : a.o b.o c.o
        ifort -o a.out a.o b.o c.o  $(LIB) lipack.a blas.a
另外一种是如下:
  LIB = -L home/dicp/lib
LIBB = -llapack -lblas
   a.out : a.o b.o c.o
        ifort -o a.out a.o b.o c.o  $(LIB) l$(LIBB)
上面两种是我常用的调用方式。
2楼2011-11-11 22:45:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jjdg(金币+1): 感谢参与 2011-11-12 00:41:44
余泽成(金币+10): 替楼主发放奖励! 2011-11-14 11:17:44
引用回帖:
2楼: Originally posted by maomao1210 at 2011-11-11 22:45:59:
首先你需要确定你安装lapack函数库了吧。假设你的lapck安装在目录home/dicp/lib中,这个目录下会有静态库lapack.a,或者liblapack.a libblas.a
假设你的代码有三个a.f90 b.f90 c.f90
调用方式有两种,  ...

首先你需要确定你安装lapack函数库了吧。假设你的lapck安装在目录home/dicp/lib中,这个目录下会有静态库lapack.a,blas.a tmglib.a 或者liblapack.a libblas.a libtmg.a
假设你的代码有三个a.f90 b.f90 c.f90
调用方式有两种,
写个简单的makefile 做个例子:
LIB = -L home/dicp/lib
a.out : a.o b.o c.o
        ifort -o a.out a.o b.o c.o  $(LIB) lipack.a blas.a
第一种意思就是 : 把静态库,lapack.a 直接和你写的程序一起编译链接,生成可执行程序a.out.

另外一种是如下:
  LIB = -L home/dicp/lib
LIBB = -llapack -lblas
   a.out : a.o b.o c.o
        ifort -o a.out a.o b.o c.o  $(LIB) $(LIBB)
第二种, -l lapack make会识别你想调用lapack库,会去相应的目录寻找前缀为liblapack.a的库文件。

如果你写makefile 上面两种都可以利用,如果不用makefile的话,建议你利用第一种。
上面两种是我常用的调用方式。
3楼2011-11-11 22:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

引用回帖:
3楼: Originally posted by maomao1210 at 2011-11-11 22:51:49:
首先你需要确定你安装lapack函数库了吧。假设你的lapck安装在目录home/dicp/lib中,这个目录下会有静态库lapack.a,blas.a tmglib.a 或者liblapack.a libblas.a libtmg.a
假设你的代码有三个a.f90 b.f90  ...

请教一下:这个还需要知道对角化的程序在lapack中的名字和用法,以及参数的类型和数目吧?你的第一种形成的a.out是不是把lapack中的函数库全部给你打包到a.out中了?还是只把你要调用的函数给打包进去了?这个可执行文件在执行的时候就可以脱离lapack的环境了吧?第二种情况呢?是不是也是脱离lapack环境呢?
弘德明志博学笃行
4楼2011-11-12 09:33:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

引用回帖:
3楼: Originally posted by maomao1210 at 2011-11-11 22:51:49:
首先你需要确定你安装lapack函数库了吧。假设你的lapck安装在目录home/dicp/lib中,这个目录下会有静态库lapack.a,blas.a tmglib.a 或者liblapack.a libblas.a libtmg.a
假设你的代码有三个a.f90 b.f90  ...

非常感谢您的详细指导!
这样的话就不需要在fortran程序开头部分应用use 语句指定应用lapack中的哪个具体函数了吧?
5楼2011-11-12 10:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)


余泽成(金币+1): 谢谢参与应助! 2011-11-13 21:39:27
引用回帖:
5楼: Originally posted by uuv2010 at 2011-11-12 10:00:41:
非常感谢您的详细指导!
这样的话就不需要在fortran程序开头部分应用use 语句指定应用lapack中的哪个具体函数了吧?

恩,不需要的,当然你是需要知道需要调用哪个lapack函数才行,比如对角化的是dicp这个函数,在你的源代码里 call dicp(*,*,*,*,*) 就可以了,大致的写法,需要你查看具体的手册。
6楼2011-11-13 11:04:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)


余泽成(金币+1): 谢谢参与应助! 2011-11-13 21:39:38
引用回帖:
4楼: Originally posted by zhangguangping at 2011-11-12 09:33:23:
请教一下:这个还需要知道对角化的程序在lapack中的名字和用法,以及参数的类型和数目吧?你的第一种形成的a.out是不是把lapack中的函数库全部给你打包到a.out中了?还是只把你要调用的函数给打包进去了?这个可 ...

那是当然的,如果你想对角化矩阵的话,当然是需要知道具体哪个是对角化的函数的嘛,对吧? 只需要把你要调用的函数 call 在你的源代码里就可以了,当然,传递变量的类型数目都需要你手册的。lapack.a,blas.a 是把编译后的二进制.o文件打包成静态库了,把这个文件放在那里都是可以用的,建议你尝试就知道了。呵呵。
7楼2011-11-13 11:08:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

引用回帖:
7楼: Originally posted by maomao1210 at 2011-11-13 11:08:19:
那是当然的,如果你想对角化矩阵的话,当然是需要知道具体哪个是对角化的函数的嘛,对吧? 只需要把你要调用的函数 call 在你的源代码里就可以了,当然,传递变量的类型数目都需要你手册的。lapack.a,blas.a 是 ...

非常感谢!
弘德明志博学笃行
8楼2011-11-13 14:32:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 uuv2010 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +8 瞬息宇宙 2026-08-31 15/750 2026-09-02 11:13 by ma0526
[基金申请] 我不理解! +18 Edward_pc 2026-08-26 26/1300 2026-09-02 11:01 by shanzhipu
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +4 瞬息宇宙 2026-09-01 8/400 2026-09-02 10:35 by 瞬息宇宙
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
信息提示
请填处理意见