d12r8.txt

Visual Basic 常用数学算法集书中收录了所有代码

Private Sub Command1_Click()
    'PROGRAM D12R8
    'Driver for routine FOURN
    NDIM = 3
    NDAT = 1024
    Dim NN(3), Data(1024)
    For I = 1 To NDIM
        NN(I) = 2 * (2 ^ I)
    Next I
    For I = 1 To NN(3)
        For J = 1 To NN(2)
            For K = 1 To NN(1)
              L = K + (J - 1) * NN(1) + (I - 1) * NN(2) * NN(1)
              LL = 2 * L - 1
              Data(LL) = LL
              Data(LL + 1) = LL + 1
            Next K
        Next J
    Next I
    ISIGN = 1
    Call FOURN(Data(), NN(), NDIM, ISIGN)
    ISIGN = -1
    Print
    Print Tab(5); "FORuble 3-dimensional Transform"
    Print
    Print Tab(5); "FORuble Trabsf.    Original Data      Ratio"
    Print
    Print Tab(5); " Real        Imag.       Real     Imag.    Real     Imag."
    Call FOURN(Data(), NN(), NDIM, ISIGN)
    For I = 1 To 4
        J = 2 * I
        K = 2 * J
        L = K + (J - 1) * NN(1) + (I - 1) * NN(2) * NN(1)
        LL = 2 * L - 1
        Print Tab(5); Format$(Data(LL), "######.#0");
        Print Tab(17); Format$(Data(LL + 1), "######.#0");
        Print Tab(29); Format$(LL, "###.#0");
        Print Tab(39); Format$(LL + 1, "###.##");
        Print Tab(47); Format$(Data(LL) / LL, "###.#0");
        Print Tab(57); Format$(Data(LL + 1) / (LL + 1), "###.#0")
    Next I
    Print
    Print Tab(5); "The product of transform lengths is: ";
    Print NN(1) * NN(2) * NN(3)
End Sub