vmath.h

Go to the documentation of this file.

 
#ifndef SPLIB_VMATH_H
#define SPLIB_VMATH_H
 
#include <sp/vector.h>
 
#ifdef __cplusplus
extern "C" {
#endif
 
#if defined(MACOS)
#pragma import on
#endif
 
  /*----@addtogroup vmathGroup----*/
 
extern void svabs(spSVector x);
extern void lvabs(spLVector x);
extern void fvabs(spFVector x);
extern void dvabs(spDVector x);
extern spSVector xsvabs(spSVector x);
extern spLVector xlvabs(spLVector x);
extern spFVector xfvabs(spFVector x);
extern spDVector xdvabs(spDVector x);
 
extern void svsqrt(spSVector x);
extern void lvsqrt(spLVector x);
extern void fvsqrt(spFVector x);
extern void dvsqrt(spDVector x);
extern spSVector xsvsqrt(spSVector x);
extern spLVector xlvsqrt(spLVector x);
extern spFVector xfvsqrt(spFVector x);
extern spDVector xdvsqrt(spDVector x);
 
extern void svsquare(spSVector x);
extern void lvsquare(spLVector x);
extern void fvsquare(spFVector x);
extern void dvsquare(spDVector x);
extern spSVector xsvsquare(spSVector x);
extern spLVector xlvsquare(spLVector x);
extern spFVector xfvsquare(spFVector x);
extern spDVector xdvsquare(spDVector x);
 
extern void svsign(spSVector x);
extern void lvsign(spLVector x);
extern void fvsign(spFVector x);
extern void dvsign(spDVector x);
extern spSVector xsvsign(spSVector x);
extern spLVector xlvsign(spLVector x);
extern spFVector xfvsign(spFVector x);
extern spDVector xdvsign(spDVector x);
 
extern void svcumsum(spSVector x);
extern void lvcumsum(spLVector x);
extern void fvcumsum(spFVector x);
extern void dvcumsum(spDVector x);
extern spSVector xsvcumsum(spSVector x);
extern spLVector xlvcumsum(spLVector x);
extern spFVector xfvcumsum(spFVector x);
extern spDVector xdvcumsum(spDVector x);
 
extern void svcumprod(spSVector x);
extern void lvcumprod(spLVector x);
extern void fvcumprod(spFVector x);
extern void dvcumprod(spDVector x);
extern spSVector xsvcumprod(spSVector x);
extern spLVector xlvcumprod(spLVector x);
extern spFVector xfvcumprod(spFVector x);
extern spDVector xdvcumprod(spDVector x);
 
extern void svcodiff(spSVector x, double coef);
extern void lvcodiff(spLVector x, double coef);
extern void fvcodiff(spFVector x, double coef);
extern void dvcodiff(spDVector x, double coef);
extern spSVector xsvcodiff(spSVector x, double coef);
extern spLVector xlvcodiff(spLVector x, double coef);
extern spFVector xfvcodiff(spFVector x, double coef);
extern spDVector xdvcodiff(spDVector x, double coef);
 
#define svdiff(x) svcodiff(x, 1.0);
#define lvdiff(x) lvcodiff(x, 1.0);
#define fvdiff(x) fvcodiff(x, 1.0);
#define dvdiff(x) dvcodiff(x, 1.0);
#define xsvdiff(x) xsvcodiff(x, 1.0);
#define xlvdiff(x) xlvcodiff(x, 1.0);
#define xfvdiff(x) xfvcodiff(x, 1.0);
#define xdvdiff(x) xdvcodiff(x, 1.0);
    
extern long svsum(spSVector x);
extern long lvsum(spLVector x);
extern float fvsum(spFVector x);
extern double dvsum(spDVector x);
 
extern long svisum(spSVector x);
extern long lvisum(spLVector x);
extern float fvisum(spFVector x);
extern double dvisum(spDVector x);
 
extern long svsqsum(spSVector x);
extern long lvsqsum(spLVector x);
extern float fvsqsum(spFVector x);
extern double dvsqsum(spDVector x);
 
extern long svabssum(spSVector x);
extern long lvabssum(spLVector x);
extern float fvabssum(spFVector x);
extern double dvabssum(spDVector x);
 
extern long svprod(spSVector x);
extern long lvprod(spLVector x);
extern float fvprod(spFVector x);
extern double dvprod(spDVector x);
 
extern long svriprod(spSVector x, long *oiprod);
extern long lvriprod(spLVector x, long *oiprod);
extern float fvriprod(spFVector x, float *oiprod);
extern double dvriprod(spDVector x, double *oiprod);
 
extern double svnorm(spSVector x, long p);
extern double lvnorm(spLVector x, long p);
extern float fvnorm(spFVector x, long p);
extern double dvnorm(spDVector x, long p);
 
extern long svdot(spSVector x, spSVector y);
extern long lvdot(spLVector x, spLVector y);
extern float fvdot(spFVector x, spFVector y);
extern double dvdot(spDVector x, spDVector y);
 
extern void svcplxdot(spSVector x, spSVector y, long *re, long *im);
extern void lvcplxdot(spLVector x, spLVector y, long *re, long *im);
extern void fvcplxdot(spFVector x, spFVector y, float *re, float *im);
extern void dvcplxdot(spDVector x, spDVector y, double *re, double *im);
 
extern double svvar(spSVector x, double *mean);
extern double lvvar(spLVector x, double *mean);
extern float fvvar(spFVector x, float *mean);
extern double dvvar(spDVector x, double *mean);
 
extern double svstd(spSVector x, double *mean);
extern double lvstd(spLVector x, double *mean);
extern float fvstd(spFVector x, float *mean);
extern double dvstd(spDVector x, double *mean);
 
#define svmean(x) ((double)svsum(x) / (double)x->length)
#define lvmean(x) ((double)lvsum(x) / (double)x->length)
#define fvmean(x) (fvsum(x) / (float)x->length)
#define dvmean(x) (dvsum(x) / (double)x->length)
 
extern short svmin(spSVector x, long *index);
extern long lvmin(spLVector x, long *index);
extern float fvmin(spFVector x, long *index);
extern double dvmin(spDVector x, long *index);
 
extern short svmax(spSVector x, long *index);
extern long lvmax(spLVector x, long *index);
extern float fvmax(spFVector x, long *index);
extern double dvmax(spDVector x, long *index);
 
extern void svscmin(spSVector x, short a);
extern void lvscmin(spLVector x, long a);
extern void fvscmin(spFVector x, float a);
extern void dvscmin(spDVector x, double a);
 
extern void svscmax(spSVector x, short a);
extern void lvscmax(spLVector x, long a);
extern void fvscmax(spFVector x, float a);
extern void dvscmax(spDVector x, double a);
 
extern void svelmin(spSVector x, spSVector y);
extern void lvelmin(spLVector x, spLVector y);
extern void fvelmin(spFVector x, spFVector y);
extern void dvelmin(spDVector x, spDVector y);
 
extern void svelmax(spSVector x, spSVector y);
extern void lvelmax(spLVector x, spLVector y);
extern void fvelmax(spFVector x, spFVector y);
extern void dvelmax(spDVector x, spDVector y);
 
#define SP_NMINMAX_OPTION_PEAK_ONLY 1
#define SP_NMINMAX_OPTION_ZERO_EDGES (1<<1)
extern spLVector xsvnmin(spSVector x, long n, long margin, unsigned long options);
extern spLVector xlvnmin(spLVector x, long n, long margin, unsigned long options);
extern spLVector xfvnmin(spFVector x, long n, long margin, unsigned long options);
extern spLVector xdvnmin(spDVector x, long n, long margin, unsigned long options);
extern spLVector xsvnmax(spSVector x, long n, long margin, unsigned long options);
extern spLVector xlvnmax(spLVector x, long n, long margin, unsigned long options);
extern spLVector xfvnmax(spFVector x, long n, long margin, unsigned long options);
extern long dvnmax(spLVector idx, spDVector x, long margin, unsigned long options);
extern spLVector xdvnmax(spDVector x, long n, long margin, unsigned long options);
    
extern void svsort(spSVector x);
extern void lvsort(spLVector x);
extern void fvsort(spFVector x);
extern void dvsort(spDVector x);
extern spSVector xsvsort(spSVector x);
extern spLVector xlvsort(spLVector x);
extern spFVector xfvsort(spFVector x);
extern spDVector xdvsort(spDVector x);
 
extern void spheapsort(void *keysptr, void *valuesptr, long length,
               int (*compare_func)(void *keysptr, long i, long j),
               void (*swap_func)(void *keysptr, void *valuesptr, long i, long j));
    
extern spLVector xsvsortidx(spSVector x);
extern spLVector xlvsortidx(spLVector x);
extern spLVector xfvsortidx(spFVector x);
extern spLVector xdvsortidx(spDVector x);
 
extern double svmedian(spSVector x);
extern double lvmedian(spLVector x);
extern float fvmedian(spFVector x);
extern double dvmedian(spDVector x);
extern spSVector xsvmedian(spSVector x, double *median);
extern spLVector xlvmedian(spLVector x, double *median);
extern spFVector xfvmedian(spFVector x, float *median);
extern spDVector xdvmedian(spDVector x, double *median);
 
extern long svzerocross(spSVector x, int slope_sign);
extern long lvzerocross(spLVector x, int slope_sign);
extern long fvzerocross(spFVector x, int slope_sign);
extern long dvzerocross(spDVector x, int slope_sign);
 
/*-- the following functions are supported for float and double only --*/
 
extern void fvangle(spFVector x);
extern spFVector xfvangle(spFVector x);
extern void dvangle(spDVector x);
extern spDVector xdvangle(spDVector x);
 
extern void fvcos(spFVector x);
extern void dvcos(spDVector x);
extern spFVector xfvcos(spFVector x);
extern spDVector xdvcos(spDVector x);
 
extern void fvsin(spFVector x);
extern void dvsin(spDVector x);
extern spFVector xfvsin(spFVector x);
extern spDVector xdvsin(spDVector x);
 
extern void fvtan(spFVector x);
extern void dvtan(spDVector x);
extern spFVector xfvtan(spFVector x);
extern spDVector xdvtan(spDVector x);
 
extern void fvcosh(spFVector x);
extern void dvcosh(spDVector x);
extern spFVector xfvcosh(spFVector x);
extern spDVector xdvcosh(spDVector x);
 
extern void fvsinh(spFVector x);
extern void dvsinh(spDVector x);
extern spFVector xfvsinh(spFVector x);
extern spDVector xdvsinh(spDVector x);
 
extern void fvtanh(spFVector x);
extern void dvtanh(spDVector x);
extern spFVector xfvtanh(spFVector x);
extern spDVector xdvtanh(spDVector x);
 
extern void fvacos(spFVector x);
extern void dvacos(spDVector x);
extern spFVector xfvacos(spFVector x);
extern spDVector xdvacos(spDVector x);
 
extern void fvasin(spFVector x);
extern void dvasin(spDVector x);
extern spFVector xfvasin(spFVector x);
extern spDVector xdvasin(spDVector x);
 
extern void fvatan(spFVector x);
extern void dvatan(spDVector x);
extern spFVector xfvatan(spFVector x);
extern spDVector xdvatan(spDVector x);
 
extern void fvatan2(spFVector y, spFVector x, spBool reverse);
extern void dvatan2(spDVector y, spDVector x, spBool reverse);
extern spFVector xfvatan2(spFVector y, spFVector x);
extern spDVector xdvatan2(spDVector y, spDVector x);
 
extern void fvexp(spFVector x);
extern void dvexp(spDVector x);
extern spFVector xfvexp(spFVector x);
extern spDVector xdvexp(spDVector x);
 
extern void fvlog(spFVector x);
extern void dvlog(spDVector x);
extern spFVector xfvlog(spFVector x);
extern spDVector xdvlog(spDVector x);
 
extern void fvlog10(spFVector x);
extern void dvlog10(spDVector x);
extern spFVector xfvlog10(spFVector x);
extern spDVector xdvlog10(spDVector x);    
 
extern void fvdecibela(spFVector x);
extern void dvdecibela(spDVector x);
extern spFVector xfvdecibela(spFVector x);
extern spDVector xdvdecibela(spDVector x);
    
extern void fvdecibelp(spFVector x);
extern void dvdecibelp(spDVector x);
extern spFVector xfvdecibelp(spFVector x);
extern spDVector xdvdecibelp(spDVector x);
 
extern void fvundecibelp(spFVector x);
extern void dvundecibelp(spDVector x);
extern spFVector xfvundecibelp(spFVector x);
extern spDVector xdvundecibelp(spDVector x);
 
extern void fvdecibel(spFVector x);
extern void dvdecibel(spDVector x);
extern spFVector xfvdecibel(spFVector x);
extern spDVector xdvdecibel(spDVector x);
 
extern void fvundecibel(spFVector x);
extern void dvundecibel(spDVector x);
extern spFVector xfvundecibel(spFVector x);
extern spDVector xdvundecibel(spDVector x);
    
extern void fvrandun(spFVector x);
extern void dvrandun(spDVector x);
 
extern spFVector xfvrandunul(spPlugin *plugin, long length, spBool unlock_flag);
extern spDVector xdvrandunul(spPlugin *plugin, long length, spBool unlock_flag);
 
extern spFVector xfvrandun(long length);
extern spDVector xdvrandun(long length);
 
extern void fvgauss(spFVector x, float mu, float sigma);
extern void dvgauss(spDVector x, double mu, double sigma);
 
extern spFVector xfvgaussul(spPlugin *plugin, long length, float mu, float sigma, spBool unlock_flag);
extern spDVector xdvgaussul(spPlugin *plugin, long length, double mu, double sigma, spBool unlock_flag);
 
extern spFVector xfvgauss(long length, float mu, float sigma);
extern spDVector xdvgauss(long length, double mu, double sigma);
 
extern void fvrandn(spFVector x);
extern void dvrandn(spDVector x);
 
extern spFVector xfvrandnul(spPlugin *plugin, long length, spBool unlock_flag);
extern spDVector xdvrandnul(spPlugin *plugin, long length, spBool unlock_flag);
 
extern spFVector xfvrandn(long length);
extern spDVector xdvrandn(long length);
 
  /*----@addtogroup vmathGroup----*/
 
extern double dvscpolyval(spDVector p, double x);
extern void dvcplxpolyval(spDVector p, double xre, double xim, double *oyre, double *oyim);
extern void dvpolyval(spDVector iox, spDVector p);
extern spDVector xdvpolyval(spDVector x, spDVector p);
 
#if defined(MACOS)
#pragma import off
#endif
 
#ifdef __cplusplus
}  /* Close scope of 'extern "C"' declaration */
#endif
 
#endif /* SPLIB_VMATH_H */