Mike Rosing on Mon, 23 Feb 2015 22:07:26 +0100 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Now what am I doing wrong? |
I have a really simple program: #include <pari/pari.h> int main() { GEN f, num, sum; int i, k; pari_init(2*1024*1024, 5*1024*512); f = factor(stoi(20765)); pari_printf("%Ps\n", f); num = gel(f, 2); sum = gen_0; k =lg(num); for(i=1; i<=k; i++) { sum = gadd(sum, gel(num, i)); pari_printf("%d %Pd\n", i, sum); } return 0; } which dies at the end of the for loop: $ ./bug-dyn [5, 1; 4153, 1] 1 1 2 2 *** bug in PARI/GP (Segmentation Fault), please report. *** Error in the PARI system. End of program.The value of sum is always correct, but exiting the for loop crashes the program. What am I doing wrong?
Thanks!