matrix_lift()

M2_bool GaussElimComputation::matrix_lift ( const Matrix * m, const Matrix ** result_remainder, const Matrix ** result_quotient )

virtual

Implements GBComputation.

Definition at line 368 of file gauss.cpp.

{
  if (m->get_ring() != R)
    {
      ERROR("encountered different rings");
      *result_remainder = nullptr;
      *result_quotient = nullptr;
      return false;
    }
  if (m->n_rows() != gens->rows()->rank())
    {
      ERROR("expected matrices to have same number of rows");
      *result_remainder = nullptr;
      *result_quotient = nullptr;
      return false;
    }
  MatrixConstructor mat_remainder(m->rows(), m->cols(), m->degree_shift());
  MatrixConstructor mat_quotient(Fsyz, m->cols(), nullptr);
  bool all_zeroes = true;
  for (int i = 0; i < m->n_cols(); i++)
    {
      vec f = R->copy_vec(m->elem(i));
      vec fsyz = nullptr;
 
      reduce(f, fsyz);
      R->negate_vec_to(fsyz);
      if (f != nullptr) all_zeroes = false;
      mat_remainder.set_column(i, f);
      mat_quotient.set_column(i, fsyz);
    }
  *result_remainder = mat_remainder.to_matrix();
  *result_quotient = mat_quotient.to_matrix();
  return all_zeroes;
}

References Matrix::cols(), Matrix::degree_shift(), Matrix::elem(), ERROR, Fsyz, gens, Matrix::get_ring(), Matrix, Matrix::n_cols(), Matrix::n_rows(), R, reduce(), Matrix::rows(), MatrixConstructor::set_column(), and MatrixConstructor::to_matrix().