saf_utility_misc.h

Go to the documentation of this file.

/*
 * Copyright 2020 Leo McCormack
 *
 * Permission to use, copy, modify, and/or distribute this software for any
 * purpose with or without fee is hereby granted, provided that the above
 * copyright notice and this permission notice appear in all copies.
 *
 * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES WITH
 * REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY
 * AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY SPECIAL, DIRECT,
 * INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM
 * LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR
 * OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR
 * PERFORMANCE OF THIS SOFTWARE.
 */

 
#ifndef SAF_MISC_H_INCLUDED
#define SAF_MISC_H_INCLUDED
 
#ifdef __cplusplus
extern "C" {
#endif /* __cplusplus */
 
#include "saf_utility_complex.h"
 
#if defined(__STDC_VERSION__) && (__STDC_VERSION__ >= 201112L) && !defined(__STDC_NO_ATOMICS__)
  #include <stdatomic.h>
  typedef _Atomic int   _Atomic_INT32;
  typedef _Atomic float _Atomic_FLOAT32;
#else
  typedef int   _Atomic_INT32;
  typedef float _Atomic_FLOAT32;
#endif
 
/* Cross-platform sleep macro (slightly modified), originally taken from:
 * https://cboard.cprogramming.com/c-programming/170381-cross-platform-wait-sleep.html */
#ifdef _WIN32
    /* For Windows (32- and 64-bit) */
# include <windows.h>
# define SAF_SLEEP(msecs) Sleep(msecs)
#elif defined(__unix) || defined(__APPLE__)
    /* For linux, OSX, and other unixes */
# ifndef _POSIX_C_SOURCE
#  define _POSIX_C_SOURCE 199309L /* or greater */
# endif
# include <time.h>
# define SAF_SLEEP(msecs) do {     \
   struct timespec ts;             \
   ts.tv_sec = msecs/1000;         \
   ts.tv_nsec = msecs%1000*1000;   \
   nanosleep(&ts, NULL);           \
} while (0)
#else
# error "Unknown system"
#endif

#define SAF_CASSERT(predicate, file) _impl_CASSERT_LINE(predicate,__LINE__,file)
#define _impl_PASTE(a,b) a##b
#define _impl_CASSERT_LINE(predicate, line, file) \
 typedef char _impl_PASTE(assertion_failed_##file##_,line)[2*!!(predicate)-1];

void convert_0_360To_m180_180(float* dirs_deg,
                              int nDirs);

int nextpow2(int numsamp);

void lagrangeWeights(int N,
                     float* x,
                     int len_x,
                     float* weights);

void findERBpartitions(/* Input Arguments */
                       float* centerFreq,
                       int nBands,
                       float maxFreqLim,
                       /* Output Arguments */
                       int** erb_idx,
                       float** erb_freqs,
                       int* nERBBands);

void randperm(/* Input Arguments */
              int len,
              /* Output Arguments */
              int* randperm_inds);

long double factorial(int n);

float matlab_fmodf(float x,  float y);

void cxcorr(float* a,
            float* b,
            float* x_ab,
            size_t la,
            size_t lb);

void rand_m1_1(float* vector,
               int length);

void rand_cmplx_m1_1(float_complex* vector,
                     int length);

void rand_0_1(float* vector,
              int length);

void convd(double* x,
           double* h,
           int len_x,
           int len_h,
           double* y);

void convz(double_complex* x,
           double_complex* h,
           int len_x,
           int len_h,
           double_complex* y);

void polyd_v(double* x,
             double* poly,
             int len_x);

void polyz_v(double_complex* x,
             double_complex* poly,
             int len_x);

void polyd_m(double* X,
             double_complex* poly,
             int size_x);

float sumf(float* values,
           int nValues);

int anyLessThanf(float* values,
                 int nValues,
                 float threshold);

void unique_i(int* input,
              int nInputs,
              int** uniqueVals,
              int** uniqueInds,
              int* nUnique);

void findCombinations(int* arrValues,
                      int nValues,
                      int nElements,
                      int** comb,
                      int* nComb);

void gexpm(float* D,
           int sizeD,
           int m1,
           float* Y);
 
 
#ifdef __cplusplus
}/* extern "C" */
#endif /* __cplusplus */
 
#endif /* SAF_MISC_H_INCLUDED */
 /* doxygen addtogroup Utilities */