version 1.5, 2004/07/04 02:31:51 |
version 1.7, 2013/10/18 01:28:52 |
|
|
/* $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.6 2008/09/19 10:55:40 iwane Exp $ */ |
|
|
#include <NTL/ZZXFactoring.h> |
#include <NTL/ZZXFactoring.h> |
#include <NTL/LLL.h> |
#include <NTL/LLL.h> |
|
|
: [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; |
|
char emss[25] = "Invalid Parameter(#)"; |
|
char emst[25] = "Invalid Parameter(type)"; |
|
|
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(emss))); |
|
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(emst)); |
|
p = oxstack_node_init(err); |
|
return (p); |
} |
} |
|
|
#if __NTL_PRINT |
#if __NTL_PRINT |
cout << "input: " << f << endl; |
std::cout << "input: " << f << std::endl; |
#endif |
#endif |
|
|
cmon_factors_t *factors = new_cmon_factors(); |
cmon_factors_t *factors = new_cmon_factors(); |
Line 52 ntl_fctr(cmo **arg, int argc) |
|
Line 60 ntl_fctr(cmo **arg, int argc) |
|
factor(*factors->cont, *factors->f, f); |
factor(*factors->cont, *factors->f, f); |
|
|
#if __NTL_PRINT |
#if __NTL_PRINT |
cout << "fctr : " << *factors->f << endl; |
std::cout << "fctr : " << *factors->f << std::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 81 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; |
|
char emss[25] = "Invalid Parameter(#)"; |
|
char emst[25] = "Invalid Parameter(type)"; |
|
|
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(emss))); |
|
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(emst))); |
|
p = oxstack_node_init(err); |
|
return (p); |
} |
} |
|
|
#if __NTL_PRINT |
#if __NTL_PRINT |
cout << "input: " << (*mat->mat) << endl; |
std::cout << "input: " << (*mat->mat) << std::endl; |
#endif |
#endif |
|
|
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; |
std::cout << "output: " << (*mat->mat) << std::endl; |
cout << U << endl; |
std::cout << U <<std:: endl; |
#endif |
#endif |
|
|
return ((cmo *)mat); |
p = oxstack_node_init(NULL); |
|
p->p = mat; |
|
return (p); |
} |
} |
|
|
|
|
Line 112 ntl_lll(cmo **arg, int argc) |
|
Line 131 ntl_lll(cmo **arg, int argc) |
|
|
|
#include <unistd.h> |
#include <unistd.h> |
#include <gc/gc.h> |
#include <gc/gc.h> |
#include <strstream> |
|
#include "gmp.h" |
|
|
|
void ntl_free(void *ptr, size_t size) {} |
void ntl_free(void *ptr, size_t size) {} |
|
|
Line 148 main(int argc, char *argv[]) |
|
Line 165 main(int argc, char *argv[]) |
|
#endif |
#endif |
|
|
{ |
{ |
cmo_zz *n; |
|
|
|
istrstream istr("[3 -3 -6]"); |
|
ZZ cont; |
ZZ cont; |
ZZX fac; |
ZZX fac; |
vec_pair_ZZX_long facs; |
vec_pair_ZZX_long facs; |
cmo_indeterminate *x = new_cmo_indeterminate((cmo *)new_cmo_string("x")); |
// char xx[2] = "x"; |
istr >> fac; |
// cmo_indeterminate *x = new_cmo_indeterminate((cmo *)new_cmo_string(xx)); |
|
|
|
// std::istringstream istr("[3 -3 -6]"); |
|
// istr >> fac; |
|
|
|
SetCoeff(fac, 0, 3); |
|
SetCoeff(fac, 1, -3); |
|
SetCoeff(fac, 2, -6); |
|
|
|
|
factor(cont, facs, fac); |
factor(cont, facs, fac); |
|
|
// mpz_clear(ppp->mpz); |
// mpz_clear(ppp->mpz); |
|
} |
|
|
|
if (0) { |
mp_set_memory_functions(GC_malloc, ntl_realloc, ntl_free); |
mp_set_memory_functions(GC_malloc, ntl_realloc, ntl_free); |
for (int i = 0;; i++) { |
for (int i = 0;; i++) { |
|
|
cmon_factors_t *mat = new_cmon_factors(cont, facs, x); |
// cmon_factors_t *mat = new_cmon_factors(cont, facs, x); |
cmo *aaa = convert_cmon((cmo *)mat); |
// cmo *aaa = convert_cmon((cmo *)mat); |
delete_cmon((cmo *)mat); |
// delete_cmon((cmo *)mat); |
|
|
if (i % 1000 == 0) { |
if (i % 1000 == 0) { |
printf("GC-counts: %d, size: %u\n", GC_gc_no, GC_get_heap_size()); |
printf("GC-counts: %lu, size: %u\n", GC_gc_no, GC_get_heap_size()); |
GC_gcollect(); |
GC_gcollect(); |
usleep(1); |
usleep(1); |
} |
} |