你好,游客 登录 注册 搜索
背景:
阅读新闻

windows下GSL的配置 - Augusdi的专栏

[日期:2013-04-14] 来源:  作者: [字体: ]

一、GSL介绍

GNU科学计算函数库GSL(GNU Scientific Library)是一个强大的C/C++数值计算函数库,它是一个自由软件,是GNU项目软件的一个部分,遵循GPL协议。GSL是一个为C和C++程序员提供的科学数值运算库。该科学计算库异常强大,函数库提供了大量的数值计算程序

 

,如随机函数、特殊函数和拟合函数等等,整个函数库大约有1000多个函数,几乎涵盖了科学计算的各个方面。提供了如下方面的支持:

Complex Numbers          Roots of Polynomials           Special Functions

Vectors and Matrices        Permutations                  Sorting

BLAS Support             Linear Algebra                 Eigensystems

Fast Fourier Transforms      Quadrature                    Random Numbers

Quasi-Random Sequences    Random Distributions           Statistics

Histograms                N-Tuples                     Monte Carlo Integration

Simulated Annealing        Differential Equations           Interpolation

Numerical Differentiation    Chebyshev Approximation       Series Acceleration

Discrete Hankel Transforms   Root-Finding                 Minimization

Least-Squares Fitting        Physical Constants             IEEE Floating-Point

Discrete Wavelet Transforms                              Basis splines

 

该函数库的主页是:http://www.gnu.org/software/gsl/gsl.html。不过遗憾的是原始GSL并不支持不支持windows平台,可所幸的是有人做了GSL在windows上的移植工作,详见http://gnuwin32.sourceforge.net/packages/gsl.htm,目前版本是1.8。

二、下载gsl

1、从http://gnuwin32.sourceforge.net/packages/gsl.htm下载Complete package, except sources和Sources两个exe文件。

 

三、安装

1、 首先安装从http://gnuwin32.sourceforge.net/packages/gsl.htm下载的两个文件gsl-1.8.exe和gsl-1.8-src.exe。

 

四、设置Visual C++ 6.0编译环境

 1、生成lib文件。发现安装目录lib下并没有libgsl.lib,libgslcblas.lib这两个文件,倒是有两个扩展名为def和a(linux下库文件包格式)的文件,因此必须进行转换。

l          开始菜单,点击运行,输入cmd。

l          进入gsl库的lib目录下依次输入以下两条语句:
    lib /machine:i386 /def:libgsl.def
    lib /machine:i386 /def:libgslcblas.def

再看lib目录下,发现有了libgsl.lib,libgslcblas.lib这两个文件。

 2、将x:/Program Files/GnuWin32l/bin中的libgsl.dll和libgslcblas.dll复制到x:/VC98/Bin;将/include整个Gsl目录复制到x:/VC98/include下;/lib目录下的所有.lib文件全部复制到x:/VC98/Lib下。

3、 新建一个工程用于测试,。然后在该项目的project-settings-link,在object/library modules中加入你用到的库文件,如libgsl.lib libgslcblas.lib,用空格隔开。

五、测试Gsl函数库

#define GSL_DLL
#include <gsl/gsl_spline.h>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <gl/glut.h>
#include <gl/gl.h>
void Display()
{
    glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);
    
const size_t n = 4;
    double x[] = {0,0.333336,0.666666,1
};
    double y[] = {0,0.5,0.9,0
};
    gsl_interp* interps[3] = 
{NULL,NULL,NULL};
    interps[0
] = gsl_interp_alloc(gsl_interp_linear,n);
    interps[1
] = gsl_interp_alloc(gsl_interp_polynomial,n);
    interps[2
] = gsl_interp_alloc(gsl_interp_cspline,n);
    gsl_interp_init(interps[0
],x,y,n);
    gsl_interp_init(interps[1
],x,y,n);
    gsl_interp_init(interps[2
],x,y,n);
    gsl_interp_accel* acc = gsl_interp_accel_alloc();
    glBegin(GL_LINE_STRIP);
        
for(double t=0.0; t<=1.025; t+=0.025)
        
{
            glColor3f(1,0,0
);
            glVertex3f(t,gsl_interp_eval(interps[0],x,y,t,acc),0.0
);
        }
    glEnd();
    glBegin(GL_LINE_STRIP);
        
for(double t=0.0; t<=1.025; t+=0.025)
        
{
            glColor3f(0,1,0
);
            glVertex3f(t,gsl_interp_eval(interps[1],x,y,t,acc),0.0
);
        }
    glEnd();
    glBegin(GL_LINE_STRIP);
        
for(double t=0.0; t<=1.025; t+=0.025)
        
{
            glColor3f(0,0,1
);
            glVertex3f(t,gsl_interp_eval(interps[2],x,y,t,acc),0.0
);
        }
    glEnd();
    gsl_interp_accel_free(acc);
    gsl_interp_free(interps[
0]);
    gsl_interp_free(interps[1
]);
    gsl_interp_free(interps[2
]);
    glutSwapBuffers();
}
int main(int argc, char** argv)
{
    glutInit(&argc,argv);
    glutInitDisplayMode(GLUT_DOUBLE|GLUT_RGBA);
    glutInitWindowSize(
512,512);
    glutCreateWindow("GSL Interpolation"
);
    glutDisplayFunc(Display);
    glClearColor(
1,1,1,1);
    glMatrixMode(GL_MODELVIEW);
    glLoadIdentity();
    glTranslatef(-
1,-1,0);
    glScalef(2,2,1
);
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
    glViewport(
0,0,512,512);
    glLineWidth(4.0
);
    glutMainLoop();
    
return 0;
}运行结果:

 

说明Gsl函数库已经可以使用了。





收藏 推荐 打印 | 录入:admin | 阅读:
相关新闻