24小时热门版块排行榜    

查看: 4474  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 范进中举一文的中心思想 +4 炎黄贵胄 2026-08-22 5/250 2026-08-24 02:09 by 流流伤
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +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 紫杉醇
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见