version 1.7, 2003/01/17 07:40:10 |
version 1.14, 2003/08/21 12:44:06 |
|
|
/* -*- mode: C; coding: euc-japan -*- */ |
/* -*- mode: C; coding: euc-japan -*- */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.6 2003/01/17 06:49:53 ohara Exp $ */ |
/* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.13 2003/06/05 21:12:07 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); |
Line 22 static char* new_string_set_cmo_zz(cmo_zz *c); |
|
Line 22 static char* new_string_set_cmo_zz(cmo_zz *c); |
|
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 63 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 112 void resize_mpz(mpz_ptr mpz, int size) |
|
Line 120 void resize_mpz(mpz_ptr mpz, int size) |
|
/* 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_ATOMIC(sizeof(cmo_null)); |
m->tag = CMO_NULL; |
m->tag = CMO_NULL; |
return m; |
return m; |
} |
} |
Line 120 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_ATOMIC(sizeof(cmo_int32)); |
c->tag = CMO_INT32; |
c->tag = CMO_INT32; |
c->i = i; |
c->i = i; |
return c; |
return c; |
} |
} |
|
|
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_ATOMIC(strlen(s)+1); |
strcpy(c->s, s); |
strcpy(c->s, s); |
}else { |
}else { |
c->s = NULL; |
c->s = NULL; |
Line 141 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 149 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 159 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 169 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; |
} |
} |
|
|
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 184 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 219 cmo_zz* new_cmo_zz_size(int size) |
|
Line 227 cmo_zz* new_cmo_zz_size(int size) |
|
|
|
cmo_zero* new_cmo_zero() |
cmo_zero* new_cmo_zero() |
{ |
{ |
cmo_zero* m = malloc(sizeof(cmo_zero)); |
cmo_zero* m = MALLOC_ATOMIC(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_ATOMIC(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 234 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_ATOMIC(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 249 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 257 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 266 cmo_distributed_polynomial* new_cmo_distributed_polyno |
|
Line 274 cmo_distributed_polynomial* new_cmo_distributed_polyno |
|
return c; |
return c; |
} |
} |
|
|
|
cmo_polynomial_in_one_variable* new_cmo_polynomial_in_one_variable(int var) |
|
{ |
|
cmo_polynomial_in_one_variable* c = MALLOC(sizeof(cmo_polynomial_in_one_variable)); |
|
c->tag = CMO_POLYNOMIAL_IN_ONE_VARIABLE; |
|
c->length = 0; |
|
c->head->next = c->head; |
|
c->head->prev = c->head; |
|
c->var = var; |
|
return c; |
|
} |
|
|
|
cmo_recursive_polynomial* new_cmo_recursive_polynomial(cmo* ringdef, cmo* coef) |
|
{ |
|
cmo_recursive_polynomial* c = MALLOC(sizeof(cmo_recursive_polynomial)); |
|
c->tag = CMO_RECURSIVE_POLYNOMIAL; |
|
c->ringdef = ringdef; |
|
c->coef = coef; |
|
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)); |
c->tag = CMO_ERROR2; |
c->tag = CMO_ERROR2; |
c->ob = ob; |
c->ob = ob; |
return c; |
return c; |
Line 293 static char *new_string_set_cmo_int32(int integer) |
|
Line 321 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 305 static char *new_string_set_cmo_list(cmo_list *m) |
|
Line 333 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 313 static char *new_string_set_cmo_list(cmo_list *m) |
|
Line 341 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 321 static char *new_string_set_cmo_list(cmo_list *m) |
|
Line 349 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 330 static char *new_string_set_cmo_double(cmo_double *m) |
|
Line 357 static char *new_string_set_cmo_double(cmo_double *m) |
|
char buff[1024]; |
char buff[1024]; |
char *s; |
char *s; |
|
|
sprintf(buff, "%lf", m->d); |
sprintf(buff, "%.20f", m->d); |
s = malloc(strlen(buff)+1); |
s = MALLOC_ATOMIC(strlen(buff)+1); |
strcpy(s, buff); |
strcpy(s, buff); |
|
|
return s; |
return s; |
Line 339 static char *new_string_set_cmo_double(cmo_double *m) |
|
Line 366 static char *new_string_set_cmo_double(cmo_double *m) |
|
|
|
char *new_string_set_cmo(cmo *m) |
char *new_string_set_cmo(cmo *m) |
{ |
{ |
symbol_t symp; |
|
switch(m->tag) { |
switch(m->tag) { |
case CMO_ZZ: |
case CMO_ZZ: |
return new_string_set_cmo_zz((cmo_zz *)m); |
return new_string_set_cmo_zz((cmo_zz *)m); |
Line 354 char *new_string_set_cmo(cmo *m) |
|
Line 380 char *new_string_set_cmo(cmo *m) |
|
case CMO_64BIT_MACHINE_DOUBLE: |
case CMO_64BIT_MACHINE_DOUBLE: |
return new_string_set_cmo_int32(m); |
return new_string_set_cmo_int32(m); |
default: |
default: |
#ifdef DEBUG |
ox_printf("unconvertible <%s>\n", get_symbol_by_tag(m->tag)); |
symp = lookup_by_tag(m->tag); |
|
ox_printf("I do not know how to convert %s to a string.\n", symp->key); |
|
#endif |
|
/* yet not implemented. */ |
/* yet not implemented. */ |
return NULL; |
return NULL; |
} |
} |