=================================================================== RCS file: /home/cvs/OpenXM/src/ox_toolkit/cmo.c,v retrieving revision 1.22 retrieving revision 1.25 diff -u -p -r1.22 -r1.25 --- OpenXM/src/ox_toolkit/cmo.c 2007/03/14 10:30:54 1.22 +++ OpenXM/src/ox_toolkit/cmo.c 2015/08/18 02:24:04 1.25 @@ -1,5 +1,5 @@ /* -*- mode: C; coding: euc-japan -*- */ -/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.21 2005/10/12 04:03:37 takayama Exp $ */ +/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.24 2015/08/17 05:18:35 noro Exp $ */ /* This module includes functions for sending/receiveng CMO's. @@ -11,6 +11,7 @@ #include #include #include +#include #include "ox_toolkit.h" #include "parse.h" @@ -253,6 +254,21 @@ cmo_qq* new_cmo_qq() return c; } +cmo_bf* new_cmo_bf() +{ + cmo_bf* c = MALLOC(sizeof(cmo_bf)); + c->tag = CMO_BIGFLOAT; + mpfr_init(c->mpfr); + return c; +} + +cmo_complex* new_cmo_complex() +{ + cmo_complex* c = MALLOC(sizeof(cmo_complex)); + c->tag = CMO_COMPLEX; + return c; +} + cmo_qq* new_cmo_qq_set_mpq(mpq_ptr q) { cmo_qq* c = new_cmo_qq(); @@ -268,6 +284,22 @@ cmo_qq* new_cmo_qq_set_mpz(mpz_ptr num, mpz_ptr den) return c; } +cmo_bf* new_cmo_bf_set_mpfr(mpfr_ptr num) +{ + cmo_bf* c = new_cmo_bf(); + mpfr_init2(c->mpfr,num->_mpfr_prec); + mpfr_set(c->mpfr,num,MPFR_RNDN); + return c; +} + +cmo_complex* new_cmo_complex_set_re_im(cmo *re,cmo *im) +{ + cmo_complex* c = new_cmo_complex(); + c->re = re; + c->im = im; + return c; +} + cmo_zero* new_cmo_zero() { cmo_zero* m = MALLOC_ATOMIC(sizeof(cmo_zero)); @@ -448,4 +480,18 @@ char *new_string_set_cmo(cmo *m) /* yet not implemented. */ return NULL; } +} + +int cmo_to_int(cmo *n) +{ + switch(n->tag) { + case CMO_ZERO: + return 0; + case CMO_INT32: + return ((cmo_int32 *)n)->i; + case CMO_ZZ: + return mpz_get_si(((cmo_zz *)n)->mpz); + default: + return 0; + } }