- 1、本文档共3页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
CLapack 在vs2008下的编译和使用说明
CLapack 在vs2008 下的静态编译和使用说明
1 CLapack 源代码下载
Lapack 全称为Linear Algebra Package,以BLAS (Basic Linear Algebra Subprograms)为基础,
用Fortran77 语言编写,可用于计算诸如求解线性代数方程、线性系统方程组的最小平方解、计
算特征值和特征向量等问题。是美国阿贡国家实验室、Kentuchy 大学与UC Davis 在1992 年共
同开发,为跨平台软件,支持Unix、Linux、OS2、Windows 等系统。当前最新版本为3.2.1,
其主页为/lapack/
CLapack 为Lapack 的C 语言接口,其主页为/clapack/
2. VS2008 下编译CLapack 3.1.1
文件的下载网址:/clapack/CLAPACK-3.1.1-VisualStudio.zip
也 可 以 使 用 CMake 编 译 库 文 件 , 方 法 请 见
/view/d286adeb172ded630b1cb644.html
这里我们直接使用VS2008 静态编译CLAPACK-3.1.1-VisualStudio 。
第一步,删除lib 目录下的所有内容,
第二步,使用vs2008 打开clapack.vc proj 文件,
选择 debug 或者 release 版进行编译,编译前应设置每个项目的属性-配置属性-c/c++-代
码生成-运行时库,根据的debug、release 选择/MTd 、/ MDd /MT 或者/MD。
第三步,依次编译libf2c 子项目,tmglib 子项目,blas 子项目和CLAPACK 子项目。
会生成四个lib 文件:libf2cd.lib,BLASd.lib,tmglibd.lib ,clapackd.lib 。
第四步,设定VC 环境
在工具-选项-VC++ 目录-包含文件或者库文件分别增加相应目录
本文本别为:
C:\clapack-3.1.1-VisualStudio\include 和C:\clapack-3.1.1-VisualStudio\lib\debug
3.在项目中使用CLAPACK
在MFC 下新建一个win32 控制台应用程序,
第一步,打开项目的属性-配置属性-c/c++-代码生成-运行时库,根据的debug 或者release
选择/MTd 或者MDd。这里与编译lib 的环境一致。
第二步,将libf2cd.lib,BLASd.lib,tmglibd.lib ,clapackd.lib 加入属性-配置属性-链接器-输
入-附加依赖项.
第三步,在项目的头文件中加入
externC
{
#include f2c.h
#include clapack.h
}
4.一个实例
广义特征值和特征向量求解。
使用函数dsygv_ ()
特征值放在*w 中,特征向量放在*a 中,其中a 为重复使用的数据结构,再输入时为a*x= (lamda)*b*x 中
的特征矩阵a,输出时为特征向量。
#include stdafx.h
#include iostream
#include string
using namespace std;
extern C{
#include f2c.h
#include blaswrap.h
#include clapack.h
};
int _tmain (int argc, _TCHAR* argv[])
{
integer itype;
itype=1;
char jobz;
jobz=V;
char uplo ;
uplo=U;
integer n;
n=3;
doublereal a[]={2,-1,-1,-1,4,-3,-1,-3,4};
integer lda;
lda=3;
doublereal b []={1,0,0,0,2,0,0,0,3};
integer ldb;
ldb=3;
doublereal *w ;
w=new doublereal[3];
doublereal *work ;
work=new doublereal[1];
work [0]=0;
//work=new doublereal[9];
//for(
文档评论(0)