version 1.15, 2003/09/15 09:31:41 |
version 1.23, 2015/08/04 05:24:44 |
|
|
/* -*- mode: C; coding: euc-japan -*- */ |
/* -*- mode: C; coding: euc-japan -*- */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.14 2003/08/21 12:44:06 ohara Exp $ */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.22 2007/03/14 10:30:54 ohara 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 165 cmo_list* new_cmo_list() |
|
Line 166 cmo_list* new_cmo_list() |
|
return c; |
return c; |
} |
} |
|
|
|
cmo_list* new_cmo_list_array(void *array[], int n) |
|
{ |
|
int i; |
|
cmo_list* c = new_cmo_list(); |
|
for(i=0; i<n; i++) { |
|
list_append(c, array[i]); |
|
} |
|
return c; |
|
} |
|
|
|
cmo_list* new_cmo_list_array_map(void *array[], int n, void *(* mapf)(void *)) |
|
{ |
|
int i; |
|
cmo_list* c = new_cmo_list(); |
|
for(i=0; i<n; i++) { |
|
list_append(c, (cmo *)mapf(array[i])); |
|
} |
|
return c; |
|
} |
|
|
cmo_monomial32* new_cmo_monomial32() |
cmo_monomial32* new_cmo_monomial32() |
{ |
{ |
cmo_monomial32* c = MALLOC(sizeof(cmo_monomial32)); |
cmo_monomial32* c = MALLOC(sizeof(cmo_monomial32)); |
Line 225 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() |
|
{ |
|
cmo_qq* c = MALLOC(sizeof(cmo_qq)); |
|
c->tag = CMO_QQ; |
|
mpq_init(c->mpq); |
|
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* c = new_cmo_qq(); |
|
mpq_set(c->mpq, q); |
|
return c; |
|
} |
|
|
|
cmo_qq* new_cmo_qq_set_mpz(mpz_ptr num, mpz_ptr den) |
|
{ |
|
cmo_qq* c = new_cmo_qq(); |
|
mpq_set_num(c->mpq, num); |
|
mpq_set_den(c->mpq, 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_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 235 cmo_zero* new_cmo_zero() |
|
Line 295 cmo_zero* new_cmo_zero() |
|
cmo_double *new_cmo_double(double d) |
cmo_double *new_cmo_double(double d) |
{ |
{ |
cmo_double* m = MALLOC_ATOMIC(sizeof(cmo_double)); |
cmo_double* m = MALLOC_ATOMIC(sizeof(cmo_double)); |
m->tag = CMO_64BIT_MACHINE_DOUBLE; |
m->tag = CMO_IEEE_DOUBLE_FLOAT; |
m->d = d; |
m->d = d; |
return m; |
return m; |
} |
} |
Line 294 cmo_recursive_polynomial* new_cmo_recursive_polynomial |
|
Line 354 cmo_recursive_polynomial* new_cmo_recursive_polynomial |
|
return c; |
return c; |
} |
} |
|
|
|
cmo_tree* new_cmo_tree(cmo_string* name, cmo_list* attributes, cmo_list* leaves) |
|
{ |
|
cmo_tree* c = MALLOC(sizeof(cmo_tree)); |
|
c->tag = CMO_TREE; |
|
c->name= name; |
|
c->attributes = attributes; |
|
c->leaves = leaves; |
|
return c; |
|
} |
|
|
|
cmo_lambda* new_cmo_lambda(cmo_list* args, cmo_tree* body) |
|
{ |
|
cmo_lambda* c = MALLOC(sizeof(cmo_lambda)); |
|
c->tag = CMO_LAMBDA; |
|
c->args = args; |
|
c->body = body; |
|
return c; |
|
} |
|
|
cmo_error2* new_cmo_error2(cmo* ob) |
cmo_error2* new_cmo_error2(cmo* ob) |
{ |
{ |
cmo_error2* c = MALLOC(sizeof(cmo_error2)); |
cmo_error2* c = MALLOC(sizeof(cmo_error2)); |
Line 347 static char *new_string_set_cmo_list(cmo_list *m) |
|
Line 426 static char *new_string_set_cmo_list(cmo_list *m) |
|
strcat(s, sp[i]); |
strcat(s, sp[i]); |
strcat(s, " , "); |
strcat(s, " , "); |
} |
} |
strcat(s, sp[len-1]); |
if (len > 0) |
|
strcat(s, sp[len-1]); |
strcat(s, " ]"); |
strcat(s, " ]"); |
return s; |
return s; |
} |
} |
Line 377 char *new_string_set_cmo(cmo *m) |
|
Line 457 char *new_string_set_cmo(cmo *m) |
|
return new_string_set_cmo_null(); |
return new_string_set_cmo_null(); |
case CMO_LIST: |
case CMO_LIST: |
return new_string_set_cmo_list((cmo_list *)m); |
return new_string_set_cmo_list((cmo_list *)m); |
case CMO_64BIT_MACHINE_DOUBLE: |
case CMO_64BIT_MACHINE_DOUBLE: |
|
case CMO_IEEE_DOUBLE_FLOAT: |
return new_string_set_cmo_double((cmo_double *)m); |
return new_string_set_cmo_double((cmo_double *)m); |
default: |
default: |
ox_printf("unconvertible <%s>\n", get_symbol_by_tag(m->tag)); |
ox_printf("unconvertible <%s>\n", get_symbol_by_tag(m->tag)); |