OpenXM/src/ox_pari/convert.c - diff

Return to convert.c CVS log

Up to [local] / OpenXM / src / ox_pari

Diff for /OpenXM/src/ox_pari/convert.c between version 1.5 and 1.7

version 1.5, 2018/03/29 02:14:30

version 1.7, 2021/03/25 07:03:21

Line 1

#include "ox_pari.h"

int with_gmp;

int gmp_check()

{

GEN z1,z2;

// z1=[*,*,1,2], z2=[*,*,2,1]

// if gmp style => z1>z2 else z1<z2

z1 = cgeti(2+2); setsigne(z1,1); setlgefint(z1,lg(z1));

z1[2] = 1; z1[3] = 2;

z2 = cgeti(2+2); setsigne(z2,1); setlgefint(z2,lg(z2));

z2[2] = 2; z2[3] = 1;

if ( cmpii(z1,z2) > 0 ) with_gmp = 1;

else with_gmp = 0;

}

GEN cmo_int32_to_GEN(cmo_int32 *c)

{

GEN z;

Line 42 GEN cmo_zz_to_GEN(cmo_zz *c)

Line 58 GEN cmo_zz_to_GEN(cmo_zz *c)

len = ABSIZ(mpz);

ptr = (long *)PTR(mpz);

z = cgeti(len+2);

for ( j = 0; j < len; j++ )

if ( with_gmp ) // least signifcant first

z[len-j+1] = ptr[j];

for ( j = 0; j < len; j++ )

z[j+2] = ptr[j];

else // most signifcant first

for ( j = 0; j < len; j++ )

z[len-j+1] = ptr[j];

setsigne(z,sgn);

setlgefint(z,lg(z));

return z;

Line 117 GEN cmo_up_to_GEN(cmo_polynomial_in_one_variable *c)

Line 137 GEN cmo_up_to_GEN(cmo_polynomial_in_one_variable *c)

z = cgetg(d+3,t_POL);

setsigne(z,1);

setvarn(z,c->var);

setlgef(z,d+3);

setlg(z,d+3);

for ( i = 2; i <= d+2; i++ )

z[i] = (long)gen_0;

for ( cell = c->head->next; cell != c->head; cell = cell->next ) {

Line 158 cmo_zz *GEN_to_cmo_string(GEN z)

Line 178 cmo_zz *GEN_to_cmo_string(GEN z)

cmo_string *c;

c = new_cmo_string(GSTR(z));

return c;

return (cmo_zz *)c;

}

cmo_zz *GEN_to_cmo_zz(GEN z)

Line 166 cmo_zz *GEN_to_cmo_zz(GEN z)

Line 186 cmo_zz *GEN_to_cmo_zz(GEN z)

cmo_zz *c;

c = new_cmo_zz();

mpz_import(c->mpz,lgef(z)-2,1,sizeof(long),0,0,&z[2]);

// -1:least significant first, 1:most significant first

mpz_import(c->mpz,lg(z)-2,with_gmp?-1:1,sizeof(long),0,0,&z[2]);

if ( signe(z) < 0 )

mpz_neg(c->mpz,c->mpz);

return c;

Line 180 cmo_qq *GEN_to_cmo_qq(GEN z)

Line 201 cmo_qq *GEN_to_cmo_qq(GEN z)

num = (GEN)z[1];

den = (GEN)z[2];

c = new_cmo_qq();

mpz_import(mpq_numref(c->mpq),lgef(num)-2,1,sizeof(long),0,0,&num[2]);

mpz_import(mpq_numref(c->mpq),lg(num)-2,1,sizeof(long),0,0,&num[2]);

mpz_import(mpq_denref(c->mpq),lgef(den)-2,1,sizeof(long),0,0,&den[2]);

mpz_import(mpq_denref(c->mpq),lg(den)-2,1,sizeof(long),0,0,&den[2]);

if ( signe(num)*signe(den) < 0 )

mpz_neg(mpq_numref(c->mpq),mpq_numref(c->mpq));

return c;

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>