我对使用Pari / GP用户定义函数和返回函数的用户定义函数有疑问。这是一个愚蠢的例子:
? f(x) = { (t) -> x + t }
%1 = (x)->(t)->x+t
? g(x) = { local(y); y = f(x); (t) -> y(t) }
%2 = (x)->local(y);y=f(x);(t)->y(t)
? h = g(2)
%3 = (t)->my(x=2);y(t)
? h(1)
*** at top-level: h(1)
*** ^----
*** in function h: y(t)
*** ^----
*** not a function in function call
*** Break loop: type 'break' to go back to GP
我期望得到h(1)=(g(2))(1)= y(1)=(f(2))(1)= 3。我是否在Pari中遇到一流的功能限制,或者我做错了什么?如果是后者,我该如何解决?
最佳答案
注意,local
关键字将变量限制在动态范围内。确保在调用y
时,来自同一个/外部作用域的那家伙h(1)
引用了预期的东西。在这种情况下,适当的修复方法如下:
?f = (x) -> { (t) -> x + t }
%1 = (x)->(t)->x+t
?g = (x) -> { local(y); (t) -> y(x)(t) }
%1 = (x)->local(y);(t)->y(x)(t)
?h = g(2)
%1 = (t)->my(x=2);y(x)(t)
?y = (x) -> f(x);
?h(1)
%1 = 3
?y = 1;
?h(1)
*** not a function in function call
另外,我强烈建议对变量使用静态(词法)作用域,以简化代码并消除引用依赖性。为此使用
my
说明符。查看下面的固定代码:?f(x) = { (t) -> x + t }
%1 = (x)->my(t);t->x+t
?g(x) = { my (y = f(x)); (t) -> y(t) }
%1 = (x)->my(y=f(x));(t)->y(t)
?h = g(2)
%1 = (t)->my(x=2,y=(t)->my(x=2);x+t);y(t)
?h(1)
%1 = 3
注意,评估最新版本的
h
花费的时间更少。另外,我建议您使用PARI / GP作者(http://pari.math.u-bordeaux1.fr/pub/pari/manuals/2.4.2/users.pdf)的代码样式关于first-class-functions - Pari/GP:用户定义的函数作为返回值,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/30318668/