transitionmixer.cc

Go to the documentation of this file.

#define TF_TRANSITIONMIXER_CC_VERSION \
  "TF_TRANSITIONMIXER_CC   V1.0"

#include "lisousi.h"
#include "functions.h"

TFourier::Tseries transitionmixer(TFourier::Tseries singlevelocityseries,
                                  TFourier::Tseries directwaveseries,
                                  const Parameters& par,
                                  const Options& opt)
{
  TFourier::Tseries series;

  // trim number of samples
  if (opt.transition1<=par.offset)
  {
    directwaveseries.setlastindex(par.nsamples+
                                  directwaveseries.f()-1);
  }
  if (opt.transition2>=par.offset)
  {
    singlevelocityseries.setlastindex(par.nsamples+
                                      singlevelocityseries.f()-1);
  }
  // report offset
  if (opt.verbose)
  {
    cout << "  offset " << par.offset << "m: ";
  }
  // enter mixing operation
  if (opt.transition1>=par.offset)
  {
    series=singlevelocityseries;
    if (opt.verbose)
    {
      cout << "100% single velocity transformation" << endl;
    }
  }
  else if (opt.transition2<=par.offset)
  {
    series=directwaveseries;
    if (opt.verbose)
    {
      cout << "100% direct wave transformation" << endl;
    }
  }
  else
  {
    const double dwf=(par.offset-opt.transition1)/
      (opt.transition2-opt.transition1);
    const double svf=1.-dwf;
    series=Tseries(par.nsamples);
    if (opt.verbose)
    {
      cout << 100.*svf << "% single velocity and "
        << 100.*dwf << "% direct wave transformation" << endl;
    }

    aff::Browser<Tseries> DWB(directwaveseries);
    aff::Browser<Tseries> SVB(singlevelocityseries);
    aff::Iterator<Tseries> SI(series);
    while (DWB.valid() && SVB.valid() && SI.valid())
    {
      (*SI) = dwf*(*DWB) + svf*(*SVB);
      ++DWB;
      ++SVB;
      ++SI;
    }
  }
  return(series);
}

/* ----- END OF transitionmixer.cc ----- */