computeLU() [3/5]

void DMatLUinPlace< M2::ARingGFFlint >::computeLU ( )

inlineprivate

Definition at line 143 of file dmat-lu-inplace.hpp.

{
  if (mIsDone) return;
  //  std::cout << "computing LU decomposition GFFlint" << std::endl;
  mp_limb_signed_t* perm = newarray_atomic(mp_limb_signed_t, mLU.numRows());
  fq_zech_mat_lu(perm, mLU.fq_zech_mat(), false, ring().flintContext());
  // Now we set mPerm:
  mPerm.clear();
  for (long i = 0; i < mLU.numRows(); i++) mPerm.push_back(perm[i]);
  mSign = (_perm_parity(perm, mLU.numRows()) == 0);
  freemem(perm);
 
  // Now we set mPivotColumns
  LUUtil<RingType>::computePivotColumns(mLU, mPivotColumns);
 
  mIsDone = true;
}