KBasis()

KBasis::KBasis ( const Matrix * bottom, const int * lo_degree, const int * hi_degree, std::vector< int > heftvec, std::vector< int > varlist, bool do_truncation, int limit )

private

Definition at line 168 of file matrix-kbasis.cpp.

    : bottom_matrix(bottom),
      mHeftVector(heftvec),
      mVariables(varlist),
      do_truncation(do_truncation0),
      weight_has_zeros(false),
      limit(limit0),
      lo_degree(lo_degree0),
      hi_degree(hi_degree0),
      kb_error(false)
{
  P = bottom->get_ring()->cast_to_PolynomialRing();
  M = P->getMonoid();
  D = P->get_degree_ring()->getMonoid();
 
  if (lo_degree == nullptr && hi_degree == nullptr) { computation_type = KB_FULL; }
  else if (mHeftVector.size() == 1) { computation_type = KB_SINGLE; }
  else { computation_type = KB_MULTI; }
 
  // Compute the (non-negative) weights of each of the variables in
  // 'mVariables'.
 
  var_wts = newarray_atomic(int, mVariables.size());
  var_degs = newarray_atomic(int, mVariables.size() * mHeftVector.size());
  int *exp =
      newarray_atomic(int, D->n_vars());  // used to hold exponent vectors
  int next = 0;
  for (int i = 0; i < mVariables.size(); i++, next += mHeftVector.size())
    {
      int v = mVariables[i];
      D->to_expvector(M->degree_of_var(v), exp);
      var_wts[i] = exponents::weight(mHeftVector.size(), exp, mHeftVector);
      if (var_wts[i] == 0) weight_has_zeros = true;
      exponents::copy(mHeftVector.size(), exp, var_degs + next);
    }
  freemem(exp);
 
  // Set the recursion variables
  kb_exp = newarray_atomic_clear(int, P->n_vars());
  kb_exp_weight = 0;
 
  if (lo_degree != nullptr)
    kb_target_lo_weight =
        exponents::weight(mHeftVector.size(), lo_degree, mHeftVector);
  if (hi_degree != nullptr)
    kb_target_hi_weight =
        exponents::weight(mHeftVector.size(), hi_degree, mHeftVector);
 
  if (lo_degree && hi_degree && mHeftVector.size() == 1 && mHeftVector[0] < 0)
    {
      int t = kb_target_lo_weight;
      kb_target_lo_weight = kb_target_hi_weight;
      kb_target_hi_weight = t;
    }
 
  kb_mon = M->make_one();
 
  mat = MatrixConstructor(bottom->rows(), 0);
  kb_exp_multidegree = D->make_one();
 
  if (mHeftVector.size() > 1 && lo_degree != nullptr)
    {
      kb_target_multidegree = D->make_one();
      exponents::copy(mHeftVector.size(), lo_degree, kb_target_multidegree);
    }
  else { kb_target_multidegree = nullptr; }
}

References bottom_matrix, computation_type, ExponentVector< int, true >::copy(), D, do_truncation, freemem(), Matrix::get_ring(), hi_degree, kb_error, kb_exp, kb_exp_multidegree, kb_exp_weight, KB_FULL, kb_mon, kb_target_hi_weight, kb_target_lo_weight, kb_target_multidegree, limit, lo_degree, M, mat, Matrix, mHeftVector, mVariables, newarray_atomic, newarray_atomic_clear, P, Matrix::rows(), var_degs, var_wts, ExponentVector< int, true >::weight(), and weight_has_zeros.

Referenced by k_basis().