HilbertStrategy.cpp

Go to the documentation of this file.

/* Frobby: Software for monomial ideal computations.
   Copyright (C) 2007 Bjarke Hammersholt Roune (www.broune.com)
 
   This program is free software; you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation; either version 2 of the License, or
   (at your option) any later version.
 
   This program is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with this program.  If not, see http://www.gnu.org/licenses/.
*/
#include "stdinc.h"
#include "HilbertStrategy.h"
 
#include "Term.h"
#include "HilbertSlice.h"
#include "Ideal.h"
#include "CoefTermConsumer.h"
#include "Projection.h"
#include "IndependenceSplitter.h"
#include "HilbertIndependenceConsumer.h"
#include "ElementDeleter.h"
 
HilbertStrategy::HilbertStrategy(CoefTermConsumer* consumer,
                                 const SplitStrategy* splitStrategy):
  SliceStrategyCommon(splitStrategy),
  _consumerCache(),
  _consumerCacheDeleter(_consumerCache),
  _consumer(consumer) {
}
 
void HilbertStrategy::run(const Ideal& ideal) {
  ASSERT(_consumer != 0);
 
  size_t varCount = ideal.getVarCount();
  Ideal sliceIdeal(varCount);
 
  if (!ideal.contains(Term(varCount))) {
    _consumer->consume(1, Term(varCount));
 
    if (ideal.getGeneratorCount() > 0) {
      Term allOnes(varCount);
      for (size_t var = 0; var < varCount; ++var)
        allOnes[var] = 1;
 
      sliceIdeal = ideal;
      sliceIdeal.product(allOnes);
    }
  }
 
  auto_ptr<Slice> slice
    (new HilbertSlice(*this, sliceIdeal, Ideal(varCount),
                      Term(varCount), _consumer));
 
  simplify(*slice);
  _tasks.addTask(slice.release());
  _tasks.runTasks();
  _consumerCacheDeleter.deleteElements();
}
 
bool HilbertStrategy::processSlice
(TaskEngine& tasks, auto_ptr<Slice> slice) {
  ASSERT(slice.get() != 0);
  ASSERT(debugIsValidSlice(slice.get()));
 
  if (slice->baseCase(getUseSimplification())) {
    freeSlice(slice);
    return true;
  }
 
  if (getUseIndependence() && _indepSplitter.analyze(*slice)) {
    independenceSplit(slice);
  } else {
    ASSERT(_split->isPivotSplit());
    pivotSplit(auto_ptr<Slice>(slice));
  }
 
  return false;
}
 
auto_ptr<HilbertSlice> HilbertStrategy::newHilbertSlice() {
  auto_ptr<Slice> slice(newSlice());
  ASSERT(debugIsValidSlice(slice.get()));
  return auto_ptr<HilbertSlice>(static_cast<HilbertSlice*>(slice.release()));
}
 
auto_ptr<Slice> HilbertStrategy::allocateSlice() {
  return auto_ptr<Slice>(new HilbertSlice(*this));
}
 
bool HilbertStrategy::debugIsValidSlice(Slice* slice) {
  ASSERT(slice != 0);
  ASSERT(dynamic_cast<HilbertSlice*>(slice) != 0);
  return true;
}
 
void HilbertStrategy::getPivot(Term& term, Slice& slice) {
  ASSERT(term.getVarCount() == slice.getVarCount());
 
  _split->getPivot(term, slice);
}
 
void HilbertStrategy::freeConsumer(auto_ptr<HilbertIndependenceConsumer>
                                   consumer) {
  ASSERT(consumer.get() != 0);
  ASSERT(std::find(_consumerCache.begin(),
                   _consumerCache.end(), consumer.get()) ==
         _consumerCache.end());
 
  consumer->clear();
  noThrowPushBack(_consumerCache, consumer);
}
 
void HilbertStrategy::independenceSplit(auto_ptr<Slice> sliceParam) {
  ASSERT(sliceParam.get() != 0);
  ASSERT(debugIsValidSlice(sliceParam.get()));
  auto_ptr<HilbertSlice> slice
    (static_cast<HilbertSlice*>(sliceParam.release()));
 
  // Construct split object.
  auto_ptr<HilbertIndependenceConsumer> autoSplit = newConsumer();
  autoSplit->reset(slice->getConsumer(), _indepSplitter, slice->getVarCount());
  HilbertIndependenceConsumer* split = autoSplit.release();
  _tasks.addTask(split); // Runs when we are done with all of this split.
 
  // Construct left slice.
  auto_ptr<HilbertSlice> leftSlice(newHilbertSlice());
  leftSlice->setToProjOf(*slice, split->getLeftProjection(),
                         split->getLeftConsumer());
  _tasks.addTask(leftSlice.release());
 
  // Construct right slice.
  auto_ptr<HilbertSlice> rightSlice(newHilbertSlice());
  rightSlice->setToProjOf(*slice, split->getRightProjection(),
                          split->getRightConsumer());
  _tasks.addTask(rightSlice.release());
 
  // Deal with slice.
  freeSlice(auto_ptr<Slice>(slice));
}
 
auto_ptr<HilbertIndependenceConsumer> HilbertStrategy::newConsumer() {
  if (_consumerCache.empty())
    return auto_ptr<HilbertIndependenceConsumer>
      (new HilbertIndependenceConsumer(this));
 
  auto_ptr<HilbertIndependenceConsumer> consumer(_consumerCache.back());
  _consumerCache.pop_back();
  return consumer;
}