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

Diff for /OpenXM/src/ox_gsl/ox_gsl.c between version 1.3 and 1.11

version 1.3, 2018/03/30 04:43:16 version 1.11, 2018/06/06 07:40:32
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/ox_gsl/ox_gsl.c,v 1.2 2018/03/29 11:52:18 takayama Exp $  /* $OpenXM: OpenXM/src/ox_gsl/ox_gsl.c,v 1.10 2018/06/04 06:39:26 ohara Exp $
 */  */
   
 #include <stdio.h>  #include <stdio.h>
Line 6 
Line 6 
 #include <setjmp.h>  #include <setjmp.h>
 #include <string.h>  #include <string.h>
 #include <unistd.h>  #include <unistd.h>
   #include <signal.h>
 #include <math.h>  #include <math.h>
 #include "ox_gsl.h"  #include "ox_gsl.h"
 #include "call_gsl.h" // need only when you bind call_gsl functions.  #include "call_gsl.h" // need only when you bind call_gsl functions.
Line 40  void gc_free(void *p,size_t size)
Line 41  void gc_free(void *p,size_t size)
 void init_gc()  void init_gc()
 { GC_INIT();  { GC_INIT();
   mp_set_memory_functions(GC_malloc,gc_realloc,gc_free);    mp_set_memory_functions(GC_malloc,gc_realloc,gc_free);
     init_dic();  // initialize ox_eval.c
 }  }
   
 void initialize_stack()  void initialize_stack()
Line 110  int sm_mathcap()
Line 112  int sm_mathcap()
 //    CMO_DISTRIBUTED_POLYNOMIAL,  //    CMO_DISTRIBUTED_POLYNOMIAL,
 //    CMO_RECURSIVE_POLYNOMIAL,  //    CMO_RECURSIVE_POLYNOMIAL,
 //    CMO_POLYNOMIAL_IN_ONE_VARIABLE,  //    CMO_POLYNOMIAL_IN_ONE_VARIABLE,
       CMO_TREE,
     CMO_ERROR2,      CMO_ERROR2,
     0};      0};
   int available_sm_command[]={    int available_sm_command[]={
Line 154  cmo *make_error2(const char *reason,const char *fname,
Line 157  cmo *make_error2(const char *reason,const char *fname,
     n = 5;      n = 5;
     argv = (cmo **) GC_malloc(sizeof(cmo *)*n);      argv = (cmo **) GC_malloc(sizeof(cmo *)*n);
     ms = (cmo *)new_cmo_string("Error"); argv[0] = ms;      ms = (cmo *)new_cmo_string("Error"); argv[0] = ms;
     if (reason != NULL) s = (char *)GC_malloc(strlen(reason)+1);      if (reason != NULL) {s = (char *)GC_malloc(strlen(reason)+1); strcpy(s,reason);
     else strcpy(s,"");      }else strcpy(s,"");
     ms = (cmo *) new_cmo_string(s); argv[1] = ms;      ms = (cmo *) new_cmo_string(s); argv[1] = ms;
     if (reason != NULL) s = (char *)GC_malloc(strlen(fname)+1);      if (fname != NULL) {s = (char *)GC_malloc(strlen(fname)+1); strcpy(s,fname);
     else strcpy(s,"");      }else strcpy(s,"");
     ms = (cmo *) new_cmo_string(s); argv[2] = ms;      ms = (cmo *) new_cmo_string(s); argv[2] = ms;
     err = (cmo *)new_cmo_int32(line); argv[3] = err;      err = (cmo *)new_cmo_int32(line); argv[3] = err;
     err = (cmo *)new_cmo_int32(code); argv[4] = err;      err = (cmo *)new_cmo_int32(code); argv[4] = err;
Line 179  int get_i()
Line 182  int get_i()
     }else if (c->tag == CMO_ZERO) {      }else if (c->tag == CMO_ZERO) {
         return(0);          return(0);
     }      }
     myhandler("get_i, not an integer",NULL,0,-1);      myhandler("get_i: not an integer",NULL,0,-1);
     return 0;      return 0;
 }  }
   
Line 208  double get_double()
Line 211  double get_double()
     }else if (c->tag == CMO_ZZ) {      }else if (c->tag == CMO_ZZ) {
        if ((mpz_cmp_si(mympz(c),(long int) 0x7fffffff)>0) ||         if ((mpz_cmp_si(mympz(c),(long int) 0x7fffffff)>0) ||
            (mpz_cmp_si(mympz(c),(long int) -0x7fffffff)<0)) {             (mpz_cmp_si(mympz(c),(long int) -0x7fffffff)<0)) {
          myhandler("get_double, out of int32",NULL,0,-1);           myhandler("get_double: out of int32",NULL,0,-1);
          return(NAN);           return(NAN);
        }         }
        return( (double) mpz_get_si(((cmo_zz *)c)->mpz));         return( (double) mpz_get_si(((cmo_zz *)c)->mpz));
Line 217  double get_double()
Line 220  double get_double()
     }else if (c->tag == CMO_ZERO) {      }else if (c->tag == CMO_ZERO) {
         return(0);          return(0);
     }      }
     myhandler("get_double, not a double",NULL,0,-1);      myhandler("get_double: not a double",NULL,0,-1);
     return(NAN);      return(NAN);
 }  }
   
