exec()

void stfinv::STFEngineNormalize::exec ( )

virtual

Start engine.

Reimplemented from stfinv::STFFourierDomainEngine.

Definition at line 79 of file stfinvnormalize.cc.

References stfinv::STFFourierDomainEngine::fftinput(), stfinv::STFFourierDomainEngine::fftoutput(), Mwaterlevel, stfinv::STFFourierDomainEngine::nfreq(), stfinv::STFBaseEngine::nreceivers(), stfinv::STFFourierDomainEngine::recordingcoeff(), stfinv::STFFourierDomainEngine::stfcoeff(), stfinv::STFFourierDomainEngine::syntheticcoeff(), stfinv::STFBaseEngine::weight(), and stfinv::STFBaseEngine::weights().

  {
    // read signals and calculate FFT
    this->fftinput();

    // cycle through frequencies
    // calculate weighted energy of all synthetic traces
    Tseries syntheticenergy(0,this->nfreq()-1);
    syntheticenergy=0.;
    double totalsyntheticenergy=0.;
    for (unsigned int i=0; i<this->nfreq(); ++i)
    {
      TAspectrum syntheticref=this->syntheticcoeff(i);
      for (unsigned int j=0; j<this->nreceivers(); ++j)
      {
        syntheticenergy(i) += this->weight(j)*this->weight(j)
          *std::norm(syntheticref(j));
      }
      totalsyntheticenergy += syntheticenergy(i);
    }

    // calculate waterlevel
    double waterlevel=Mwaterlevel*totalsyntheticenergy/this->nfreq();

    double weightsum=aff::func::sqrsum(this->weights());

    for (unsigned int i=0; i<this->nfreq(); ++i)
    {
      double recordingenergy=0.;
      TAspectrum recordingref=this->recordingcoeff(i);
      TAspectrum syntheticref=this->syntheticcoeff(i);
      for (unsigned int j=0; j<this->nreceivers(); ++j)
      {
        recordingenergy += this->weight(j)*this->weight(j)
          *std::norm(recordingref(j));
      }

      double modulus
        =std::sqrt(recordingenergy/(syntheticenergy(i)+waterlevel));

      double phase = 0.;
      for (unsigned int j=0; j<this->nreceivers(); ++j)
      {
        phase += this->weight(j)*
          (std::arg(recordingref(j))-std::arg(syntheticref(j)));
      }
      phase /= weightsum;

      this->stfcoeff(i)=std::polar(modulus, phase);
    } // for (unsigned int i=0; i<this->nfreq(); ++i)

    // provide results to user
    this->fftoutput();
  } // void STFEngineNormalize::exec()

Here is the call graph for this function: