Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
fast accurate Kneser sexp algorithm
#34
(09/05/2011, 05:05 AM)sheldonison Wrote:
(09/04/2011, 08:16 PM)JmsNxn Wrote: Sheldon, I'm wondering if I'm doing something wrong here, but whenever I try to plot anything with this function:

Code:
V(a,b,t) = {
  if (real(t)<1,  return(exp2(t,exp2(-t,a)+exp2(-t,b))));
  if (real(t)>=1, return(exp2(t,exp2(-t,a)+exp2(-1,b))));
}

exp2(t, a) = {
  if (a==2, return(2));
  if (a == 4, return(4));
  if (real(a)>2, return(superf(isuperf(a)+t)), return(superf2(isuperf2(a)+t)));
}

I repeatedly get a cryptic Win32 error:
Code:
gp 4852 child_copy: linked dll data write copy failed

This is my proposed "natural" extension to the Ackermann function, I'm confused as to why it won't plot.

I appreciate any help you can give me, thanks.
I'll take a look.... you're probably using base sqrt(2), which has two different superfunction solutions and fixed points of 2 and 4. I'll post some more after I do some experiments for base sqrt(2), along with your routine.
- Shel
I wasn't able to reproduce the error you saw, but I think I know what the problem is. sqrt(2)<eta, which is a special case for the program. The problem I suspect you're having is that I changed the code in this last round of optimizations, so that the superf/isuperf routines now use a power series, but the power series for the superfunction and inverse superfunction only gets calculated when you calculate the kneser mapping, which didn't used to be the case. I forgot that you were using the program the code only for the superfunctions, and not to generate a kneser mapping. I realize you have no interest in the esoteric kneser mapping I calculate for the sqrt(2) base, but it only takes a few seconds to caulcate. The fix for your problem is probably to execute the "loop;" command after init(sqr(2)). Then the graphs I generated for exp2(-2,t) and exp2(+2,t), worked fine.
- Shel
Code:
\r kneser.gp
init(sqrt(2));
loop;
exp2(t, a) = {
  if (a==2, return(2));
  if (a == 4, return(4));
  if (real(a)>2, return(superf(isuperf(a)+t)), return(superf2(isuperf2(a)+t)));
}

ploth(t=1.5,5,z=exp2(2,t);[real(z),imag(z)]);
ploth(t=1.5,5,z=exp2(-2,t);[real(z),imag(z)]);
For the most recent code version: go to the Nov 21st, 2011 thread.
Reply


Messages In This Thread
The pari-GP code - by sheldonison - 08/07/2010, 09:17 PM
updated kneser.gp code - by sheldonison - 08/19/2010, 02:35 AM
RE: updated kneser.gp code - by nuninho1980 - 08/19/2010, 12:08 PM
RE: updated kneser.gp code - by sheldonison - 08/20/2010, 01:05 AM
update to support B<eta - by sheldonison - 11/15/2010, 02:53 PM
RE: update to support B<eta - by nuninho1980 - 11/15/2010, 03:26 PM
another new version - by sheldonison - 11/17/2010, 06:52 PM
RE: fast accurate Kneser sexp algorithm - by sheldonison - 09/05/2011, 01:07 PM

Possibly Related Threads...
Thread Author Replies Views Last Post
  "Kneser"/Riemann mapping method code for *complex* bases mike3 2 6,775 08/15/2011, 03:14 PM
Last Post: Gottfried
  Attempt to make own implementation of "Kneser" algorithm: trouble mike3 9 15,991 06/16/2011, 11:48 AM
Last Post: mike3
  Numerical algorithm for Fourier continuum sum tetration theory mike3 12 19,944 09/18/2010, 04:12 AM
Last Post: mike3
  regular sexp: curve near h=-2 (h=-2 + eps*I) Gottfried 2 6,030 03/10/2010, 07:52 AM
Last Post: Gottfried
  Attempting to compute the kslog numerically (i.e., Kneser's construction) jaydfox 11 18,682 10/26/2009, 05:56 PM
Last Post: bo198214
  regular sexp:different fixpoints Gottfried 6 11,934 08/11/2009, 06:47 PM
Last Post: jaydfox
  sexp(strip) is winding around the fixed points Kouznetsov 8 13,201 06/29/2009, 10:05 AM
Last Post: bo198214
  sexp and slog at a microcalculator Kouznetsov 0 3,200 01/08/2009, 08:51 AM
Last Post: Kouznetsov



Users browsing this thread: 1 Guest(s)