best_subsets.hh

Go to the documentation of this file.

#ifndef _BEST_SUBSETS_HH_
#define _BEST_SUBSETS_HH_

#include <openmodeller/om.hh>

#include "bs_algorithm_factory.hh"
#include "garp_run.hh"

#include <string>
using std::string;

/****************************************************************/
/************************* GARP Algorithm ***********************/

class BestSubsets : public AlgorithmImpl, public BSAlgorithmFactory
{
public:

  BestSubsets(AlgMetadata * metadata);
  virtual ~BestSubsets();

  int initialize();
  
  int iterate();
  
  int done() const;

  float getProgress() const;

  //
  // Methods used to project the model
  //
  
  Scalar getValue( const Sample& x ) const;
  
  int getConvergence( Scalar * const val ) const;

  /*
   */
  void printListOfRuns(string msg, AlgorithmRun ** runs, int numOfRuns);

protected:

  void _getConfiguration( ConfigurationPtr& ) const;
  void _setConfiguration( const ConstConfigurationPtr & );

  int _nparam;
  AlgParameter * _alg_params;

private:

  virtual AlgorithmImpl * getBSAlgorithm() = 0;
  virtual int transferParametersToAlgorithm() = 0;

  int numActiveThreads();
  int earlyTerminationConditionMet();
  int calculateBestSubset();
  void sortRuns(AlgorithmRun ** runList, int nelements, int errorType);


  //
  // best subsets parameters
  //

  double _trainProp;
  int _totalRuns;
  double _omissionThreshold;
  int _modelsUnderOmission;
  double _commissionThreshold;
  int _commissionSampleSize;
  int _maxThreads;

  bool _softOmissionThreshold;
  int _currentModelsUnderOmissionThreshold;

  //
  // Internal data structures for Best Subsets
  //
  AlgorithmRun ** _finishedRun;
  AlgorithmRun ** _activeRun;
  AlgorithmRun ** _bestRun;

  int _numFinishedRuns;
  int _numActiveRuns;
  int _numBestRuns;

  int _done;

  mutable float _maxProgress;
};


#endif