(04/18/2022, 05:55 PM)Ember Edison Wrote: I suggest to set init_OFF's g(|u|)=log(abs(u)), fine-tuning this function to other values (like log(z)^n) makes little sense, the adjustable range of n is very small.
init_OFF(1, k) log(u)^n
10.0^±24 n>0.5
10.0^±25 n>0.5
10.0^±26 n>0.9
10.0^±10^2 0.97 ≤ n ≤ 1.4
log(z) can to work under
Code:/* 256GB, x86-64 Pari-GP Stack limit */
default(parisize, a1=274877906944); init_OFF(1, 10.0^10^9)
There is also a small problem: Is_Julia() does not work in scenarios where init(1, 10.0^-10) is used, let alone scenarios where init_OFF is needed
I'll look into it, I have a list of things to program and fix in the next iteration of this program. I thought I had already set it to log(abs(u)), but I must've forgotten to change it--I know you pointed this out before.
Is_Julia is pretty wonky for extreme values; it's difficult to implement this better. But I have added some changes already which helped with a quite a few inaccuracies. I'm still fiddling with the program to make it more stream lined.
I've also kind of depreciated and retired the protocol init_off, because I figured out how to do an automatic flag so that all you need is init. Still in a beta stage though, haven't had enough time to work on the program.
