do_spairs()

void F4GB::do_spairs ( )

private

Definition at line 1008 of file f4.cpp.

{
  if (hilbert && hilbert->nRemainingExpected() == 0)
    {
      mSPairSet.discardSPairsInCurrentDegree();
      if (M2_gbTrace >= 1)
        fprintf(stderr,
                "-- skipping degree...no elements expected in this degree\n");
      return;
    }
  reset_matrix();
  clock_t begin_time = clock();
 
  n_lcmdups = 0;
  make_matrix();
 
  if (M2_gbTrace >= 6)
    {
      fprintf(stderr, "---------\n");
      show_matrix();
      fprintf(stderr, "---------\n");
    }
 
  clock_t end_time = clock();
  clock_make_matrix += end_time - begin_time;
  double nsecs = static_cast<double>(end_time - begin_time);
  nsecs /= CLOCKS_PER_SEC;
  if (M2_gbTrace >= 2) fprintf(stderr, " make matrix time = %f\n", nsecs);
 
  if (M2_gbTrace >= 3) mMonomialHashTable.dump();
 
  double oldParallelGauss = mParallelGaussTime;
  double oldSerialGauss = mSerialGaussTime;
  double oldTailReduce = mTailReduceTime;
  double oldNewSPair = mNewSPairTime;
  double oldInsertNewGB = mInsertGBTime;
 
  begin_time = clock();
  gauss_reduce(true);
  end_time = clock();
  clock_gauss += end_time - begin_time;
 
  //  fprintf(stderr, "---------\n");
  //  show_matrix();
  //  fprintf(stderr, "---------\n");
 
  nsecs = static_cast<double>(end_time - begin_time);
  nsecs /= CLOCKS_PER_SEC;
  if (M2_gbTrace >= 2)
    {
      fprintf(stderr, " gauss time          = %f\n", nsecs);
      fprintf(stderr, " parallel gauss time          = %g\n", mParallelGaussTime - oldParallelGauss);
      fprintf(stderr, " serial gauss time            = %g\n", mSerialGaussTime - oldSerialGauss);
      fprintf(stderr, " tail reduce time             = %g\n", mTailReduceTime - oldTailReduce);
 
      fprintf(stderr, " lcm dups            = %ld\n", n_lcmdups);
      if (M2_gbTrace >= 6)
        {
          fprintf(stderr, "---------\n");
          show_matrix();
          fprintf(stderr, "---------\n");
        }
    }
  new_GB_elements();
  if (M2_gbTrace >= 2)
    {
      fprintf(stderr, " finding new spair time             = %g\n",  mNewSPairTime - oldNewSPair);
      fprintf(stderr, " number of spairs in queue          = %zu\n", mSPairSet.numberOfSPairs());
      fprintf(stderr, " insert new gb time                 = %g\n",  mInsertGBTime - oldInsertNewGB - (mNewSPairTime - oldNewSPair));
    }
  
  int ngb = INTSIZE(mGroebnerBasis);
  if (M2_gbTrace >= 1)
    {
      fprintf(stderr, " # GB elements   = %d\n", ngb);
      if (M2_gbTrace >= 5) {
        show_gb_array(mGroebnerBasis);
        mSPairSet.display();
      }
    }
 
  clear_matrix();
}

References clear_matrix(), clock_gauss, clock_make_matrix, gauss_reduce(), hilbert, INTSIZE, M2_gbTrace, make_matrix(), mGroebnerBasis, mInsertGBTime, mMonomialHashTable, mNewSPairTime, mParallelGaussTime, mSerialGaussTime, mSPairSet, mTailReduceTime, n_lcmdups, new_GB_elements(), reset_matrix(), show_gb_array(), and show_matrix().

Referenced by start_computation().