stfinvidentity.cc

Go to the documentation of this file.

#define STFINV_STFINVIDENTITY_CC_VERSION \
  "STFINV_STFINVIDENTITY_CC   V1.3"

#include <stfinv/stfinvidentity.h>
#include <stfinv/stfinvidentity_summary_usage.h>
#include <stfinv/stfinvidentity_description_usage.h>
#include <aff/functions/sqrsum.h>
#include <aff/seriesoperators.h>
#include <cmath>
#include <stfinv/tools.h>

namespace stfinv {

  const char* const STFEngineIdentity::ID="ident";

  const char* const STFEngineIdentity::description
    ="scale with amplitude factor";

  /*----------------------------------------------------------------------*/

  void STFEngineIdentity::help(std::ostream& os) const
  {
    STFEngineIdentity::classhelp(os);
  } // void STFEngineIdentity::help(std::ostream& os) const

  /*----------------------------------------------------------------------*/

  void STFEngineIdentity::usage(std::ostream& os) const
  {
    STFEngineIdentity::classusage(os);
  } // void STFEngineIdentity::usage(std::ostream& os) const

  /*----------------------------------------------------------------------*/

  const char* STFEngineIdentity::name() const
  {
    return("STFEngineIdentity");
  } //  const char const* STFEngineIdentity::name() const

  /*----------------------------------------------------------------------*/

  void STFEngineIdentity::initialize() 
  {
    // scale energy
    Mscaleenergy=(this->parameter("scaleenergy","false")=="true");
  } // void STFEngineIdentity::initialize()

  /*----------------------------------------------------------------------*/

  void STFEngineIdentity::exec() 
  {
    // effective amplitude factor
    double fac=1.;

    // scale to reproduce average energy if requested
    if (Mscaleenergy)
    {
      double recording_sqrsum=0.;;
      double synthetic_sqrsum=0.;;
      for (unsigned int i=0; i<this->nreceivers(); ++i)
      {
        synthetic_sqrsum 
          += aff::func::sqrsum(this->synthetic(i)) 
             * this->weight(i) * this->weight(i);
        recording_sqrsum 
          += aff::func::sqrsum(this->recording(i)) 
             * this->weight(i) * this->weight(i);
        fac = std::sqrt(recording_sqrsum/synthetic_sqrsum);
      }
    } // if (Mscaleenergy)

    Tseries stf=this->stf();
    stf=0.;
    stf(0)=fac/this->dt();
    for (unsigned int i=0; i<this->nreceivers(); ++i)
    {
      Tseries::Tcoc synthetic=this->synthetic(i);
      Tseries convolvedsynthetic=this->convolvedsynthetic(i);
      convolvedsynthetic.copyin(synthetic);
      if (Mscaleenergy) { convolvedsynthetic *= fac; }
    }
    for (unsigned int i=0; i<this->npairs(); ++i)
    {
      Tseries::Tcoc series=this->series(i);
      Tseries convolvedseries=this->convolvedseries(i);
      convolvedseries.copyin(series);
      if (Mscaleenergy) { convolvedseries *= fac; }
    }
  } // void STFEngineIdentity::exec()

  /*----------------------------------------------------------------------*/

  void STFEngineIdentity::classhelp(std::ostream& os)
  {
    os << stfinvidentity_summary_usage;
    os << std::endl;
    stfinv::tools::report_engine_ID<STFEngineIdentity>(os);
  } // void STFEngineIdentity::classhelp(std::ostream& os)

  /*----------------------------------------------------------------------*/

  void STFEngineIdentity::classusage(std::ostream& os)
  {
    os << stfinvidentity_description_usage;
    os << std::endl;
    Tbase::classusage(os);
    os << std::endl;
    stfinv::tools::report_engine_ID<STFEngineIdentity>(os);
  } // void STFEngineIdentity::classusage(std::ostream& os)

} // namespace stfinv

/* ----- END OF stfinvidentity.cc ----- */