BinaryTranscribe.cpp

Go to the documentation of this file.

#include <cmath>
#include "BinaryTranscribe.h"
#include "FatalException.h"
#include <boost/numeric/conversion/cast.hpp>

namespace gplib
  {
    BinaryTranscribe::BinaryTranscribe(const ttranscribed &base,
        const ttranscribed &step, const tsizev &gene) :
      basevalues(base), stepsizes(step), genesizes(gene)
      {
        //check that all three vectors have the same size
        if ((basevalues.size() != stepsizes.size()) || (basevalues.size()
            != genesizes.size()))
          throw FatalException(
              "Basevalues, Stepsizes and Genesizes do not have equal length !");
        //they all have the same size, so we only check whether one is empty
        if (basevalues.empty())
          throw FatalException("Specified an empty base value vector for binary transcribe !");

      }

    BinaryTranscribe::~BinaryTranscribe()
      {
      }

    ttranscribed BinaryTranscribe::GetValues(const tpopmember &member)
      {
        unsigned int currentindex = 0;
        double currentvalue;

        ttranscribed ReturnValues(genesizes.size());
        //for each model parameter
        for (unsigned int i = 0; i < genesizes.size(); ++i)
          {
            currentvalue = 0;
            //scan the corresponding section of the genetic string
            for (unsigned int j = currentindex; j < currentindex + genesizes(i); ++j)
              {
                //calculate the numeric value from the bits
                //some compilers have trouble using pow with some combinations of numeric types
                //so we explicitly use double and int which should work everywhere
                currentvalue += member(j) * std::pow(2.0, boost::numeric_cast<int>(j - currentindex));
              }
            //currentvalue now holds the raw integer value that is encoded in the bit string
            //so we multiply by the step size and add the basevalue to get the parameter value
            currentvalue *= stepsizes(i);
            currentvalue += basevalues(i);
            currentindex += genesizes(i);
            ReturnValues(i) = currentvalue; //store in the right component
          }
        return ReturnValues;
      }

    BinaryTranscribe::BinaryTranscribe(const BinaryTranscribe &Old) :
      basevalues(Old.basevalues), stepsizes(Old.stepsizes), genesizes(
          Old.genesizes)
      {
        // The body is empty
      }

  }