Phat Tran on Thu, 18 Mar 2010 07:08:18 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: Degree overflow in pow_monome - PARI


Sorry for the inconvenience if this email was sent to you mail box. You can unsubscribe it.

Back to the problem, I found that the function gpowgs() had caused the error. I replaced it with powgi(), but it gave me a wrong answer. (I got a correct answer by using GP). I assigned m, n, b and d 30-digit numbers.
Any help ?

Thanks,

Phat

On Wed, Mar 17, 2010 at 9:32 PM, Phat Tran <minhphat@gmail.com> wrote:
***   degree overflow in pow_monome.
  ***   Error in the PARI system. End of program.

I got the error message above while running the codes in C using PARI/GP 2.3.5 under cygwin.

    GEN m = readseq("123456789");
    GEN n = readseq("1234567");
   
    GEN a = readseq("7");
    GEN c = readseq("8");
   
    GEN b = readseq("10358");
    GEN d = readseq("49850");
   
   
    GEN x = pol_x[fetch_user_var("x")];
    GEN y = pol_x[fetch_user_var("y")];
    GEN z = pol_x[fetch_user_var("a")];
       
    z = gpowgs(gmodulo(gadd(x, gmodulo(gadd(y, gmodulss(1, n)), gsub(gpowgs(y, c), d))), gsub(gpowgs(x, a), b)), m);

I wonder if the number m is too big for PARI to be handled. 

Thanks,

Minhphat