Return to kerntest.c CVS log

Up to [local] / OpenXM_contrib / pari-2.2 / src / test

Annotation of OpenXM_contrib/pari-2.2/src/test/kerntest.c, Revision 1.1

1.1     ! noro        1: #include <pari.h>
        !             2: #include <anal.h>
        !             3:
        !             4: GEN   gzero, gun, gdeux;
        !             5: ulong top, bot, avma, memused = 0;
        !             6: long  DEBUGLEVEL,DEBUGMEM = 0;
        !             7:
        !             8: void specinit()
        !             9: {
        !            10:   long size = 100000L;
        !            11:   bot = (long)malloc(size);
        !            12:   top = avma = bot + size;
        !            13:   gzero = malloc(2);
        !            14:   gzero[0] = evaltyp(t_INT)|evallg(2);
        !            15:   gzero[1] = evallgefint(2);
        !            16:   gun   = stoi(1);
        !            17:   gdeux = stoi(2);
        !            18: }
        !            19:
        !            20: void pari_err(long x, ...) { exit (0); }
        !            21: GEN gerepileuptoint(long av, GEN q){ return q; }
        !            22: void gerepilemanysp(long av, long tetpil, GEN* gptr[], long n){}
        !            23: void gerepilemany(long av, GEN* gptr[], long n){}
        !            24:
        !            25: void sorstring(char* b, long x)
        !            26: {
        !            27: #ifdef LONG_IS_64BIT
        !            28:   printf(b,(ulong)x>>32,x & MAXHALFULONG);
        !            29: #else
        !            30:   printf(b,x);
        !            31: #endif
        !            32: }
        !            33:
        !            34: void _voir(GEN x)
        !            35: {
        !            36:   long tx = typ(x), lx = lg(x),i;
        !            37:   sorstring(VOIR_STRING2,x[0]);
        !            38:   if (tx==t_INT) lx = lgefint(x);
        !            39:   for (i=1; i < lx; i++) sorstring(VOIR_STRING2,x[i]);
        !            40:   printf("\n");
        !            41: }
        !            42:
        !            43: int main()
        !            44: {
        !            45:   GEN x,y,r,z, xr,yr;
        !            46:   specinit();
        !            47:   x = stoi(187654321L);
        !            48:   y = stoi(-12345678L);
        !            49:   printf("INT: %ld\n", itos(x));
        !            50:   printf("conv:"); _voir(x);
        !            51:   printf("+:"); _voir(addii(x,y));
        !            52:   printf("-:"); _voir(subii(x,y));
        !            53:   printf("*:"); _voir(mulii(x,y));
        !            54:   printf("/:"); _voir(dvmdii(x,y, &z));
        !            55:   printf("rem:"); _voir(z);
        !            56:   printf("pow:\n");
        !            57:   z = mulii(x,x); _voir(z);
        !            58:   z = mulii(z,z); _voir(z);
        !            59:   z = mulii(z,z); _voir(z);
        !            60:   z = mulii(z,z); _voir(z);
        !            61:   z = mulii(z,z); _voir(z);
        !            62:   printf("invmod:"); invmod(y,z,&r); _voir(r);
        !            63:   xr = cgetr(DEFAULTPREC); affir(x, xr);
        !            64:   yr = cgetr(DEFAULTPREC); affir(y, yr);
        !            65:   printf("\nREAL: %f\n", rtodbl(xr));
        !            66:   printf("conv1:"); _voir(xr);
        !            67:   printf("conv2:"); _voir(dbltor(rtodbl(xr)));
        !            68:   printf("+:"); _voir(addrr(xr,yr));
        !            69:   printf("-:"); _voir(subrr(xr,yr));
        !            70:   printf("*:"); _voir(mulrr(xr,yr));
        !            71:   printf("/:"); _voir(divrr(xr,yr));
        !            72:   printf("gcc bug?:"); _voir(divrs(dbltor(3.),2));
        !            73:   return 0;
        !            74: }