MTRoughness.cpp

Go to the documentation of this file.

#include "MTRoughness.h"
#include <gsl/gsl_math.h>

namespace gplib
  {
    MTRoughness::MTRoughness()
      {
      }

    MTRoughness::~MTRoughness()
      {
      }

    void MTRoughness::SafeParallel(const ttranscribed &member)
      {
        const unsigned int length = member.size() / 3; //we have 3 parameters in the model, so size/3 layers
        double roughness = 0; // init returnvalue
        const int fitexp = GetFitExponent();

        for (unsigned int i = 1; i < length; ++i) // for all layers except the top
          {
            roughness += gsl_pow_int((member(i) - member(i - 1)), fitexp); // add the squared difference of the values divided by reference
          }
        SetRMS(roughness);
      }

    double MTRoughness::PostParallel(const ttranscribed &member)
      {
        return GetRMS();
      }

    MTRoughness::MTRoughness(const MTRoughness &Old) :
      GeneralObjective(Old)
      {

      }
    MTRoughness& MTRoughness::operator=(const MTRoughness& source)
      {
        if (this == &source)
          return *this;
        GeneralObjective::operator=(source);
        return *this;
      }
  }