[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.2 and 1.3

version 1.2, 2003/11/08 12:34:00 version 1.3, 2003/11/15 09:06:20
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/ox_ntl/ntl.cpp,v 1.1 2003/11/03 03:11:21 iwane Exp $ */  /* $OpenXM: OpenXM/src/ox_ntl/ntl.cpp,v 1.2 2003/11/08 12:34:00 iwane Exp $ */
   
 #include <NTL/ZZXFactoring.h>  #include <NTL/ZZXFactoring.h>
   #include <NTL/LLL.h>
   
 #include "ox_toolkit.h"  #include "ox_toolkit.h"
 #include "ntl.h"  #include "ntl.h"
Line 9 
Line 10 
 #define __NTL_PRINT (1)  #define __NTL_PRINT (1)
 #endif  #endif
   
   #define __NTL_PRINT (1)
   #define DPRINTF(x)      printf x; fflush(stdout)
   
 /****************************************************************************  /****************************************************************************
  *   *
  * Factorize polynomial over the integers.   * Factorize polynomial over the integers.
Line 27  ntl_fctr(cmo **arg, int argc)
Line 31  ntl_fctr(cmo **arg, int argc)
         cmo *poly = arg[0];          cmo *poly = arg[0];
         cmo_indeterminate *x;          cmo_indeterminate *x;
         ZZX f;          ZZX f;
         ZZ c;  
   
         int ret;          int ret;
   
         vec_pair_ZZX_long factors;  
   
         cmo_list *ans;  
         cmo *fcts;  
   
         if (argc != 1) {          if (argc != 1) {
                 return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(#)")));                  return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(#)")));
         }          }
Line 50  ntl_fctr(cmo **arg, int argc)
Line 47  ntl_fctr(cmo **arg, int argc)
         cout << "input: " << f << endl;          cout << "input: " << f << endl;
 #endif  #endif
   
         factor(c, factors, f);          cmon_factors_t *factors = new_cmon_factors();
           factors->x = x;
   
           factor(*factors->cont, *factors->f, f);
   
 #if __NTL_PRINT  #if __NTL_PRINT
         cout << "fctr : " << factors << endl;          cout << "fctr : " << *factors->f << endl;
 #endif  #endif
   
         cmo_zz *zz = ZZ_to_cmo_zz(c);  
   
         ans = new_cmo_list();          return ((cmo *)factors);
   }
   
         fcts = (cmo *)vec_pair_ZZX_long_to_cmo(factors, x);  /****************************************************************************
    *
    * Factorize polynomial over the integers.
    *
    * PARAM : I : arg  : polynomial
    *       : I : argc :
    * RETURN: [num,[[factor1,multiplicity1],[factor2,multiplicity2],...]]
    *
    * EX    : ntl_fctr([2*x^11+4*x^8-2*x^6+2*x^5-4*x^3-2],1) ==>
            : [2,[[x-1,1],[x^4+x^3+x^2+x+1,1],[x+1,2],[x^2-x+1,2]]]
    *
    ****************************************************************************/
   cmo *
   ntl_lll(cmo **arg, int argc)
   {
           cmo *poly = arg[0];
           cmo_indeterminate *x;
           ZZX f;
           int ret;
           cmon_mat_zz_t *mat;
   
   DPRINTF(("lll start\n"));
   
         list_appendl(ans, zz, (cmo *)fcts, NULL);          if (argc != 1) {
   DPRINTF(("invalid argc %d\n", argc));
                   return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(#)")));
           }
   
         return ((cmo *)ans);          mat = new_cmon_mat_zz();
           ret = cmo_to_mat_zz(*mat->mat, arg[0]);
   DPRINTF(("ret = %d, convert\n", ret));
           if (ret != NTL_SUCCESS) {
   DPRINTF(("convert failed\n", ret));
                   delete_cmon_mat_zz(mat);
                   /* format error */
                   return ((cmo *)new_cmo_error2((cmo *)new_cmo_string("Invalid Parameter(type)")));
           }
   
   DPRINTF(("convert success\n", ret));
   
   #if __NTL_PRINT
   DPRINTF(("convert success\n", ret));
           cout << "input: " << (*mat->mat) << endl;
   #endif
   
   
   DPRINTF(("lll start\n", ret));
           ZZ det2;
           long rd = LLL(det2, *mat->mat);
   
   
   #if __NTL_PRINT
           cout << "output: " << (*mat->mat) << endl;
   #endif
   
           return ((cmo *)mat);
 }  }
   
   
 #if __NTL_DEBUG  
   #if __NTL_DEBUG /*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@*/
   
 int  int
 main(int argc, char *argv[])  main(int argc, char *argv[])
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>