GSL的Linux安装和使用教程（小白教程）

TAICHIFEI

已于 2024-01-19 12:24:42 修改

阅读量7.2k

点赞数 4

文章标签：人工智能 c++

于 2023-07-13 17:55:54 首次发布

本文链接：https://blog.csdn.net/weixin_42035282/article/details/131708094

版权

————————————————————————————————

官网：GSL - GNU Scientific Library - GNU Project - Free Software Foundation
教程文档：GNU Scientific Library — GSL 2.7 documentation
下载文档：https://www.gnu.org/software/gsl/doc/latex/gsl-ref.pdf
源码：https://github.com/ampl/gsl

————————————————————————————————

一、GSL简介

GSL （GNU Scientific Library）是面向 C 和 C++ 程序员的数值库。专门为应用数学和科学技术领域的数值计算提供支持的软件库。

GSL 使用C语言编写，同时也为其他语言做了相应的封装。该库易于编译，并且不依赖于其他包。

GSL 在GNU通用公共许可下是免费的。该函数库提供了广泛的数学算法的实现函数，包括随机数生成器，特殊函数和最小二乘拟合等等。目前该函数库提供有超过1000个函数，这些函数包含的范围有：复数计算、多项式求根、特殊函数、向量和矩阵运算、排列、组合、排序、线性代数、特征值和特征向量、快速傅里叶变换(FFT)、数值积分、随机数生成、随机数分布、统计、蒙特卡洛积分、模拟退火、常微分方程组、插值、数值微分、方程求根、最小二乘拟合、小波变换等。

二、GSL下载

GSL是在 GNU/Linux 系统上开发的，可以使用 gcc 命令来查看当前系统上安装的 GCC 版本。

打开终端，输入以下命令：

gcc --version

最新版下载链接：Index of /gnu/gsl

gsl-2.7.1.tar.gz

2021-12-01 00:33

7.2M

cd 下载目录

1. 解压压缩包：$ tar -zxvf gsl-2.7.1.tar.gz

2. 创建安装目录(自定义):

mkdir /home/tx-deepocean/Library

cd /home/tx-deepocean/Library

# 创建安装目录gsl_gcc

mkdir gsl_gcc

3. 进入源码目录:$ cd 下载目录/gsl-2.7.1

4. 配置configure文件：执行

$ ./configure CC="gcc" --prefix=/home/tx-deepocean/Library/gsl_gcc

# 参数含义
# ./configure -h: 可以得到配置的帮助信息
# CC:指定编译工具链
# --prefix: 指定库的安装目录
# --host: 指定安装平台,可选，交叉编译要制定安装平台如在arm上编译 --host=arm

5.使用make进行编译

$ make

6.使用make install进行安装

$ make install

7.使用make installcheck检查安装是否成功

$ make installcheck

8. 编辑.bashrc文件
$ vi ~/.bashrc

在文件末尾添加如下代码：

export C_INCLUDE_PATH=$C_INCLUDE_PATH:/home/tx-deepocean/Library/gsl_gcc/include
export CPLUS_INCLUDE_PATH=$CPLUS_INCLUDE_PATH:/home/tx-deepocean/Library/gsl_gcc/include
export LD_LIBRARY_PATH=$LD_LIBRARY_PATH::/home/tx-deepocean/Library/gsl_gcc/lib
export LIBRARY_PATH=$LIBRARY_PATH::/home/tx-deepocean/Library/gsl_gcc/lib

保存该文件，立刻加载修改后的设置，使之生效：source ~/.bashrc

三、GSL使用

1.新建工程文件夹

mkdir /home/tx-deepocean/Library/gsl_test

2.新建makefile文件并编辑

#进入gsl_test文件夹并新建makefile文件
$ touch makefile
# 编辑makefile文件
$ vim makefile

makefile文件内容（如果出现错误注意tab的缩进，gcc前要求一个tab的缩进，最好自己手敲一遍，直接粘贴可能出问题）

IPATH=-I/home/ccchenji/Desktop/gsl_gcc/include
LPATH=-L/home/ccchenji/Desktop/gsl_gcc/lib
#这里直接链接的静态库文件
main:main.c
	gcc main.c $(IPATH)  $(LPATH) -l:libgsl.a -l:libgslcblas.a -o main

.PHONY:clean
clean:
	rm -rf main

3.新建main.c文件并编辑

#进入gsl_test文件夹并新建main.c文件
$ touch main.c
# 编辑makefile文件
$ vim main.c

#include <stdio.h>
#include <gsl/gsl_matrix.h>
int main (void)
{
        int i, j;
        gsl_matrix * m = gsl_matrix_alloc (5, 3);
        for (i = 0; i < 5; i++)
                for (j = 0; j < 3; j++)
                        gsl_matrix_set (m, i, j, 100*i + j);
        for (i = 0; i < 5; i++)
                for (j = 0; j < 3; j++)
                        printf ("m(%d,%d) = %g\n", i, j,
                                        gsl_matrix_get (m, i, j));
        gsl_matrix_free (m);
        return 0;
}

4.使用make命令编译main.c文件

$ make

#编译后生成main文件

#运行main文件

$ ./main

运行结果：

四、GSL常用函数--矩阵运算

1. 加法、减法

#include <stdio.h>
#include <gsl/gsl_matrix.h>

void print_matrix(const gsl_matrix *m)
{
    int i, j;
    for (i = 0; i < m->size1; i++)
    {
        for (j = 0; j < m->size2; j++)
        {
            printf("%g ", gsl_matrix_get(m, i, j));
        }
        printf("\n");
    }
    printf("\n");
}

