Re: DEP: deprecate complex dtypes in PchipInterpolator and Akima1DInterpolator

Evgeni Burovski <[email protected]> Mon, 25 Dec 2023 23:22:39 +0300
Newsgroups gmane.comp.python.scientific.devel
Message-ID <CAMRo0iu-AgZHqaqDzQRVKBg16xhNb5Sei5wUvhPZX3Tgyz3oEg@mail.gmail.com>
--===============1073320683723828586==
Content-Type: multipart/alternative; boundary="000000000000824a3b060d5b5413"

--000000000000824a3b060d5b5413
Content-Type: text/plain; charset="UTF-8"

>
>
>
> Is there any reason not to fix the underlying problem instead?
>
>
The underlying problem is not fixable. Since pchip interpolator is
nonlinear in data, one cannot expect that

pchip(x, re + i*im).c  == pchip(x, re).c + i*pchip(x, im).c

cf
https://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.PchipInterpolator.html


> On Mon, Dec 25, 2023, 10:36 Lucas Colley <[email protected]> wrote:
>
>> Hi all,
>>
>> https://github.com/scipy/scipy/pull/19752 deprecates complex dtypes in
>> PchipInterpolator and Akima1DInterpolator. It also clarifies the
>> documentation for these interpolators and RegularGridInterpolator w.r.t.
>> complex data.
>>
>> In Evgeni(@ev-br)'s words, these interpolators return "nonsense" for
>> complex data. This deprecation moves us towards raising an exception on the
>> suggestion of Marten(@mhvk), as "people may be getting hit with bad
>> results, better that they know!".
>>
>> If anyone objects to the deprecation, please let us know!
>>
>> Cheers,
>> Lucas
>> _______________________________________________
>> SciPy-Dev mailing list -- [email protected]
>> To unsubscribe send an email to [email protected]
>> https://mail.python.org/mailman3/lists/scipy-dev.python.org/
>> Member address: [email protected]
>>
> _______________________________________________
> SciPy-Dev mailing list -- [email protected]
> To unsubscribe send an email to [email protected]
> https://mail.python.org/mailman3/lists/scipy-dev.python.org/
> Member address: [email protected]
>

--000000000000824a3b060d5b5413
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr"><div class=3D"gmail_quote"><blockquote class=3D"gmail_quot=
e" style=3D"margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204)=
;padding-left:1ex"><div dir=3D"auto"><br><br>Is there any reason not to fix=
 the underlying problem instead?<br><br></div></blockquote><div><br></div><=
div>The underlying problem is not fixable. Since pchip interpolator is nonl=
inear in data, one cannot expect that</div><div><br></div><div>pchip(x, re =
+ i*im).c=C2=A0 =3D=3D pchip(x, re).c=C2=A0+ i*pchip(x, im).c</div><div><br=
></div><div>cf=C2=A0<a href=3D"https://docs.scipy.org/doc/scipy/reference/g=
enerated/scipy.interpolate.PchipInterpolator.html">https://docs.scipy.org/d=
oc/scipy/reference/generated/scipy.interpolate.PchipInterpolator.html</a></=
div><div>=C2=A0</div><blockquote class=3D"gmail_quote" style=3D"margin:0px =
0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div=
 dir=3D"auto"><div><div class=3D"gmail_quote"><div dir=3D"ltr" class=3D"gma=
il_attr">On Mon, Dec 25, 2023, 10:36 Lucas Colley &lt;<a href=3D"mailto:luc=
[email protected]" target=3D"_blank">[email protected]</a>&gt; wro=
te:<br></div><blockquote class=3D"gmail_quote" style=3D"margin:0px 0px 0px =
0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi all,<br>
<br>
<a href=3D"https://github.com/scipy/scipy/pull/19752" rel=3D"noreferrer nor=
eferrer" target=3D"_blank">https://github.com/scipy/scipy/pull/19752</a> de=
precates complex dtypes in PchipInterpolator and Akima1DInterpolator. It al=
so clarifies the documentation for these interpolators and RegularGridInter=
polator w.r.t. complex data.<br>
<br>
In Evgeni(@ev-br)&#39;s words, these interpolators return &quot;nonsense&qu=
ot; for complex data. This deprecation moves us towards raising an exceptio=
n on the suggestion of Marten(@mhvk), as &quot;people may be getting hit wi=
th bad results, better that they know!&quot;.<br>
<br>
If anyone objects to the deprecation, please let us know!<br>
<br>
Cheers,<br>
Lucas<br>
_______________________________________________<br>
SciPy-Dev mailing list -- <a href=3D"mailto:[email protected]" rel=3D"no=
referrer" target=3D"_blank">[email protected]</a><br>
To unsubscribe send an email to <a href=3D"mailto:[email protected]=
g" rel=3D"noreferrer" target=3D"_blank">[email protected]</a><br>
<a href=3D"https://mail.python.org/mailman3/lists/scipy-dev.python.org/" re=
l=3D"noreferrer noreferrer" target=3D"_blank">https://mail.python.org/mailm=
an3/lists/scipy-dev.python.org/</a><br>
Member address: <a href=3D"mailto:[email protected]" rel=3D"noreferr=
er" target=3D"_blank">[email protected]</a><br>
</blockquote></div></div></div>
_______________________________________________<br>
SciPy-Dev mailing list -- <a href=3D"mailto:[email protected]" target=3D=
"_blank">[email protected]</a><br>
To unsubscribe send an email to <a href=3D"mailto:[email protected]=
g" target=3D"_blank">[email protected]</a><br>
<a href=3D"https://mail.python.org/mailman3/lists/scipy-dev.python.org/" re=
l=3D"noreferrer" target=3D"_blank">https://mail.python.org/mailman3/lists/s=
cipy-dev.python.org/</a><br>
Member address: <a href=3D"mailto:[email protected]" target=3D"_bl=
ank">[email protected]</a><br>
</blockquote></div></div>

--000000000000824a3b060d5b5413--

--===============1073320683723828586==
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

_______________________________________________
SciPy-Dev mailing list -- [email protected]
To unsubscribe send an email to [email protected]
https://mail.python.org/mailman3/lists/scipy-dev.python.org/
Member address: [email protected]

--===============1073320683723828586==--