Open NAThompson opened 3 months ago
Negative exponents is currently not supported in this package, e.g., it will make printing fail but it's probably not too hard to add support for it. Then you could easily do sum(o[i, j] * x[i] * y[j] for i in axes(o, 1), j in axes(o, 2))
In Polynomials.jl, we can use an
OffsetArray
to easily construct Laurent polynomials:However, I cannot do the same for an
OffsetMatrix
withMultivariatePolynomials
. Reading through the documentation, it appears the this library is more focuses on very sophisticated manipulations on polynomials that you could type out by hand, rather than large polynomials whose coefficients are stored in a matrix-so hopefully this issue is not too annoying. However, I felt like it might be natural to type something like the following, as well as increase compatibility with the univariate case: