[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.5 and 1.6

version 1.5, 2004/07/04 02:31:51 version 1.6, 2008/09/19 10:55:40
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/ox_ntl/ntl.cpp,v 1.4 2003/11/17 12:04:20 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 24 
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 55  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 71  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)
 {  {
           oxstack_node *p = arg[0];
           cmo *err;
         ZZX f;          ZZX f;
         int ret;          int ret;
         cmon_mat_zz_t *mat;          cmon_mat_zz_t *mat;
   
         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);
         }          }
   
         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);
         if (ret != NTL_SUCCESS) {          if (ret != NTL_SUCCESS) {
                 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);
         }          }
   
 #if __NTL_PRINT  #if __NTL_PRINT
Line 96  ntl_lll(cmo **arg, int argc)
Line 110  ntl_lll(cmo **arg, int argc)
   
         ZZ det2;          ZZ det2;
         mat_ZZ U;          mat_ZZ U;
         long rd = LLL(det2, *mat->mat, U);  
   
 #if __NTL_PRINT  #if __NTL_PRINT
         cout << "output: " << (*mat->mat) << endl;          cout << "output: " << (*mat->mat) << endl;
         cout << U << endl;          cout << U << endl;
 #endif  #endif
   
         return ((cmo *)mat);          p = oxstack_node_init(NULL);
           p->p = mat;
           return (p);
 }  }
   
   
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>