Augusdi的专栏

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (int argc, char **argv){ double a, b; size_t n; if (argc != 4) {

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ double data[] = { 1.0 , 1/2.0, 1/3.0, 1/4.0, 1/2

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")double f (double x, void * params) { double alpha = *(double *) params; double f = lo

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ int i, j; gsl_matrix * m = gsl_matrix_alloc (10, 3); for (i = 0; i < 10

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ int i, j, k = 0; gsl_matrix * m = gsl_matrix_alloc (100, 100); gsl_matrix

#include #include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intfunc (double t, const double y[], double f[], void *params){ doubl

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ double P, Q; double x = 2.0; P = gsl_cdf_ugaussian_P (x); printf ("prob(

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")doublef (double x, void *p){ if (x < 0.5) return 0.25; else return 0.75;}

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")int main (void) { gsl_combination * c; size_t i; printf ("All subsets of {0,1,2,3} by size

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ double c = GSL_CONST_MKSA_SPEED_OF_LIGHT; double au = GSL_CONST_MKSA_ASTRONO

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.dll")intmain (void){ gsl_block * b = gsl_block_alloc (100); printf ("length of block = %u\n", b

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ double a[] = { 0.11, 0.12, 0.13, 0.21, 0.22, 0.23 }; double

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")struct data{ double x; double y; double z;};int sel_func (void *ntuple_data, vo

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ double x = 1, oldsum = 0, sum = 0; int i = 0; gsl_ieee_env_setu

#define GSL_DLL#include #include #pragma comment(lib, "libgsl_dll_d.lib")#pragma comment(lib, "libgslcblas_dll_d.lib")intmain (void){ double x = 5.0; double y = gsl_sf_bessel_J0 (x); pr

#define GSL_DLL#include #include #include #include #include #include #pragma comment(lib, "libgsl_dll_d.lib")#pragma comment(lib, "libgslcblas_dll_d.lib")/* Computation of the integral,

#include #include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")#define REAL(z,i) ((z)[2*(i)])#define IMAG(z,i) ((z)[2*(i)+1])intmain (void

#include #include #include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ int i; const int n = 100;// double data[n];

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ const gsl_rng_type * T; gsl_rng * r; gsl_histogram2d * h = gsl_h

#include #include #include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ int i; double xi, yi, x[10], y[10]; printf (

#include #include #include #include #include #include #include #include "expfit.c"#define N 40void print_state (size_t iter, gsl_multifit_fdfsolver * s);#pragma comment(lib, "libgsl_d.li

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")double f (double x, void * params){ return pow (x, 1.5);}intmain (void){ gsl_fu

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ double x; const gsl_rng_type * T; gsl_rng * r; gsl_rng_env_s

#include #include #include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ int N = 4; double x[4] = {0.00, 0.10, 0.27, 0

#define GSL_DLL#include #include #pragma comment(lib, "libgsl_dll_d.lib")#pragma comment(lib, "libgslcblas_dll_d.lib")intmain (void){ double x = 5.0; double expected = -0.177596771314338

#include #include #include #pragma comment(lib, "libgsl_dll_d.lib")#pragma comment(lib, "libgslcblas_dll_d.lib")intmain (void){ double x = 5.0; gsl_sf_result result; double expected =

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain(void){ double data[5] = {17.2, 18.1, 16.5, 18.3, 12.6}; double mean, variance, larges

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")#define N 20intmain (void){ double t[N]; double sum_accel, err; double sum = 0;

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ int i; gsl_vector * v = gsl_vector_alloc (3); for (i = 0; i < 3; i++)

几个有用的数学库，GMP，MPIR，MPFR，MPC，CGAL，GSL。前两个是大整数库，MPFR是浮点库，MPC是复数库，CGAL是计算几何库，GSL是gnu的科学计算库，C语言写的。GSL有些问题，比如阶乘，实际上指提供到170的阶乘，超过就溢出了，同时，为了提高速度这个是硬编码进去的，实际上没有计算。用了一个数组，保存了0~170阶乘的结果，用的时候直接查找。所以一些使用它的软件，例如oct

1.介绍2.库的使用 本章描述如何编译使用GSL的程序，介绍GSL的一般用法。2.1例子程序 2.2编译和链接 2.3共享库2.4与ANSI C兼容性2.5inline函数 2.6长双精度long double2.7可移植性函数 2.8可优化函数2.9对不同的数据类型的支持2.10与C++兼容性2.11数组别名 2.12线程安全 2.13遗弃的函数 2.14代码重用3.出错处理 本章描述GSL函数

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")/* set up parameters for this simulated annealing run *//* how many points do we try be

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain(void){ double data[5] = {17.2, 18.1, 16.5, 18.3, 12.6}; double median, uppe

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ const gsl_rng_type * T; gsl_rng * r; size_t i, k = 5, N = 100000; doubl

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ const gsl_rng_type * T; gsl_rng * r; int i, n = 10; gsl_rng_env_setup()

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")gsl_rng * r; /* global generator */intmain (void){ const gsl_rng_type * T; gsl_rng_env_se

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void){ double a_data[] = { 0.18, 0.60, 0.57, 0.96, 0.41, 0.24, 0

#include #include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")double fn1 (double x, void * params){ return cos(x) + 1.0;}intmain (void

#include #include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")struct data{ double x; double y; double z;};intmain (void){ const gsl_rng_t

#include #include #pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (void) { gsl_permutation * p = gsl_permutation_alloc (3); gsl_permutation_init (p);