algutils.cc

Gambit 是一个游戏库理论软件

//
// $Source: /home/gambit/CVS/gambit/sources/nash/algutils.cc,v $
// $Date: 2002/08/27 18:29:36 $
// $Revision: 1.3 $
//
// DESCRIPTION:
// Implementation of utility classes for algorithms
//
// This file is part of Gambit
// Copyright (c) 2002, The Gambit Project
//
// 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, write to the Free Software
// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
//

#include "algutils.h"
#include <assert.h>

AlgParams::AlgParams(void)
  : m_accuracy(0), trace(0), stopAfter(0), precision(precDOUBLE) 
{ }

AlgParams::~AlgParams()
{ }

void AlgParams::SetAccuracy(const gNumber &p_accuracy)
{ m_accuracy = p_accuracy; }


FuncMinParams::FuncMinParams(void)
  : maxits1(100), maxitsN(20), tol1(2.0e-10), tolN(1.0e-10) 
{ }

FuncMinParams::~FuncMinParams()
{ }

void FuncMinParams::SetAccuracy(const gNumber &p_accuracy)
{ 
  m_accuracy = p_accuracy;
  maxits1 = maxitsN = 200;
  tol1 = tolN = p_accuracy * p_accuracy;
}

//#include "algutils.imp"

//------------------------------------------------------------------------
//                      class gFact<T>: Member functions
//------------------------------------------------------------------------

template <class T> gFact<T>::gFact(void) 
  : m_answer(0)
{ }

template <class T> gFact<T>::gFact(const gFact<T> &p_fact) 
  : m_answer(0)
{ 
  if (p_fact.Checked()) {
    Set(p_fact.Answer());
  }
}

template <class T> gFact<T>::~gFact() 
{ 
  if (m_answer)  {
    delete m_answer;
  }
}


template <class T> gFact<T> &gFact<T>::operator=(const gFact<T> &p_fact)
{ 
  if (this != &p_fact) {
    if (m_answer) {
      delete m_answer;
      m_answer = 0;
    }

    if (p_fact.Checked()) { 
      Set(p_fact.Answer());
    }
  } 
  return *this;
}

template <class T> const T &gFact<T>::Answer(void) const
{
  if (!m_answer) { 
    throw Unknown();
  }
  return *m_answer;
}

template <class T> void gFact<T>::Set(const T &p_state)
{
  if (m_answer && *m_answer != p_state) { 
    throw Contradiction();
  }
  if (!m_answer) { 
    m_answer = new T(p_state);
  }
}

template <class T> gFact<T>::Contradiction::~Contradiction()
{ }

template <class T> gText gFact<T>::Contradiction::Description(void) const 
{
  return "Contradictory information in gFact";
}

template <class T> gFact<T>::Unknown::~Unknown()
{ }

template <class T> gText gFact<T>::Unknown::Description(void) const 
{
  return "Unknown information in gFact";
}

template class gFact<gNumber>;
template class gFact<gTriState>;

#include "math/gpvector.h"
template class gFact<gPVector<gNumber> >;

// Following is an intitial, more complicated implementation of the above
// I am checking it in (commented out), so I don't forget how to 
// pass pointers to member functions of a class, since this may be useful 
// elsewhere.  

/*
template <class C, class T> gFact<C, T>::gFact(const C &m,FactFinder f) 
  : ms(m), answer(0), gofer(f)
{ }

template <class C, class T> gFact<C, T>::gFact(const gFact<C, T> &f) 
  : ms(f.ms), answer(f.answer), gofer(f.gofer)
{
  // Don't use this constructor.  
  assert(0);  
}

template <class C, class T> gFact<C, T>::~gFact(void) 
{ 
  if(answer) delete answer;
}


template <class C, class T> gFact<C, T> &gFact<C, T>::operator=(const gFact<C, T> &b)
{ 
  if (this != &b) {
    if(b.Checked()) 
      Set(b.Answer());
  } 
  return *this;
}

template <class C, class T> const T &gFact<C, T>::Answer(void) const
{
  if(!answer) 
    answer = new T((ms.*gofer)());
  return *answer;
}

template <class C, class T> void gFact<C, T>::Set(T state)
{
  if(!answer) 
    answer = new T(state);
  if(*answer!=state ) 
    throw Contradiction();
}

template <class C, class T> gFact<C, T>::Contradiction::~Contradiction()
{ }

template <class C, class T> gText gFact<C, T>::Contradiction::Description(void) const 
{
  return "Contradictory information in gFact";
}

class MixedSolution;
class BehavSolution;
template class gFact<MixedSolution, gNumber>;
template class gFact<MixedSolution, gTriState>;

template class gFact<BehavSolution, gNumber>;
template class gFact<BehavSolution, gTriState>;
#include "gpvector.h"
template class gFact<BehavSolution, gPVector<gNumber> >;

*/