So, after fighting it with a bit, I realized that the reason Prem wasn't working properly was because sometimes I would write xy instead of x*y or x y, which I'm assuming it interpreted as a whole new variable. It was really weird, because it kept saying things like r was 0, when my program was ensuring a gcd of 1 the line before.

tl;dr version: Prem is working for me now, and it was my own error that prevented it from working properly before. Thank you very much for your help!

So, after fighting it with a bit, I realized that the reason Prem wasn't working properly was because sometimes I would write xy instead of x*y or x y, which I'm assuming it interpreted as a whole new variable. It was really weird, because it kept saying things like r was 0, when my program was ensuring a gcd of 1 the line before.

tl;dr version: Prem is working for me now, and it was my own error that prevented it from working properly before. Thank you very much for your help!

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!

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!

Thanks very much! I had looked through that thread a bit, but apparently did not follow it through long enough - sorry about that!

Thanks very much! I had looked through that thread a bit, but apparently did not follow it through long enough - sorry about that!