[BACK]Return to cmo.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / ox_toolkit

Diff for /OpenXM/src/ox_toolkit/cmo.c between version 1.20 and 1.26

version 1.20, 2005/07/26 12:52:04 version 1.26, 2016/06/29 05:07:23
Line 1 
Line 1 
 /* -*- 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.25 2015/08/18 02:24:04 noro Exp $ */
   
 /*  /*
    This module includes functions for sending/receiveng CMO's.     This module includes functions for sending/receiveng CMO's.
Line 11 
Line 11 
 #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_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* 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_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* 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 480  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;
     }
 }  }
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.26
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>