Line 291  char *get_string() {
Line 294  char *get_string() {
   return(NULL);    return(NULL);
 }  }
   
   void test_ox_eval() {
     cmo *c;
     double d=0;
     pop();
     c=pop();
     if (Debug) {
       ox_printf("cmo *c="); print_cmo(c); ox_printf("\n");
     }
     init_dic();
     register_entry("x",1.25);
     if (eval_cmo(c,&d) == 0) myhandler("eval_cmo failed",NULL,0,-1);
     push((cmo *)new_cmo_double(d));
   }
   
 int sm_executeFunction()  int sm_executeFunction()
 {  {
     cmo_string *func = (cmo_string *)pop();      cmo_string *func = (cmo_string *)pop();
Line 298  int sm_executeFunction()
Line 315  int sm_executeFunction()
         push(make_error2("sm_executeFunction, not CMO_STRING",NULL,0,-1));          push(make_error2("sm_executeFunction, not CMO_STRING",NULL,0,-1));
         return -1;          return -1;
     }      }
       init_dic();
     // Test functions      // Test functions
     if (strcmp(func->s, "add_int32") == 0) {      if (strcmp(func->s, "add_int32") == 0) {
         my_add_int32();          my_add_int32();
Line 307  int sm_executeFunction()
Line 325  int sm_executeFunction()
         show_double_list();          show_double_list();
     }else if (strcmp(func->s,"restart")==0) {      }else if (strcmp(func->s,"restart")==0) {
         pop(); restart();          pop(); restart();
       }else if (strcmp(func->s,"test_ox_eval")==0) {
           test_ox_eval();
     // The following functions are defined in call_gsl.c      // The following functions are defined in call_gsl.c
     }else if (strcmp(func->s,"gsl_sf_lngamma_complex_e")==0) {      }else if (strcmp(func->s,"gsl_sf_lngamma_complex_e")==0) {
         call_gsl_sf_lngamma_complex_e();          call_gsl_sf_lngamma_complex_e();
       }else if (strcmp(func->s,"gsl_integration_qags")==0) {
           call_gsl_integration_qags();
       }else if (strcmp(func->s,"gsl_monte_plain_integrate")==0) {
           call_gsl_monte_plain_integrate();
     }else {      }else {
         push(make_error2("sm_executeFunction, unknown function",NULL,0,-1));          push(make_error2("sm_executeFunction, unknown function",NULL,0,-1));
         return -1;          return -1;
Line 382  void myhandler(const char *reason,const char *file,int
Line 406  void myhandler(const char *reason,const char *file,int
   fprintf(fp,"%d\n",line);    fprintf(fp,"%d\n",line);
   if (file != NULL) fprintf(fp,"%s\n",file); else fprintf(fp,"file?\n");    if (file != NULL) fprintf(fp,"%s\n",file); else fprintf(fp,"file?\n");
   if (reason != NULL) fprintf(fp,"%s\n",reason); else fprintf(fp,"reason?\n");    if (reason != NULL) fprintf(fp,"%s\n",reason); else fprintf(fp,"reason?\n");
   fflush(NULL); fclose(fp);  // BUG. the contents are deleted after it is closed.    fflush(NULL); fclose(fp);
   // m = make_error2(reason,file,line,gsl_errno);    // m = make_error2(reason,file,line,gsl_errno);
   //  send_ox_cmo(fd_rw, m);  ox_flush(fd_rw);    //  send_ox_cmo(fd_rw, m);  ox_flush(fd_rw);
   // send error packet even it is not asked. Todo, OK? --> no    // send error packet even it is not asked. Todo, OK? --> no
   restart();    restart();
 }  }
   
 void push_error_from_file() {  void push_error_from_file() {
   FILE *fp;    FILE *fp;
 #define BUF_SIZE 1024  #define BUF_SIZE 1024
Line 400  void push_error_from_file() {
Line 423  void push_error_from_file() {
   cmo *m;    cmo *m;
   fprintf(stderr,"push_error_from_file()\n");    fprintf(stderr,"push_error_from_file()\n");
   sprintf(logname,"/tmp/ox_gsl-%d.txt",(int) getpid());    sprintf(logname,"/tmp/ox_gsl-%d.txt",(int) getpid());
   fp = fopen(logname,"w");    fp = fopen(logname,"r");
   if (fp == NULL) return;    if (fp == NULL) {
       fprintf(stderr,"open %s is failed\n",logname); return;
     }
   fgets(cmd,BUF_SIZE-2,fp); sscanf(cmd,"%d",&gsl_errno);    fgets(cmd,BUF_SIZE-2,fp); sscanf(cmd,"%d",&gsl_errno);
   fgets(cmd,BUF_SIZE-2,fp); sscanf(cmd,"%d",&line);    fgets(cmd,BUF_SIZE-2,fp); sscanf(cmd,"%d",&line);
   fgets(file,BUF_SIZE-2,fp);  #define remove_newline(s) {char *tmp_pos; if ((tmp_pos=strchr(s,'\n')) != NULL) *tmp_pos = '\0';}
   fgets(reason,BUF_SIZE-2,fp);    fgets(file,BUF_SIZE-2,fp);  remove_newline(file);
     fgets(reason,BUF_SIZE-2,fp); remove_newline(reason);
   fclose(fp);    fclose(fp);
   m = make_error2(reason,file,line,gsl_errno);    m = make_error2(reason,file,line,gsl_errno);
   push(m);    push(m);
Line 431  int main()
Line 457  int main()
     fd_rw = oxf_open(3);      fd_rw = oxf_open(3);
     oxf_determine_byteorder_server(fd_rw);      oxf_determine_byteorder_server(fd_rw);
   }    }
   #if defined(__CYGWIN__)
     void *mysignal(int sig,void (*handler)(int m));
     mysignal(SIGUSR1,usr1_handler);
   #else
   signal(SIGUSR1,usr1_handler);    signal(SIGUSR1,usr1_handler);
   #endif
   
   while(1) {    while(1) {
     receive();      receive();
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.11
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>