JuliaApproximation / JuliaApproximation/ApproxFun.jl
Error solving nonlinear boundary value problem / worked in a previous version of Julia 111 and Approxfun
Nobody has claimed this yet.
- Dominant language
- Julia
- Stars
- 559
- Forks
- 71
- PR merge metrics
- No merged PRs in 30d
Description
I try to solve the following nonlinear Differential Equation/Boundary Value Problem:
y'''' + 4.0*y'''/x = 0
with the following boundary conditions
y(0) = y'(0) = 0 y''(1)=0 y'''(1)=1
in a previous version the following solution worked:
using ApproxFun
x = Fun(identity,0..5)
u0 = 0.0x
N = u->[u(0), u'(0), u''(1), u'''(1)-1, u''''+4*u^3/x]
u = newton(N,u0)
Warning: Maximum numer of iterations 15 reached, with approximate accuracy of 1.439e-13
the results were ok
if I try to solve the same problem with the latest version
I rceive the following error messages:
ERROR: AssertionError: abs(first(g)) < tol
Stacktrace:
if I change the statement N = u-> (because of singularity at x=0) to
N = u->[u(0), u'(0), u''(1), u'''(1)-1,
if x<1e-4 u''''+12*y'y^2 else u''''+4u^3/x end]
u = newton(N,u0)
I receive the same error message
ERROR: AssertionError: abs(first(g)) < tol
Stacktrace:
Can anybody help me
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by reproducing the nonlinear boundary value problem with ApproxFun's newton entry point and inspect the omitted assertion stacktrace. Compare the previous and latest versions around the singular expression at x=0; done means the supplied problem solves without AssertionError: abs(first(g)) < tol.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- julia
- Domain
- tooling
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100