public static int fnFIE_fft_fwd_DDtoDD_in(
	FHANDLE hFFT,
	DOUBLE_PTR dpRe,
	DOUBLE_PTR dpIm,
	int iOrder
)

Public Shared Function fnFIE_fft_fwd_DDtoDD_in ( 
	hFFT As FHANDLE,
	dpRe As DOUBLE_PTR,
	dpIm As DOUBLE_PTR,
	iOrder As Integer
) As Integer

//    $Revision: 1.1 $

using System;
using System.Collections.Generic;
using System.Text;
using fvalgcli;

namespace TC
{
    public partial class FIE
    {
        /// <summary>
        /// splitted format の複素数データの１次元順方向フーリエ変換（インプレース処理).
        /// </summary>
        [FvPluginExecute]
        public void fnFIE_fft_fwd_DDtoDD_in()
        {
            int status = (int)f_err.F_ERR_NONE;

            int iOrder = 5;                                // データ長の次数.
            int length = 1 << iOrder;                    // 入力データの長さ.
            FHANDLE hFFT = FHANDLE.Zero;                // FFTハンドル.
            DOUBLE_PTR pDataRe = DOUBLE_PTR.Zero;        // 入力データの実部(＝出力データの実部).
            DOUBLE_PTR pDataIm = DOUBLE_PTR.Zero;        // 入力データの虚部(＝出力データの虚部).
            DOUBLE_PTR pDataReOrg = DOUBLE_PTR.Zero;    // 結果比較用データの実部.
            DOUBLE_PTR pDataImOrg = DOUBLE_PTR.Zero;    // 結果比較用データの虚部.

            try
            {
                // データの作成.                
                pDataRe = DOUBLE_PTR.alloc(length);
                pDataIm = DOUBLE_PTR.alloc(length);
                pDataReOrg = DOUBLE_PTR.alloc(length);
                pDataImOrg = DOUBLE_PTR.alloc(length);
                for (int i = 0; i < length; i++)
                {
                    double image = (double)i * 10 + 0.1;
                    double real = (double)i * 10 + 0.2;
                    pDataRe[i] = real;
                    pDataIm[i] = image;
                    pDataReOrg[i] = real;
                    pDataImOrg[i] = image;
                }

                // FFTハンドルの初期化.
                status = api.fnFIE_fft_1D_alloc(ref hFFT, iOrder, f_fft_direction.F_FFT_INVERSE, f_fft_normalize_type.F_FFT_DIV_BY_SQRTN, f_fft_data_type.F_1D_FFT_FCOMPLEX, f_fft_data_type.F_1D_FFT_FCOMPLEX);
                Assert.IsTrue(status == (int)f_err.F_ERR_NONE, "エラーが発生しました。({0})", (f_err)status);

                // 処理の実行.
                status = api.fnFIE_fft_fwd_DDtoDD_in(hFFT, pDataRe, pDataIm, iOrder);
                Assert.IsTrue(status == (int)f_err.F_ERR_NONE, "エラーが発生しました。({0})", (f_err)status);

                // 結果比較.
                bool result = false;
                for (int i = 0; i < length; i++)
                {
                    if (pDataRe[i] != pDataReOrg[i]
                        || pDataIm[i] != pDataImOrg[i])
                    {
                        result = true;
                        break;
                    }
                }
                if (!result)
                {
                    Assert.Fail("処理前と値が変わっていません");
                }
            }
            finally
            {
                // オブジェクトの開放.
                hFFT.Dispose();
                pDataRe.Dispose();
                pDataIm.Dispose();
                pDataReOrg.Dispose();
                pDataImOrg.Dispose();
            }
        }
    }
}

'    $Revision: 1.1 $

Imports System.Collections.Generic
Imports System.Text
Imports fvalgcli

Public Partial Class FIE
    ''' <summary>
    ''' splitted format の複素数データの１次元順方向フーリエ変換（インプレース処理).
    ''' </summary>
    <FvPluginExecute> _
    Public Sub fnFIE_fft_fwd_DDtoDD_in()
        Dim status As Integer = CInt(f_err.F_ERR_NONE)

        Dim iOrder As Integer = 5
        ' データ長の次数.
        Dim length As Integer = 1 << iOrder
        ' 入力データの長さ.
        Dim hFFT As FHANDLE = FHANDLE.Zero
        ' FFTハンドル.
        Dim pDataRe As DOUBLE_PTR = DOUBLE_PTR.Zero
        ' 入力データの実部(＝出力データの実部).
        Dim pDataIm As DOUBLE_PTR = DOUBLE_PTR.Zero
        ' 入力データの虚部(＝出力データの虚部).
        Dim pDataReOrg As DOUBLE_PTR = DOUBLE_PTR.Zero
        ' 結果比較用データの実部.
        Dim pDataImOrg As DOUBLE_PTR = DOUBLE_PTR.Zero
        ' 結果比較用データの虚部.
        Try
            ' データの作成.                
            pDataRe = DOUBLE_PTR.alloc(length)
            pDataIm = DOUBLE_PTR.alloc(length)
            pDataReOrg = DOUBLE_PTR.alloc(length)
            pDataImOrg = DOUBLE_PTR.alloc(length)
            For i As Integer = 0 To length - 1
                Dim image As Double = CDbl(i) * 10 + 0.1
                Dim real As Double = CDbl(i) * 10 + 0.2
                pDataRe(i) = real
                pDataIm(i) = image
                pDataReOrg(i) = real
                pDataImOrg(i) = image
            Next

            ' FFTハンドルの初期化.
            status = api.fnFIE_fft_1D_alloc(hFFT, iOrder, f_fft_direction.F_FFT_INVERSE, f_fft_normalize_type.F_FFT_DIV_BY_SQRTN, f_fft_data_type.F_1D_FFT_FCOMPLEX, f_fft_data_type.F_1D_FFT_FCOMPLEX)
            Assert.IsTrue(status = CInt(f_err.F_ERR_NONE), "エラーが発生しました。({0})", CType(status, f_err))

            ' 処理の実行.
            status = api.fnFIE_fft_fwd_DDtoDD_in(hFFT, pDataRe, pDataIm, iOrder)
            Assert.IsTrue(status = CInt(f_err.F_ERR_NONE), "エラーが発生しました。({0})", CType(status, f_err))

            ' 結果比較.
            Dim result As Boolean = False
            For i As Integer = 0 To length - 1
                If pDataRe(i) <> pDataReOrg(i) OrElse pDataIm(i) <> pDataImOrg(i) Then
                    result = True
                    Exit For
                End If
            Next
            If Not result Then
                Assert.Fail("処理前と値が変わっていません")
            End If
        Finally
            ' オブジェクトの開放.
            hFFT.Dispose()
            pDataRe.Dispose()
            pDataIm.Dispose()
            pDataReOrg.Dispose()
            pDataImOrg.Dispose()
        End Try
    End Sub
End Class