05/12/2022, 11:58 AM
I forgot to mention
To avoid the use of constants for the gaussian method we can " anchor " the point 0 :
We use this modification
f(s) = exp( R(s-1) f(s-1) )
Now R(s) = t(s) erf(s)^2 = (1 + erf(s))/2 * erf(s)^2.
So that R(s) = 0 when s = 0.
Therefore
f(1) = exp(0 * ..) = 1.
This helps numerically and theoretically , however R(s) is slightly more complicated.
But R(s) behaves similar on the complex plane as t(s) so it works.
Extending this idea is worth considering.
But at least
sexp(s-1) = lim ln^[n] f(s + n)
where sexp(0) = 1.
regards
tommy1729
Tom Marcel Raes
To avoid the use of constants for the gaussian method we can " anchor " the point 0 :
We use this modification
f(s) = exp( R(s-1) f(s-1) )
Now R(s) = t(s) erf(s)^2 = (1 + erf(s))/2 * erf(s)^2.
So that R(s) = 0 when s = 0.
Therefore
f(1) = exp(0 * ..) = 1.
This helps numerically and theoretically , however R(s) is slightly more complicated.
But R(s) behaves similar on the complex plane as t(s) so it works.
Extending this idea is worth considering.
But at least
sexp(s-1) = lim ln^[n] f(s + n)
where sexp(0) = 1.
regards
tommy1729
Tom Marcel Raes

