JuliaApproximation / JuliaApproximation/ApproxFun.jl

Error solving nonlinear boundary value problem / worked in a previous version of Julia 111 and Approxfun

Open
#747 3 comments 0 reactions 0 assignees View on GitHub

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

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. 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

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.