compute()

void KBasis::compute ( )

private

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

{
  if (limit == 0) return;
  std::vector<int> zero_vars;
  if (weight_has_zeros)
    {
      for (int i = 0; i < mVariables.size(); i++)
        {
          if (var_wts[i] == 0) { zero_vars.push_back(i); }
        }
    }
 
  for (int i = 0; i < bottom_matrix->n_rows(); i++)
    {
      if (system_interrupted()) return;
      kb_comp = i;
 
      // Make the monomial ideal: this should contain only
      // monomials involving 'mVariables'.
      kb_monideal = bottom_matrix->make_monideal(i, true);
      // the true means: over ZZ, don't consider monomials with non-unit lead
      // coeffs
 
      if (kb_monideal->is_one()) continue;
      if (hi_degree == nullptr)
        {
          // check here that kb_monideal is 0-dimensional
          // (at least for the variables being used):
          if (!all_have_pure_powers(kb_monideal, mVariables))
            {
              kb_error = true;
              ERROR("module given is not finite over the base");
              return;
            }
        }
      else if (zero_vars.size() > 0)
        {
          // if we have any variables with zero degrees, then kb_monideal needs
          // to be 0-dimensional in those variables.
          if (!all_have_pure_powers(kb_monideal, zero_vars))
            {
              kb_error = true;
              ERROR(
                  "module given is not finite over the zero-degree variables");
              return;
            }
        }
 
      const int *component_degree = bottom_matrix->rows()->degree(i);
      D->to_expvector(component_degree, kb_exp_multidegree);
      kb_exp_weight =
          exponents::weight(mHeftVector.size(), kb_exp_multidegree, mHeftVector);
 
      // Do the recursion
      switch (computation_type)
        {
          case KB_FULL:
            basis0_full();
            break;
          case KB_SINGLE:
            basis0_singly_graded();
            break;
          case KB_MULTI:
            basis0_multi_graded();
            break;
        }
    }
}

References all_have_pure_powers(), basis0_full(), basis0_multi_graded(), basis0_singly_graded(), bottom_matrix, computation_type, D, ERROR, hi_degree, kb_comp, kb_error, kb_exp_multidegree, kb_exp_weight, KB_FULL, kb_monideal, KB_MULTI, KB_SINGLE, limit, mHeftVector, mVariables, system_interrupted(), var_wts, ExponentVector< int, true >::weight(), and weight_has_zeros.

Referenced by k_basis().