|
39 | 39 |
|
40 | 40 |
|
41 | 41 | function coefficients(f::AbstractVector{T}, a::Jacobi, b::Chebyshev) where T |
42 | | - if domain(a) == domain(b) && (!isapproxinteger(a.a-0.5) || !isapproxinteger(a.b-0.5)) |
| 42 | + if domainscompatible(a, b) && !(isapproxinteger_addhalf(a.a) && isapproxinteger_addhalf(a.b)) |
43 | 43 | jac2cheb(f, strictconvert(T,a.a), strictconvert(T,a.b)) |
44 | 44 | else |
45 | 45 | defaultcoefficients(f,a,b) |
46 | 46 | end |
47 | 47 | end |
48 | 48 | function coefficients(f::AbstractVector{T}, a::Chebyshev, b::Jacobi) where T |
49 | 49 | isempty(f) && return f |
50 | | - if domain(a) == domain(b) && (!isapproxinteger(b.a-0.5) || !isapproxinteger(b.b-0.5)) |
| 50 | + if domainscompatible(a, b) && !(isapproxinteger_addhalf(b.a) && isapproxinteger_addhalf(b.b)) |
51 | 51 | cheb2jac(f, strictconvert(T,b.a), strictconvert(T,b.b)) |
52 | 52 | else |
53 | 53 | defaultcoefficients(f,a,b) |
54 | 54 | end |
55 | 55 | end |
56 | 56 |
|
57 | 57 | function coefficients(f::AbstractVector,a::Jacobi,b::Jacobi) |
58 | | - if domain(a) == domain(b) && (!isapproxinteger(a.a-b.a) || !isapproxinteger(a.b-b.b)) |
| 58 | + if domainscompatible(a, b) && !(isapproxinteger(a.a-b.a) && isapproxinteger(a.b-b.b)) |
59 | 59 | jac2jac(f,a.a,a.b,b.a,b.b) |
60 | 60 | else |
61 | 61 | defaultcoefficients(f,a,b) |
|
0 commit comments