CoCoA » CoCoALib

Try the following program; and then again but making P a NewPolyRing_DMPI.
The times I get are: clean=20s, DMPI=25s. There must be something wrong.

  const long N = 10000000;
  const long nvars = 12;
  cout << ShortDescription << endl;
  ring Q = RingQQ();
  SparsePolyRing P = NewPolyRing(Q, nvars);
  clog << "P = " << P << endl;
  double StartTime = CpuTime();
  vector<RingElem> v; v.reserve(N);

  RandomSeqBigInt rnd(-power(2,130),power(2,130));
  RandomSeqLong RndExp(0,100);
  vector<long> expv(nvars);
  for (long i=0; i < N; ++i)
  {
    RingElem coeff = RingElem(Q, BigRat(NextValue(rnd),NextValue(rnd)));
    for (long j=0; j < nvars; ++j)
      expv[j] = NextValue(RndExp);
    v.push_back(monomial(P, coeff, expv));
  }
  clog << "vector filled after time " << CpuTime()-StartTime << "\n";