convolve()

void stfinv::STFFourierDomainEngine::convolve ( )

private

convolve synthetics with stf

Convolve all synthetics with the source time function correction filter.

The filter response is available in terms of the Fourier coefficients of the correction filter as are returned by STFFourierDomainEngine::stfspec(). Convolution takes place in the Fourier domain simply by a multiplication of Fourier coefficients. Since all Fourier coefficients are expected to be scaled appropriately to represent samples for the Fourier integral transform, no scaling has to be applied here.

Definition at line 503 of file stfinvfourier.cc.

References stfinv::STFBaseEngine::convolvedseries(), stfinv::STFBaseEngine::convolvedsynthetic(), Mfftengineinput, Mfftengineoutput, nfreq(), stfinv::STFBaseEngine::npairs(), stfinv::STFBaseEngine::nreceivers(), stfcoeff(), stfspec(), stfinv::STFBaseEngine::synthetic(), and syntheticspec().

Referenced by fftoutput().

  {
    TAspectrum synthetic=this->syntheticspec();
    TAspectrum convolvedsynthetic
      =aff::subarray(Mfftengineoutput.spectrum())()(0,this->nreceivers()-1);
    convolvedsynthetic.shape().setfirst(0,0);
    convolvedsynthetic.shape().setfirst(1,0);
    TAspectrum stfcoeff=this->stfspec();

    // cycle through all receivers
    for (unsigned int i=0; i<this->nreceivers(); ++i)
    {
      for (unsigned int j=0; j<this->nfreq(); ++j)
      {
        convolvedsynthetic(j,i)=synthetic(j,i)*stfcoeff(j);
      } // for (unsigned int j=0; j<this->nfreq(); ++j)
    } // for (unsigned int i=0; i<this->nreceivers(); ++i)

    // copy signals for additional pairs
    if (this->npairs()>0)
    {
      // reference to input series Fourier coefficients of additional pairs
      TAspectrum inspecarray=Mfftengineinput.spectrum();
      TAspectrum seriesspecarray
        =aff::subarray(inspecarray)()(2*this->nreceivers(),
                                         (2*this->nreceivers()
                                          +this->npairs()-1));
      seriesspecarray.shape().setfirst(0,0);
      seriesspecarray.shape().setfirst(1,0);

      // reference to convolved series Fourier coefficients of additional
      // pairs
      TAspectrum convolvedseries
        =aff::subarray(Mfftengineoutput.spectrum())()(this->nreceivers()+1,
                                                      this->nreceivers()
                                                      +this->npairs());
      convolvedseries.shape().setfirst(0,0);
      convolvedseries.shape().setfirst(1,0);
      // cycle through all pairs
      for (unsigned int i=0; i<this->npairs(); ++i)
      {
        for (unsigned int j=0; j<this->nfreq(); ++j)
        {
          convolvedseries(j,i)=seriesspecarray(j,i)*stfcoeff(j);
        } // for (unsigned int j=0; j<this->nfreq(); ++j)
      } // for (unsigned int i=0; i<this->npairs(); ++i)
    } // if (this->npairs()>0)

  } // void STFFourierDomainEngine::convolve() 

Here is the call graph for this function:

Here is the caller graph for this function: