|
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