-
Notifications
You must be signed in to change notification settings - Fork 71
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
Newton for simple PDE gives NaN result #887
Comments
dstahlke
added a commit
to dstahlke/ApproxFun.jl
that referenced
this issue
Jun 12, 2023
Test case (requires also JuliaApproximation/ApproxFunBase.jl#479): using ApproxFun xdom = Chebyshev(-1..1) ydom = Chebyshev(-1..1) domain = xdom * ydom x,y = Fun(identity, domain) Dx = Derivative(Chebyshev()^2, [1,0]) Dy = Derivative(Chebyshev()^2, [0,1]) N(u, v) = [ 2*u - x; 3*v + y ] u0 = one(x) * one(y) v0 = one(x) * one(y) u, v = newton(N, [u0, v0]) Bug report: JuliaApproximation#887
Both these test cases are fixed by the following PRs: |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
This simple test gives a NaN result:
Result:
And I don't know if it's related, but with a two-dimensional domain it gives a different error:
Result:
The text was updated successfully, but these errors were encountered: