Homec4science
Diffusion µSpectre prototype implementation (add_eigenstrain) Browse

muSpectre_proto/tests/test_fftw_c2c.cc
54cf92443a69add_eigenstrain

Last Change
Raw File

test_fftw_c2c.cc
View Options

/**
* file test_fftw_c2c.cc
*
* @author Till Junge <till.junge@epfl.ch>
*
* @date 11 May 2017
*
* @brief test fftw redundant class
*
* @section LICENCE
*
* Copyright (C) 2017 Till Junge
*
* µSpectre 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 3, or (at
* your option) any later version.
*
* µSpectre 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 GNU Emacs; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*/
#include <boost/test/unit_test.hpp>
#include <iostream>
#include "complex"
#include "system/fftw_engine_c2c.hh"
#include "tests.hh"
namespace muSpectre {
template <Dim_t DimS, Dim_t DimM>
struct fftw_c2c_fixture: public FFTW_EngineC2C<DimS,DimM> {
fftw_c2c_fixture()
:FFTW_EngineC2C<DimS,DimM>(std_nb_pix){}
const static typename FFTW_EngineC2C<DimS,DimM>::Index std_nb_pix;
};
template<> const Index_t<1> fftw_c2c_fixture<1,1>::std_nb_pix{3};
template<> const Index_t<2> fftw_c2c_fixture<2,2>::std_nb_pix{5,3};
template<> const Index_t<3> fftw_c2c_fixture<3,3>::std_nb_pix{7,5,3};
using fix2 = fftw_c2c_fixture<2, 2>;
BOOST_FIXTURE_TEST_CASE(numpy_comparison, fix2) {
const size_t dim = 2;
init(FFT_PlanFlags::estimate);
Real in_arr[tot_nb_pixels*dim*dim];
for (size_t i = 0; i < tot_nb_pixels; ++i) {
for (size_t d = 0; d < dim*dim; ++d) {
in_arr[i*dim*dim+d] = i;
}
}
Real Gdata[dim*dim*dim*dim*tot_nb_pixels] = {0};
Real * out_arr = convolve(Gdata, in_arr);
std::complex<Real>* cptr = reinterpret_cast<std::complex<Real>*>(out_arr);
}
} // muSpectre
c4science · Help