24小时热门版块排行榜    

查看: 4476  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 号外号外 +8 GZUZM 2026-08-24 8/400 2026-08-24 19:22 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 人气不行了 +9 fansofjerry 2026-08-21 9/450 2026-08-24 17:48 by zhaozhengxin
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见