exec()

void stfinv::STFEngineFDLeastSquares::exec ( )

virtual

Start engine.

Reimplemented from stfinv::STFFourierDomainEngine.

Definition at line 111 of file stfinvfdleastsquares.cc.

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

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

    // calculate waterlevel
    double waterlevel=0.;
    for (unsigned int i=0; i<this->nreceivers(); ++i)
    {
      waterlevel += this->weight(i)*this->weight(i)
        *aff::func::sqrsum(this->synthetic(i));
    } // for (unsigned int i=0; i<this->nreceivers(); ++i)
    waterlevel *= Mwaterlevel;

    STFINV_debug(Mdebug&2, "STFEngineFDLeastSquares::exec()",
                 "waterlevel=" << waterlevel);

    STFINV_debug(Mdebug&2, "STFEngineFDLeastSquares::exec()",
                 "nfreq=" << this->nfreq());
    // cycle through frequencies
    for (unsigned int i=0; i<this->nfreq(); ++i)
    {
      STFINV_debug(Mdebug&4, "STFEngineFDLeastSquares::exec()",
                   "frequency=" << i);
      TAspectrum syntheticref=this->syntheticcoeff(i);
      TAspectrum recordingref=this->recordingcoeff(i);
      STFINV_debug(Mdebug&4, "STFEngineFDLeastSquares::exec()",
                   "syntheticref.first(0)=" << syntheticref.first(0) << " " <<
                   "syntheticref.last(0)=" << syntheticref.last(0) << " "
                   "syntheticref.first(1)=" << syntheticref.first(1) << " " <<
                   "syntheticref.last(1)=" << syntheticref.last(1) << " ");
      TAspectrum::Tvalue numerator, denominator;
      denominator=waterlevel;
      numerator=0.;
      for (unsigned int j=0; j<this->nreceivers(); ++j)
      {
        numerator += this->weight(j)*this->weight(j)
          *conj(syntheticref(j))*recordingref(j);
        denominator += this->weight(j)*this->weight(j)
          *conj(syntheticref(j))*syntheticref(j);
      }
      this->stfcoeff(i)=numerator/denominator;
    } // for (unsigned int i=0; i<this->nfreq(); ++i)

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

Here is the call graph for this function: