/[marvin]/marvin/src/libsip/algo/fft.cc
ViewVC logotype

Diff of /marvin/src/libsip/algo/fft.cc

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.6 by exos, Fri Sep 12 03:06:18 2003 UTC revision 1.7 by vishnou, Fri Sep 19 18:03:15 2003 UTC
# Line 53  fft::fft ( DiscSignal &signal, bool inv) Line 53  fft::fft ( DiscSignal &signal, bool inv)
53    double angle = -2.0 * PI;    double angle = -2.0 * PI;
54    unsigned nbi;    unsigned nbi;
55    unsigned ns = signal.length();    unsigned ns = signal.length();
56    double *result_r_tmp = (double *)malloc(sizeof(double) * ns);    double *result_r_tmp = new double[ns];
57    double *result_i_tmp = (double *)malloc(sizeof(double) * ns);    double *result_i_tmp = new double[ns];
58    DiscSignal &tmpsig = signal;    DiscSignal &tmpsig = signal;
59        
60    if ( !puiss_de_2(ns) )    if ( !puiss_de_2(ns) )
61      {      {
62        unsigned l = get_puiss(ns);        unsigned l = get_puiss(ns);
63        double *tmp = (double *)malloc(sizeof(double) * ns);        double *tmp = new double[ns];
64        int i;        int i;
65        for (i = 0; i < ns; i++)        for (i = 0; i < ns; i++)
66          tmp[i] = signal[i];          tmp[i] = signal[i];
# Line 148  DiscSignal &fft::get_magnitude() Line 148  DiscSignal &fft::get_magnitude()
148  {  {
149    DiscSignal *result;    DiscSignal *result;
150    
151    result = new DiscSignal((double *)malloc(sizeof(double) * result_i->length()), result_i->length());    result = new DiscSignal(new double[result_i->length()], result_i->length());
152    for (int i = 0; i < result_i->length(); ++i)    for (int i = 0; i < result_i->length(); ++i)
153      (*result)[i] = sqrt ((*result_i)[i] * (*result_i)[i] + (*result_r)[i] * (*result_r)[i]);      (*result)[i] = sqrt ((*result_i)[i] * (*result_i)[i] + (*result_r)[i] * (*result_r)[i]);
154    

Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26