/*** FAAD2 - Freeware Advanced Audio (AAC) Decoder including SBR decoding** Copyright (C) 2003-2005 M. Bakker, Nero AG, http://www.nero.com**** This program is free software; you can redistribute it and/or modify** it under the terms of the GNU General Public License as published by** the Free Software Foundation; either version 2 of the License, or** (at your option) any later version.**** This program is distributed in the hope that it will be useful,** but WITHOUT ANY WARRANTY; without even the implied warranty of** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the** GNU General Public License for more details.**** You should have received a copy of the GNU General Public License** along with this program; if not, write to the Free Software** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.**** Any non-GPL usage of this software or parts of this software is strictly** forbidden.**** The "appropriate copyright message" mentioned in section 2c of the GPLv2** must read: "Code from FAAD2 is copyright (c) Nero AG, www.nero.com"**** Commercial non-GPL licensing of this software is possible.** For more info contact Nero AG through Mpeg4AAClicense@nero.com.**** $Id: cfft.c,v 1.35 2007/11/01 12:33:29 menno Exp $**/ /* * Algorithmically based on Fortran-77 FFTPACK * by Paul N. Swarztrauber(Version 4, 1985). * * Does even sized fft only */ /* isign is +1 for backward and -1 for forward transforms */ #include "common.h"#include "structs.h" #include "cfft.h"#include "cfft_tab.h"  /* static function declarations */static void passf2pos(const uint16_t ido, const uint16_t l1, const complex_t *cc,                      complex_t *ch, const complex_t *wa);static void passf2neg(const uint16_t ido, const uint16_t l1, const complex_t *cc,                      complex_t *ch, const complex_t *wa);static void passf3(const uint16_t ido, const uint16_t l1, const complex_t *cc,                   complex_t *ch, const complex_t *wa1, const complex_t *wa2, const int8_t isign);static void passf4pos(const uint16_t ido, const uint16_t l1, const complex_t *cc, complex_t *ch,                      const complex_t *wa1, const complex_t *wa2, const complex_t *wa3);static void passf4neg(const uint16_t ido, const uint16_t l1, const complex_t *cc, complex_t *ch,                      const complex_t *wa1, const complex_t *wa2, const complex_t *wa3);static void passf5(const uint16_t ido, const uint16_t l1, const complex_t *cc, complex_t *ch,                   const complex_t *wa1, const complex_t *wa2, const complex_t *wa3,                   const complex_t *wa4, const int8_t isign);static void cffti1(uint16_t n, complex_t *wa, uint16_t *ifac);  /*----------------------------------------------------------------------   passf2, passf3, passf4, passf5. Complex FFT passes fwd and bwd.  ----------------------------------------------------------------------*/ static void passf2pos(const uint16_t ido, const uint16_t l1, const complex_t *cc,                      complex_t *ch, const complex_t *wa){    uint16_t i, k, ah, ac;     if (ido == 1)    {#ifdef FUZZING_BUILD_MODE_UNSAFE_FOR_PRODUCTION        // TODO: remove this code once fuzzer proves it is totally unreahable        // For supported frame lengths that have odd number of factor 2 it is        // never the last factor; consequently `ido` should never be 1.        __builtin_trap();        /*#endif        for (k = 0; k < l1; k++)        {            ah = 2*k;            ac = 4*k;             RE(ch[ah])    = RE(cc[ac]) + RE(cc[ac+1]);            RE(ch[ah+l1]) = RE(cc[ac]) - RE(cc[ac+1]);            IM(ch[ah])    = IM(cc[ac]) + IM(cc[ac+1]);            IM(ch[ah+l1]) = IM(cc[ac]) - IM(cc[ac+1]);        }#ifdef FUZZING_BUILD_MODE_UNSAFE_FOR_PRODUCTION        */#endif    } else {        for (k = 0; k < l1; k++)        {            ah = k*ido;            ac = 2*k*ido;             for (i = 0; i < ido; i++)            {                complex_t t2;                 RE(ch[ah+i]) = RE(cc[ac+i]) + RE(cc[ac+i+ido]);                RE(t2)       = RE(cc[ac+i]) - RE(cc[ac+i+ido]);                 IM(ch[ah+i]) = IM(cc[ac+i]) + IM(cc[ac+i+ido]);                IM(t2)       = IM(cc[ac+i]) - IM(cc[ac+i+ido]); #if 1                ComplexMult(&IM(ch[ah+i+l1*ido]), &RE(ch[ah+i+l1*ido]),                    IM(t2), RE(t2), RE(wa[i]), IM(wa[i]));#else                ComplexMult(&RE(ch[ah+i+l1*ido]), &IM(ch[ah+i+l1*ido]),                    RE(t2), IM(t2), RE(wa[i]), IM(wa[i]));#endif            }        }    }} static void passf2neg(const uint16_t ido, const uint16_t l1, const complex_t *cc,                      complex_t *ch, const complex_t *wa){    uint16_t i, k, ah, ac;     if (ido == 1)    {#ifdef FUZZING_BUILD_MODE_UNSAFE_FOR_PRODUCTION        // TODO: remove this code once fuzzer proves it is totally unreahable        // For supported frame lengths that have odd number of factor 2 it is        // never the last factor; consequently `ido` should never be 1.        __builtin_trap();        /*#endif        for (k = 0; k < l1; k++)        {            ah = 2*k;            ac = 4*k;             RE(ch[ah])    = RE(cc[ac]) + RE(cc[ac+1]);            RE(ch[ah+l1]) = RE(cc[ac]) - RE(cc[ac+1]);            IM(ch[ah])    = IM(cc[ac]) + IM(cc[ac+1]);            IM(ch[ah+l1]) = IM(cc[ac]) - IM(cc[ac+1]);        }#ifdef FUZZING_BUILD_MODE_UNSAFE_FOR_PRODUCTION        */#endif    } else {        for (k = 0; k < l1; k++)        {            ah = k*ido;            ac = 2*k*ido;             for (i = 0; i < ido; i++)            {                complex_t t2;                 RE(ch[ah+i]) = RE(cc[ac+i]) + RE(cc[ac+i+ido]);                RE(t2)       = RE(cc[ac+i]) - RE(cc[ac+i+ido]);                 IM(ch[ah+i]) = IM(cc[ac+i]) + IM(cc[ac+i+ido]);                IM(t2)       = IM(cc[ac+i]) - IM(cc[ac+i+ido]); #if 1                ComplexMult(&RE(ch[ah+i+l1*ido]), &IM(ch[ah+i+l1*ido]),                    RE(t2), IM(t2), RE(wa[i]), IM(wa[i]));#else                ComplexMult(&IM(ch[ah+i+l1*ido]), &RE(ch[ah+i+l1*ido]),                    IM(t2), RE(t2), RE(wa[i]), IM(wa[i]));#endif            }        }    }}  static void passf3(const uint16_t ido, const uint16_t l1, const complex_t *cc,                   complex_t *ch, const complex_t *wa1, const complex_t *wa2,                   const int8_t isign){    static real_t taur = FRAC_CONST(-0.5);    static real_t taui = FRAC_CONST(0.866025403784439);    uint16_t i, k, ac, ah;    complex_t c2, c3, d2, d3, t2;     if (ido == 1)    {#ifdef FUZZING_BUILD_MODE_UNSAFE_FOR_PRODUCTION        // TODO: remove this code once fuzzer proves it is totally unreahable        // 3 is never the the biggest factor for supported frame lengths;        // consequently `ido` should never be 1.        __builtin_trap();        /*#endif        if (isign == 1)        {            for (k = 0; k < l1; k++)            {                ac = 3*k+1;                ah = k;                 RE(t2) = RE(cc[ac]) + RE(cc[ac+1]);                IM(t2) = IM(cc[ac]) + IM(cc[ac+1]);                RE(c2) = RE(cc[ac-1]) + MUL_F(RE(t2),taur);                IM(c2) = IM(cc[ac-1]) + MUL_F(IM(t2),taur);                 RE(ch[ah]) = RE(cc[ac-1]) + RE(t2);                IM(ch[ah]) = IM(cc[ac-1]) + IM(t2);                 RE(c3) = MUL_F((RE(cc[ac]) - RE(cc[ac+1])), taui);                IM(c3) = MUL_F((IM(cc[ac]) - IM(cc[ac+1])), taui);                 RE(ch[ah+l1]) = RE(c2) - IM(c3);                IM(ch[ah+l1]) = IM(c2) + RE(c3);                RE(ch[ah+2*l1]) = RE(c2) + IM(c3);                IM(ch[ah+2*l1]) = IM(c2) - RE(c3);            }        } else {            for (k = 0; k < l1; k++)            {                ac = 3*k+1;                ah = k;                 RE(t2) = RE(cc[ac]) + RE(cc[ac+1]);                IM(t2) = IM(cc[ac]) + IM(cc[ac+1]);                RE(c2) = RE(cc[ac-1]) + MUL_F(RE(t2),taur);                IM(c2) = IM(cc[ac-1]) + MUL_F(IM(t2),taur);                 RE(ch[ah]) = RE(cc[ac-1]) + RE(t2);                IM(ch[ah]) = IM(cc[ac-1]) + IM(t2);                 RE(c3) = MUL_F((RE(cc[ac]) - RE(cc[ac+1])), taui);                IM(c3) = MUL_F((IM(cc[ac]) - IM(cc[ac+1])), taui);                 RE(ch[ah+l1]) = RE(c2) + IM(c3);                IM(ch[ah+l1]) = IM(c2) - RE(c3);                RE(ch[ah+2*l1]) = RE(c2) - IM(c3);                IM(ch[ah+2*l1]) = IM(c2) + RE(c3);            }        }#ifdef FUZZING_BUILD_MODE_UNSAFE_FOR_PRODUCTION        */#endif    } else {        if (isign == 1)        {            for (k = 0; k < l1; k++)            {                for (i = 0; i < ido; i++)                {                    ac = i + (3*k+1)*ido;                    ah = i + k * ido;                     RE(t2) = RE(cc[ac]) + RE(cc[ac+ido]);                    RE(c2) = RE(cc[ac-ido]) + MUL_F(RE(t2),taur);                    IM(t2) = IM(cc[ac]) + IM(cc[ac+ido]);                    IM(c2) = IM(cc[ac-ido]) + MUL_F(IM(t2),taur);                     RE(ch[ah]) = RE(cc[ac-ido]) + RE(t2);                    IM(ch[ah]) = IM(cc[ac-ido]) + IM(t2);                     RE(c3) = MUL_F((RE(cc[ac]) - RE(cc[ac+ido])), taui);                    IM(c3) = MUL_F((IM(cc[ac]) - IM(cc[ac+ido])), taui);                     RE(d2) = RE(c2) - IM(c3);                    IM(d3) = IM(c2) - RE(c3);                    RE(d3) = RE(c2) + IM(c3);                    IM(d2) = IM(c2) + RE(c3); #if 1                    ComplexMult(&IM(ch[ah+l1*ido]), &RE(ch[ah+l1*ido]),                        IM(d2), RE(d2), RE(wa1[i]), IM(wa1[i]));                    ComplexMult(&IM(ch[ah+2*l1*ido]), &RE(ch[ah+2*l1*ido]),                        IM(d3), RE(d3), RE(wa2[i]), IM(wa2[i]));#else                    ComplexMult(&RE(ch[ah+l1*ido]), &IM(ch[ah+l1*ido]),                        RE(d2), IM(d2), RE(wa1[i]), IM(wa1[i]));                    ComplexMult(&RE(ch[ah+2*l1*ido]), &IM(ch[ah+2*l1*ido]),                        RE(d3), IM(d3), RE(wa2[i]), IM(wa2[i]));#endif                }            }        } else {            for (k = 0; k < l1; k++)            {                for (i = 0; i < ido; i++)                {                    ac = i + (3*k+1)*ido;                    ah = i + k * ido;                     RE(t2) = RE(cc[ac]) + RE(cc[ac+ido]);                    RE(c2) = RE(cc[ac-ido]) + MUL_F(RE(t2),taur);                    IM(t2) = IM(cc[ac]) + IM(cc[ac+ido]);                    IM(c2) = IM(cc[ac-ido]) + MUL_F(IM(t2),taur);                     RE(ch[ah]) = RE(cc[ac-ido]) + RE(t2);                    IM(ch[ah]) = IM(cc[ac-ido]) + IM(t2);                     RE(c3) = MUL_F((RE(cc[ac]) - RE(cc[ac+ido])), taui);                    IM(c3) = MUL_F((IM(cc[ac]) - IM(cc[ac+ido])), taui);                     RE(d2) = RE(c2) + IM(c3);                    IM(d3) = IM(c2) + RE(c3);                    RE(d3) = RE(c2) - IM(c3);                    IM(d2) = IM(c2) - RE(c3); #if 1                    ComplexMult(&RE(ch[ah+l1*ido]), &IM(ch[ah+l1*ido]),                        RE(d2), IM(d2), RE(wa1[i]), IM(wa1[i]));                    ComplexMult(&RE(ch[ah+2*l1*ido]), &IM(ch[ah+2*l1*ido]),                        RE(d3), IM(d3), RE(wa2[i]), IM(wa2[i]));#else                    ComplexMult(&IM(ch[ah+l1*ido]), &RE(ch[ah+l1*ido]),                        IM(d2), RE(d2), RE(wa1[i]), IM(wa1[i]));                    ComplexMult(&IM(ch[ah+2*l1*ido]), &RE(ch[ah+2*l1*ido]),                        IM(d3), RE(d3), RE(wa2[i]), IM(wa2[i]));#endif                }            }        }    }}  static void passf4pos(const uint16_t ido, const uint16_t l1, const complex_t *cc,                      complex_t *ch, const complex_t *wa1, const complex_t *wa2,                      const complex_t *wa3){    uint16_t i, k, ac, ah;     if (ido == 1)    {        for (k = 0; k < l1; k++)        {            complex_t t1, t2, t3, t4;             ac = 4*k;            ah = k;             RE(t2) = RE(cc[ac])   + RE(cc[ac+2]);            RE(t1) = RE(cc[ac])   - RE(cc[ac+2]);            IM(t2) = IM(cc[ac])   + IM(cc[ac+2]);            IM(t1) = IM(cc[ac])   - IM(cc[ac+2]);            RE(t3) = RE(cc[ac+1]) + RE(cc[ac+3]);            IM(t4) = RE(cc[ac+1]) - RE(cc[ac+3]);            IM(t3) = IM(cc[ac+3]) + IM(cc[ac+1]);            RE(t4) = IM(cc[ac+3]) - IM(cc[ac+1]);             RE(ch[ah])      = RE(t2) + RE(t3);            RE(ch[ah+2*l1]) = RE(t2) - RE(t3);             IM(ch[ah])      = IM(t2) + IM(t3);            IM(ch[ah+2*l1]) = IM(t2) - IM(t3);             RE(ch[ah+l1])   = RE(t1) + RE(t4);            RE(ch[ah+3*l1]) = RE(t1) - RE(t4);             IM(ch[ah+l1])   = IM(t1) + IM(t4);            IM(ch[ah+3*l1]) = IM(t1) - IM(t4);        }    } else {        for (k = 0; k < l1; k++)        {            ac = 4*k*ido;            ah = k*ido;             for (i = 0; i < ido; i++)            {                complex_t c2, c3, c4, t1, t2, t3, t4;                 RE(t2) = RE(cc[ac+i]) + RE(cc[ac+i+2*ido]);                RE(t1) = RE(cc[ac+i]) - RE(cc[ac+i+2*ido]);                IM(t2) = IM(cc[ac+i]) + IM(cc[ac+i+2*ido]);                IM(t1) = IM(cc[ac+i]) - IM(cc[ac+i+2*ido]);                RE(t3) = RE(cc[ac+i+ido]) + RE(cc[ac+i+3*ido]);                IM(t4) = RE(cc[ac+i+ido]) - RE(cc[ac+i+3*ido]);                IM(t3) = IM(cc[ac+i+3*ido]) + IM(cc[ac+i+ido]);                RE(t4) = IM(cc[ac+i+3*ido]) - IM(cc[ac+i+ido]);                 RE(c2) = RE(t1) + RE(t4);                RE(c4) = RE(t1) - RE(t4);                 IM(c2) = IM(t1) + IM(t4);                IM(c4) = IM(t1) - IM(t4);                 RE(ch[ah+i]) = RE(t2) + RE(t3);                RE(c3)       = RE(t2) - RE(t3);                 IM(ch[ah+i]) = IM(t2) + IM(t3);                IM(c3)       = IM(t2) - IM(t3); #if 1                ComplexMult(&IM(ch[ah+i+l1*ido]), &RE(ch[ah+i+l1*ido]),                    IM(c2), RE(c2), RE(wa1[i]), IM(wa1[i]));                ComplexMult(&IM(ch[ah+i+2*l1*ido]), &RE(ch[ah+i+2*l1*ido]),                    IM(c3), RE(c3), RE(wa2[i]), IM(wa2[i]));                ComplexMult(&IM(ch[ah+i+3*l1*ido]), &RE(ch[ah+i+3*l1*ido]),                    IM(c4), RE(c4), RE(wa3[i]), IM(wa3[i]));#else                ComplexMult(&RE(ch[ah+i+l1*ido]), &IM(ch[ah+i+l1*ido]),                    RE(c2), IM(c2), RE(wa1[i]), IM(wa1[i]));                ComplexMult(&RE(ch[ah+i+2*l1*ido]), &IM(ch[ah+i+2*l1*ido]),                    RE(c3), IM(c3), RE(wa2[i]), IM(wa2[i]));                ComplexMult(&RE(ch[ah+i+3*l1*ido]), &IM(ch[ah+i+3*l1*ido]),                    RE(c4), IM(c4), RE(wa3[i]), IM(wa3[i]));#endif            }        }    }} static void passf4neg(const uint16_t ido, const uint16_t l1, const complex_t *cc,                      complex_t *ch, const complex_t *wa1, const complex_t *wa2,                      const complex_t *wa3){    uint16_t i, k, ac, ah;     if (ido == 1)    {        for (k = 0; k < l1; k++)        {            complex_t t1, t2, t3, t4;             ac = 4*k;            ah = k;             RE(t2) = RE(cc[ac])   + RE(cc[ac+2]);            RE(t1) = RE(cc[ac])   - RE(cc[ac+2]);            IM(t2) = IM(cc[ac])   + IM(cc[ac+2]);            IM(t1) = IM(cc[ac])   - IM(cc[ac+2]);            RE(t3) = RE(cc[ac+1]) + RE(cc[ac+3]);            IM(t4) = RE(cc[ac+1]) - RE(cc[ac+3]);            IM(t3) = IM(cc[ac+3]) + IM(cc[ac+1]);            RE(t4) = IM(cc[ac+3]) - IM(cc[ac+1]);             RE(ch[ah])      = RE(t2) + RE(t3);            RE(ch[ah+2*l1]) = RE(t2) - RE(t3);             IM(ch[ah])      = IM(t2) + IM(t3);            IM(ch[ah+2*l1]) = IM(t2) - IM(t3);             RE(ch[ah+l1])   = RE(t1) - RE(t4);            RE(ch[ah+3*l1]) = RE(t1) + RE(t4);             IM(ch[ah+l1])   = IM(t1) - IM(t4);            IM(ch[ah+3*l1]) = IM(t1) + IM(t4);        }    } else {        for (k = 0; k < l1; k++)        {            ac = 4*k*ido;            ah = k*ido;             for (i = 0; i < ido; i++)            {                complex_t c2, c3, c4, t1, t2, t3, t4;                 RE(t2) = RE(cc[ac+i]) + RE(cc[ac+i+2*ido]);                RE(t1) = RE(cc[ac+i]) - RE(cc[ac+i+2*ido]);                IM(t2) = IM(cc[ac+i]) + IM(cc[ac+i+2*ido]);                IM(t1) = IM(cc[ac+i]) - IM(cc[ac+i+2*ido]);                RE(t3) = RE(cc[ac+i+ido]) + RE(cc[ac+i+3*ido]);                IM(t4) = RE(cc[ac+i+ido]) - RE(cc[ac+i+3*ido]);                IM(t3) = IM(cc[ac+i+3*ido]) + IM(cc[ac+i+ido]);                RE(t4) = IM(cc[ac+i+3*ido]) - IM(cc[ac+i+ido]);                 RE(c2) = RE(t1) - RE(t4);                RE(c4) = RE(t1) + RE(t4);                 IM(c2) = IM(t1) - IM(t4);                IM(c4) = IM(t1) + IM(t4);                 RE(ch[ah+i]) = RE(t2) + RE(t3);                RE(c3)       = RE(t2) - RE(t3);                 IM(ch[ah+i]) = IM(t2) + IM(t3);                IM(c3)       = IM(t2) - IM(t3); #if 1                ComplexMult(&RE(ch[ah+i+l1*ido]), &IM(ch[ah+i+l1*ido]),                    RE(c2), IM(c2), RE(wa1[i]), IM(wa1[i]));                ComplexMult(&RE(ch[ah+i+2*l1*ido]), &IM(ch[ah+i+2*l1*ido]),                    RE(c3), IM(c3), RE(wa2[i]), IM(wa2[i]));                ComplexMult(&RE(ch[ah+i+3*l1*ido]), &IM(ch[ah+i+3*l1*ido]),                    RE(c4), IM(c4), RE(wa3[i]), IM(wa3[i]));#else                ComplexMult(&IM(ch[ah+i+l1*ido]), &RE(ch[ah+i+l1*ido]),                    IM(c2), RE(c2), RE(wa1[i]), IM(wa1[i]));                ComplexMult(&IM(ch[ah+i+2*l1*ido]), &RE(ch[ah+i+2*l1*ido]),                    IM(c3), RE(c3), RE(wa2[i]), IM(wa2[i]));                ComplexMult(&IM(ch[ah+i+3*l1*ido]), &RE(ch[ah+i+3*l1*ido]),                    IM(c4), RE(c4), RE(wa3[i]), IM(wa3[i]));#endif            }        }    }} 

Lines 1–500 of 1048. Load more