Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fatou.jl #1

Open
4 tasks
bennn opened this issue Mar 21, 2023 · 1 comment
Open
4 tasks

Fatou.jl #1

bennn opened this issue Mar 21, 2023 · 1 comment

Comments

@bennn
Copy link
Member

bennn commented Mar 21, 2023

https://github.com/chakravala/Fatou.jl

  • try running some examples from Fatou.jl and generating pictures
  • do they run with TrackedFloats?
  • are there parameters we can tweak to introduce NaNs & change the picture?
  • what does the TF log say about the NaNs?
@bennn bennn mentioned this issue Mar 21, 2023
5 tasks
@ashton314
Copy link
Member

Tried running Fatou with FloatTracker. Adding support for complex numbers works well, but Fatou explicitly checks for something of type Complex{Float64}, which we cannot handle. Problematic code is in Fatou.jl:327:

323 function orbit(K::Define{FT,QT,CT,M,N,P,D,im},z0::Complex{Float64}) where {FT,QT,CT,M,N,P,D}
324     M ? (z = K.seed) : (z = P ? plane(z0) : z0)
325     zn = 0x0000
326     while (N ? (K.Q(z,z0)::Float64>K.ϵ)::Bool : (K.Q(z,z0)::Float64<K.ϵ))::Bool && K.N>zn
327         z = K.F(z,z0)::Complex{Float64}
328         zn+=0x0001
329     end; #end
330     # return the normalized argument of z or iteration count
331     return (zn::UInt16,(D ? disk(z) : z)::Complex{Float64})
332 end

I'm going to give up on this for now. It could be an interesting library to fuzz, but we'd have to work out how to make it past this block.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants