OpenXM/src/ox_toolkit/cmo.c - diff

Return to cmo.c CVS log

Up to [local] / OpenXM / src / ox_toolkit

Diff for /OpenXM/src/ox_toolkit/cmo.c between version 1.17 and 1.27

-version 1.17, 2003/11/02 03:09:16
+version 1.27, 2016/06/30 01:14:00
 Line 1
 Line 1
 Line 1
  /* -*- mode: C; coding: euc-japan -*- */
- /* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.16 2003/09/18 12:46:08 ohara Exp $ */
+ /* $OpenXM: OpenXM/src/ox_toolkit/cmo.c,v 1.26 2016/06/29 05:07:23 ohara Exp $ */
  /*
     This module includes functions for sending/receiveng CMO's.
-    Some commnets is written in Japanese by the EUC-JP coded
-    character set.
  */
  #include <stdio.h>
  #include <stdlib.h>
  #include <stdarg.h>
  #include <string.h>
+ #include <mpfr.h>
  #include "ox_toolkit.h"
  #include "parse.h"
-Line 165  cmo_list* new_cmo_list()
+Line 164  cmo_list* new_cmo_list()
 Line 165  cmo_list* new_cmo_list()
 Line 164  cmo_list* new_cmo_list()
      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* c = MALLOC(sizeof(cmo_monomial32));
-Line 225  cmo_zz* new_cmo_zz_size(int size)
+Line 244  cmo_zz* new_cmo_zz_size(int size)
 Line 225  cmo_zz* new_cmo_zz_size(int size)
 Line 244  cmo_zz* new_cmo_zz_size(int size)
      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_BIGFLOAT32;
+     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();
+     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_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));
-Line 235  cmo_zero* new_cmo_zero()
+Line 308  cmo_zero* new_cmo_zero()
 Line 235  cmo_zero* new_cmo_zero()
 Line 308  cmo_zero* new_cmo_zero()
  cmo_double *new_cmo_double(double d)
  {
      cmo_double* m = MALLOC_ATOMIC(sizeof(cmo_double));
-     m->tag = CMO_64BIT_MACHINE_DOUBLE;
+     m->tag = CMO_IEEE_DOUBLE_FLOAT;
      m->d = d;
      return m;
  }
-Line 397  char *new_string_set_cmo(cmo *m)
+Line 470  char *new_string_set_cmo(cmo *m)
 Line 397  char *new_string_set_cmo(cmo *m)
 Line 470  char *new_string_set_cmo(cmo *m)
          return new_string_set_cmo_null();
      case CMO_LIST:
          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);
      default:
          ox_printf("unconvertible <%s>\n", get_symbol_by_tag(m->tag));
          /* 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;
+   }
  }

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>