to_M2_vec()

vec ResF4toM2Interface::to_M2_vec ( const ResPolyRing & R, const ResPolynomial & f, const FreeModule * F )

static

Definition at line 94 of file res-f4-m2-interface.cpp.

{
  const PolynomialRing* origR = F->get_ring()->cast_to_PolynomialRing();
  const Monoid* M = origR->getMonoid();
 
  monomial m1 = M->make_one();
 
  Nterm** comps = newarray(Nterm*, F->rank());
  Nterm** last = newarray(Nterm*, F->rank());
  for (int i = 0; i < F->rank(); i++)
    {
      comps[i] = nullptr;
      last[i] = nullptr;
    }
 
  int* exp = new int[M->n_vars()];
 
  const res_monomial_word* w = f.monoms.data();
  for (int i = 0; i < f.len; i++)
    {
      component_index comp;
      R.monoid().to_expvector(w, exp, comp);
      w = w + R.monoid().monomial_size(w);
      M->from_expvector(exp, m1);
      ring_elem a =
          R.vectorArithmetic().ringElemFromElementArray(f.coeffs,i);
      Nterm* g = origR->make_flat_term(a, m1);
      g->next = nullptr;
      if (last[comp] == nullptr)
        {
          comps[comp] = g;
          last[comp] = g;
        }
      else
        {
          last[comp]->next = g;
          last[comp] = g;
        }
    }
  vec result = nullptr;
  for (int i = 0; i < F->rank(); i++)
    {
      if (comps[i] != nullptr)
        {
          vec v = origR->make_vec(i, comps[i]);
          origR->add_vec_to(result, v);
          comps[i] = nullptr;
          last[i] = nullptr;
        }
    }
 
  delete[] exp;
  return result;
}

References Ring::add_vec_to(), Ring::cast_to_PolynomialRing(), ResPolynomial::coeffs, Monoid::from_expvector(), FreeModule::get_ring(), PolynomialRing::getMonoid(), ResPolynomial::len, PolynomialRing::make_flat_term(), Monoid::make_one(), Ring::make_vec(), ResPolyRing::monoid(), monomial, ResMonoidDense::monomial_size(), ResPolynomial::monoms, Monoid::n_vars(), newarray, Nterm::next, FreeModule::rank(), result(), VectorArithmetic::ringElemFromElementArray(), ResMonoidDense::to_expvector(), and ResPolyRing::vectorArithmetic().

Referenced by to_M2_matrix().