insectC的专栏

      由于实验室里课题的需要，寒假里开始学习GNU Scientific Library(GSL)。这个库是一系列科学数值计算子程序的集合。包含了线性代数、多项式求根、特殊函数、统计、数值微分与积分等等子程序，感觉相当实用。      当然，按我刨根问底的个性，自然要结合着源代码来学习。这个寒假要痛苦并快乐着了~

     我是在linux下使用gsl库的。如果要在windows之下使用，则可以参考下文：     http://www.foelin.com/2009/04/windows-%E7%8E%AF%E5%A2%83%E4%B8%8B%E4%BD%BF%E7%94%A8gsl.html     首先，我们可以从ftp://ftp.gnu.org/gnu/gsl/下载最新版本的gsl库，我下载

    作为一个数值计算库，gsl主要关注浮点数的运算。且完全继承了c语言的浮点数体系，即分为单精度float型（32-bit），双精度double型（64-bit）以及扩展精度浮点数long double型（80-bit）。    对应不同精度的需要，gsl中一些结构体、函数也有double,float与long double三套。其中，double型为这些结构体与函数的默认数据类型，而fl

    gsl中的向量与矩阵是基于一种底层的数据类型，即数据块（block）实现的。因此我们的分析，就从数据块开始。    与复数一样，gsl中也有各种不同数据类型的数据块，详细列表如下：gsl_blockdoublegsl_block_floatfloatgsl_block_long_doublelong doublegsl_b

       如果在程序中需要调用gsl中的矩阵类型，首先需要声明#include       与向量一样，gsl中的矩阵同样也构建在数据块之上。其声明如下：［gsl_matrix_double.h］......typedef struct { size_t size1;//矩阵的行数 size_t size2;//矩阵的列数 size_t

      IEEE754标准是现在最通用的二进制算术运算标准，它被用于Sun，DEC，HP，IBM工作站以及所有的个人计算机。IEEE算术运算包括两类浮点数：单精度（float型，32位字长）和双精度（double型，64位字长）。除此之外，还有一个可选的扩展精度格式(long double，80位字长）。      规范化的浮点数可以表示为                     

     在gsl中，由于有float、double、long double多种数据精度，因此如果针对每一种数据精度都分别实现诸如内存分配等函数，代码量将骤增数倍。因此各种数据类型共用这些代码就显得相当重要。gsl利用c语言的预编译命令实现了这种代码的重用功能。      这里我们以最为基础的gsl_block_alloc函数模板为例，剖析gsl代码重用的原理。     下面给出该函数的模