[BACK]Return to ntl.cpp CVS log [TXT][DIR] Up to [local] / OpenXM / src / ox_ntl

Diff for /OpenXM/src/ox_ntl/ntl.cpp between version 1.3 and 1.6

version 1.3, 2003/11/15 09:06:20 version 1.6, 2008/09/19 10:55:40
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/ox_ntl/ntl.cpp,v 1.2 2003/11/08 12:34:00 iwane Exp $ */  /* $OpenXM: OpenXM/src/ox_ntl/ntl.cpp,v 1.5 2004/07/04 02:31:51 iwane Exp $ */
   
 #include <NTL/ZZXFactoring.h>  #include <NTL/ZZXFactoring.h>
 #include <NTL/LLL.h>  #include <NTL/LLL.h>
Line 10 
Line 10 
 #define __NTL_PRINT (1)  #define __NTL_PRINT (1)
 #endif  #endif
   
 #define __NTL_PRINT (1)  
 #define DPRINTF(x)      printf x; fflush(stdout)  #define DPRINTF(x)      printf x; fflush(stdout)
   
 /****************************************************************************  /****************************************************************************
Line 25 
Line 24 
          : [2,[[x-1,1],[x^4+x^3+x^2+x+1,1],[x+1,2],[x^2-x+1,2]]]           : [2,[[x-1,1],[x^4+x^3+x^2+x+1,1],[x+1,2],[x^2-x+1,2]]]
  *   *
  ****************************************************************************/   ****************************************************************************/
 cmo *  oxstack_node *
 ntl_fctr(cmo **arg, int argc)  ntl_fctr(oxstack_node **arg, int argc)
 {  {
         cmo *poly = arg[0];          oxstack_node *p = arg[0];
           cmo *poly = p->c;
         cmo_indeterminate *x;          cmo_indeterminate *x;
           cmo *err;
         ZZX f;          ZZX f;
         int ret;          int ret;
   
         if (argc != 1) {          if (argc != 1) {
                 return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(#)")));                  err = ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(#)")));
                   p = oxstack_node_init(err);
                   return (p);
         }          }
   
         ret = cmo_to_ZZX(f, poly, x);          ret = cmo_to_ZZX(f, poly, x);
         if (ret != NTL_SUCCESS) {          if (ret != NTL_SUCCESS) {
                 /* format error */                  /* format error */
                 return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(type)")));                  err = (cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(type)"));
                   p = oxstack_node_init(err);
                   return (p);
         }          }
   
 #if __NTL_PRINT  #if __NTL_PRINT
Line 56  ntl_fctr(cmo **arg, int argc)
Line 61  ntl_fctr(cmo **arg, int argc)
         cout << "fctr : " << *factors->f << endl;          cout << "fctr : " << *factors->f << endl;
 #endif  #endif
   
           p = oxstack_node_init(NULL);
           p->p = factors;
   
         return ((cmo *)factors);          return (p);
 }  }
   
 /****************************************************************************  /****************************************************************************
Line 72  ntl_fctr(cmo **arg, int argc)
Line 79  ntl_fctr(cmo **arg, int argc)
          : [2,[[x-1,1],[x^4+x^3+x^2+x+1,1],[x+1,2],[x^2-x+1,2]]]           : [2,[[x-1,1],[x^4+x^3+x^2+x+1,1],[x+1,2],[x^2-x+1,2]]]
  *   *
  ****************************************************************************/   ****************************************************************************/
 cmo *  oxstack_node *
 ntl_lll(cmo **arg, int argc)  ntl_lll(oxstack_node **arg, int argc)
 {  {
         cmo *poly = arg[0];          oxstack_node *p = arg[0];
         cmo_indeterminate *x;          cmo *err;
         ZZX f;          ZZX f;
         int ret;          int ret;
         cmon_mat_zz_t *mat;          cmon_mat_zz_t *mat;
   
 DPRINTF(("lll start\n"));  
   
         if (argc != 1) {          if (argc != 1) {
 DPRINTF(("invalid argc %d\n", argc));                  err = ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(#)")));
                 return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(#)")));                  p = oxstack_node_init(err);
                   return (p);
         }          }
   
         mat = new_cmon_mat_zz();          mat = new_cmon_mat_zz();
         ret = cmo_to_mat_zz(*mat->mat, arg[0]);          ret = cmo_to_mat_zz(*mat->mat, p->c);
 DPRINTF(("ret = %d, convert\n", ret));  
         if (ret != NTL_SUCCESS) {          if (ret != NTL_SUCCESS) {
 DPRINTF(("convert failed\n", ret));  
                 delete_cmon_mat_zz(mat);                  delete_cmon_mat_zz(mat);
                 /* format error */                  /* format error */
                 return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(type)")));                  err = ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(type)")));
                   p = oxstack_node_init(err);
                   return (p);
         }          }
   
 DPRINTF(("convert success\n", ret));  
   
 #if __NTL_PRINT  #if __NTL_PRINT
 DPRINTF(("convert success\n", ret));  
         cout << "input: " << (*mat->mat) << endl;          cout << "input: " << (*mat->mat) << endl;
 #endif  #endif
   
   
 DPRINTF(("lll start\n", ret));  
         ZZ det2;          ZZ det2;
         long rd = LLL(det2, *mat->mat);          mat_ZZ U;
   
   
 #if __NTL_PRINT  #if __NTL_PRINT
         cout << "output: " << (*mat->mat) << endl;          cout << "output: " << (*mat->mat) << endl;
           cout << U << endl;
 #endif  #endif
   
         return ((cmo *)mat);          p = oxstack_node_init(NULL);
           p->p = mat;
           return (p);
 }  }
   
   
   
 #if __NTL_DEBUG /*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@*/  #if __NTL_DEBUG /*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@*/
   
   #include <unistd.h>
   #include <gc/gc.h>
   #include <strstream>
   #include "gmp.h"
   
   void ntl_free(void *ptr, size_t size) {}
   
   void *ntl_realloc(void *org, size_t old, size_t size)
   {
           void *ptr = GC_realloc(org, size);
           return (ptr);
   }
   
   
 int  int
 main(int argc, char *argv[])  main(int argc, char *argv[])
 {  {
   #if 0
   {
         ZZX f;          ZZX f;
         const char *str = "[1 -1 -1 1]";          const char *str = "12345";
         int num = 1;          int num = 1;
         char *var = "x";          char *var = "x";
         cmo *c;          cmo *c;
Line 140  main(int argc, char *argv[])
Line 159  main(int argc, char *argv[])
         arg[0] = (cmo *)poly;          arg[0] = (cmo *)poly;
   
         c = ntl_fctr(arg, num);          c = ntl_fctr(arg, num);
   }
   #endif
   
   {
           cmo_zz *n;
   
           istrstream istr("[3 -3 -6]");
           ZZ cont;
           ZZX fac;
           vec_pair_ZZX_long facs;
           cmo_indeterminate *x = new_cmo_indeterminate((cmo *)new_cmo_string("x"));
           istr >> fac;
   
           factor(cont, facs, fac);
   
           //      mpz_clear(ppp->mpz);
           mp_set_memory_functions(GC_malloc, ntl_realloc, ntl_free);
           for (int i = 0;; i++) {
   
                   cmon_factors_t *mat = new_cmon_factors(cont, facs, x);
                   cmo *aaa = convert_cmon((cmo *)mat);
                   delete_cmon((cmo *)mat);
   
                   if (i % 1000 == 0) {
                           printf("GC-counts: %d,   size: %u\n", GC_gc_no, GC_get_heap_size());
                           GC_gcollect();
                           usleep(1);
                   }
           }
   }
   
   
         return (0);          return (0);
 }  }
   
 #endif  #endif /* __NTL_DEBUG @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@*/
   
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>