<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On 

<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Sun, 8 Apr 2018 08:14:47 +1000</span>, 

<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Virginia Frey </span>wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">The current implementation of Slepians in scipy's signal module only<br>
offers the zeroth order, so it is not possible to use for e.g. the<br>
multitaper spectral analysis. In addition to that, because of the way<br>
the eigensystem is solved, it breaks with a MemoryError when using high<br>
numbers of points (i.e. a few thousand). The request for high numbers of<br>
points is not unreasonable: people do take long time-traces and the FFT<br>
makes it possible to compute Fourier transforms on those.<br></blockquote><div><br></div><div>In SciPy 1.1.0, <span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">windows</span>.slepian will be deprecated and replaced by <span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">windows.dpss</span>

:</div><div><br></div><div><a href="https://scipy.github.io/devdocs/generated/scipy.signal.windows.dpss.html" target="_blank">https://scipy.github.io/<wbr>devdocs/generated/scipy.<wbr>signal.windows.dpss.html</a></div><div><a href="https://github.com/scipy/scipy/pull/7802" target="_blank">https://github.com/scipy/<wbr>scipy/pull/7802</a><br></div><div><a href="https://github.com/scipy/scipy/issues/4354" target="_blank">https://github.com/scipy/<wbr>scipy/issues/4354</a><br></div><div><br></div><div>Does that do what you want?</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I have working code for this, which is based on the implementation<br>
suggested in [5] but is in pure Python.</blockquote><div><br></div><div>Aren't there copyright/license issues with using code from Numerical Recipes?</div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Once we have the full set of Slepians available, from there we (or I)<br>
could also implement the multitaper spectral analysis technique<br></blockquote><div><br></div><div>That would be useful!</div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">So, what do you guys think of my proposal? If you like it, please advise<br>
me on how to proceed as this is the first time I'm planning to<br>
contribute. From what I took of the web page, I'll start by downloading<br>
the repo - and then should I make a new branch? I'm also happy to send<br>
some code through via email for your review first.</blockquote><div><br></div><div>Can you read through 

<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><a href="https://github.com/scipy/scipy/pull/7802" target="_blank">https://github.com/scipy/<wbr>scipy/pull/7802</a></span> and comment if you see a problem with it, or it doesn't do all that you would want it to do?  </div><div><br></div><div>Also we were a little confused about how best to normalize the amplitude of the functions.  In your work, do you use the L2 norm?</div><div><br></div></div></div></div>