24小时热门版块排行榜    

查看: 4495  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 8 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +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 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见