Recent changes to 2682: Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/Recent changes to 2682: Function zeta fails numerically for large numbers that aren't even integersenTue, 28 Jan 2014 23:15:18 -0000#2682 Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/?limit=25#18ef<div class="markdown_content"><p>I believe Euler-Maclaurin series is valid over the entire plane, but I'd have to go back read over everything again.</p>
<p>But as I mentioned zetam1 is a feature request.</p></div>Raymond ToyTue, 28 Jan 2014 23:15:18 -0000https://sourceforge.net5971474cbee69b62920437403ccd662ff912a6ab#2682 Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/?limit=25#f386<div class="markdown_content"><p>Great! So on what region should we be using that method? Using it, you can calculate zeta by <em>just adding one at the end</em>. (which is why I was/am slightly sceptical)</p>
<p>I note that the paper that describes the series we're using starts by talking about Euler-Maclaurin summation of the standard series and explains that their method is faster. But maybe this is something about real vs. complex arguments?</p></div>Rupert SwarbrickTue, 28 Jan 2014 08:22:48 -0000https://sourceforge.net84703c22d7a9026de87165ef6a1bc892e1cc9786#2682 Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/?limit=25#a01d<div class="markdown_content"><p>I'm sure there are better ways, but you can do sum(1/n^s, n, 1, inf) via a Euler-Maclaurin series to get zetam1(s). The Euler-Maclaurin series converges quite a bit faster than the original series.</p></div>Raymond ToyTue, 28 Jan 2014 02:19:28 -0000https://sourceforge.net68573e3229c789cadf890ac7a34779017ba8de44#2682 Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/?limit=25#543f<div class="markdown_content"><p>Maybe I've misunderstood, but doesn't that only help if you've got a way to calculate zeta(s)-1? I mean, yes, you can calculate zeta(s) with higher floating precision then subtract 1 and return the result in just a double, but that seems pretty niche...</p>
<p>I haven't seen any papers that say "Ahah, near the real line, you can compute zeta(x)-1 with the following series" (and then don't do something silly like adding and subtracting terms of order one), but maybe I've not been looking in the right place?</p></div>Rupert SwarbrickSun, 26 Jan 2014 21:21:26 -0000https://sourceforge.net1b5f21081df5647174686103f2d775ff664e2ca5#2682 Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/?limit=25#90ed<div class="markdown_content"><p>I think zetam1 is useful, but that's at least a feature request, not a bug.<br />
If you're computing zeta - 1, a zetam1 function could allow you do the work with doubles instead of bigfloats, or less precise bigfloats than would otherwise would be needed to compute zeta - 1 to decent accuracy.</p></div>Raymond ToySun, 26 Jan 2014 21:15:55 -0000https://sourceforge.net2632ed370f90822889abd2a16be4ebd62965716b#2682 Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/?limit=25#ad10<div class="markdown_content"><ul>
<li><strong>status</strong>: open --> closed</li>
</ul></div>Rupert SwarbrickSun, 26 Jan 2014 18:50:18 -0000https://sourceforge.net3d0dff4c8de8a5cbdca21fc1f04c9c77f9d17ecfFunction zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/<div class="markdown_content"><p>The zeta function provides incorrect numerical output for numbers that aren't even integers greater than some number between 33.44 and 33.45. The output for 33.44 is 1.000000000128721, and for 33.45 is 0.0. It should be clear from the definition ( <a href="http://dlmf.nist.gov/25.2" rel="nofollow">http://dlmf.nist.gov/25.2</a> ) that for finite positive real values zeta should be > 1 (>= 1.0 given the accuracy limitations of floating point). For values this large, direct calculation using the definition should suffice.</p>
<p>Another potentially useful thing to do would be to provide a function zetam1(s) that is formally equal to zeta(s) - 1 that can work for large values of s where the difference between zeta(s) and 1 would be lost in floating point cancellation.</p></div>Sean LakeThu, 23 Jan 2014 23:38:25 -0000https://sourceforge.net721e2147d8babd75ae5d32ee89d7acaf9743a693Function zeta fails numerically for large numbers that aren't even integershttps://sourceforge.net/p/maxima/bugs/2682/<div class="markdown_content"><p>Ticket 2682 has been modified: Function zeta fails numerically for large numbers that aren't even integers<br />
Edited By: Rupert Swarbrick (rswarbrick)<br />
Status updated: u'open' => u'closed'</p></div>Sean LakeThu, 23 Jan 2014 23:38:25 -0000https://sourceforge.netc12f5813c9a4fd4d8c25b9590dda24c4c4c838e9