basic_algebra.h

Go to the documentation of this file.

#ifndef BASICALGEBRAH

#include "defines.h"
#include "boolean.h"

#include <stdio.h>

#if !defined(_LAPACK) && !defined(_GSL)
  #error "CuikSuite requires LAPACK or GSL (LAPACK preferred)"
#endif

#ifdef _CBLAS
  #ifdef _GSL_CBLAS
    /* Use the CBLAS included in the GSL */
    #include <gsl/gsl_cblas.h>
  #else
    /* Use the system CBLAS */
    #include <cblas.h>
  #endif

  #define CBLAS_INLINE inline
#else

  #define CBLAS_INLINE
#endif 

#ifdef _GSL

  #define ROW_MAJOR 1
#endif

#ifdef _LAPACK

  #define ROW_MAJOR 0
#endif

#define RC2INDEX(i,j,nr,nc) (ROW_MAJOR?(i)*(nc)+(j):(j)*(nr)+(i))



/*******************************************************************************/
/* The functions below have a naif implementation and a more efficient one     */
/* based on blas                                                               */

double GeneralDotProduct(unsigned int s,double *v1,double *v2);

void ScaleVector(double f,unsigned int s,double *v);

void ScaleVector2(double f,unsigned int s,double *v,double *vout);

void AccumulateVector(unsigned int s,double *v1,double *v2);

void SumVector(unsigned int s,double *v1,double *v2,double *v);

void SumVectorScale(unsigned int s,double *v1,double w,double *v2,double *v);

void CosVector(unsigned int s,double *v,double *co);

void SinVector(unsigned int s,double *v,double *si);

unsigned int MaxVectorElement(unsigned int m,double *v);

double MaxVector(unsigned int m,double *v);

unsigned int MinVectorElement(unsigned int m,double *v);

double MinVector(unsigned int m,double *v);

void Vector2Range(double l,double u,double mode,
                  unsigned int m,double *v);

void SubtractVector(unsigned int s,double *v1,double *v2);

void DifferenceVector(unsigned int s,double *v1,double *v2,double *v);

void DifferenceVectorSubset(unsigned int s,
                            boolean *subset,
                            double *v1,double *v2,double *v);
void DifferenceVectorTopology(unsigned int s,unsigned int *tp,
                              double *v1,double *v2,double *v);

void DifferenceVectorTopologySubset(unsigned int s,unsigned int *tp,
                                    boolean *subset,
                                    double *v1,double *v2,double *v);

double Norm(unsigned int s,double *v);

double NormWithStride(unsigned int s,unsigned int st,double *v);

double SquaredDistance(unsigned int s,double *v1,double *v2);

double SquaredDistanceSubset(unsigned int s,boolean *subset,
                             double *v1,double *v2);

double Distance(unsigned int s,double *v1,double *v2);

double SquaredDistanceTopology(unsigned int s,unsigned int *tp,double *v1,double *v2);

double SquaredDistanceTopologyMin(double t2,unsigned int s,unsigned int *tp,
                                  double *v1,double *v2);

double DistanceTopology(unsigned int s,unsigned int *tp,
                        double *v1,double *v2);

double DistanceTopologySubset(unsigned int s,unsigned int *tp,
                              boolean *subset,
                              double *v1,double *v2);

double DistanceTopologyMin(double t,unsigned int s,unsigned int *tp,
                           double *v1,double *v2);

boolean CrossTopologyBorder(unsigned int s,unsigned int *tp,double *v1,double *v2);

void Normalize(unsigned int s,double *v);

double Mean(unsigned int s,double *v);

double StdDev(unsigned int s,double m,double *v);

void ArrayPi2Pi(unsigned int n,unsigned int *t,double *a);

void GetRow(unsigned int k,unsigned int r,unsigned int c,double *m,double *v);

void GetColumn(unsigned int k,unsigned int r,unsigned int c,double *m,double *v);

void SetRow(double *v,unsigned int k,unsigned int r,unsigned int c,double *m);

void SetColumn(double *v,unsigned int k,unsigned int r,unsigned int c,double *m);

double RowSquaredNorm(unsigned int k,unsigned int r,unsigned int c,double *m);

double ColumnSquaredNorm(unsigned int k,unsigned int r,unsigned int c,double *m);

void MatrixVectorProduct(unsigned int r,unsigned int c,double *A,double *b,double *o);

void TMatrixVectorProduct(unsigned int r,unsigned int c,double *A,double *b,double *o);

void TMatrixVectorStrideProduct(unsigned int r,unsigned int c,double *A,unsigned int s,double *b,double *o);

void MatrixMatrixProduct(unsigned int ra,unsigned int ca,double *A,
                         unsigned int cb,double *B,double *C);
void TMatrixMatrixProduct(unsigned int ra,unsigned int ca,double *A,
                          unsigned int cb,double *B,double *C);

double MinCosinusBetweenSubSpaces(unsigned int m,unsigned int k,double *T1,double *T2);

void SubMatrixFromMatrix(unsigned int nr1,unsigned int nc1,double *m1,
                         unsigned int nri,unsigned int nci,
                         unsigned int nr,unsigned int nc,double *m);

void SubMatrixFromTMatrix(unsigned int nr1,unsigned int nc1,double *m1,
                          unsigned int nri,unsigned int nci,
                          unsigned int nr,unsigned int nc,double *m);

void PrintVector(FILE *f,char *label,unsigned int n,double *v);

void SaveVector(FILE *f,unsigned int n,double *v);

void PrintMatrix(FILE *f,char *label,unsigned int r,unsigned int c,double *m);

void PrintTMatrix(FILE *f,char *label,unsigned int r,unsigned int c,double *m);

#endif