-
Notifications
You must be signed in to change notification settings - Fork 7
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
Fourier{BigFloat} #155
Comments
GenericFFT.jl try loading that it might just work |
I thought GenericFFT was a mirror of the Fourier functionality in FastTransforms? |
Did you try loading julia> using ApproxFun, GenericFFT
julia> f = (r,θ)-> sin(100*r*cos(θ))
#19 (generic function with 1 method)
julia> S = Chebyshev(big"0.1"..big"1.0")⊗Fourier(big"0.0"..big"2.0"*π)
Chebyshev(0.1000000000000000000000000000000000000000000000000000000000000000000000000000002 .. 1.0) ⊗ Fourier(【0.0,6.283185307179586476925286766559005768394338798750211641949889184615632812572396❫)
julia> F = ProductFun(LowRankFun(f, S; gridx = 256, gridy = 256))
ProductFun on Chebyshev(0.1000000000000000000000000000000000000000000000000000000000000000000000000000002 .. 1.0) ⊗ Fourier(【0.0,6.283185307179586476925286766559005768394338798750211641949889184615632812572396❫)
julia> coefficients(F); # get the array of coefficients
julia> f(big"0.123", big"0.456")
-0.9988661226402205869689991698845071591411740075338145723676043325127213580466635
julia> F(big"0.123", big"0.456")
-0.998866122640220586968999169884507159141174007533814572367604332512721358046603
|
|
Odd as ApproxFunFourier.jl uses the exact same transform ( So I'm confused how the ApproxFun code is working |
Ah I see we special cased the transform for |
Is there a way to work with Fourier and BigFloat? It goes via FFTW which is causing some issues.
The text was updated successfully, but these errors were encountered: