jens
f66d39c11c
- use blaze matrix library instead of Eigen git-svn-id: http://moon:8086/svn/software/trunk/libsrc/cpp@943 b431acfa-c32f-4a4a-93f1-934dc6c82436
115 lines
2.1 KiB
C++
115 lines
2.1 KiB
C++
/*
|
|
* File: Vector.hpp
|
|
* Author: jens
|
|
*
|
|
* Created on 19. Dezember 2014, 18:24
|
|
*/
|
|
|
|
#ifndef _RADIO_VECTOR_HPP_
|
|
#define _RADIO_VECTOR_HPP_
|
|
|
|
#define VECTOR_USE_BLAZE
|
|
|
|
#include <complex>
|
|
#include <radio/cpx.h>
|
|
#include <radio/real.h>
|
|
|
|
#ifndef radio_float_t
|
|
#define radio_float_t float
|
|
#endif
|
|
|
|
#ifdef VECTOR_USE_EIGEN
|
|
#include <Eigen/Dense>
|
|
#endif
|
|
|
|
#ifdef VECTOR_USE_BLAZE
|
|
#include <blaze/Blaze.h>
|
|
#endif
|
|
|
|
#ifdef VECTOR_USE_ARMADILLO
|
|
#include <armadillo>
|
|
#endif
|
|
|
|
namespace Radio
|
|
{
|
|
|
|
typedef radio_float_t RealScalar;
|
|
typedef std::complex<radio_float_t> ComplexScalar;
|
|
|
|
#ifdef VECTOR_USE_EIGEN
|
|
|
|
typedef Eigen::Matrix<RealScalar, Eigen::Dynamic, 1> RVec;
|
|
typedef Eigen::Matrix<ComplexScalar, Eigen::Dynamic, 1> CVec;
|
|
|
|
typedef Eigen::Matrix<RealScalar, Eigen::Dynamic, Eigen::Dynamic> RMat;
|
|
typedef Eigen::Matrix<ComplexScalar, Eigen::Dynamic, Eigen::Dynamic> CMat;
|
|
|
|
#endif
|
|
|
|
#ifdef VECTOR_USE_BLAZE
|
|
typedef blaze::DynamicVector<RealScalar> RVec;
|
|
typedef blaze::DynamicVector<ComplexScalar> CVec;
|
|
|
|
typedef blaze::DynamicMatrix<RealScalar> RMat;
|
|
typedef blaze::DynamicMatrix<ComplexScalar> CMat;
|
|
|
|
inline void ror(auto &vec)
|
|
{
|
|
vec = elements(vec, [size=vec.size()](size_t i){ return (i-1) % size; }, vec.size());
|
|
}
|
|
|
|
inline void rol(auto &vec)
|
|
{
|
|
vec = elements(vec, [size=vec.size()](size_t i){ return (i+1) % size; }, vec.size());
|
|
}
|
|
|
|
inline void insert_left(auto &vec, auto const &x)
|
|
{
|
|
ror(vec);
|
|
vec[0] = x;
|
|
}
|
|
|
|
inline void insert_right(auto &vec, auto const &x)
|
|
{
|
|
rol(vec);
|
|
vec[vec.size()-1] = x;
|
|
}
|
|
|
|
#endif
|
|
|
|
#ifdef VECTOR_USE_ARMADILLO
|
|
#error "Vector: Implement armadillo support!"
|
|
#endif
|
|
|
|
static ComplexScalar toComplexScalar(cpx_t v)
|
|
{
|
|
return ComplexScalar(v.real, v.imag);
|
|
}
|
|
|
|
static cpx_t toCpx(ComplexScalar v)
|
|
{
|
|
return Cpx(v.real(), v.imag());
|
|
}
|
|
|
|
static void toComplexVector(CVec &dst, cpx_t const *pSrc, uint32_t len)
|
|
{
|
|
for (uint32_t i=0; i < len; i++)
|
|
{
|
|
dst[i] = ComplexScalar(pSrc[i].real, pSrc[i].imag);
|
|
}
|
|
}
|
|
|
|
static void toCpxVector(cpx_t *pDst, CVec const &src, uint32_t len)
|
|
{
|
|
for (uint32_t i=0; i < len; i++)
|
|
{
|
|
pDst[i].real = src[i].real();
|
|
pDst[i].imag = src[i].imag();
|
|
}
|
|
}
|
|
|
|
} // ::Radio
|
|
|
|
#endif /* _RADIO_VECTOR_HPP_ */
|
|
|