I am looking for q and r as rational functions of y and z, and I am implementing the code for any finite field. I was told earlier on this site, though, that ?Domains was left unfinished, and is sort of a skeleton.
Basically, I am working with elements of F[x.y.z], and I have constructed them using alias(a = RootOf('irreducible polynomial mod p, of degree n')), and then using Expand() mod p and collect() to perform basic operations. For this part of the code, though, I want to treat them as being polynomials in x, and perform polynomial division. Is there an easier way to do this using Domains? I've been trying PRem, and it doesn't seem to be doing what I would like in most situations.
Thanks for your help!