Question: Bug in PolynomialToPDE?

Has anyone else run into the following bug in the PolynomialTools package?

> PolynomialTools[PolynomialToPDE]([z*f], [z], [f]);

                            d
                           [-- f(z)]
                            dz

> PolynomialTools[PolynomialToPDE]([z[1]*f], [z[1]], [f]);

                           [f(z[1])]

I believe the first behavior is correct, but when I use subscripted variables, I get the second behavior for linear terms.  (Nonlinear terms appear to work correctly:

>  PolynomialTools[PolynomialToPDE]([z[1]^2*f], [z[1]], [f]);
                                2
                               d
                             [------ f(z[1])]
                                  2
                              dz[1]

I have observed this in Maple 2024.2 and previous versions.  Does anyone know if it's fixed in Maple 2025?  (I don't have access yet.)  Does anyone know a good workaround for polynomials with both linear and nonlinear terms?

Thanks!

----Josh