/usr/include/Cauchy/Eigen3/FFT.h is in cauchy-dev 0.9.0-0ubuntu1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 | /*
* Copyright (c) 2010 Cyrille Berger <cberger@cberger.net>
*
* This library 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 2, or (at your option) any later version of the License.
*
* This library 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
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this library; see the file COPYING.RUNTIME. If not, write to
* the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
* Boston, MA 02110-1301, USA.
*
* The file COPYING.RUNTIME.EXCEPTION contains an exception that allow
* to use this file with other open source license
*/
#ifndef _CAUCHY_FFT_H_
#define _CAUCHY_FFT_H_
#include <unsupported/Eigen/FFT>
namespace Cauchy {
Eigen::MatrixXcd fft(const Eigen::MatrixXd& m)
{
Eigen::FFT<double> fft;
Eigen::MatrixXcd r(m.rows(), m.cols());
for(int i = 0; i < m.cols(); ++i)
{
Eigen::VectorXcd rv;
fft.fwd(rv, Eigen::VectorXd(m.col(i)));
r.col(i) = rv;
}
return r;
}
Eigen::MatrixXcd ifft(const Eigen::MatrixXcd& m)
{
Eigen::FFT<double> fft;
Eigen::MatrixXcd r(m.rows(), m.cols());
for(int i = 0; i < m.cols(); ++i)
{
Eigen::VectorXcd rv;
fft.inv(rv, Eigen::VectorXcd(m.col(i)));
r.col(i) = rv;
}
return r;
}
}
#endif
|