Add Support for Arbitrary Precision Arithmetic with BigFloat

Is your feature request related to a problem? Please describe. I tried running ‘pysr’ on a 1,000 row array with 4 integer input variables and one integer output variable - a Goedel Number.

From Mathematica:

GoedelNumber[l_List] := Times @@ MapIndexed[Prime[First[#2]]^#1 &, l]

E.g.

Data file:
# 7	1	5	8	6917761200000

julia> 2^7*3^1*5^5*7^8
6917761200000

The model returned:

Complexity  Loss       Score     Equation
1           Inf       NaN       0.22984365

I am just learning ‘pysr’ and maybe it’s just ‘user error’. However, Inf and Nan suggest that Goedel numbers may exceed Float64.

Describe the solution you’d like Not sure what happened, because the largest Goedel number in the input is: 1.6679880978201e+23

Additional context I didn’t see any parameters to set ‘verbose’ mode or ‘debugging’ information.

GoedelTableFourParameters.txt

[dbl001]

btw - I showed ‘PYSR’ yesterday to a UCSC Pure Math PhD student. He commented 'this is going to put me out of a job. 😉

[MilesCranmer]

Recall that you edited your local copy to convert to BigFloats. Install pysr with pip install -U pysr to get the normal one back.