博客访问： 171215
博文数量： 66
博客积分： 1410
博客等级： 上尉
技术积分： 730
注册时间：
APP发帖 享双倍积分
IT168企业级官微
微信号：IT168qiye
系统架构师大会
微信号：SACC2013
简单的关于这两个库的介绍*****************************************************blas: Basic Linear Algebra SubprogramsThe BLAS (Basic Linear Algebra Subprograms) are high quality "building block" routines for performing basic vector and matrix operations. Level 1 BLAS do vector-vector operations, Level 2 BLAS do matrix-vector operations, and Level 3 BLAS do matrix-matrix operations. Because the BLAS are efficient, portable, and widely available, they're commonly used in the development of high quality linear algebra software, LINPACK and LAPACK for example.A Fortran77 reference implementation of the BLAS is located in the blas directory of Netlib.blas 是许多数值计算软件库的核心, 一般是用 Fortran77 实现的, 但也有一些其它的包装, 如 cblas 就是 C 语言, 也有 C++ 的包装, boost/ublas 是 C++ template class 的实现; 另外还有一些特别的实现, 如 intel MKL, AMD core math libraryblas就是做向量、矩阵的基本运算，如加、减、乘等操作 ********************************************************lapack: Linear Algebra PACKageLAPACK is written in Fortran77 and provides routines for solving systems of simultaneous linear equations, least-squares solutions of linear systems of equations, eigenvalue problems, and singular value problems. The associated matrix factorizations (LU, Cholesky, QR, SVD, Schur, generalized Schur) are also provided, as are related computations such as reordering of the Schur factorizations and estimating condition numbers. Dense and banded matrices are handled, but not general sparse matrices. In all areas, similar functionality is provided for real and complex matrices, in both single and double precision.lapack 用 blas 做底层运算许多高层的数学库都用 blas 和 lapack 做底层, 如 cvmlib, MTL 等等lapack是做线性运算的，最常用的就是求解线性方程组，AX=b, X=A, 还有如 矩阵分解、求逆，求矩阵特征值、奇异值等 *********************************************************
阅读(5005) | 评论(1) | 转发(0) |
相关热门文章
给主人留下些什么吧！~~
你能不能把你编译lapack的过程写写啊，具体的什么环境之类，谢谢啊
或者发到我的邮箱chenxuelong_
请登录后评论。LAPACK函数介绍-1
LAPACK函数介绍_工学_高等教育_教育专区。今日推荐 88份文档 2014全国高考状元联手分享状元笔记 衡水中学文科学霸高中数学笔记 清华附中文科学霸高中政治笔记 东北师大...
blas函数介绍_计算机软件及应用_IT/计算机_专业资料。blas: Basic Linear Algebra...lapack 用 blas 做底层运算 许多高层的数学库都用 blas 和 lapack 做底层, ...
cpplapack 由于 BLAS 和 LAPACK 的每个函数的参数都很多,传递起来比较麻烦,有很多包都实现了 对 LAPACK 和 BLAS 的封装,cpplapack 是其中一个,编译的时候需要 ...
Matlab 中的 很多数希疏矩阵函数的原型都可以在其中找到。性能非常出色。 ? TNT(Template Numerical Toolkit) C++的数值计算库,目标很大,据说还要集成 LAPACK++,...
数学库_理学_高等教育_教育专区。值得推荐的 数学库 相关函数1 ACML(AMD Core...基本线性代数库) 3 GOTO 和 ATLAS 库 4 LAPACK Linear Algebra Package, 线性...
google 与 Pfizer 的员工也介绍了 R 在自己公司中...另一方面,它必须具有良好的可 扩展性,即它...以下只写一个简单的 R 函数作为 说明,在 R 目录...
s_n=dblquad(@(x,y)x.^y,0,1,1,2) %匿名函数表示 s_n=dblquad('...(saxpy 是 LAPACK 术语) %---用二重循环体实现矩阵乘法---C2=zeros(size(A...
(saxpy 是 LAPACK 术语) %---用二重循环体实现...非奇异时, 线性代数教科书常介绍的线性方程 Ax ? ...%匿名函数写成的 ode45 所需得导数函数 [tt,yy]...
然后设置需要编译的库函数如下 #lib: lapacklib tmglib lib: blaslib lapacklib tmglib 默认的情况是不编译 blas 库的。(除非事先编译好了 blas 库,否则应一...
关键词:BLASLAPACK 1/2 相关文档推荐 突出应用背景知识介绍 彰显... 暂无评价...IMSL:BLAS库函数简表,Fo... 5页 免费 2009年7月高等教育自学考试... 3页...君，已阅读到文档的结尾了呢~~
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
LAPACK函数介绍不完全.doc
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口成功使用MKL中的LAPACK95
昨天和今天在折腾intel composer XE 2013的数学核函数库MKL里面的LAPACK
以及LAPACK95。终于弄通了，我猜想其他的库，比如BLAS， BLAS95, FFTW都是类似。
以求矩阵的逆为例。（虽然在Intel& Math Kernel
Library Reference Manual（Document Number: US, MKL
11.2）中已经明确提醒: 在求解Ax=b时，不要直接以x=inv(A）*b
这样求解，以提高稳定性和精度。但是还是有极少情况下需要对大型矩阵的求逆，比如在对随机变量进行非参数估计其均值的时候，就需要求自协方差矩阵的逆。）
========================================
第一个例子： 测试LAPACK
========================================
LAPACK, 也就是LAPACK 的 FORTRAN
77版本。在源代码中不需要添加任何use语句，直接调用?gettf() +
?getri()就可了。就像调用?abs()等内部函数一样，？代表传递参数的类型,包括s, d,
c和z，它们分别代表单精度实型，双精度实型，单精度复实型，双精度复实型。
相比于LAPACK95， LAPACK
FORTRAN77版本的函数的形参个数比较多，调用还得区分清楚实数类型。
方法一：使用command
prompt下的ifort编译命令：&
test_lapack77.f90&
使用ifort编译时，必须加上/Qmkl选项，以告诉编译器要使用mkl库。
方法二：在Microsoft Visual studio * 的library属性里把use
Intel Math Kernel Library
加上就可以了，在其下拉菜单中随便选一个就行，除了NO。
*******************************************************************
& &Fortran 代码
*******************************************************************
**************************************************************************
! Build this test program by Intel(R) Fortran Compiler in
Window7 OS :
! & & ifort
test_lapack77.f90 /Qmkl
--------------------------------------------------------------------------
! Reference:
! Intel? Math Kernel Library Reference Manual, Document
Number: US, MKL 11.2
! &pages: 411~419, 596~598
!**************************************************************************
& & PROGRAM
test_lapack77
& & IMPLICIT NONE
& & INTEGER,PARAMETER :: N =
& & REAL,DIMENSION(N,N) ::
toeplitz, a
& & integer,dimension(N) ::
& & real,dimension(N) ::
& & INTEGER :: i, j,
-------------------------------
toeplitz=reshape((/3,-1,-1,4,-2,-1,-3,2,1/),(/N,N/))
toeplitz=transpose(toeplitz)
& & a = toeplitz
-------------------------------
write(*,'(/,2x,a)')'a='
write(*,10)((a(i,j),j=1,N),i=1,N)
&10 format(3x,f6.2)
-------------------------------
call sgetrf(N,N,a,N,ipiv,info)
call sgetri(N,a,N,ipiv,work,N,info)
-------------------------------
write(*,'(/,2x,a)')'Inv(a)='
write(*,10)((a(i,j),j=1,N),i=1,N)
-------------------------------
& & write(*,*)'The pivot
& & write(*,*)ipiv
write(*,*)'workspace:'
& & write(*,*)work
& & write(*,*)'CALL
& & write(*,*)info
& & END PROGRAM
test_lapack77
========================================
第二个例子： 测试LAPACK95
========================================
相比于LAPACK
FORTRAN77版本， LAPACK95的函数的形参个数非常简洁，调用时不需要区分清楚实数类型。因为LAPACK95已经用pure
procedure(纯过程)把不同实数类型的函数封装在一个interface里面了。
方法一：&使用command
prompt下的ifort编译命令：&
& ifort&/Qmkl&mkl_lapack95_lp64.lib
test_lapack95.f90
使用ifort编译时，除了要用到/Qmkl&，还需要用到mkl_lapack95_lp64.lib文件，它可以在intel的安装目录下找到，比如我的D:\Program
Files (x86)\Intel\Composer XE
2013\mkl\lib\intel64。把它拷贝到与test_lapack95.f90同一个文件夹下。
另外，在源代码中要使用&“
F95_PRECISION 和 &USE&LAPACK95 ”或者
& “USE&MKL95_PRECISION 和
& &USE&MKL95_LAPACK”，二者选一种即可。
& & 除了在Microsoft Visual studio * 的library属性里把use Intel
Math Kernel Library 加上，还要在其Linker
-& input中的Additional Dependencies加上mkl_lapack95_lp64.lib
*******************************************************************
&Fortran 代码
*******************************************************************
!**************************************************************************
! Build this test program by Intel(R) Fortran Compiler in
Window7 OS :&
! & &ifort /Qmkl
mkl_lapack95_lp64.lib test_lapack95.f90 &
!**************************************************************************
& & PROGRAM
test_lapack95
! & &USE F95_PRECISION
! & &USE LAPACK95
& & USE MKL95_PRECISION
& & USE MKL95_LAPACK
& & IMPLICIT NONE
& & INTEGER,PARAMETER :: N =
& & REAL,DIMENSION(N,N) ::
toeplitz, a
& & integer,dimension(N) ::
& & real,dimension(N) ::
& & INTEGER :: i, j,
-------------------------------
toeplitz=reshape((/3,-1,-1,4,-2,-1,-3,2,1/),(/N,N/))
toeplitz=transpose(toeplitz)
& & a = toeplitz
-------------------------------
write(*,'(/,2x,a)')'a='
write(*,10)((a(i,j),j=1,N),i=1,N)
&10 format(3x,f6.2)
-------------------------------
LAPACK F95
-------------------------------
& & call getrf(a,ipiv,info)
& & call getri(a,ipiv,info)
-------------------------------
write(*,'(/,2x,a)')'Inv(a)='
write(*,10)((a(i,j),j=1,N),i=1,N)
& & END PROGRAM
test_lapack95
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。人人网 - 抱歉
哦，抱歉，好像看不到了
现在你可以:
看看其它好友写了什么
北京千橡网景科技发展有限公司：
文网文[号··京公网安备号·甲测资字
文化部监督电子邮箱：wlwh@··
文明办网文明上网举报电话： 举报邮箱：&&&&&&&&&&&&