OpenXM/src/ox_toolkit/ox.c - diff

Return to ox.c CVS log

Up to [local] / OpenXM / src / ox_toolkit

Diff for /OpenXM/src/ox_toolkit/ox.c between version 1.38 and 1.48

version 1.38, 2015/08/04 05:24:44

version 1.48, 2016/06/30 01:14:00

Line 1

/* -*- mode: C; coding: euc-japan -*- */

/* $OpenXM: OpenXM/src/ox_toolkit/ox.c,v 1.37 2014/04/07 04:00:10 iwane Exp $ */

/* $OpenXM: OpenXM/src/ox_toolkit/ox.c,v 1.47 2016/06/29 05:07:23 ohara Exp $ */

This module includes functions for sending/receiveng CMO's.

Some commnets is written in Japanese by the EUC-JP coded

character set.

#include <stdio.h>

#include <stdlib.h>

#include <stdarg.h>

#include <string.h>

#include <unistd.h>

#include <fcntl.h>

#include <sys/file.h>

#include <time.h>

#include <limits.h>

#if !defined(_MSC_VER)

#include <unistd.h>

#include <sys/file.h>

#endif

#include <mpfr.h>

/* XXX : defined in mpfr-impl.h */

#define MPFR_PREC(x) ((x)->_mpfr_prec)

#define MPFR_EXP(x) ((x)->_mpfr_exp)

#define MPFR_MANT(x) ((x)->_mpfr_d)

#define MPFR_LAST_LIMB(x) ((MPFR_PREC (x) - 1) / GMP_NUMB_BITS)

#define RAT_CEIL(nm,dn) (((nm)+(dn)-1)/((dn)))

#define MPFR_LIMB_SIZE(x) (MPFR_LAST_LIMB (x) + 1)

#define MPFR_LIMB_SIZE_REAL(x) (RAT_CEIL(MPFR_PREC((x)),sizeof(mp_limb_t)*CHAR_BIT) * (sizeof(mp_limb_t)/sizeof(int)) )

#define MPFR_LIMB_SIZE_BODY(x) (RAT_CEIL(MPFR_PREC((x)),sizeof(unsigned int)*CHAR_BIT))

#if SIZEOF_LONG==4

typedef long long L64;

Line 70 static cmo_zz* receive_cmo_zz(OXFILE *oxfp);

Line 72 static cmo_zz* receive_cmo_zz(OXFILE *oxfp);

static void receive_mpz(OXFILE *oxfp, mpz_ptr mpz);

static int send_cmo_zz(OXFILE *oxfp, cmo_zz* c);

static int send_mpz(OXFILE *oxfp, mpz_ptr mpz);

static cmo_bf* receive_cmo_bf(OXFILE *oxfp);

static void receive_mpfr(OXFILE *oxfp, mpfr_ptr mpfr);

static int send_cmo_bf(OXFILE *oxfp, cmo_bf* c);

static int send_mpfr(OXFILE *oxfp, mpfr_ptr mpfr);

/* hook functions. (yet not implemented) */

Line 244 static cmo_bf* receive_cmo_bf(OXFILE *oxfp)

Line 248 static cmo_bf* receive_cmo_bf(OXFILE *oxfp)

return new_cmo_bf_set_mpfr(num);

}

static cmo_complex* receive_cmo_complex(OXFILE *oxfp)

{

cmo *re, *im;

re = receive_cmo(oxfp);

im = receive_cmo(oxfp);

return new_cmo_complex_set_re_im(re,im);

}

static cmo_zero* receive_cmo_zero(OXFILE *oxfp)

{

return new_cmo_zero();

Line 368 cmo *receive_cmo_tag(OXFILE *oxfp, int tag)

Line 382 cmo *receive_cmo_tag(OXFILE *oxfp, int tag)

case CMO_QQ:

m = (cmo *)receive_cmo_qq(oxfp);

break;

case CMO_BIGFLOAT:

case CMO_BIGFLOAT32:

m = (cmo *)receive_cmo_bf(oxfp);

break;

case CMO_COMPLEX:

m = (cmo *)receive_cmo_complex(oxfp);

break;

case CMO_ZERO:

m = (cmo *)receive_cmo_zero(oxfp);

break;

Line 656 static int send_cmo_bf(OXFILE *oxfp, cmo_bf* c)

Line 673 static int send_cmo_bf(OXFILE *oxfp, cmo_bf* c)

return 0;

}

static int send_cmo_complex(OXFILE *oxfp, cmo_complex* c)

{

send_cmo(oxfp, c->re);

send_cmo(oxfp, c->im);

return 0;

}

