rules_range.cpp

Go to the documentation of this file.

#include <openmodeller/Random.hh>
#include <math.h>
#include <string.h>

#include <openmodeller/Log.hh>
#include <openmodeller/Sample.hh>

#include "rules_range.hh"
#include "bioclim_histogram.hh"


// ==========================================================================
//  RangeRule implelentation
// ==========================================================================
RangeRule::RangeRule() : 
  GarpRule() 
{}

RangeRule::RangeRule(int numGenes) : 
  GarpRule(numGenes) 
{ }

RangeRule::RangeRule(Scalar prediction, int numGenes, 
                     const Sample& chrom1, const Sample& chrom2, 
                     const double * performances) : 
  GarpRule(prediction, numGenes, chrom1, chrom2, performances) 
{}

// ==========================================================================
RangeRule::~RangeRule() { }

// ==========================================================================
void RangeRule::initialize(const BioclimHistogram& histogram)
{
  int i, j;
  Random rnd;

  // loop iterates through variables
  for(i = 0; i < _numGenes; i++)
    {
      j = rnd.get(_numGenes);
      
      Scalar a = 0 , b = 0;
      histogram.getBioclimRange(_prediction, j, a, b);
      _chrom1[j] = a;
      _chrom2[j] = b;
    }
}

// ==========================================================================
bool RangeRule::applies(const Sample& sample) const
{
  // visit each of the genes
  for (int i = 0; i < _numGenes; i++)
    {
      if (!(equalEps(_chrom1[i], -1.0) && equalEps(_chrom2[i], +1.0)))
        if (!between(sample[i], _chrom1[i], _chrom2[i]))
          return false;
    }

  return true;
}

// ==========================================================================
int RangeRule::getStrength(const Sample& sample) const
{
  //printf("GetStrength(%+7.3f)\n", _prediction);
  for (int i = 0; i < _numGenes; i++)
    {
      if (!membership(_chrom1[i], _chrom2[i], sample[i])) 
  { 
    //printf("Strength = 0\n"); 
    return 0; 
  }
    } 
  
  //printf("Strength = 1\n");
  return 1;
}

// ==========================================================================
void RangeRule::log()
{
  Log::instance()->info( "Range: " );
  GarpRule::log();
}

// ==========================================================================