int main(void)
{
    gsl_matrix *a = gsl_matrix_alloc(2, 2);
    gsl_matrix *b = gsl_matrix_alloc(2, 2);

    gsl_matrix_set(a, 0, 0, 1.0);
    gsl_matrix_set(a, 0, 1, 2.0);
    gsl_matrix_set(a, 1, 0, 3.0);
    gsl_matrix_set(a, 1, 1, 4.0);

    gsl_matrix_set(b, 0, 0, 0.5);
    gsl_matrix_set(b, 0, 1, 1.0);
    gsl_matrix_set(b, 1, 0, 1.5);
    gsl_matrix_set(b, 1, 1, 2.0);

    printf("Matrix a:\n");
    print_matrix(a);

    printf("Matrix b:\n");
    print_matrix(b);

    gsl_matrix_add(a, b);
    printf("Matrix a after addition:\n");
    print_matrix(a);

    gsl_matrix_sub(a, b);
    printf("Matrix a after subtraction:\n");
    print_matrix(a);

    gsl_matrix_free(a);
    gsl_matrix_free(b);

    return 0;
}

运行结果：

2.矩阵乘法

要使用用gcc编译gsl_blas_dgemm 函数必须使用动态库，静态库链接会报错(交叉编译工具链默认使用静态库，对交叉编译工具链，可以使用静态库)

2.1.配置动态链接库路径

#编辑.bashrc文件
$ vi ~/.bashrc
#在文件末尾加入以下代码，以下代码的作用是将动态连接库加入到gcc默认搜索路径
export LD_LIBRARY_PATH=/home/tx-deepocean/Library/gsl_gcc/lib:$LD_LIBRARY_PATH

2.2.重写makefile文件

#注意：ccchenji是我系统的用户名，这里要换成你的系统的用户名。
IPATH=-I/home/tx-deepocean/Library/gsl_gcc/include
LPATH=-L/home/tx-deepocean/Library/gsl_gcc/lib
#这里直接链接的静态库文件
main:main.c
	gcc main.c $(IPATH)  $(LPATH) -l:libgsl.a -l:libgslcblas.a -o main

.PHONY:clean
clean:
rm -rf main

3.gsl_blas_dgemm()函数说明

/*
参数说明
1.CBLAS_TRANSPOSE_t 可取值: 
  CblasNoTrans：矩阵不需要转置
  CblasTrans: 矩阵需要转置
2.alpha、beta数乘系数
3.计算表达式： C = alpha * op( A ) * op( B ) + beta * C
*/

//函数原型
int gsl_blas_dgemm( CBLAS_TRANSPOSE_t TransA, 
					CBLAS_TRANSPOSE_t TransB, 
					double alpha,
					const gsl_matrix *A,
					const gsl_matrix *B, 
					double beta, 
					gsl_matrix *C)

4.gsl_blas_dgemm()函数使用例程

#include <stdio.h>
#include <gsl/gsl_blas.h>
#include <gsl/gsl_matrix.h>

int main(void)
{
    double a[] = {0.11, 0.12, 0.13, 0.21, 0.22, 0.23};
    double b[] = {1011, 1012, 1021, 1022, 1031, 1032};
    double c[] = {0.00, 0.00, 0.00, 0.00};

    gsl_matrix_view A = gsl_matrix_view_array(a, 2, 3);
    gsl_matrix_view B = gsl_matrix_view_array(b, 3, 2);

    gsl_matrix *C = gsl_matrix_alloc(2, 2);

    /* Compute C = A B */
    gsl_blas_dgemm(CblasNoTrans, CblasNoTrans,
                   1.0, &A.matrix, &B.matrix,
                   0.0, C);

    printf("[ %g, %g\n", gsl_matrix_get(C, 0, 0), gsl_matrix_get(C, 0, 1));
    printf("  %g, %g ]\n", gsl_matrix_get(C, 1, 0), gsl_matrix_get(C, 1, 1));

    gsl_matrix_free(C);

    return 0;
}

运行结果：

3.矩阵求逆

#include <stdio.h>
#include <gsl/gsl_linalg.h>

int main (void)
{
        double a_data[] = { 1.0, 0.6, 0.0,
                0.0, 1.5, 1.0,
                0.0, 1.0, 1.0 };
        /*
         * Inverse is
         *    1  -1.2   1.2
         *    0   2.0  -2.0
         *    0  -2.0   3.0
         */
        double inva[9];
        int s, i, j;

        gsl_matrix_view m
                = gsl_matrix_view_array(a_data, 3, 3);
        gsl_matrix_view inv
                = gsl_matrix_view_array(inva,3,3);
        gsl_permutation * p = gsl_permutation_alloc (3);

        printf("The matrix is\n");
        for (i = 0; i < 3; ++i)
                for (j = 0; j < 3; ++j)
                        printf(j==2?"%6.3f\n":"%6.3f ", gsl_matrix_get(&m.matrix,i,j));

        gsl_linalg_LU_decomp (&m.matrix, p, &s);    
        gsl_linalg_LU_invert (&m.matrix, p, &inv.matrix);

        printf("The inverse is\n");
        for (i = 0; i < 3; ++i)
                for (j = 0; j < 3; ++j)
                        printf(j==2?"%6.3f\n":"%6.3f ",gsl_matrix_get(&inv.matrix,i,j));
        gsl_permutation_free (p);
        return 0;
}