SWAnisoRoughness.cpp

Go to the documentation of this file.

#include "SWAnisoRoughness.h"

namespace gplib
  {
    SWAnisoRoughness::SWAnisoRoughness() :
      veldiffweight(0.0), anisovelweight(1.0), strikediffweight(0.0),
          deltastrikediffweight(0.0)
      {
      }

    SWAnisoRoughness::~SWAnisoRoughness()
      {
      }

    void SWAnisoRoughness::SafeParallel(const ttranscribed &member)
      {
        const unsigned int length = member.size() / 6; //we have 6 parameters in the model, so size/6 layers
        double roughness = std::abs(member(2 * length)); // the loop starts at 1, so we skip this value in the 2nd summation
        const int fitexp = GetFitExponent();
        const double anglescale = 90.; //value for phi from 0 to 180
        for (unsigned int i = 1; i < length; ++i) // for all layers except the top
          {
            roughness += veldiffweight * std::pow((member(i + 1 * length)
                - member(i - 1 + 1 * length)), fitexp); //minimize difference in S velocity between layers
            roughness += anisovelweight * std::pow(std::abs(member(i + 5
                * length)), fitexp); //minimize anisovel factor in each layer
            roughness += strikediffweight * std::pow((member(i + 2 * length))
                / anglescale, fitexp); //minimize strike within each layer
            roughness += deltastrikediffweight * std::pow(std::abs(member(i + 3
                * length) / 30), fitexp); //minimize deltaangle in each layer (for a range :-60/ 60)
          }
        SetRMS(std::pow(roughness, 1. / fitexp));
      }

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

    SWAnisoRoughness::SWAnisoRoughness(const SWAnisoRoughness &Old) :
      GeneralObjective(Old), veldiffweight(Old.veldiffweight), anisovelweight(
          Old.anisovelweight), strikediffweight(Old.strikediffweight),
          deltastrikediffweight(Old.deltastrikediffweight)
      {

      }
    SWAnisoRoughness& SWAnisoRoughness::operator=(
        const SWAnisoRoughness& source)
      {
        if (this == &source)
          return *this;
        GeneralObjective::operator=(source);
        veldiffweight = source.veldiffweight;
        anisovelweight = source.anisovelweight;
        strikediffweight = source.strikediffweight;
        deltastrikediffweight = source.deltastrikediffweight;
        return *this;
      }
  }