putoutput()

void stfinv::STFFourierDomainEngine::putoutput ( )

private

copy workspace time series for convolved synthetics and stf to user memory

Read out time series data from Mfftengineoutput. The time series containers in Mfftengineoutput are expected to be larger than those used for the API of the library (because of padding).

Definition at line 455 of file stfinvfourier.cc.

References stfinv::STFBaseEngine::convolvedseries(), stfinv::STFBaseEngine::convolvedsynthetic(), stfinv::STFBaseEngine::dt(), Mfftengineoutput, stfinv::STFBaseEngine::npairs(), stfinv::STFBaseEngine::nreceivers(), and stfinv::STFBaseEngine::stf().

Referenced by fftoutput().

  {
    // scaling factor for source correction filter
    double scalingstf=Mfftengineoutput.scale_series(this->dt());
    // scaling factor for convolved synthetics
    double scalingcs=Mfftengineoutput.scale_spectrum(this->dt())
      *scalingstf;

    // cycle through receivers
    for (unsigned int i=0; i<this->nreceivers(); ++i)
    {
      // get references to time series 
      stfinv::Tseries convolvedsynthetics=this->convolvedsynthetic(i);
      // copyin function copies as many elements as possible
      convolvedsynthetics.copyin(Mfftengineoutput.series(i)*scalingcs);
    } // for (unsigned int i=0; i<this->nreceivers(); ++i)
      
    // copy stf too; get reference to series and use copyin
    stfinv::Tseries stf=this->stf();
    stf.copyin(Mfftengineoutput.series(this->nreceivers())
               *scalingstf);

    // cycle through additional time series pairs
    for (unsigned int i=0; i<this->npairs(); ++i)
    {
      // get references to time series 
      stfinv::Tseries convolvedseries=this->convolvedseries(i);
      // copyin function copies as many elements as possible
      convolvedseries.copyin(Mfftengineoutput.series(1+i+this->nreceivers())
                             *scalingcs);
    } // for (unsigned int i=0; i<this->npairs(); ++i)

  } // void STFFourierDomainEngine::putoutput() 

Here is the call graph for this function:

Here is the caller graph for this function: