Ejemplos de HNF en C++ (Cpp)

Ejemplo n.º 1

Archivo: MatrixTest.c Proyecto: JamesHirschorn/QFCL

NTL_CLIENT

int main()
{
   mat_ZZ B, X;
   vec_ZZ v, w;

   cin >> B;
   cin >> v;

   ZZ d;

   double t;
   cerr << "matrix inverse...";
   t = GetTime();
   inv(d, X, B);
   cerr << (GetTime()-t) << "\n";

   cout << d << "\n";
   cout << X << "\n";

   cout << "\n\n\n";

   cerr << "hensel solve...";
   t = GetTime();
   HenselSolve1(d, w, B, v);
   cerr << (GetTime()-t) << "\n";

   cout << d << "\n";
   cout << w << "\n";

   cout << "\n\n\n";

   ZZX f;

   cerr << "char poly...";
   t = GetTime();
   CharPoly(f, B);
   cerr << (GetTime()-t) << "\n";

   cout << f << "\n";

   cout << "\n\n\n";

   cerr << "HNF...";
   t = GetTime();
   HNF(X, B, d);
   cerr << (GetTime()-t) << "\n";

   cout << X;

   return 0;
}

Ejemplo n.º 2

Archivo: ntlwrap.cpp Proyecto: saraedum/sage-renamed

static CYTHON_INLINE mat_ZZ* mat_ZZ_HNF(const mat_ZZ* A, const struct ZZ* D)
{
    mat_ZZ* W = new mat_ZZ();
    HNF(*W, *A, *D);
    return W;
}