Преобразование Фурье на ассемблере

так ты просто присвоишь указателю адрес своего массива, без всякого выравнивания памяти.

double *signal = new double[size_data];// = {7,8,11,15,17,14,9,5,2,-1,-4};

for(int i=0;i<size_data;i++)
        signal[i] = sin(2*pi*12800*i/250000);

double *smb_signal = (double*)_aligned_malloc(sizeof(double)*(size_data<<1), 16);//new double[size_data<<1];
double *temp = smb_signal;
 
for(int i=0; i<size_data; i++)
{
    *(temp++) = signal[i];
    *(temp++) = 0;
}

//---------------------------------------------------------------------------
 
#include <vcl.h>
#include <math.h>
#include <stdio.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma link "Chart"
#pragma link "TeEngine"
#pragma link "TeeProcs"
#pragma link "Series"
#pragma resource "*.dfm"
TForm1 *Form1;
#define _USE_MATH_DEFINES
 
 
double *FFT( int *dIn, int nn )
{
    int i, j, n, m, mmax, istep;
    double tempr, tempi, wtemp, theta, wpr, wpi, wr, wi;
 
    int isign = -1;
    double *data = new double [ nn * 2 + 1 ];
 
    for( i = 0; i < nn; i++ )
    {
        data[ i * 2 ] = 0;
        data[ i * 2 + 1 ] = dIn[ i ];
    }
 
    n = nn << 1;
    j = 1;
    i = 1;
    while( i < n )
    {
        if( j > i )
        {
            tempr = data[ i ]; data[ i ] = data[ j ]; data[ j ] = tempr;
            tempr = data[ i + 1 ]; data[ i + 1 ] = data[ j + 1 ]; data[ j + 1 ] = tempr;
        }
        m = n >> 1;
        while( ( m >= 2 ) && ( j > m ) )
        {
            j = j - m;
            m = m >> 1;
        }
        j = j + m;
        i = i + 2;
    }
    mmax = 2;
    while( n > mmax )
    {
        istep = 2 * mmax;
        theta = 2.0 * M_PI / ( isign * mmax );
        wtemp = sin( 0.5 * theta );
        wpr = -2.0 * wtemp * wtemp;
        wpi = sin( theta );
        wr = 1.0;
        wi = 0.0;
        m = 1;
        while( m < mmax )
        {
            i = m;
            while( i < n )
            {
                j = i + mmax;
                tempr = wr * data[ j ] - wi * data[ j + 1 ];
                tempi = wr * data[ j + 1 ] + wi * data[ j ];
                data[ j ] = data[ i ] - tempr;
                data[ j + 1 ] = data[ i + 1 ] - tempi;
                data[ i ] = data[ i ] + tempr;
                data[ i + 1 ] = data[ i + 1 ] + tempi;
                i = i + istep;
            }
            wtemp = wr;
            wr = wtemp * wpr - wi * wpi + wr;
            wi = wi * wpr + wtemp * wpi + wi;
            m = m + 2;
        }
        mmax = istep;
    }
    double *dOut = new double [ nn / 2 ];
 
    for( i = 0; i < ( nn / 2 ); i++ )
    {
        dOut[ i ] = sqrt( data[ i * 2 ] * data[ i * 2 ] + data[ i * 2 + 1 ] * data[ i * 2 + 1 ] );
    }
 
    delete []data;
    return dOut;
}
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
    : TForm(Owner)
{
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button2Click(TObject *Sender)
{
 int *dsin = new int [ 1801 ], *pdsin = dsin;
 
    for( double x = 0.; x <= 10. * M_PI; x += M_PI / 180.0 )
    {
        *pdsin++ =  sin( x ) * 1000. + cos( 0.5 * x ) * 1000.;
        Chart1->Series[0]->AddXY(x, sin( x ) * 1000. + cos( 0.5 * x ) * 1000.);
    }
 
    double *dfourier = FFT( dsin, 1024 );
    delete []dsin;
 
    for( int i = 0; i < 512; i++ )
    {
        printf( "%d\t%lf\r\n", i, dfourier[ i ] );
    }
 
    for( int i = 0; i < 512; i++ )
    {
    Chart2->Series[0]->AddXY(i,dfourier[ i ]);
    }
 
    delete []dfourier;
 
}
//---------------------------------------------------------------------------

#include <malloc.h>
#include <math.h>
#include <stdio.h>
 
extern "C" void FFT_SSE2( double *in_data, double *out_data, int fftFrameSize);
 
const double pi = 3.1415926535897932384626433832795;
int size_data = 1 << 10;    // 1024
 
int main(int argc, char* argv[])
{
    double *in_signal = (double*) _aligned_malloc(sizeof(double)*(size_data<<1), 16);
    double *out_signal = (double*) _aligned_malloc(sizeof(double)*(size_data<<1), 16);
 
    for(int i=0; i<size_data*2; i+=2)
    {
        in_signal[i] = sin(2*pi*1000*i/44100);  
        in_signal[i+1] = .0;
    }
 
    FFT_SSE2(in_signal, out_signal, size_data);
 
    FILE *f = fopen("spec.txt", "w+");
    for(int i=0; i<size_data*2; i+=2)
        fprintf(f, "%.15e%c", sqrt(out_signal[i]*out_signal[i]+out_signal[i+1]*out_signal[i+1]), ' '); 
    fclose(f);
 
    _aligned_free(in_signal);
    _aligned_free(out_signal);
    return 0;
}