version 1.20, 2005/07/26 12:52:04 |
version 1.24, 2015/08/17 05:18:35 |
|
|
/* -*- mode: C; coding: euc-japan -*- */ |
/* -*- mode: C; coding: euc-japan -*- */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.19 2005/07/20 17:48:03 ohara Exp $ */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.23 2015/08/04 05:24:44 noro Exp $ */ |
|
|
/* |
/* |
This module includes functions for sending/receiveng CMO's. |
This module includes functions for sending/receiveng CMO's. |
|
|
#include <stdlib.h> |
#include <stdlib.h> |
#include <stdarg.h> |
#include <stdarg.h> |
#include <string.h> |
#include <string.h> |
|
#include <mpfr.h> |
#include "ox_toolkit.h" |
#include "ox_toolkit.h" |
#include "parse.h" |
#include "parse.h" |
|
|
Line 110 cmo *list_nth(cmo_list* this, int n) |
|
Line 111 cmo *list_nth(cmo_list* this, int n) |
|
return NULL; |
return NULL; |
} |
} |
|
|
void **list_to_array(cmo_list *c) |
|
{ |
|
int i; |
|
int n = list_length(c); |
|
void **array = MALLOC(sizeof(void *)*(n+1)); |
|
for(i=0; i<n; i++) { |
|
array[i] = list_nth(c, i); |
|
} |
|
array[n] = NULL; |
|
return array; |
|
} |
|
|
|
void **list_to_array_map(cmo_list *c, void *(* mapf)(void *)) |
|
{ |
|
int i; |
|
int n = list_length(c); |
|
void **array = MALLOC(sizeof(void *)*(n+1)); |
|
for(i=0; i<n; i++) { |
|
array[i] = mapf(list_nth(c, i)); |
|
} |
|
array[n] = NULL; |
|
return array; |
|
} |
|
|
|
/* for GNU mpz */ |
/* for GNU mpz */ |
void resize_mpz(mpz_ptr mpz, int size) |
void resize_mpz(mpz_ptr mpz, int size) |
{ |
{ |
Line 189 cmo_list* new_cmo_list() |
|
Line 166 cmo_list* new_cmo_list() |
|
return c; |
return c; |
} |
} |
|
|
cmo_list* new_cmo_list_set(void *array[], int n) |
cmo_list* new_cmo_list_array(void *array[], int n) |
{ |
{ |
int i; |
int i; |
cmo_list* c = new_cmo_list(); |
cmo_list* c = new_cmo_list(); |
Line 199 cmo_list* new_cmo_list_set(void *array[], int n) |
|
Line 176 cmo_list* new_cmo_list_set(void *array[], int n) |
|
return c; |
return c; |
} |
} |
|
|
cmo_list* new_cmo_list_map(void *array[], int n, void *(* mapf)(void *)) |
cmo_list* new_cmo_list_array_map(void *array[], int n, void *(* mapf)(void *)) |
{ |
{ |
int i; |
int i; |
cmo_list* c = new_cmo_list(); |
cmo_list* c = new_cmo_list(); |
Line 269 cmo_zz* new_cmo_zz_size(int size) |
|
Line 246 cmo_zz* new_cmo_zz_size(int size) |
|
return c; |
return c; |
} |
} |
|
|
cmo_qq* new_cmo_qq_noinit() |
cmo_qq* new_cmo_qq() |
{ |
{ |
cmo_qq* c = MALLOC(sizeof(cmo_qq)); |
cmo_qq* c = MALLOC(sizeof(cmo_qq)); |
c->tag = CMO_QQ; |
c->tag = CMO_QQ; |
c->num = c->den = NULL; |
mpq_init(c->mpq); |
return c; |
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_qq* new_cmo_qq_set_mpq(mpq_ptr q) |
cmo_qq* new_cmo_qq_set_mpq(mpq_ptr q) |
{ |
{ |
cmo_qq* c = new_cmo_qq_noinit(); |
cmo_qq* c = new_cmo_qq(); |
c->num = new_cmo_zz_set_mpz(mpq_numref(q)); |
mpq_set(c->mpq, q); |
c->den = new_cmo_zz_set_mpz(mpq_denref(q)); |
|
return c; |
return c; |
} |
} |
|
|
cmo_qq* new_cmo_qq_set_mpz(mpz_ptr num, mpz_ptr den) |
cmo_qq* new_cmo_qq_set_mpz(mpz_ptr num, mpz_ptr den) |
{ |
{ |
cmo_qq* c = new_cmo_qq_noinit(); |
cmo_qq* c = new_cmo_qq(); |
c->num = new_cmo_zz_set_mpz(num); |
mpq_set_num(c->mpq, num); |
c->den = new_cmo_zz_set_mpz(den); |
mpq_set_den(c->mpq, den); |
return c; |
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_zero* new_cmo_zero() |
cmo_zero* new_cmo_zero() |
{ |
{ |
cmo_zero* m = MALLOC_ATOMIC(sizeof(cmo_zero)); |
cmo_zero* m = MALLOC_ATOMIC(sizeof(cmo_zero)); |
Line 473 char *new_string_set_cmo(cmo *m) |
|
Line 465 char *new_string_set_cmo(cmo *m) |
|
/* yet not implemented. */ |
/* yet not implemented. */ |
return NULL; |
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; |
|
} |
} |
} |