static int send_cmo_recursive_polynomial(OXFILE *oxfp, cmo_recursive_polynomial* c)

{

send_cmo(oxfp, (cmo *)c->ringdef);

Line 722 void send_cmo(OXFILE *oxfp, cmo* c)

Line 746 void send_cmo(OXFILE *oxfp, cmo* c)

case CMO_QQ:

send_cmo_qq(oxfp, (cmo_qq *)c);

break;

case CMO_BIGFLOAT:

case CMO_BIGFLOAT32:

send_cmo_bf(oxfp, (cmo_bf *)c);

break;

case CMO_COMPLEX:

send_cmo_complex(oxfp, (cmo_complex *)c);

break;

case CMO_DISTRIBUTED_POLYNOMIAL:

send_cmo_distributed_polynomial(oxfp, (cmo_distributed_polynomial *)c);

break;

Line 756 static int send_mpz(OXFILE *oxfp, mpz_ptr mpz)

Line 783 static int send_mpz(OXFILE *oxfp, mpz_ptr mpz)

int len = abs(mpz->_mp_size) * n;

int *ptr = (int *)mpz->_mp_d;

int size;

send_int32(oxfp, mpz->_mp_size * n);

#if 0

if (len > 0 && ptr[len-1] == 0) {

len--;

}

size = mpz->_mp_size < 0 ? -len : len;

send_int32(oxfp, size);

#endif

for(i=0; i<len; i++) {

send_int32(oxfp, ptr[i]);

}

return 0;

}

void send_int64(OXFILE *oxfp,UL64 a)

int send_int64(OXFILE *oxfp,UL64 a)

{

send_int32(oxfp, a>>32);

return oxfp->send_double(oxfp,((double *)&a)[0]);

send_int32(oxfp, a&0xffffffff);

}

UL64 receive_int64(OXFILE *oxfp)

{

UL64 u,l;

double d = receive_double(oxfp);

return ((UL64 *)&d)[0];

u = receive_int32(oxfp);

l = receive_int32(oxfp);

return (u<<32)|l;

}

static void receive_mpfr(OXFILE *oxfp, mpfr_ptr mpfr)

{

int sgn,prec,len,i;

int sgn,prec,exp,len_r,len,i;

long *ptr;

unsigned int *ptr;

L64 exp;

sgn = receive_int32(oxfp);

prec = receive_int32(oxfp);

sgn = receive_int32(oxfp);

exp = receive_int64(oxfp);

exp = receive_int32(oxfp);

/* len = length as an int array (int = 4bytes) */

len = receive_int32(oxfp);

mpfr_init2(mpfr,prec); /* initialized by NaN */

MPFR_SIGN(mpfr) = sgn;

mpfr_init2(mpfr,prec);

MPFR_EXP(mpfr) = exp;

MPFR_SIGN(mpfr) = sgn;

*(MPFR_MANT(mpfr))=0; /* initialized by 0 */

ptr = mpfr->_mpfr_d;

ptr = (unsigned int *)MPFR_MANT(mpfr);

mpfr->_mpfr_exp = exp;

len_r = MPFR_LIMB_SIZE_REAL(mpfr);

#if SIZEOF_LONG==4

for(i=(len_r-len); i<len_r; i++) {

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

ptr[i] = receive_int32(oxfp);

}

#else

len >>= 1;

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

ptr[i] = receive_int64(oxfp);

#endif

}

static int send_mpfr(OXFILE *oxfp, mpfr_ptr mpfr)

{

int i,len_r,len;

unsigned int *ptr;

int i,len;

send_int32(oxfp, MPFR_PREC(mpfr));

long *ptr;

send_int32(oxfp, MPFR_SIGN(mpfr));

send_int32(oxfp, MPFR_EXP(mpfr));

send_int32(oxfp, MPFR_SIGN(mpfr));

len_r = MPFR_LIMB_SIZE_REAL(mpfr);

send_int32(oxfp, MPFR_PREC(mpfr));

len = MPFR_LIMB_SIZE_BODY(mpfr);

send_int64(oxfp, MPFR_EXP(mpfr));

ptr = (unsigned int *)MPFR_MANT(mpfr);

len = MPFR_LIMB_SIZE(mpfr);

send_int32(oxfp, len);

ptr = MPFR_MANT(mpfr);

for(i=(len_r-len); i<len_r; i++) {

#if SIZEOF_LONG==4

send_int32(oxfp, ptr[i]);

send_int32(oxfp, len);

}

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

send_int32(oxfp,ptr[i]);

#else /* SIZEOF_LONG==8 */

send_int32(oxfp, 2*len);

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

send_int64(oxfp,ptr[i]);

#endif

return 0;

}

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