使用GSL-GNU科学计算库

2021-09-12 09:31:54 阅读：237 来源： 互联网

标签：GNU 科学计算 poly GSL complex workspace double gsl

这里使用的自带vc工程的1.8版本，地址 http://gnuwin32.sourceforge.net/packages/gsl.htm

这个网页里面有GSL1.8版本的，里面有目录VC8，下面有libgsl.sln, GSL的官方文档在

http://www.gnu.org/software/gsl/doc/latex/gsl-ref.pdf GSL都是C的API和一些结构体，我们

试了下多项式求根的例子，多项式是x^5+1，配置好编译好的gsl的lib, dll和header

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

int main(void)
{
    int i;
    /* coefficients of P(x) = -1 + x^5 */
    double a[6] = { -1, 0, 0, 0, 0, 1 };
    double z[10];
    gsl_poly_complex_workspace* w
        = gsl_poly_complex_workspace_alloc(6);
    gsl_poly_complex_solve(a, 6, w, z);
    gsl_poly_complex_workspace_free(w);
    for (i = 0; i < 5; i++)
    {
        printf("z%d = %+.18f %+.18f\n",
            i, z[2 * i], z[2 * i + 1]);
    }
    return 0;
}

最后控制台输出：

z0 = -0.809016994374947673 +0.587785252292473359
z1 = -0.809016994374947673 -0.587785252292473359
z2 = +0.309016994374947507 +0.951056516295152976
z3 = +0.309016994374947507 -0.951056516295152976
z4 = +0.999999999999999889 +0.000000000000000000

把实根和复根都求出来了，换一下系数，double a[6] = { -1, 0.1, 0.2, 0.3, 0.4, 1 };

多项式是-1+0.1*x+0.2*x^2+0.3*x^3+0.4*x^4+x^5

输出：

z0 = +0.835460440971803497 +0.000000000000000000
z1 = +0.247415004045392517 +1.007728810338902381
z2 = +0.247415004045392517 -1.007728810338902381
z3 = -0.865145224531293833 +0.602636011854109310
z4 = -0.865145224531293833 -0.602636011854109310

标签：GNU,科学计算,poly,GSL,complex,workspace,double,gsl
来源： https://www.cnblogs.com/abcstar/p/15257398.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

使用GSL-GNU科学计算库