24小时热门版块排行榜    

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

uuv2010

荣誉版主 (职业作家)

优秀版主

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

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

» 猜你喜欢

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

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

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的回帖
信息提示
请填处理意见