[Numpy-discussion] Polynomial evaluation inconsistencies

Eric Wieser wieser.eric+numpy at gmail.com
Mon Jul 2 02:57:24 EDT 2018


I think the `x` is just noise there, especially if it's ignored (that is,
`T[0](x*2)` doesn't do anything reasonable).

    Chebyshev.literal(lambda T: 1*T[0] + 2*T[1] + 3*T[2])

Would work, but honestly I don't think that provides much clarity. I think
the value here is mainly for "simple" polynomials.

On Sun, 1 Jul 2018 at 23:42 Maxwell Aifer <maifer at haverford.edu> wrote:

> Say we add a constructor to the polynomial base class that looks something
> like this:
>
>
> -------------------------------------------------------------------------------------------
>    @classmethod
>     def literal(cls, f):
>         def basis_function_getter(self, deg):
>             coefs = [0]*deg + [1]
>             return lambda _: cls(coefs)
>         basis = type('',(object,),{'__getitem__': basis_function_getter})()
>         return f(basis, None)
>
> -------------------------------------------------------------------------------------------
>
>
> Then the repr for, say, a Chebyshev polynomial could look like this:
>
> >>> Chebyshev.literal(lambda T,x: 1*T[0](x) + 2*T[1](x) + 3*T[2](x))
>
> Does this sound like a good idea to anyone?
>
> Max
>
>
> On Sat, Jun 30, 2018 at 6:47 PM, Charles R Harris <
> charlesr.harris at gmail.com> wrote:
>
>>
>>
>> On Sat, Jun 30, 2018 at 4:42 PM, Charles R Harris <
>> charlesr.harris at gmail.com> wrote:
>>
>>>
>>>
>>> On Sat, Jun 30, 2018 at 3:41 PM, Eric Wieser <
>>> wieser.eric+numpy at gmail.com> wrote:
>>>
>>>> Since the one of the arguments for the decreasing order seems to just
>>>> be textual representation - do we want to tweak the repr to something like
>>>>
>>>> Polynomial(lambda x: 2*x**3 + 3*x**2 + x + 0)
>>>>
>>>> (And add a constructor that calls the lambda with Polynomial(1))
>>>>
>>>> Eric
>>>>
>>>
>>> IIRC there was a proposal for that. There is the possibility of adding
>>> renderers for latex and html that could be used by Jupyter, and I think the
>>> ordering was an option.
>>>
>>
>> See https://github.com/numpy/numpy/issues/8893 for the proposal. BTW, if
>> someone would like to work on this, go for it.
>>
>> Chuck
>>
>>>>>>>
>>>
>>
>> _______________________________________________
>> NumPy-Discussion mailing list
>> NumPy-Discussion at python.org
>> https://mail.python.org/mailman/listinfo/numpy-discussion
>>
>> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion at python.org
> https://mail.python.org/mailman/listinfo/numpy-discussion
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20180701/ef1a643e/attachment-0001.html>


More information about the NumPy-Discussion mailing list