tests pass

Author: dlfivefifty

src/classical/jacobi.jl

@@ -83,6 +83,8 @@ jacobi(a,b, d::AbstractInterval{T}) where T = Jacobi{float(promote_type(eltype(a
 
 Jacobi(P::Legendre{T}) where T = Jacobi(zero(T), zero(T))
 
+_basis(w::JacobiWeight) = Weighted(Jacobi(w.a, w.b))
+
 """
      jacobip(n, a, b, z)
 

src/classical/legendre.jl

@@ -178,4 +178,10 @@ function __dot(::ExpansionLayout, ::ExpansionLayout, a, b)
     Q,d = basis(b),coefficients(b)
     c' * (P'Q) * d
 end
-__dot(_, _, a, b) = dot(expand(a), expand(b))
+__dot(_, _, a, b) = dot(expand(a), expand(b))
+__dot(::ZerosLayout, _, a, b) = __dot(ZerosLayout(), ZerosLayout(), a, b)
+__dot(_, ::ZerosLayout, a, b) = __dot(ZerosLayout(), ZerosLayout(), a, b)
+function __dot(::ZerosLayout, ::ZerosLayout, a, b)
+    axes(a) == axes(b) || throw(DimensionMismatch("axes must match"))
+    zero(promote_type(eltype(a), eltype(b)))
+end