rawZZfactor()

gmp_arrayZZ rawZZfactor

(

gmp_ZZ

x

)

Definition at line 46 of file flint.cpp.

{
  fmpz_t n;
  fmpz_init(n);
  fmpz_set_mpz(n, x);
  //std::cout << "factoring fmpz " << static_cast<void*>(n) << std::endl;
  fmpz_factor_t factor;
  fmpz_factor_init(factor);
  fmpz_factor(factor,n);
  int len = factor->num;
  gmp_arrayZZ result = getmemarraytype(gmp_arrayZZ,2*len+1);
  result->len = 2*len+1;
  __mpz_struct *tmp;
  // The sign is the first element of the result.
  tmp = newitem(__mpz_struct);
  mpz_init(tmp);
  mpz_set_si(tmp, factor->sign);
  mpz_reallocate_limbs(tmp);
  result->array[0] = tmp;
  for (int i=0; i<len; i++) {
    // Get the i-th factor
    tmp = newitem(__mpz_struct);
    mpz_init(tmp);
    fmpz_get_mpz(tmp,factor->p + i);
    mpz_reallocate_limbs(tmp);
    result->array[2*i+1] = tmp;
    // Get the i-th factor and its exponent
    tmp = newitem(__mpz_struct);
    mpz_init(tmp);
    fmpz_get_mpz(tmp,(fmpz *)(factor->exp + i));
    mpz_reallocate_limbs(tmp);
    result->array[2*i+2] = tmp;
  }
  fmpz_factor_clear(factor);
  fmpz_clear(n);
  return result;
}

References getmemarraytype, mpz_reallocate_limbs(), newitem, result(), and x.