fp_LLL()

bool fp_LLL	(	MutableMatrix *	M,
		MutableMatrix *	U,
		int	strategy )

Definition at line 9 of file fplll-interface.cpp.

{
  (void) M;
  (void) U;
  (void) strategy;
#ifndef HAVE_FPLLL
  ERROR("fplll is not available (configure M2 with fplll!)");
  return 0;
#else
  assert(U == NULL);
  double delta = 0.99;
  double eta = 0.51;
  fplll::LLLMethod method = fplll::LM_WRAPPER;
  fplll::FloatType floatType = fplll::FT_DEFAULT;
  int precision = 0;
  int flags = fplll::LLL_DEFAULT;
 
  int ncols = static_cast<int>(M->n_rows());
  int nrows = static_cast<int>(M->n_cols());
 
  fplll::ZZ_mat<mpz_t> mat(nrows, ncols);
 
  for (int i = 0; i < nrows; i++)
    for (int j = 0; j < ncols; j++)
      {
        ring_elem a;
        if (M->get_entry(j, i, a))
          {
            mpz_set(mat[i][j].get_data(), a.get_mpz());
          }
      }
 
  int result =
      lll_reduction(mat, delta, eta, method, floatType, precision, flags);
 
  switch (result)
    {
      case fplll::RED_SUCCESS:
        break;
      case fplll::RED_BABAI_FAILURE:
        ERROR("Error in fpLLL");
        return 0;
      case fplll::RED_LLL_FAILURE:
        ERROR("infinite loop in LLL");
        return 0;
      default:
        ERROR("unknown error in fpLLL");
        return 0;
    }
 
  /* Put this back into M */
  mpz_t a;
  mpz_init(a);
 
  for (int j = 0; j < ncols; j++)
    for (int i = 0; i < nrows; i++)
      {
        mpz_set(a, mat[i][j].get_data());
        ring_elem b = globalZZ->from_int(a);
        M->set_entry(j, i, b);
      }
  mpz_clear(a);
  return true;
#endif
}

References ERROR, MutableMatrix::get_entry(), ring_elem::get_mpz(), globalZZ, MutableMatrix::n_cols(), MutableMatrix::n_rows(), result(), and MutableMatrix::set_entry().

Referenced by rawLLL().