dft()

void fourier::adft_func::dft ( double  w, cplx &  result, cplx &  err  )

inline

Computes the cubically corrected DFT.

Purpose

Computes the cubically corrected DFT by integrating over the individual intervals.

Parameters

w	frequency \(\omega\) in the Fourier integral
result	value of the Fourier integral
err	estimate of the error

Definition at line 317 of file fourier.hpp.

References fourier::dft_cplx().

Referenced by cntr::green_equilibrium_mat(), and cntr::smooth_box::set_norm().

                                            {
      int i,n_i;
      double a,b;
      cplx res_i,err_i,*f_i;
      if(n_<=0) {std::cerr << "adft_cplx: no interval in f" << std::endl; abort();}
      result=0;
      err=0;
  
      for(i=0;i<n_;i++){
         f_i=interval_data(i);
         a=a_[i];
         b=b_[i];
         n_i=len_[i];
         fourier::dft_cplx(w,n_i,a,b,f_i,res_i,err_i);
         result += res_i;
         err += err_i;
      }
   }

Here is the call graph for this function:

Here is the caller graph for this function: