version 1.9, 2003/03/23 20:17:34 |
version 1.12, 2003/06/03 16:06:48 |
|
|
/* -*- mode: C; coding: euc-japan -*- */ |
/* -*- mode: C; coding: euc-japan -*- */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.8 2003/02/03 23:13:23 ohara Exp $ */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.11 2003/06/02 10:25:56 ohara Exp $ */ |
|
|
/* |
/* |
This module includes functions for sending/receiveng CMO's. |
This module includes functions for sending/receiveng CMO's. |
|
|
#include "ox_toolkit.h" |
#include "ox_toolkit.h" |
#include "parse.h" |
#include "parse.h" |
|
|
static cell* new_cell(); |
static cell* new_cell(cmo *ob, int e); |
static char* new_string_set_cmo_null(); |
static char* new_string_set_cmo_null(); |
static char* new_string_set_cmo_int32(int integer); |
static char* new_string_set_cmo_int32(int integer); |
static char* new_string_set_cmo_list(cmo_list *c); |
static char* new_string_set_cmo_list(cmo_list *c); |
#if defined(WITH_GMP) |
|
static char* new_string_set_cmo_zz(cmo_zz *c); |
static char* new_string_set_cmo_zz(cmo_zz *c); |
#endif /* WITH_GMP */ |
|
static char* new_string_set_cmo_double(cmo_double *m); |
static char* new_string_set_cmo_double(cmo_double *m); |
|
|
/* functions for a cmo_list */ |
/* functions for a cmo_list */ |
static cell* new_cell(cmo *ob) |
static cell* new_cell(cmo *ob, int e) |
{ |
{ |
cell* h = malloc(sizeof(cell)); |
cell* h = MALLOC(sizeof(cell)); |
h->next = NULL; |
h->next = NULL; |
h->prev = NULL; |
h->prev = NULL; |
h->cmo = ob; |
h->cmo = ob; |
|
h->exp = e; |
return h; |
return h; |
} |
} |
|
|
Line 65 static void list_cons(cell *head, cell *new) |
|
Line 64 static void list_cons(cell *head, cell *new) |
|
|
|
cmo_list *list_append(cmo_list* this, cmo* ob) |
cmo_list *list_append(cmo_list* this, cmo* ob) |
{ |
{ |
list_cons(this->head, new_cell(ob)); |
list_cons(this->head, new_cell(ob, 0)); |
this->length++; |
this->length++; |
return this; |
return this; |
} |
} |
|
|
|
cmo_list *list_append_monomial(cmo_list* this, cmo* coef, int exp) |
|
{ |
|
list_cons(this->head, new_cell(coef, exp)); |
|
this->length++; |
|
return this; |
|
} |
|
|
/* call as list_appendl(List, ob1, ob2, ob3, NULL) */ |
/* call as list_appendl(List, ob1, ob2, ob3, NULL) */ |
cmo_list *list_appendl(cmo_list* this, ...) |
cmo_list *list_appendl(cmo_list* this, ...) |
{ |
{ |
Line 104 cmo *list_nth(cmo_list* this, int n) |
|
Line 110 cmo *list_nth(cmo_list* this, int n) |
|
return NULL; |
return NULL; |
} |
} |
|
|
#if defined(WITH_GMP) |
|
/* for GNU mpz */ |
/* for GNU mpz */ |
void resize_mpz(mpz_ptr mpz, int size) |
void resize_mpz(mpz_ptr mpz, int size) |
{ |
{ |
_mpz_realloc(mpz, abs(size)); |
_mpz_realloc(mpz, abs(size)); |
mpz->_mp_size = size; |
mpz->_mp_size = size; |
} |
} |
#endif /* WITH_GMP */ |
|
|
|
/* functions named new_cmo_*. */ |
/* functions named new_cmo_*. */ |
cmo_null* new_cmo_null() |
cmo_null* new_cmo_null() |
{ |
{ |
cmo_null* m = malloc(sizeof(cmo_null)); |
cmo_null* m = MALLOC(sizeof(cmo_null)); |
m->tag = CMO_NULL; |
m->tag = CMO_NULL; |
return m; |
return m; |
} |
} |
Line 124 cmo_null* new_cmo_null() |
|
Line 128 cmo_null* new_cmo_null() |
|
cmo_int32* new_cmo_int32(int i) |
cmo_int32* new_cmo_int32(int i) |
{ |
{ |
cmo_int32* c; |
cmo_int32* c; |
c = malloc(sizeof(cmo_int32)); |
c = MALLOC(sizeof(cmo_int32)); |
c->tag = CMO_INT32; |
c->tag = CMO_INT32; |
c->i = i; |
c->i = i; |
return c; |
return c; |
Line 132 cmo_int32* new_cmo_int32(int i) |
|
Line 136 cmo_int32* new_cmo_int32(int i) |
|
|
|
cmo_string* new_cmo_string(char* s) |
cmo_string* new_cmo_string(char* s) |
{ |
{ |
cmo_string* c = malloc(sizeof(cmo_string)); |
cmo_string* c = MALLOC(sizeof(cmo_string)); |
c->tag = CMO_STRING; |
c->tag = CMO_STRING; |
if (s != NULL) { |
if (s != NULL) { |
c->s = malloc(strlen(s)+1); |
c->s = MALLOC(strlen(s)+1); |
strcpy(c->s, s); |
strcpy(c->s, s); |
}else { |
}else { |
c->s = NULL; |
c->s = NULL; |
Line 145 cmo_string* new_cmo_string(char* s) |
|
Line 149 cmo_string* new_cmo_string(char* s) |
|
|
|
cmo_mathcap* new_cmo_mathcap(cmo* ob) |
cmo_mathcap* new_cmo_mathcap(cmo* ob) |
{ |
{ |
cmo_mathcap* c = malloc(sizeof(cmo_mathcap)); |
cmo_mathcap* c = MALLOC(sizeof(cmo_mathcap)); |
c->tag = CMO_MATHCAP; |
c->tag = CMO_MATHCAP; |
c->ob = ob; |
c->ob = ob; |
return c; |
return c; |
Line 153 cmo_mathcap* new_cmo_mathcap(cmo* ob) |
|
Line 157 cmo_mathcap* new_cmo_mathcap(cmo* ob) |
|
|
|
cmo_list* new_cmo_list() |
cmo_list* new_cmo_list() |
{ |
{ |
cmo_list* c = malloc(sizeof(cmo_list)); |
cmo_list* c = MALLOC(sizeof(cmo_list)); |
c->tag = CMO_LIST; |
c->tag = CMO_LIST; |
c->length = 0; |
c->length = 0; |
c->head->next = c->head; |
c->head->next = c->head; |
Line 163 cmo_list* new_cmo_list() |
|
Line 167 cmo_list* new_cmo_list() |
|
|
|
cmo_monomial32* new_cmo_monomial32() |
cmo_monomial32* new_cmo_monomial32() |
{ |
{ |
cmo_monomial32* c = malloc(sizeof(cmo_monomial32)); |
cmo_monomial32* c = MALLOC(sizeof(cmo_monomial32)); |
c->tag = CMO_MONOMIAL32; |
c->tag = CMO_MONOMIAL32; |
return c; |
return c; |
} |
} |
Line 173 cmo_monomial32* new_cmo_monomial32_size(int size) |
|
Line 177 cmo_monomial32* new_cmo_monomial32_size(int size) |
|
cmo_monomial32* c = new_cmo_monomial32(); |
cmo_monomial32* c = new_cmo_monomial32(); |
if (size>0) { |
if (size>0) { |
c->length = size; |
c->length = size; |
c->exps = malloc(sizeof(int)*size); |
c->exps = MALLOC(sizeof(int)*size); |
} |
} |
return c; |
return c; |
} |
} |
|
|
#if defined(WITH_GMP) |
|
cmo_zz* new_cmo_zz() |
cmo_zz* new_cmo_zz() |
{ |
{ |
cmo_zz* c = malloc(sizeof(cmo_zz)); |
cmo_zz* c = MALLOC(sizeof(cmo_zz)); |
c->tag = CMO_ZZ; |
c->tag = CMO_ZZ; |
mpz_init(c->mpz); |
mpz_init(c->mpz); |
return c; |
return c; |
Line 189 cmo_zz* new_cmo_zz() |
|
Line 192 cmo_zz* new_cmo_zz() |
|
|
|
cmo_zz* new_cmo_zz_noinit() |
cmo_zz* new_cmo_zz_noinit() |
{ |
{ |
cmo_zz* c = malloc(sizeof(cmo_zz)); |
cmo_zz* c = MALLOC(sizeof(cmo_zz)); |
c->tag = CMO_ZZ; |
c->tag = CMO_ZZ; |
return c; |
return c; |
} |
} |
Line 221 cmo_zz* new_cmo_zz_size(int size) |
|
Line 224 cmo_zz* new_cmo_zz_size(int size) |
|
resize_mpz(c->mpz, size); |
resize_mpz(c->mpz, size); |
return c; |
return c; |
} |
} |
#endif /* WITH_GMP */ |
|
|
|
cmo_zero* new_cmo_zero() |
cmo_zero* new_cmo_zero() |
{ |
{ |
cmo_zero* m = malloc(sizeof(cmo_zero)); |
cmo_zero* m = MALLOC(sizeof(cmo_zero)); |
m->tag = CMO_ZERO; |
m->tag = CMO_ZERO; |
return m; |
return m; |
} |
} |
|
|
cmo_double *new_cmo_double(double d) |
cmo_double *new_cmo_double(double d) |
{ |
{ |
cmo_double* m = malloc(sizeof(cmo_double)); |
cmo_double* m = MALLOC(sizeof(cmo_double)); |
m->tag = CMO_64BIT_MACHINE_DOUBLE; |
m->tag = CMO_64BIT_MACHINE_DOUBLE; |
m->d = d; |
m->d = d; |
return m; |
return m; |
Line 240 cmo_double *new_cmo_double(double d) |
|
Line 242 cmo_double *new_cmo_double(double d) |
|
|
|
cmo_dms_generic* new_cmo_dms_generic() |
cmo_dms_generic* new_cmo_dms_generic() |
{ |
{ |
cmo_dms_generic* m = malloc(sizeof(cmo_dms_generic)); |
cmo_dms_generic* m = MALLOC(sizeof(cmo_dms_generic)); |
m->tag = CMO_DMS_GENERIC; |
m->tag = CMO_DMS_GENERIC; |
return m; |
return m; |
} |
} |
|
|
cmo_ring_by_name* new_cmo_ring_by_name(cmo* ob) |
cmo_ring_by_name* new_cmo_ring_by_name(cmo* ob) |
{ |
{ |
cmo_ring_by_name* c = malloc(sizeof(cmo_ring_by_name)); |
cmo_ring_by_name* c = MALLOC(sizeof(cmo_ring_by_name)); |
c->tag = CMO_RING_BY_NAME; |
c->tag = CMO_RING_BY_NAME; |
c->ob = ob; |
c->ob = ob; |
return c; |
return c; |
Line 255 cmo_ring_by_name* new_cmo_ring_by_name(cmo* ob) |
|
Line 257 cmo_ring_by_name* new_cmo_ring_by_name(cmo* ob) |
|
|
|
cmo_indeterminate* new_cmo_indeterminate(cmo* ob) |
cmo_indeterminate* new_cmo_indeterminate(cmo* ob) |
{ |
{ |
cmo_indeterminate* c = malloc(sizeof(cmo_indeterminate)); |
cmo_indeterminate* c = MALLOC(sizeof(cmo_indeterminate)); |
c->tag = CMO_INDETERMINATE; |
c->tag = CMO_INDETERMINATE; |
c->ob = ob; |
c->ob = ob; |
return c; |
return c; |
Line 263 cmo_indeterminate* new_cmo_indeterminate(cmo* ob) |
|
Line 265 cmo_indeterminate* new_cmo_indeterminate(cmo* ob) |
|
|
|
cmo_distributed_polynomial* new_cmo_distributed_polynomial() |
cmo_distributed_polynomial* new_cmo_distributed_polynomial() |
{ |
{ |
cmo_distributed_polynomial* c = malloc(sizeof(cmo_distributed_polynomial)); |
cmo_distributed_polynomial* c = MALLOC(sizeof(cmo_distributed_polynomial)); |
c->tag = CMO_DISTRIBUTED_POLYNOMIAL; |
c->tag = CMO_DISTRIBUTED_POLYNOMIAL; |
c->length = 0; |
c->length = 0; |
c->head->next = c->head; |
c->head->next = c->head; |
Line 274 cmo_distributed_polynomial* new_cmo_distributed_polyno |
|
Line 276 cmo_distributed_polynomial* new_cmo_distributed_polyno |
|
|
|
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)); |
c->tag = CMO_ERROR2; |
c->tag = CMO_ERROR2; |
c->ob = ob; |
c->ob = ob; |
return c; |
return c; |
} |
} |
|
|
|
|
#if defined(WITH_GMP) |
|
/* Following functions translate cmo's to (asciiz) strings. */ |
/* Following functions translate cmo's to (asciiz) strings. */ |
static char *new_string_set_cmo_zz(cmo_zz *c) |
static char *new_string_set_cmo_zz(cmo_zz *c) |
{ |
{ |
return mpz_get_str(NULL, 10, c->mpz); |
return mpz_get_str(NULL, 10, c->mpz); |
} |
} |
#endif /* WITH_GMP */ |
|
|
|
static char *new_string_set_cmo_null() |
static char *new_string_set_cmo_null() |
{ |
{ |
Line 301 static char *new_string_set_cmo_int32(int integer) |
|
Line 301 static char *new_string_set_cmo_int32(int integer) |
|
char *s; |
char *s; |
|
|
sprintf(buff, "%d", integer); |
sprintf(buff, "%d", integer); |
s = malloc(strlen(buff)+1); |
s = MALLOC(strlen(buff)+1); |
strcpy(s, buff); |
strcpy(s, buff); |
|
|
return s; |
return s; |
Line 313 static char *new_string_set_cmo_list(cmo_list *m) |
|
Line 313 static char *new_string_set_cmo_list(cmo_list *m) |
|
int i; |
int i; |
int size = 0; |
int size = 0; |
int len = list_length(m); |
int len = list_length(m); |
char **sp = malloc(len*sizeof(cmo *)); |
char **sp = ALLOCA(len*sizeof(cmo *)); |
|
|
cell* cp = list_first(m); |
cell* cp = list_first(m); |
for(i = 0; i < len; i++) { |
for(i = 0; i < len; i++) { |
Line 321 static char *new_string_set_cmo_list(cmo_list *m) |
|
Line 321 static char *new_string_set_cmo_list(cmo_list *m) |
|
size += strlen(sp[i]) + 3; |
size += strlen(sp[i]) + 3; |
cp = list_next(cp); |
cp = list_next(cp); |
} |
} |
s = malloc(size+2); |
s = MALLOC(size+2); |
strcpy(s, "[ "); |
strcpy(s, "[ "); |
for(i = 0; i < len - 1; i++) { |
for(i = 0; i < len - 1; i++) { |
strcat(s, sp[i]); |
strcat(s, sp[i]); |
Line 329 static char *new_string_set_cmo_list(cmo_list *m) |
|
Line 329 static char *new_string_set_cmo_list(cmo_list *m) |
|
} |
} |
strcat(s, sp[len-1]); |
strcat(s, sp[len-1]); |
strcat(s, " ]"); |
strcat(s, " ]"); |
free(sp); |
|
return s; |
return s; |
} |
} |
|
|
Line 339 static char *new_string_set_cmo_double(cmo_double *m) |
|
Line 338 static char *new_string_set_cmo_double(cmo_double *m) |
|
char *s; |
char *s; |
|
|
sprintf(buff, "%lf", m->d); |
sprintf(buff, "%lf", m->d); |
s = malloc(strlen(buff)+1); |
s = MALLOC(strlen(buff)+1); |
strcpy(s, buff); |
strcpy(s, buff); |
|
|
return s; |
return s; |
Line 348 static char *new_string_set_cmo_double(cmo_double *m) |
|
Line 347 static char *new_string_set_cmo_double(cmo_double *m) |
|
char *new_string_set_cmo(cmo *m) |
char *new_string_set_cmo(cmo *m) |
{ |
{ |
switch(m->tag) { |
switch(m->tag) { |
#if defined(WITH_GMP) |
|
case CMO_ZZ: |
case CMO_ZZ: |
return new_string_set_cmo_zz((cmo_zz *)m); |
return new_string_set_cmo_zz((cmo_zz *)m); |
#endif /* WITH_GMP */ |
|
case CMO_INT32: |
case CMO_INT32: |
return new_string_set_cmo_int32(((cmo_int32 *)m)->i); |
return new_string_set_cmo_int32(((cmo_int32 *)m)->i); |
case CMO_STRING: |
case CMO_STRING: |