win32-gas 175 ms for FPU sqrt, sum= 66977004.90197226 332 ms for the proxy, sum= 66977016.92774981 win32-gcc -O 2 105 ms for FPU sqrt, sum= 66977004.90197226 135 ms for the proxy, sum= 66977016.91725335 win64-gas64 102 ms for FPU sqrt, sum= 66977004.90197226 280 ms for the proxy, sum= 66977016.9277...

@jj2007
if you compile with optimization using gcc the non-asm version will be completely eliminated because the calculations are not used outside the timing loop, that's why I include some kind of sum in the loop and print it out after

nice :-)

Hi dodicat
minimax approximations are better, I presented my scheme as a novelty, I don't think anyone else would be crazy enough to come up with it

just for fun to see if it was possible to approximate functions other than the gamma function using the Lanczos scheme to illustrate I will use a system of linear equations of degree 5, this will not give us much precision but it will show the algorithm lets approximate sqr between 1 and 2, here I w...

hi caseih
I just tried 2^2^3 on Maple 2016 but it complains with "syntax error, ambiguous use of `^`, please use parentheses"
pari/gp evaluates 2^2^3 to 256, FB and C++ give 64
just noticed that Ed's Eval can easily be changed to associativity = rightassoc

thanks Ed :-)

hello jaskin your example is an unfortunate bomb, it seems to overwrites a critical part of memory causing havoc try this example from the FB manual, I added storing to string '' test.bas '' compile with fbc -exx -lang fblite test.bas #lang "fblite" Sub Generate_Error dim as string s On Er...

@Ed Davis
here's challenge for you, how would you add functions of more than one argument like hypot and atan2 ?

nice :-)

congrats meucat :-)

where is the FreeBASIC-1.07.3-win64-gcc-5.2.0 folder located? is it in restricted area like "C:\Program Files" ? you are only supposed to run the batch file by drag-and-drop the FreeBASIC-1.07.3-win64-gcc-5.2.0 folder onto it anyway, all the batch file does is copy the lib and bi files to ...

meucat wrote:(by the way... I still use UBASIC in my old windows XP machine to run arbitrary precision so want to move these ubasic programs to FreeBasic environment)

if you need any trig or exponential functions then I recommend that you use mpfr instead of gmp

hi meucat I am at a loss as to why it doesn't work for you, may I ask you that you rename your present FreeBasic folder and install the latest official release from the news section of this forum? also, re-download the gmp+mpfr archive and after extracting it drag-and-drop your new FreeBasic folder ...

but.. when I run it only prints almost instantly "automatic casting from mpq 1/3 to mpf" and then program stops with no errors or warnings to check not sure what you mean, but do this, at the end of the program put: Sleep Sleep will cause the program to wait until a key is pressed, this p...