math/fourier.h

Go to the documentation of this file.

 
#ifndef FOURIER_H
#define FOURIER_H
 
#include <complex.h>
#include <stdbool.h>
#include <stddef.h>
#include <stdint.h>
 
#include "oakfield/field.h"
 
#ifndef SIM_FOURIER_ASSUME_VALID
#define SIM_FOURIER_ASSUME_VALID 0
#endif
 
#ifdef __cplusplus
extern "C" {
#endif
 
#ifndef SIM_HAVE_SIMCOMPLEXDOUBLE
typedef struct SimComplexDouble {
    double re; 
    double im; 
} SimComplexDouble;
#endif
 
#ifndef SIM_FOURIER_MINIBLEP_SPAN
#define SIM_FOURIER_MINIBLEP_SPAN 4.0
#endif
 
#ifndef SIM_FOURIER_MINIBLEP_MAX_WINDOW
#define SIM_FOURIER_MINIBLEP_MAX_WINDOW 0.5
#endif
 
#ifndef SIM_FOURIER_INTEGRATOR_LEAK
#define SIM_FOURIER_INTEGRATOR_LEAK 1e-9
#endif
 
/*==============================================================================
 *                           BLIT (Dirichlet kernel)
 *============================================================================*/
 
double sim_fourier_dirichlet(double phase_radians, int harmonic_count);
 
SimComplexDouble sim_fourier_dirichlet_complex(SimComplexDouble phase_radians, int harmonic_count);
 
double sim_fourier_blit(double phase_radians, int harmonic_count);
 
SimComplexDouble sim_fourier_blit_complex(SimComplexDouble phase_radians, int harmonic_count);
 
double sim_fourier_saw_blit(double phase_radians, double phase_increment_radians,
                            int harmonic_count, double *state);
 
SimComplexDouble sim_fourier_saw_blit_complex(SimComplexDouble phase_radians,
                                              double phase_increment_radians, int harmonic_count,
                                              SimComplexDouble *state);
 
double sim_fourier_square_blit(double phase_radians, double phase_increment_radians,
                               int harmonic_count, double duty, double *state);
 
SimComplexDouble sim_fourier_square_blit_complex(SimComplexDouble phase_radians,
                                                 double phase_increment_radians, int harmonic_count,
                                                 double duty, SimComplexDouble *state);
 
double sim_fourier_triangle_blit(double phase_radians, double phase_increment_radians,
                                 int harmonic_count, double *velocity_state,
                                 double *position_state);
 
SimComplexDouble sim_fourier_triangle_blit_complex(SimComplexDouble phase_radians,
                                                   double phase_increment_radians,
                                                   int harmonic_count,
                                                   SimComplexDouble *velocity_state,
                                                   SimComplexDouble *position_state);
 
/*==============================================================================
 *                           PolyBLEP / PolyBLAMP
 *============================================================================*/
 
double sim_fourier_polyblep(double phase, double dphase);
 
SimComplexDouble sim_fourier_polyblep_complex(double phase, double dphase);
 
double sim_fourier_polyblamp(double phase, double dphase);
 
SimComplexDouble sim_fourier_polyblamp_complex(double phase, double dphase);
 
/*==============================================================================
 *                           miniBLEP / miniBLAMP
 *============================================================================*/
 
double sim_fourier_miniblep(double phase, double dphase);
 
SimComplexDouble sim_fourier_miniblep_complex(double phase, double dphase);
 
double sim_fourier_miniblamp(double phase, double dphase);
 
SimComplexDouble sim_fourier_miniblamp_complex(double phase, double dphase);
 
#ifdef __cplusplus
}
#endif
 
#endif /* FOURIER_H */