I have a FEDVR package
But that is a bit limited at the moment, in that it only supports Dirichlet1 boundary conditions, since I've yet to fix the issue with restriction matrices (i.e. dropping the first and last basis functions)
I guess I don't really understand the support for Complex numbers in ApproxFun. If I try to evaluate a complex argument to a Fun(cos) then it returns zero if the imaginary part is non zero. Curiously, it will return a complex number with a correct real part if the imaginary part of the input is zero. Also, I can easily create Fun's with complex coefficients that will return a complex number when evaluating a real number, but again will return zero if the imaginary part is non-zero. Here are a few examples. (I believe all of these work in Chebfun, btw.)
julia (v1.2)> f = Fun(cos);
julia (v1.2)> f(.1)
0.9950041652780257
julia (v1.2)> f(.1 + .1im)
0.0
julia (v1.2)> f(.1 + .0im)
0.9950041652780257 + 0.0im
julia (v1.2)> g = Fun(Chebyshev(), randn(Complex{Float64}, 20));
julia (v1.2)> g(.1)
2.1013219596855888 - 0.11052912191219874im
julia (v1.2)> g(.1 + .1im)
0.0 + 0.0im
julia (v1.2)> cos(.1 + .1im)
0.9999833333373015 - 0.00999998888888977im
If you want to evaluate on a domain in the complex plane you can use say a
Circle or Segment
I’d be shocked if those work in Chebfun unless it is doing extrapolation which is numerically dangerous. But you can use
extrapolate(f,x) in ApproxFun if you wish
Thanks. I see your point. It does look, however, that Chebfun takes more of a "let the buyer beware" type of view. It looks like it will happily evaluate many complex args. As long as you are inside of the Bernstein ellipse, life is good. It will even plot the Bernstein ellipse as well. As you say, I can use extrapolate for those cases, but that does beg the question, is returning zero the best way to flag that error?
If it were redesigned it would probably throw an error: it doesn’t really make sense to automatically do analytic continuation while also supporting piecewise functions
using ApproxFun, Plots, LinearAlgebra
const c₀ = 3.0*10^8
d = Circle(0.0, 1.0, true) # Defines a circle
Δ = Laplacian(d) # Represent the Laplacian
f = ones(∂(d)) # one at the boundary
ω = Interval(0.0..10.0^9)
k = ω/c₀
u = \([Dirichlet(d); Δ+k*I], [f;0.]; # Solve the PDE
tolerance=1E-5)
surface(u) # Surface plotI get the error 'Implement Laplacian(Laurant(:clock:), 1)'
When looking in to the code for ApproxFun, i can find no reference to the laplacian, only a confusing macro for implementing general derivative operators.
Where would I go to work out how to implement new Operators?
I was also wondering if it was yet possible to build spaces of arbitrary dimension to solve 3D PDEs, possibly with extra non spatial dimensions to represent free parameters
do you actually want a
Disk()? That sort of existed at one point but the code is crusty. We do have a version in MultivariateOrthogonalPolynomials.jl which uses @MikaelSlevinsky’s awesome FastTransform to do function approximation on the disk, but there’s no build script for that
For 3D, we’re starting to implement that in MultivariateOrthogonalPolynomials.jl for the cone (random yes, it’s because I’m writing a paper about the cone with Yuan Xu)
But honestly ApproxFun needs a significant refactor before 3D can be tackled seriously: it’s possible to get it working but requires fighting hacky code the whole way, and is not amendable to high performance
There is work behind the scenes on ContinuumArrays.jl to provide a more sound underpinning that scales to 3D, but it won’t be worked on seriously in the next 6 months or so
If you are interested in helping with that I’m happy to talk more.
Hi. I'm working with the Gross-Pitaevskii Equation in two dimensions, with periodic boundary conditions, and I'm trying out the periodic parts of ApproxFun. Could someone please remake the ApproxFun manual? Apparently it's been 10 months since
PeriodicInterval was updated to PeriodicSegment in the source tree, but the Latest version at github.io is that out of date.
I've noticed that ApproxFun doesn't reexport DomainSets. So the manual describes ProductDomain and UnionDomain, but these aren't exported by ApproxFun.
Is this a bug or a documentation error? I.e., was the intention to export those symbols from ApproxFun, or should the manual say that you need a "using DomainSets" to access them?
Also, there doesn't seem to be any documentation for boundary. I'll raise an issue for that.
Is there a reasonable way to represent functions with discontinuities (i.e., the heaviside function) in this package? I was thinking about trying to solve an boundary value problem with an RHS of that form, but for something like
H=Fun(x->(0<x)), the package returns ┌ Warning: Maximum number of coefficients 1048577 reached in constructing Fun.
└ @ ApproxFunBase /Users/guardian/.julia/packages/ApproxFunBase/RLdLv/src/constructors.jl:138
there are obvious Gibbs oscillations at the jump discontinuity
I am trying to solve the Lane Emden equation. It's working with n as an integer but I get an error for real values
Is there any workaround to have the code work solving non linear pde equation with real value exponent?
MethodError: no method matching ^(::ApproxFun.DualFun{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},ConstantOperator{Float64,Chebyshev{Interval{:closed,:closed,Int64},Float64}}}, ::Float64)
Closest candidates are:
^(!Matched::Missing, ::Number) at missing.jl:115
^(!Matched::Float64, ::Float64) at math.jl:870
^(!Matched::Irrational{:ℯ}, ::Number) at mathconstants.jl:91
...
Stacktrace:
[1] (::var"#65#66")(::ApproxFun.DualFun{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},ConstantOperator{Float64,Chebyshev{Interval{:closed,:closed,Int64},Float64}}}) at ./In[49]:6
[2] newton(::var"#65#66", ::Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}}; maxiterations::Int64, tolerance::Float64) at /projects/d1d068b5-3a9d-4e16-9252-52734dea2c8e/.julia/packages/ApproxFun/IC3js/src/Extras/autodifferentiation.jl:98
[3] newton(::Function, ::Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}}) at /projects/d1d068b5-3a9d-4e16-9252-52734dea2c8e/.julia/packages/ApproxFun/IC3js/src/Extras/autodifferentiation.jl:95
[4] top-level scope at In[49]:7x = Fun(identity, 0..10)
n = 2.3
u₀ = 0.0x # initial guess
N = u -> [u'(0), u(0)-1, x*u'' + 2*u' + x*u^n]
u = newton(N, u₀) # perform Newton iteration in function space
plot(u)
Hello.
Could someone please comment on the ease of using ApproxFun.jl for the simulation of a 1D and 2D Fokker-Planck equations (modeling e.g. particle distributions)? In essence, I want to simulate an advection-diffusion equation with variable advection and diffusion coefficients (nonconstant over the space domain).
Thank you!
PS: sorry if this is not the right place for such questions.
Could someone please comment on the ease of using ApproxFun.jl for the simulation of a 1D and 2D Fokker-Planck equations (modeling e.g. particle distributions)? In essence, I want to simulate an advection-diffusion equation with variable advection and diffusion coefficients (nonconstant over the space domain).
Thank you!
PS: sorry if this is not the right place for such questions.
Hello ApproxFun devs; I am solving a nonlinear ODE problem over a domain that could include a singularity for the function. The goal is to stop the solution before it gets there, but its location is not known a priori. Just by extending the endpoint of the domain, I can go from getting a converged solution to the following error:
ERROR: MethodError: no method matching setdomain(::ApproxFunBase.PiecewiseSpace{Tuple{ApproxFunOrthogonalPolynomials.Chebyshev{DomainSets.ChebyshevInterval{Float64},Float64},ApproxFunSingularities.JacobiWeight{ApproxFunOrthogonalPolynomials.Chebyshev{IntervalSets.Interval{:closed,:closed,Float64},Float64},IntervalSets.Interval{:closed,:closed,Float64},Float64,Int64},ApproxFunSingularities.JacobiWeight{ApproxFunOrthogonalPolynomials.Chebyshev{IntervalSets.Interval{:closed,:closed,Float64},Float64},IntervalSets.Interval{:closed,:closed,Float64},Float64,Int64}},DomainSets.UnionDomain{Float64,Tuple{DomainSets.ChebyshevInterval{Float64},IntervalSets.Interval{:closed,:closed,Float64},IntervalSets.Interval{:closed,:closed,Float64}}},Float64}, ::IntervalSets.Interval{:closed,:closed,Float64}and if I extend it further, it results in memory exhaustion. All three of the domain spans (working, MethodError, memory exhaustion) include the singularity. I have tried providing a coarse value to the tolerances argument to newton() and it doesn't change the outcomes. This is on ApproxFun 0.12.5. Ideally the whole thing would error out consistently and gracefully so it could be caught; certainly the MethodError and memory exhaustion shouldn't be occurring.
I'd prefer not to post my code publicly, but would be happy to make available to any ApproxFun devs.
Hello people, can someone guide me as to how can we incorporate neumann boundary conditions while solving a heat equations using ApproxFun.jl? Currently I am using ApproxFun for space discretization and DifferentialEquations.jl for time integration (solving ode system). I tried using the Dirichlet(S,1) where S is the Chebyshev space for applying the neumann boundary condition, but I do not know how to set up the ODE system. Any help would be appreciated