Skip to content

Conversation

dstahlke
Copy link

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:
#887

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
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant