Quickest way to compute the Abel function on the Shell-Thron boundary
#1
I'm currently writing a protocol to evaluate the modified bennet operators:

\[
\alpha <s>_{\varphi} y = \exp_{y^{1/y}}^{\circ s}\left(\log^{\circ s}_{y^{1/y}}(\alpha) + y + \varphi\right)\\
\]

For \(\varphi\) a complex number--I'm mostly just dealing with real positive values at the moment. The goal is to evaluate the function \(\varphi(\alpha,y,s)\), such that these operators will satisfy Goodstein's equation:

\[
\alpha <s> \left(\alpha <s+1> y\right) = \alpha <s+1> (y+1)\\
\]

But for the moment, I'm just concerning myself with calculating the first function. Everything works great so far, but I'm scratching my head for when \(y^{1/y} = \partial \mathfrak{S}\)-- when it's on the boundary of the Shell-Thron region (equivalently \(|\log(y)| = 1\)). Now I know we can construct a repelling and attracting Abel function about these points--and I know all the theory. But I just realized, I've never actually seen a program that constructs it. I know Sheldon has a program for handling it, but I really don't want to go digging through all the matrix theory. I just want a quick formula. I know if you make a variable change that it becomes pretty elementary.

So for the moment, I can construct \(\alpha <s>_{\varphi} y\), for pretty much the entire complex plane in \(y\) (excluding branch cuts), excluding where \(y^{1/y} \in \mathfrak{S}\). This is primarily because I don't know a fast way to get both abel functions... I could program in a way, but I think it's going to be way too slow. This program is already pretty slow as it is (we have to consistently reinitialize to account for varying bases of the exponential). I don't want to slow it down any more. I was wondering if there's anywhere on this forum that has an easy to read program I can adapt for this.

...I just hope I don't have to write too much just to handle the case \(|\log(y)| = 1\) -_-....

Edit: I thought I'd add that I know how to write in the neutral case but it just nukes the speed of my code. I know how to program in the \(\eta\) case, but I'm wondering what the current fastest way is. For the moment, I'm just returning \(0\) anywhere on the boundary, because it just nukes my code and makes everything so fucking slow for these values of \(y\).
Reply
#2
The function is so much of a Riemann sum generator that the quick and easy formula would be the one with a logarithm at the top of your post but multiplied by a P function variable. Multiplying the top gets you past where I can see you are not getting how to add Riemann sums (specifically after a zeration or two. That's what the P does, multiplies the parameters to get two sums to retrieve a variable
Reply


Possibly Related Threads…
Thread Author Replies Views Last Post
  FractionalIteration now in Wolfram Function Repository Daniel 0 795 07/10/2023, 07:16 AM
Last Post: Daniel
  Terse Schroeder & Abel function code Daniel 1 1,421 10/16/2022, 07:03 AM
Last Post: Daniel
  Revisting my accelerated slog solution using Abel matrix inversion jaydfox 22 49,682 05/16/2021, 11:51 AM
Last Post: Gottfried
  An incremental method to compute (Abel) matrix inverses bo198214 3 16,804 07/20/2010, 12:13 PM
Last Post: Gottfried
  Attempting to compute the kslog numerically (i.e., Kneser's construction) jaydfox 11 38,797 10/26/2009, 05:56 PM
Last Post: bo198214
  Computing Abel function at a given center jaydfox 10 27,953 11/30/2007, 06:44 PM
Last Post: andydude



Users browsing this thread: 1 Guest(s)