matrix_fft.h

Go to the documentation of this file.

/*
 * This work is licensed under a Creative Commons 
 * Attribution-Noncommercial-Share Alike 3.0 United States License.
 * 
 *    http://creativecommons.org/licenses/by-nc-sa/3.0/us/
 * 
 * You are free:
 * 
 *    to Share - to copy, distribute, display, and perform the work
 *    to Remix - to make derivative works
 * 
 * Under the following conditions:
 * 
 *    Attribution. You must attribute the work in the manner specified by the
 *    author or licensor (but not in any way that suggests that they endorse you
 *    or your use of the work).
 * 
 *    Noncommercial. You may not use this work for commercial purposes.
 * 
 *    Share Alike. If you alter, transform, or build upon this work, you may
 *    distribute the resulting work only under the same or similar license to
 *    this one.
 * 
 * For any reuse or distribution, you must make clear to others the license
 * terms of this work. The best way to do this is by including this header.
 * 
 * Any of the above conditions can be waived if you get permission from the
 * copyright holder.
 * 
 * Apart from the remix rights granted under this license, nothing in this
 * license impairs or restricts the author's moral rights.
 */


#ifndef MATRIX_FFT_H
#define MATRIX_FFT_H


#include <jwsc/config.h>

#include <stdlib.h>

#include "jwsc/matrix/matrix.h"


#ifdef __cplusplus
namespace jwsc {
extern "C" {
#endif


#ifdef JWSC_HAVE_FFTW3F
void real_to_complex_matrix_dft_f
(
    Matrix_cf**     m_out,
    const Matrix_f* m_in
);
#endif

#ifdef JWSC_HAVE_FFTW3
void real_to_complex_matrix_dft_d
(
    Matrix_cd**     m_out,
    const Matrix_d* m_in
);
#endif

#ifdef JWSC_HAVE_FFTW3F
void complex_to_real_matrix_dft_f
(
    Matrix_f**       m_out,
    const Matrix_cf* m_in,
    uint8_t          ncols_even
);
#endif

#ifdef JWSC_HAVE_FFTW3
void complex_to_real_matrix_dft_d
(
    Matrix_d**       m_out,
    const Matrix_cd* m_in,
    uint8_t          ncols_even
);
#endif

void expand_symmetric_matrix_dft_f
(
    Matrix_cf**      m_out, 
    const Matrix_cf* m_in, 
    uint8_t          ncols_even
);

void expand_symmetric_matrix_dft_d
(
    Matrix_cd**      m_out, 
    const Matrix_cd* m_in, 
    uint8_t          ncols_even
);

Error* shrink_symmetric_matrix_dft_f(Matrix_cf** m_out, const Matrix_cf* m_in);

Error* shrink_symmetric_matrix_dft_d(Matrix_cd** m_out, const Matrix_cd* m_in);

void scale_matrix_dft_f(Matrix_f** m_out, const Matrix_f* m_in);

void scale_matrix_dft_d(Matrix_d** m_out, const Matrix_d* m_in);

void shift_matrix_f
(
    Matrix_f**       m_out,
    const Matrix_f*  m_in
);

void shift_matrix_d
(
    Matrix_d**       m_out,
    const Matrix_d*  m_in
);

void shift_matrix_cf
(
    Matrix_cf**       m_out,
    const Matrix_cf*  m_in
);

void shift_matrix_cd
(
    Matrix_cd**       m_out,
    const Matrix_cd*  m_in
);

void unshift_matrix_f
(
    Matrix_f**       m_out,
    const Matrix_f*  m_in
);

void unshift_matrix_d
(
    Matrix_d**       m_out,
    const Matrix_d*  m_in
);

void unshift_matrix_cf
(
    Matrix_cf**       m_out,
    const Matrix_cf*  m_in
);

void unshift_matrix_cd
(
    Matrix_cd**       m_out,
    const Matrix_cd*  m_in
);

#ifdef __cplusplus
}
}
#endif


#endif