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 <<a href=3D"mailto:luc= [email protected]" target=3D"_blank">[email protected]</a>> 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)'s words, these interpolators return "nonsense&qu= ot; for complex data. This deprecation moves us towards raising an exceptio= n on the suggestion of Marten(@mhvk), as "people may be getting hit wi= th bad results, better that they know!".<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==--