-
Notifications
You must be signed in to change notification settings - Fork 74
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Incorrect dot product fallback #397
Comments
At one time, we had |
So you are not implementing the dot product as the integral |
No, right now we don't have it implemented for the |
This should probably get the |
Is it a bug? Or does it need a special method. Right now the fall back is to the So would having
and
suffice? |
The result is incorrect or throws an error. Clearly a bug.
Em ter., 2 de abr. de 2024, 12:27, john verzani ***@***.***>
escreveu:
… Is it a bug? Or does it need a special method. Right now the fall back is
to the dot method for two iterables which will error for polynomials of
different degree. For orthogonal polynomials, we could use the integral
over [-1,1], but in general the integral over (-oo, oo) will always yield
an infinite answer, so that isn't the best.
So would having
Polynomials.dot(p::AbstractPolynomial, q::AbstractPolynomial) = throw(ArgumentError("..."))
and
SpecialPolynomials.dot(p::AbstractOrthogonalPolynomial, q::AbstractOrthogonalPolynomial) = integrate(p * conj(q), domain(p)...)
suffice?
—
Reply to this email directly, view it on GitHub
<#397 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AAZQW3IBATYDM5KZWWH4QH3Y3LE7PAVCNFSM5OXVHNE2U5DIOJSWCZC7NNSXTN2JONZXKZKDN5WW2ZLOOQ5TEMBTGIZTMMZXHE3A>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
closed by #563 |
The fallback is producing results that are not mathematically correct:
The text was updated successfully, but these errors were encountered: