I tried using the fast_callable function with a List and keep getting the "Unable to convert...to an integer" error. I am using the list to simply pass parameters to the function.
Are there any special rules or techniques for passing Lists or Arrays to fast-callable ?
Thanks.
The error you encounter is not with fast_callable, but with defining a symbolic function.
In general, it is recommended to write one instruction per line.
This would let you see that while your first example works:
sage: k = 1.8
sage: func(m) = bessel_J(m^2, k)
sage: fmax = fast_callable(func, vars=[m])
sage: print func(1.)
0.581516951731165
sage: print fmax(1.)
0.581516951731165
the second example fails at the stage of defining func(m):
sage: k = 1.8
sage: v = [0.0, 1.0, 4.0, 9.0]
sage: func(m) = bessel_J(v[m], k)
Traceback (most recent call last)
...
TypeError: unable to convert m to an integer
A simpler example of this failure is as follows:
sage: l = [1, 2, 3]
sage: f(m) = l[m]
Traceback (most recent call last)
...
TypeError: unable to convert m to an integer
Thanks for your comment!
It appears that the problem isn't actually with using a List with fast_callable but when I try to use the List index as a variable. I was able to do a workaround by defining a new function by summing all of the terms which refer to this List in the function using a simple loop. I am then able to use fast_callable with this new function. The speedup is almost 10-fold !
I personally suspect you just can't do it without extending fast_callable from its current signature. See the doc:
Currently, x can be an expression object, an element of SR, or a (univariate or multivariate) polynomial; this list will probably be extended soon.
However, hopefully someone else will provide a workaround, perhaps involving creating a list function or something.
Asked: 2016-07-14 09:20:41 +0200
Seen: 509 times
Last updated: Jul 17 '16
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
Hi! Can you give a minimal working (or not working, I guess) code example of exactly what you are trying to do?
Hi, thanks for your comment - I should've included a working snippet example of what I am trying to do.
0.581516951731165; 0.581516951731165
Traceback (click to the left of this block for traceback) ... TypeError: unable to convert m to an integer
Even better, would you edit your question and add the example there?
To display blocks of code, either indent them with 4 spaces, or select the corresponding lines and click the "code" button (the icon with '101 010'). Can you edit your question to do that?
To display inline code, surround it within "backticks" or "backquotes" `.