sim_noise_source.h

Go to the documentation of this file.

 
#ifndef OAKFIELD_SIM_NOISE_SOURCE_H
#define OAKFIELD_SIM_NOISE_SOURCE_H
 
#include "field.h"
 
#include <stdbool.h>
#include <stddef.h>
#include <stdint.h>
 
#ifdef __cplusplus
extern "C" {
#endif
 
#define SIM_NOISE_SOURCE_STREAM_STIMULUS_WHITE 0x85EBCA77C2B2AE63ULL
#define SIM_NOISE_SOURCE_STREAM_STOCHASTIC 0x000000000BAD5EEDULL
#define SIM_NOISE_SOURCE_STREAM_ORNSTEIN 0xD8D12A6C9B4E41B7ULL
 
typedef struct SimNoiseSourceRng {
    uint64_t state;        
    uint64_t inc;          
    bool normal_has_spare; 
    double normal_spare;   
} SimNoiseSourceRng;
 
void sim_noise_source_seed(SimNoiseSourceRng *rng, uint64_t initstate, uint64_t initseq);
 
double sim_noise_source_uniform(SimNoiseSourceRng *rng);
 
void sim_noise_source_normal_pair(SimNoiseSourceRng *rng, double *n0, double *n1);
 
double sim_noise_source_normal(SimNoiseSourceRng *rng);
 
double sim_noise_source_gaussian_gain(double sigma, double dt, bool scale_by_dt);
 
SimResult sim_noise_source_ensure_capacity(double **real_state, SimComplexDouble **complex_state,
                                           size_t *capacity, size_t count, bool is_complex);
 
void sim_noise_source_apply_white_real(double *dst, size_t count, double mean, double gain,
                                       SimNoiseSourceRng *rng);
 
void sim_noise_source_apply_white_complex(SimComplexDouble *dst, size_t count, double mean_re,
                                          double mean_im, double gain_re, double gain_im,
                                          SimNoiseSourceRng *rng);
 
void sim_noise_source_temporal_params(double sigma, double tau, double alpha, double dt,
                                      double *out_decay, double *out_variance);
 
void sim_noise_source_apply_temporal_real(double *dst, double *noise_state, size_t count,
                                          double output_scale, double decay, double variance,
                                          SimNoiseSourceRng *rng);
 
void sim_noise_source_apply_temporal_complex(SimComplexDouble *dst, SimComplexDouble *noise_state,
                                             size_t count, double output_scale, double decay,
                                             double variance, SimNoiseSourceRng *rng);
 
void sim_noise_source_apply_mean_reverting_real(double *dst, size_t count, double mean,
                                                double decay, double variance,
                                                SimNoiseSourceRng *rng);
 
void sim_noise_source_apply_mean_reverting_complex(SimComplexDouble *dst, size_t count, double mean,
                                                   double decay, double variance,
                                                   SimNoiseSourceRng *rng);
 
#ifdef __cplusplus
}
#endif
 
#endif /* OAKFIELD_SIM_NOISE_SOURCE_H */