exec()

void stfinv::STFEngineIdentity::exec ( )

virtual

Start engine.

Reimplemented from stfinv::STFBaseEngine.

Definition at line 88 of file stfinvidentity.cc.

References stfinv::STFBaseEngine::convolvedseries(), stfinv::STFBaseEngine::convolvedsynthetic(), stfinv::STFBaseEngine::dt(), Mscaleenergy, stfinv::STFBaseEngine::npairs(), stfinv::STFBaseEngine::nreceivers(), stfinv::STFBaseEngine::recording(), stfinv::STFBaseEngine::series(), stfinv::STFBaseEngine::stf(), stfinv::STFBaseEngine::synthetic(), and stfinv::STFBaseEngine::weight().

  {
    // 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()

Here is the call graph for this function: