{"id":260,"date":"2015-08-24T00:24:31","date_gmt":"2015-08-23T16:24:31","guid":{"rendered":"http:\/\/sqrt-1.me\/?p=260"},"modified":"2016-11-19T14:41:57","modified_gmt":"2016-11-19T06:41:57","slug":"%e7%bb%98%e5%88%b6%e6%95%b4%e7%b3%bb%e6%95%b0%e5%a4%9a%e9%a1%b9%e5%bc%8f%e5%a4%8d%e6%95%b0%e6%a0%b9%e7%9a%84%e5%88%86%e5%b8%83","status":"publish","type":"post","link":"https:\/\/sqrt-1.me\/?p=260","title":{"rendered":"\u7ed8\u5236\u6574\u7cfb\u6570\u591a\u9879\u5f0f\u590d\u6570\u6839\u7684\u5206\u5e03"},"content":{"rendered":"<p>\u6211\u5728\u767e\u5ea6\u8d34\u5427Mathematica\u5427\u7684<a href=\"http:\/\/tieba.baidu.com\/p\/3622255435\">\u8fd9\u4e2a<\/a>\u5e16\u5b50\u770b\u5230\u6709\u4eba\u5728\u7814\u7a76\u5982\u4e0b\u95ee\u9898\uff1a<\/p>\n<p><strong>\u6c42\u89e3\u6240\u6709\u7cfb\u6570\u4e3a1\u6216-1\u7684n\u6b21\u591a\u9879\u5f0f\u7684\u6839\uff0c\u7136\u540e\u628a\u6240\u6709\u8fd9\u6837\u7684\u6839\u7ed8\u5236\u5230\u590d\u5e73\u9762\u4e0a\u3002<\/strong><\/p>\n<p>\u6211\u60f3\u8d77\u4ee5\u524d\u5728<a href=\"http:\/\/www.matrix67.com\/blog\/archives\/2615\">Matrix67\u7684\u535a\u5ba2<\/a>\u4e0a\u4e5f\u770b\u5230\u8fc7\u8fd9\u4e2a\u95ee\u9898\uff0c\u91cc\u9762\u4ecb\u7ecd\u4e86Sam Derbyshire\u7528Mathematica\u8dd1\u4e86\u56db\u5929\u56db\u591c\uff08\u4e5f\u6709\u8bf4\u6cd5\u662f\u4e09\u5929\u4e09\u591c\uff09\u751f\u6210\u4e86n=24\u7684\u6240\u6709\u6839\uff0c\u7136\u540e\u7528Java\u7a0b\u5e8f\u7ed8\u56fe\uff0c\u5f97\u5230\u4e86\u4e0b\u56fe\uff1a<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/\" alt=\"\" \/><a href=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots1.png\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-262\" src=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots1.png\" alt=\"roots1\" width=\"620\" height=\"438\" srcset=\"https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots1.png 620w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots1-300x212.png 300w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 984px) 61vw, (max-width: 1362px) 45vw, 600px\" \/><\/a><\/p>\n<p>\u8d34\u5427\u4e2d\u7684\u201c\u4e09\u5206\u949f\u590d\u73b0\u201d\u662f\u7b97\u5230n=17\u7684\u60c5\u51b5\u3002\u7531\u4e8en\u589e\u52a01\u5bfc\u81f4\u9700\u8981\u6c42\u89e3\u65b9\u7a0b\u7684\u6570\u91cf\u7ffb\u500d\uff0c\u6bcf\u4e2a\u65b9\u7a0b\u4e2d\u6839\u7684\u4e2a\u6570\u4e5f\u589e\u52a01\uff0c\u6240\u4ee5\u6211\u4f30\u8ba1\u4e86\u4e00\u4e0b\uff0cn=24\u65f6\u5927\u7ea6\u9700\u8981\u6c42\u51fa24*2^24=402M\u4e2a\u89e3\u3002\uff0824\u6b21\u591a\u9879\u5f0f\u5305\u62ec\u5e38\u6570\u9879\u670925\u4e2a\u7cfb\u6570\uff0c\u4f46\u662f\u7531\u4e8e\u5bf9\u79f0\u6027\u53ef\u4ee5\u7701\u53bb\u4e00\u534a\u8ba1\u7b97\u3002\uff09\u8fd9\u4e9b\u8ba1\u7b97\u663e\u7136\u4e0d\u9700\u8981\u4ee5\u5929\u4e3a\u6570\u91cf\u7ea7\u7684\u65f6\u95f4\u3002<\/p>\n<p>\u6211\u51b3\u5b9a\u7528C\u8bed\u8a00\u89e3\u51b3\u8fd9\u4e2a\u95ee\u9898\uff0c\u540c\u65f6\u6539\u5584\u4e00\u4e0b\u663e\u793a\u7684\u6548\u679c\uff0c\u5e76\u5c1d\u8bd5\u7ed8\u5236\u5206\u8fa8\u7387\u66f4\u9ad8\u7684\u56fe\u3002<\/p>\n<p>\u89e3\u51b3\u8fd9\u4e2a\u95ee\u9898\u6709\u4e24\u4e2a\u6bd4\u8f83\u9ebb\u70e6\u7684\u4e8b\uff1a<\/p>\n<ul>\n<li>\u5982\u4f55\u6c42\u89e3\u591a\u9879\u5f0f\u7684\u6839<\/li>\n<li>\u56fe\u4e0a\u7684\u989c\u8272\u6309\u7167\u4ec0\u4e48\u7ed8\u5236<\/li>\n<\/ul>\n<p><!--more--><\/p>\n<p>\u7b2c\u4e00\u4e2a\u95ee\u9898\u6211\u67e5\u627e\u4e86\u4e00\u4e9b\u65b9\u6848\uff0c\u6700\u540e\u89c9\u5f97<a href=\"https:\/\/ring0.me\/\">boj<\/a>\u7684\u5efa\u8bae\u6bd4\u8f83\u597d\uff1a\u4f7f\u7528<a href=\"http:\/\/www.gnu.org\/software\/gsl\/gsl.html\">GNU Scientific Library<\/a>\u5e93\u3002\u8fd9\u4e2a\u5e93\u91cc\u9762\u6709\u5f88\u591a\u79d1\u5b66\u8ba1\u7b97\u7684\u51fd\u6570\uff0c\u5176\u4e2d\u7684<a href=\"http:\/\/www.gnu.org\/software\/gsl\/manual\/html_node\/Roots-of-Polynomials-Examples.html#Roots-of-Polynomials-Examples\">gsl_poly_complex_solve<\/a>\u5c31\u53ef\u4ee5\u6c42\u89e3\u591a\u9879\u5f0f\u7684\u6839\u3002<\/p>\n<p>\u4e8e\u662f\u6211\u5199\u51fa\u6c42\u89e3\u591a\u9879\u5f0f\u7684\u7a0b\u5e8f\uff0c\u8fd9\u4e2a\u7a0b\u5e8f\u628a\u6240\u6709\u7684\u6839\u90fd\u4ee5\u4e8c\u8fdb\u5236\u5f62\u5f0f\u76f4\u63a5\u4fdd\u5b58\u5728\u6587\u4ef6\u4e2d\u3002fork\u51fa8\u4e2a\u8fdb\u7a0b\u53ef\u4ee5\u628aCPU\u8dd1\u6ee1\u3002<\/p>\n<p>\u6ce8\uff1a\u8bf7\u5728cygwin\u6216linux\u4e0b\u7f16\u8bd1\uff0c\u7f16\u8bd1\u51fa\u9519\u8bf7\u5148\u68c0\u67e5\u662f\u5426\u6b63\u786e\u5b89\u88c5GSL\u5e93\uff0c\u7136\u540e\u5c1d\u8bd5\u5728\u7f16\u8bd1\u547d\u4ee4\u4e2d\u52a0 <span class=\"lang:sh decode:true  crayon-inline \" >-lgsl -lgslcblas -lm<\/span> \u53c2\u6570<\/p>\n<pre class=\"lang:c decode:true \" >\/*\u7a0b\u5e8f\u529f\u80fd\uff1a\u6c42\u89e3\u6240\u6709\u6307\u5b9a\u6b21\u6570\u7684\u3001\u7cfb\u6570\u4e3a-1\u30011\u7684\u591a\u9879\u5f0f\uff0c\u5e76\u628a\u590d\u6570\u89e3\u4fdd\u5b58\u5230r[0-7].data\u4e2d\r\n\u8f93\u51fa\u683c\u5f0f\uff1a\u4e8c\u8fdb\u5236\u6587\u4ef6\uff0c\u987a\u5e8f\u6392\u5217\u6bcf\u4e2a\u89e3(double)\u7684\u5b9e\u90e8\u548c\u865a\u90e8*\/\r\n\r\n#include &lt;stdio.h&gt;\r\n#include &lt;gsl\/gsl_poly.h&gt; \/\/GNU science library \u591a\u9879\u5f0f\r\n#include &lt;sys\/wait.h&gt;\r\n#include &lt;unistd.h&gt;\r\n\r\n#define deg 24 \/\/\u591a\u9879\u5f0f\u6b21\u6570\r\n\r\nint main(){\r\n  int i,j,p,total;\r\n  double a[deg+1],z[deg*2]; \/\/a:\u591a\u9879\u5f0f\u7cfb\u6570 z:\u6c42\u89e3\u7ed3\u679c\r\n  FILE *fp;\r\n  char fn[10];\r\n  a[deg]=1.0; \/\/\u51cf\u5c11\u4e00\u534a\u8fd0\u7b97\r\n  \r\n  for(p=0;p&lt;8;p++){\r\n    if(fork()==0){\r\n      sprintf(fn,\"r%d.data\",p);\r\n      fp=fopen(fn,\"wb\");\r\n      for(i=0;i&lt;3;i++) \/\/\u8be5\u8fdb\u7a0b\u7684\u4e09\u4e2a\u786e\u5b9a\u7cfb\u6570\r\n        a[deg-1-i]=(p&amp;(1&lt;&lt;i))?1.0:-1.0;\r\n      total=(1&lt;&lt;(deg-3));\r\n      for(i=0;i&lt;total;i++){\r\n        for(j=0;j&lt;deg-3;j++)\r\n          a[j]=(i&amp;(1&lt;&lt;j))?1.0:-1.0;\r\n        gsl_poly_complex_workspace *w=gsl_poly_complex_workspace_alloc(deg+1);\r\n        gsl_poly_complex_solve(a,deg+1,w,z); \/\/\u6c42\u89e3\r\n        gsl_poly_complex_workspace_free(w);\r\n        fwrite(z,sizeof(double),deg*2,fp);\r\n\r\n        if(i%10000==0)printf(\"Process%d : %.2f%%\\n\",p,i*100.0\/total);\r\n      }\r\n      printf(\"Process%d : Completed\\n\",p);\r\n      exit(0);\r\n    }\r\n  }\r\n  while(wait(NULL)&gt;0); \/\/\u7b49\u5f85\u5176\u4ed6\u8fdb\u7a0b\r\n  return 0;\r\n}<\/pre>\n<p>\u7a0b\u5e8f\u51e0\u5206\u949f\u5c31\u8dd1\u5b8c\u4e86\uff0c\u751f\u6210\u4e866G\u7684\u6570\u636e\u3002<\/p>\n<p>\u7136\u540e\u60f3\u529e\u6cd5\u89e3\u51b3\u7b2c\u4e8c\u4e2a\u95ee\u9898\u3002<\/p>\n<p>\u9996\u5148\u8981\u5f04\u660e\u767d\u56fe\u4e2d\u7684\u989c\u8272\u8868\u793a\u4ec0\u4e48\u3002<strong>\u4e00\u4e2a\u50cf\u7d20\u7684\u989c\u8272\u8868\u793a\u90a3\u4e2a\u50cf\u7d20\u8303\u56f4\u5185\u6839\u7684\u5bc6\u5ea6\u3002<\/strong>\u6240\u4ee5\u7b80\u5355\u7684\u529e\u6cd5\u5c31\u662f\u6309\u7167\u671f\u671b\u7684\u5206\u8fa8\u7387\uff08\u4f8b\u59822000*2000\uff09\u5efa\u7acb\u4e00\u4e2a\u4e8c\u7ef4\u6570\u7ec4\uff0c\u7136\u540e\u6839\u636e\u6bcf\u4e2a\u6839\u7684\u4f4d\u7f6e\u7edf\u8ba1\u6bcf\u4e2a\u50cf\u7d20\u8303\u56f4\u5185\u6839\u7684\u4e2a\u6570\u3002<\/p>\n<p>\u90a3\u4e48\u53c8\u5982\u4f55\u6839\u636e\u6839\u7684\u4e2a\u6570\u7b97\u51fa\u989c\u8272\uff1f\u989c\u8272\u503c\u662f\u6709\u8303\u56f4\u9650\u5236\u7684\uff0c\u4f8b\u59820~255\uff0c\u6240\u4ee5\u6211\u60f3\u5230\u8ba1\u7b97\u6839\u4e2a\u6570\u7684\u6700\u5927\u503c\u3002\u4f46\u662f\u7ecf\u8fc7\u6d4b\u8bd5\u53d1\u73b0\uff0c\u753b\u51fa\u7684\u56fe\u5f88\u6697\uff0c\u51e0\u4e4e\u770b\u4e0d\u6e05\uff0c\u53ea\u6709\u67d0\u4e9b\u70b9\u5904\u662f\u4eae\u7684\uff0c\u8fd9\u8bf4\u660e\u6839\u805a\u96c6\u5728\u67d0\u4e9b\u70b9\u5904\uff0c\u8fd9\u4e9b\u70b9\u9644\u8fd1\u6839\u7684\u4e2a\u6570\u548c\u5176\u4ed6\u5e73\u51e1\u70b9\u9644\u8fd1\u6839\u7684\u4e2a\u6570\u4e4b\u6bd4\u53ef\u80fd\u662f\u53d1\u6563\u7684\uff0c\u4e0d\u80fd\u901a\u8fc7\u7edf\u8ba1\u6700\u5927\u503c\u6765\u51b3\u5b9a\u989c\u8272\u3002<\/p>\n<p>\u6211\u76f4\u63a5\u624b\u5de5\u8bbe\u5b9a\u4e86\u4e00\u4e2a\u4eae\u5ea6\u7684\u5fae\u8c03\u503c\uff0c\u8c03\u8282\u5230\u6bd4\u8f83\u6ee1\u610f\u7684\u4f4d\u7f6e\u3002\u6839\u7684\u4e2a\u6570\u4ece\u5c11\u5230\u591a\uff0c\u989c\u8272\u4ece\u9ed1\u8272\u6e10\u53d8\u5230\u6a59\u8272\uff08RGB=255,127,0\uff09\uff0c\u8d85\u8fc7\u8bbe\u5b9a\u503c\u7684\u5c11\u91cf\u90e8\u5206\u518d\u6e10\u53d8\u5230\u767d\u8272\u3002\u6d4b\u8bd5\u53d1\u73b0\u6697\u90e8\u4e0d\u6e05\u695a\uff0c\u4e8e\u662f\u7528\u5e73\u65b9\u6839\u51fd\u6570\u628a\u4eae\u5ea6\u5f80\u4e0a\u8c03\u4e86\u4e00\u4e0b\uff08\u548c\u5f00\u6839\u53f7\u518d\u4e58\u4ee510\u7684\u8c03\u5206\u65b9\u6cd5\u4e00\u4e2a\u9053\u7406\uff09\u3002<\/p>\n<p>\u6839\u7684\u5bc6\u5ea6\u548c\u989c\u8272\u7684\u5173\u7cfb\u5982\u4e0b\u56fe\u6240\u793a\uff1a<\/p>\n<p><a href=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots2.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots2.png\" alt=\"roots2\" width=\"587\" height=\"416\" class=\"alignnone size-full wp-image-267\" srcset=\"https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots2.png 587w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots2-300x213.png 300w\" sizes=\"auto, (max-width: 587px) 85vw, 587px\" \/><\/a><\/p>\n<p>\u7a0b\u5e8f\u5982\u4e0b\uff1a<\/p>\n<pre class=\"lang:c decode:true \" >\/*\u7a0b\u5e8f\u529f\u80fd\uff1a\u6839\u636er[0-7].data\u4e2d\u7684\u590d\u6570\u6570\u636e\uff0c\u751f\u6210\u5bc6\u5ea6\u56fe*\/\r\n#include &lt;stdio.h&gt;\r\n#include &lt;stdlib.h&gt;\r\n#include &lt;math.h&gt;\r\n\r\n#define DIM 1024 \/\/\u56fe\u7247\u5206\u8fa8\u7387 DIM*DIM\r\n#define RANGE 2.0 \/\/\u7ed8\u5236\u8303\u56f4 x,y\u2208[-RANGE,RANGE]\r\n#define B 12.0 \/\/\u4eae\u5ea6\u5fae\u8c03\u53c2\u6570\uff0c\u8d8a\u5927\u8d8a\u6697\r\n#define BUFFSIZE ((1&lt;&lt;24)*24\/8*2) \/\/\u8bfb\u6587\u4ef6\u7f13\u51b2\u533a\u5927\u5c0f\uff0c24\u6b21\u591a\u9879\u5f0f\u7ea6800M\r\n\r\ndouble *data; \/\/\u8bfb\u6587\u4ef6\u7f13\u51b2\u533a\r\nint (*den)[DIM]; \/\/\u50cf\u7d20\u70b9\u8303\u56f4\u5185\u6839\u7684\u4e2a\u6570\r\nunsigned char (*bitmap)[DIM][3]; \/\/\u7528\u4e8e\u751f\u6210\u4f4d\u56fe\r\n\r\nint total=0,max;\r\n\r\nvoid countroot(int cnt){ \/\/\u6309\u7167\u50cf\u7d20\u533a\u57df\u5bf9\u65b9\u7a0b\u7684\u6839\u8ba1\u6570\r\n  double x,y;\r\n  int i,px,py;\r\n  for(i=0;i&lt;cnt;i++){\r\n    x=data[i*2];\r\n    y=data[i*2+1];\r\n    px=DIM\/2*x\/RANGE+DIM\/2; \/\/\u8ba1\u7b97\u5bf9\u5e94\u7684\u50cf\u7d20\u5750\u6807\r\n    py=DIM\/2*y\/RANGE+DIM\/2;\r\n    if(px&gt;0&amp;&amp;px&lt;DIM&amp;&amp;py&gt;0&amp;&amp;py&lt;DIM)\r\n      den[px][py]++;\r\n  }\r\n  total+=cnt; \/\/\u603b\u70b9\u6570\u8ba1\u6570\r\n}\r\n\r\nvoid colorfunc(int d,unsigned char *color){ \/\/\u6839\u5bc6\u5ea6\u4e0e\u989c\u8272\u7684\u6620\u5c04\u51fd\u6570\r\n  int t;\r\n  t=(long long)255*sqrt((double)d\/max);\r\n  if(t&lt;256){ \/\/\u9ed1\u81f3\u6a59\u6e10\u53d8\r\n    color[0]=t;\r\n    color[1]=t\/2;\r\n    color[2]=0;\r\n  }else if(t&lt;512){ \/\/\u6a59\u81f3\u767d\r\n    color[0]=255;\r\n    color[1]=128+(t-256)\/2;\r\n    color[2]=t-256;\r\n  }else{ \/\/\u767d\r\n    color[0]=color[1]=color[2]=255;\r\n  }\r\n}\r\n\r\nvoid genpic(){ \/\/\u6839\u636e\u5bc6\u5ea6\u6570\u636e\u751f\u6210\u4f4d\u56fe\r\n  int i,j;\r\n  printf(&quot;Generating pic...\\n&quot;);\r\n  max=B*total\/DIM\/DIM; \/\/\u989c\u8272255\u5bf9\u5e94\u7684\u6839\u5bc6\u5ea6\r\n  for(j=0;j&lt;DIM;j++)\r\n    for(i=0;i&lt;DIM;i++)\r\n      colorfunc(den[i][j],bitmap[j][i]);\r\n}\r\n\r\nvoid savefile(){ \/\/\u4fdd\u5b58\u6587\u4ef6\r\n  char fn[20];\r\n  sprintf(fn,&quot;roots_%d.ppm&quot;,DIM);\r\n  FILE *fp=fopen(fn,&quot;wb&quot;);\r\n  printf(&quot;Saving to %s...\\n&quot;,fn);\r\n  fprintf(fp,&quot;P6\\n%d %d\\n255\\n&quot;,DIM,DIM);\r\n  fwrite(bitmap,1,DIM*DIM*3,fp);\r\n  fclose(fp);\r\n  printf(&quot;Completed!\\n&quot;);\r\n}\r\n\r\nint main(){\r\n  int p,cnt;\r\n  char fn[10];\r\n  FILE *fp;\r\n  data=malloc(BUFFSIZE*sizeof(double));\r\n  den=malloc(DIM*DIM*sizeof(int));\r\n  bitmap=malloc(DIM*DIM*3);\r\n  if(!data||!den||!bitmap){\r\n    printf(&quot;Out of memory!\\n&quot;);\r\n    exit(-1);\r\n  }\r\n  for(p=0;p&lt;8;p++){\r\n    sprintf(fn,&quot;r%d.data&quot;,p);\r\n    printf(&quot;Loading : %s\\n&quot;,fn);\r\n    fp=fopen(fn,&quot;rb&quot;);\r\n    fseek(fp,0,SEEK_END); \/\/\u5b9a\u4f4d\u5230\u6587\u4ef6\u5c3e\r\n    cnt=ftell(fp)\/sizeof(double)\/2; \/\/\u6839\u636e\u6587\u4ef6\u5927\u5c0f\u8ba1\u7b97\u6570\u636e\u91cf\r\n    fseek(fp,0,SEEK_SET); \/\/\u5b9a\u4f4d\u5230\u6587\u4ef6\u5934\r\n    fread(data,sizeof(double),2*cnt,fp);\r\n    fclose(fp);\r\n    printf(&quot;Loaded, count=%d\\n&quot;,cnt);\r\n    printf(&quot;Processing : %s\\n&quot;,fn);\r\n    countroot(cnt);\r\n  }\r\n  genpic();\r\n  savefile();\r\n  return 0;\r\n}\r\n<\/pre>\n<p>\u8fd0\u884c\u7ed3\u679c\uff1a<\/p>\n<p><a href=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots3.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots3.png\" alt=\"roots3\" width=\"1024\" height=\"1024\" class=\"alignnone size-full wp-image-269\" srcset=\"https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots3.png 1024w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots3-150x150.png 150w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots3-300x300.png 300w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 1362px) 62vw, 840px\" \/><\/a><\/p>\n<p>\u6ce8\uff1a\u5206\u8fa8\u7387\u592a\u9ad8\u4e4b\u540e\u6709\u4e9b\u56fe\u7247\u67e5\u770b\u5668\u53ef\u80fd\u4e0d\u8bc6\u522b\uff0c\u6211\u6d4b\u8bd5\u547d\u4ee4\u884c\u7248\u7684ffmpeg\u53ef\u4ee5\u6253\u5f00\uff0c\u5546\u4e1a\u8f6f\u4ef6Photoshop\u66f4\u65b9\u4fbf\u4e00\u4e9b\u3002<\/p>\n<p>\u5982\u679c\u4e0d\u8003\u8651\u78c1\u76d8\u8bfb\u5199\uff08\u4e2d\u95f4\u6587\u4ef6\u6211\u4fdd\u5b58\u5728\u4e86ramdisk\uff09\uff0c\u751f\u6210\u51e0\u4e07\u4e58\u51e0\u4e07\u5206\u8fa8\u7387\u7684\u56fe\u4e5f\u4e0d\u8fc7\u51e0\u5206\u949f\u3002<\/p>\n<p>\u540e\u6765\u6211\u53c8\u8dd1\u4e86n=26\u7684\u60c5\u51b5\uff0c\u53d1\u73b0\u9664\u4e86\u7ec6\u8282\u653e\u5927\u540e\u6bd4\u8f83\u5e73\u6ed1\u4e4b\u5916\u6ca1\u6709\u663e\u8457\u533a\u522b\uff0c\u6545\u4e0d\u5355\u72ec\u8d34\u56fe\u3002\u4f46\u662f\u8fd0\u884c\u523098%\u65f6\u9047\u5230\u4e86\u4e00\u4e2a\u9519\u8bef\uff1a<\/p>\n<p><code>gsl: zsolve.c:78: ERROR: root solving qr method failed to converge<br \/>\nDefault GSL error handler invoked.<\/code><\/p>\n<p>\u5e76\u4e0d\u662f\u5185\u5b58\u4e0d\u8db3\u6216\u6ea2\u51fa\u5bfc\u81f4\uff0c\u5e94\u8be5\u662fGSL\u89e3\u67d0\u4e2a\u65b9\u7a0b\u65f6\u51fa\u4e86\u95ee\u9898\u3002<\/p>\n<p>\u6709\u4eba\u95ee\u6211\u4e3a\u4ec0\u4e48\u4e0d\u662f\u8ba1\u7b97\u51fa\u6839\u5206\u5e03\u540e\u76f4\u63a5\u7edf\u8ba1\u5e76\u751f\u6210\u56fe\u7247\uff0c\u800c\u662f\u4fdd\u5b58\u4e2d\u95f4\u6587\u4ef6\u3002\u6211\u89c9\u5f97\u7a0b\u5e8f\u5728\u8fd0\u884c\u65f6\u7ecf\u5e38\u9700\u8981\u7528\u540c\u4e00\u7ec4\u6570\u636e\u751f\u6210\u4e0d\u540c\u5206\u8fa8\u7387\u7684\u56fe\u7247\uff0c\u6240\u4ee5\u524d\u8005\u53cd\u800c\u4e0d\u65b9\u4fbf\u3002<\/p>\n<p>\u6700\u540e\u8d34\u51e0\u5f20\u7ec6\u8282\u56fe\uff0c\u611f\u53d7\u4e00\u4e0b\u6570\u5b66\u4e4b\u7f8e\uff1a<\/p>\n<p><a href=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots5.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots5.png\" alt=\"roots5\" width=\"1024\" height=\"1024\" class=\"alignnone size-full wp-image-272\" srcset=\"https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots5.png 1024w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots5-150x150.png 150w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots5-300x300.png 300w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 1362px) 62vw, 840px\" \/><\/a><\/p>\n<p><a href=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots6.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots6.png\" alt=\"roots6\" width=\"1024\" height=\"1024\" class=\"alignnone size-full wp-image-273\" srcset=\"https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots6.png 1024w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots6-150x150.png 150w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots6-300x300.png 300w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 1362px) 62vw, 840px\" \/><\/a><\/p>\n<p><a href=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots4.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots4.png\" alt=\"roots4\" width=\"1024\" height=\"1024\" class=\"alignnone size-full wp-image-271\" srcset=\"https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots4.png 1024w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots4-150x150.png 150w, https:\/\/sqrt-1.me\/wp-content\/uploads\/2015\/08\/roots4-300x300.png 300w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 1362px) 62vw, 840px\" \/><\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u6211\u5728\u767e\u5ea6\u8d34\u5427Mathematica\u5427\u7684\u8fd9\u4e2a\u5e16\u5b50\u770b\u5230\u6709\u4eba\u5728\u7814\u7a76\u5982\u4e0b\u95ee\u9898\uff1a \u6c42\u89e3\u6240\u6709\u7cfb\u6570\u4e3a1\u6216-1\u7684n\u6b21\u591a\u9879\u5f0f\u7684\u6839 &hellip; <a href=\"https:\/\/sqrt-1.me\/?p=260\" class=\"more-link\">\u7ee7\u7eed\u9605\u8bfb<span class=\"screen-reader-text\">\u201c\u7ed8\u5236\u6574\u7cfb\u6570\u591a\u9879\u5f0f\u590d\u6570\u6839\u7684\u5206\u5e03\u201d<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"footnotes":"","_links_to":"","_links_to_target":""},"categories":[2,3],"tags":[],"class_list":["post-260","post","type-post","status-publish","format-standard","hentry","category-c-language","category-math"],"_links":{"self":[{"href":"https:\/\/sqrt-1.me\/index.php?rest_route=\/wp\/v2\/posts\/260","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sqrt-1.me\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/sqrt-1.me\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/sqrt-1.me\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/sqrt-1.me\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=260"}],"version-history":[{"count":11,"href":"https:\/\/sqrt-1.me\/index.php?rest_route=\/wp\/v2\/posts\/260\/revisions"}],"predecessor-version":[{"id":334,"href":"https:\/\/sqrt-1.me\/index.php?rest_route=\/wp\/v2\/posts\/260\/revisions\/334"}],"wp:attachment":[{"href":"https:\/\/sqrt-1.me\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=260"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sqrt-1.me\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=260"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sqrt-1.me\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=260"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}