// This file is part of Eigen, a lightweight C++ template library
// for linear algebra.
//
// Copyright (C) 2014 Jianwei Cui <thucjw@gmail.com>
//
// This Source Code Form is subject to the terms of the Mozilla
// Public License v. 2.0. If a copy of the MPL was not distributed
// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
#include "main.h"
#include <complex>
#include <cmath>
#include <Eigen/CXX11/Tensor>
using Eigen::Tensor;
template <
int DataLayout>
static void test_1D_fft_ifft_invariant(
int sequence_length) {
Tensor<
double ,
1 , DataLayout> tensor(sequence_length);
tensor.setRandom();
array<
int ,
1 > fft;
fft[
0 ] =
0 ;
Tensor<std::complex<
double >,
1 , DataLayout> tensor_after_fft;
Tensor<std::complex<
double >,
1 , DataLayout> tensor_after_fft_ifft;
tensor_after_fft = tensor.
template fft<Eigen::BothParts, Eigen::FFT_FORWARD>(fft);
tensor_after_fft_ifft = tensor_after_fft.
template fft<Eigen::BothParts, Eigen::FFT_RE
VERSE>(fft);0 ), sequence_length);0 ), sequence_length);for (int i = 0 ; i < sequence_length; ++i) {static_cast <float >(tensor(i)), static_cast <float >(std::real(tensor_after_fft_ifft(i))));template <int DataLayout>static void test_2D_fft_ifft_invariant(int dim0, int dim1) {double , 2 , DataLayout> tensor(dim0, dim1);int , 2 > fft;0 ] = 0 ;1 ] = 1 ;double >, 2 , DataLayout> tensor_after_fft;double >, 2 , DataLayout> tensor_after_fft_ifft;template fft<Eigen::BothParts, Eigen::FFT_FORWARD>(fft);template fft<Eigen::BothParts, Eigen::FFT_REVERSE>(fft);0 ), dim0);1 ), dim1);0 ), dim0);1 ), dim1);for (int i = 0 ; i < dim0; ++i) {for (int j = 0 ; j < dim1; ++j) {//std::cout << "[" << i << "][" << j << "]" << " Original data: " << tensor(i,j) << " Transformed data:" << tensor_after_fft_ifft(i,j) << std::endl; static_cast <float >(tensor(i,j)), static_cast <float >(std::real(tensor_after_fft_ifft(i,j))));template <int DataLayout>static void test_3D_fft_ifft_invariant(int dim0, int dim1, int dim2) {double , 3 , DataLayout> tensor(dim0, dim1, dim2);int , 3 > fft;0 ] = 0 ;1 ] = 1 ;2 ] = 2 ;double >, 3 , DataLayout> tensor_after_fft;double >, 3 , DataLayout> tensor_after_fft_ifft;template fft<Eigen::BothParts, Eigen::FFT_FORWARD>(fft);template fft<Eigen::BothParts, Eigen::FFT_REVERSE>(fft);0 ), dim0);1 ), dim1);2 ), dim2);0 ), dim0);1 ), dim1);2 ), dim2);for (int i = 0 ; i < dim0; ++i) {for (int j = 0 ; j < dim1; ++j) {for (int k = 0 ; k < dim2; ++k) {static_cast <float >(tensor(i,j,k)), static_cast <float >(std::real(tensor_after_fft_ifft(i,j,k))));template <int DataLayout>static void test_sub_fft_ifft_invariant(int dim0, int dim1, int dim2, int dim3) {double , 4 , DataLayout> tensor(dim0, dim1, dim2, dim3);int , 2 > fft;0 ] = 2 ;1 ] = 0 ;double >, 4 , DataLayout> tensor_after_fft;double , 4 , DataLayout> tensor_after_fft_ifft;template fft<Eigen::BothParts, Eigen::FFT_FORWARD>(fft);template fft<Eigen::RealPart, Eigen::FFT_REVERSE>(fft);0 ), dim0);1 ), dim1);2 ), dim2);3 ), dim3);0 ), dim0);1 ), dim1);2 ), dim2);3 ), dim3);for (int i = 0 ; i < dim0; ++i) {for (int j = 0 ; j < dim1; ++j) {for (int k = 0 ; k < dim2; ++k) {for (int l = 0 ; l < dim3; ++l) {static_cast <float >(tensor(i,j,k,l)), static_cast <float >(tensor_after_fft_ifft(i,j,k,l)));4 ));16 ));32 ));1024 *1024 ));4 ,4 ));8 ,16 ));16 ,32 ));1024 ,1024 ));4 ,4 ,4 ));8 ,16 ,32 ));16 ,4 ,8 ));256 ,256 ,256 ));4 ,4 ,4 ,4 ));8 ,16 ,32 ,64 ));16 ,4 ,8 ,12 ));64 ,64 ,64 ,64 )); Messung V0.5 in Prozent C=91 H=92 G=91
¤ Dauer der Verarbeitung: 0.14 Sekunden
(vorverarbeitet am 2026-06-06)
¤ *© Formatika GbR, Deutschland
