blocks|key|2038121|text|通过向tx和Fixed添加额外的元素修复了此问题。我猜这就是你让我做的，@ViniciusMiranda。代码现在显示为|type|unstyled|depth|inlineStyleRanges|entityRanges|data|2038122|++++++++++++++++double+tx[2*NSTEPS%2B1]+=+{0},+theight[2*NSTEPS%2B1]+=+{0};
++++++++++++++++double+terrain[(int)+(2*NSTEPS*100)]+=+{0};
++++++++++++++++double+terrainsl[(int)+(2*NSTEPS*100)]+=+{0};

++++++++++++for+(j+=+0;+j+<+2*NSTEPS%2B1;+j%2B%2B)
++++++++++++{
++++++++++++++++tx[j]+=+(double)+j*space;
++++++++++++++++theight[j]+=+randomground();
++++++++++++}

++++++++++++gsl_interp_accel+*acc+=+gsl_interp_accel_alloc();
++++++++++++gsl_spline+*spline+=+gsl_spline_alloc(gsl_interp_akima,+2*NSTEPS%2B1);
++++++++++++gsl_spline_init(spline,+tx,+theight,+2*NSTEPS%2B1);

++++++++++++for+(q+=+0;+q<+2*NSTEPS*100;+q%2B%2B)
++++++++++++{
++++++++++++++++terrain[q]+=+gsl_spline_eval(spline,xi,acc);
++++++++++++++++terrainsl[q]+=+gsl_spline_eval_deriv(spline,xi,acc);
++++++++++++++++xi+=+xi%2B(double)+space/refine;
++++++++++++}|code-block|syntax|javascript|2038123|我仍然不明白为什么需要这个修复。|2038124|entityMap^0|0|0|0^^$0|@$1|2|3|4|5|6|7|K|8|@]|9|@]|A|$]]|$1|B|3|C|5|D|7|L|8|@]|9|@]|A|$E|F]]|$1|G|3|H|5|6|7|M|8|@]|9|@]|A|$]]|$1|I|3|-4|5|6|7|N|8|@]|9|@]|A|$]]]|J|$]]

Fixed the problem by adding an extra element to tx and theight. I'm guessing this is what you had asked me to do, @ViniciusMiranda. The code now reads

<pre><code> double tx[2*NSTEPS+1] = {0}, theight[2*NSTEPS+1] = {0};
 double terrain[(int) (2*NSTEPS*100)] = {0};
 double terrainsl[(int) (2*NSTEPS*100)] = {0};

 for (j = 0; j &lt; 2*NSTEPS+1; j++)
 {
 tx[j] = (double) j*space;
 theight[j] = randomground();
 }

 gsl_interp_accel *acc = gsl_interp_accel_alloc();
 gsl_spline *spline = gsl_spline_alloc(gsl_interp_akima, 2*NSTEPS+1);
 gsl_spline_init(spline, tx, theight, 2*NSTEPS+1);

 for (q = 0; q&lt; 2*NSTEPS*100; q++)
 {
 terrain[q] = gsl_spline_eval(spline,xi,acc);
 terrainsl[q] = gsl_spline_eval_deriv(spline,xi,acc);
 xi = xi+(double) space/refine;
 }
</code></pre>

I still don't understand why this fix was required though.

blocks|key|4061742|text|Akima样条是一种局部的、子样条插值。对于函数f(x)，如果您尝试获取x_i+<=+x+<=+x_i%2B1的x的值，则Akima需要用于j+=+i-2,+i-1,+i,+i%2B1,+i%2B2,+i%2B3.的(x_j,+f_j)对|type|unstyled|depth|inlineStyleRanges|offset|length|style|CODE|entityRanges|data|4061743|entityMap^0|O|4|10|H|1I|1|1W|V|2S|A|0^^$0|@$1|2|3|4|5|6|7|H|8|@$9|I|A|J|B|C]|$9|K|A|L|B|C]|$9|M|A|N|B|C]|$9|O|A|P|B|C]|$9|Q|A|R|B|C]]|D|@]|E|$]]|$1|F|3|-4|5|6|7|S|8|@]|D|@]|E|$]]]|G|$]]

Akima spline is a local, sub-spline interpolation. For a function <code>f(x)</code>, if you are trying to obtain the value of <code>x</code> for <code>x_i &lt;= x &lt;= x_i+1</code>, Akima spline needs the pairs <code>(x_j, f_j)</code> for <code>j = i-2, i-1, i, i+1, i+2, i+3.</code>

blocks|key|2025541|text|我在我的代码zigzag.sourceforge.net中发现了同样的错误，并通过对gsl+library+source中的代码行进行注释、编译和重新安装来解决。|type|unstyled|depth|inlineStyleRanges|offset|length|style|CODE|entityRanges|data|2025542|在版本1.14之前的|2025543|中，interp.c中有行150|blockquote|2025544|2025545|2025546|//+GSL_ERROR_VAL("interpolation+error",+GSL_EDOM,+GSL_NAN);|code-block|syntax|javascript|2025547|我用评论来解决这个问题！|2025548|entityMap^0|6|M|16|I|0|0|0|0|0|0|0^^$0|@$1|2|3|4|5|6|7|V|8|@$9|W|A|X|B|C]|$9|Y|A|Z|B|C]]|D|@]|E|$]]|$1|F|3|G|5|6|7|10|8|@]|D|@]|E|$]]|$1|H|3|I|5|J|7|11|8|@]|D|@]|E|$]]|$1|K|3|-4|5|6|7|12|8|@]|D|@]|E|$]]|$1|L|3|-4|5|6|7|13|8|@]|D|@]|E|$]]|$1|M|3|N|5|O|7|14|8|@]|D|@]|E|$P|Q]]|$1|R|3|S|5|6|7|15|8|@]|D|@]|E|$]]|$1|T|3|-4|5|6|7|16|8|@]|D|@]|E|$]]]|U|$]]

I found the same error in my code <code>zigzag.sourceforge.net</code> and resolved by commenting the line in the <code>gsl library source</code>, compiling it and reinstalling.

<blockquote>
 Before version 1.14 there is in interp.c the line 150 
</blockquote>

<pre><code>// GSL_ERROR_VAL("interpolation error", GSL_EDOM, GSL_NAN);
</code></pre>

that I comment to solve the break !

I'm getting a 'gsl: interp.c:150: ERROR: interpolation error' with the following code. Some googling says that this error occurs when you try to extrapolate using the interp function but I don't see how that is happening here. Help would be greatly appreciated. Thanks.

The function <code>randomground()</code> just returns a random number (double).

<pre><code>#define NSTEPS 100 

int main()
{
 int j, q, space = 1, refine = 100;
 double xi = 0.0, tx[2*NSTEPS] = {0}, theight[2*NSTEPS] = {0};

 double terrain[(int) (2*NSTEPS*100)] = {0};
 double terrainsl[(int) (2*NSTEPS*100)] = {0};

 for (j = 0; j &lt; 2*NSTEPS; j++)
 {
 tx[j] = (double) j*space;
 theight[j] = randomground();
 }

 gsl_interp_accel *acc = gsl_interp_accel_alloc();
 gsl_spline *spline = gsl_spline_alloc(gsl_interp_akima, 2*NSTEPS);
 gsl_spline_init(spline, tx, theight, 2*NSTEPS);

 for (q = 0; q&lt; 2*NSTEPS*100; q++)
 {
 terrain[q] = gsl_spline_eval(spline,xi,acc);
 terrainsl[q] = gsl_spline_eval_deriv(spline,xi,acc);
 xi = xi+(double) space/refine;
 }
return 0;
}
</code></pre>

Interpolation error with gsl akima spline

翻译质量差，导致语言生硬或混乱。

没有提供实际的解决方法或示例。

解答不清晰，无法理解或解决问题。

页面排版不美观，阅读体验差。

文章

问答

视频

学习中心

腾讯云实验室

直播

竞赛

腾讯云代码分析专区

腾讯iOA零信任安全管理系统专区

腾讯云架构师技术同盟交流圈

腾讯云数据库专区

腾讯云顾问专区

腾讯云原生专区

腾讯混元专区

腾讯云TCE专区

腾讯云Lighthouse专区

腾讯云HAI专区

腾讯云Edgeone专区

腾讯云存储专区

腾讯云智能专区

腾讯轻联专区 

腾讯云开发专区

TAPD专区

腾讯轻量云游戏服专区

腾讯云最具价值专家

腾讯云架构师技术同盟

腾讯云创作之星

腾讯云开发者先锋

腾讯云代码助手

云原生构建

TAPD 敏捷项目管理

Cloud Studio

SDK中心

API中心

命令行工具

涵盖代码开发、场景应用、自动测试全流程，助你从零构建专属AI助手

一站式MCP教程库，解锁AI应用新玩法

我得到一个'gsl: interp.c:150: ERROR:插值错误‘，代码如下。一些谷歌搜索表明，当您尝试使用interp函数进行外推时，会出现此错误，但我不明白这是如何发生的。我们将非常感谢您的帮助。谢谢。函数randomground()只返回一个随机数(double)。#define NSTEPS 100   ...

问gsl akima样条插值误差
EN

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问gsl akima样条插值误差EN