The "r8brain-free-src" library namespace. More...

Classes
class	CDSPBlockConvolver
	Single-block overlap-save convolution processing class. More...

class	CDSPFIRFilter
	Calculation and storage class for FIR filters. More...

class	CDSPFIRFilterCache
	FIR filter cache class. More...

class	CDSPFracDelayFilterBank
	Sinc function-based fractional delay filter bank class. More...

class	CDSPFracDelayFilterBankCache
	Fractional delay filter cache class. More...

class	CDSPFracInterpolator
	Fractional delay filter bank-based interpolator class. More...

class	CDSPHBDownsampler
	Half-band downsampler class. More...

class	CDSPHBUpsampler
	Half-band upsampling class. More...

class	CDSPProcessor
	The base virtual class for DSP processing algorithms. More...

class	CDSPRealFFT
	Real-valued FFT transform class. More...

class	CDSPRealFFTKeeper
	A "keeper" class for real-valued FFT transform objects. More...

class	CDSPRealFFTSetup
	Real-valued FFT setup class. More...

class	CDSPResampler
	The master sample rate converter (resampler) class. More...

class	CDSPResampler16
	The resampler class for 16-bit resampling. More...

class	CDSPResampler16IR
	The resampler class for 16-bit impulse response resampling. More...

class	CDSPResampler24
	The resampler class for 24-bit resampling. More...

class	CDSPSincFilterGen
	Sinc function-based FIR filter generator class. More...

class	CFixedBuffer
	Templated memory buffer class for element buffers of fixed capacity. More...

class	CPtrKeeper
	Pointer-to-object "keeper" class with automatic deletion. More...

class	CRefKeeper
	Reference "keeper" class with automatic unref() call. More...

class	CSineGen
	Sine signal generator class. More...

class	CSinglyLinkedListItem
	Class for items of a singly-linked list. More...

class	CStdClassAllocator
	The default base class for objects created on heap. More...

class	CStdMemAllocator
	The default base class for objects that allocate blocks of memory. More...

class	CSyncKeeper
	A "keeper" class for CSyncObject-based synchronization. More...

class	CSyncObject
	Multi-threaded synchronization object class. More...

Typedefs
typedef double	realfft_t
	Forward transform's native type.

Enumerations
enum	EDSPFilterPhaseResponse { fprLinearPhase = 0 , fprMinPhase }
	Enumeration of filter's phase responses. More...

Functions
template<typename T>
T *	align_ptr (T *const ptr, const uintptr_t align)
	Forces the provided `ptr` pointer to be aligned to `align` bytes.

double	besselI0 (const double x)
	1st kind, 0th order modified Bessel function of a value.

double	calcFIRFilterGroupDelay (const double *const flt, const int fltlen, const double th)
	FIR filter's group delay calculation function.

void	calcFIRFilterResponse (const double *flt, int fltlen, const double th, double &re0, double &im0, const int fltlat=0)
	FIR filter's frequency response calculation.

void	calcMinPhaseTransform (double const Kernel, const int KernelLen, const int LenMult=2, const bool DoFinalMul=true, double const DCGroupDelay=R8B_NULL)
	Calculates the minimum-phase transform of the filter kernel, using a discrete Hilbert transform in cepstrum domain.

void	calcSpline2p8Coeffs (double *const c, const double xm3, const double xm2, const double xm1, const double x0, const double x1, const double x2, const double x3, const double x4)
	Calculates 2nd order spline coefficients, using 8 points.

void	calcSpline3p4Coeffs (double const c, const double const y)
	Calculates 3rd order spline coefficients, using 4 points.

void	calcSpline3p6Coeffs (double const c, const double const y)
	Calculates 3rd order spline coefficients, using 6 points.

void	calcSpline3p8Coeffs (double *const c, const double xm3, const double xm2, const double xm1, const double x0, const double x1, const double x2, const double x3, const double x4)
	Calculates 3rd order spline coefficients, using 8 points.

double	clampr (const double Value, const double minv, const double maxv)
	Clamps a value to be within the specified min-max range.

template<typename T>
T *	construct_ptr (void *const ptr, const size_t c)
	Performs placement `new` to turn a block of unoccupied memory into a "constructed" array of elements of type `T`.

bool	findGCD (double l, double s, double &GCD)
	Interatively searches for a greatest common denominator (GCD) of 2 numbers.

int	getBitOccupancy (const int v)
	Calculate the exact number of bits a value needs for representation.

bool	getWholeStepping (const double SSampleRate, const double DSampleRate, int &ResInStep, int &ResOutStep)
	Evaluates source and destination sample rate ratio and returns the required input and output stepping.

template<typename T>
T	max (const T &v1, const T &v2)
	Returns maximum of two values.

template<typename T>
T	min (const T &v1, const T &v2)
	Returns minimum of two values.

void	normalizeFIRFilter (double *const p, const int l, const double DCGain, const int pstep=1)
	FIR filter's gain normalization.

double	pow_a (const double v, const double p)
	Power of an absolute value.

double	sqr (const double x)
	Returns square ot a value.

Variables
R8B_CONST double	R8B_2PI = 6.28318530717958648
	Equals `2*pi`.

R8B_CONST double	R8B_3PI = 9.42477796076937972
	Equals `3*pi`.

R8B_CONST double	R8B_PI = 3.14159265358979324
	Equals `pi`.

R8B_CONST double	R8B_PId2 = 1.57079632679489662
	Equals `0.5*pi`.

Detailed Description

The "r8brain-free-src" library namespace.

The "r8brain-free-src" sample rate converter library namespace.

Enumeration Type Documentation

◆ EDSPFilterPhaseResponse

enum r8b::EDSPFilterPhaseResponse

Enumeration of filter's phase responses.

Enumerator

fprLinearPhase

Linear-phase response. Features a linear-phase, high-latency response, with the latency expressed as an integer value.

fprMinPhase

Minimum-phase response. Features a minimal-latency response, but the response's phase is non-linear. The latency is usually expressed as a non-integer value, and is usually small, but is never equal to zero. The minimum-phase filter is obtained from a linear-phase filter. Note that since in the context of r8brain-free-src other filters (interpolation, half-band) remain linear-phase, the resulting phase will be "intermediate". The minimum-phase transformation has precision limits: this may skew both the -3 dB point and attenuation of the filter being transformed: as it was measured, the skew happens purely at random, and in most cases is within tolerable range. In a small (1%) random subset of cases the skew is bigger and cannot be predicted. Minimum-phase transform requires 64-bit floating-point FFT; results with 32-bit float FFT are far from optimal.

Function Documentation

◆ align_ptr()

template<typename T>

T * r8b::align_ptr	(	T *const	ptr,
		const uintptr_t	align )

Forces the provided ptr pointer to be aligned to align bytes.

Works with power-of-2 alignments only.

Parameters

ptr	Pointer to align.
align	Alignment, in bytes, power-of-2.

Template Parameters

T	Pointer's element type.

Returns: Aligned pointer.

◆ besselI0()

double r8b::besselI0 ( const double x )

1st kind, 0th order modified Bessel function of a value.

Parameters

x	Input value.

Returns: Calculated zero-th order modified Bessel function of the first kind of the input value. Approximate value.

◆ calcFIRFilterGroupDelay()

double r8b::calcFIRFilterGroupDelay	(	const double *const	flt,
		const int	fltlen,
		const double	th )

FIR filter's group delay calculation function.

Function calculates group delay of the specified FIR filter at the specified circular frequency. The group delay is calculated by evaluating the filter's response at close side-band frequencies of "th".

Parameters

flt	FIR filter's coefficients.
fltlen	Number of coefficients (taps) in the filter.
th	Circular frequency [0; pi].

Returns: Resulting group delay at the specified frequency, in samples.

◆ calcFIRFilterResponse()

void r8b::calcFIRFilterResponse	(	const double *	flt,
		int	fltlen,
		const double	th,
		double &	re0,
		double &	im0,
		const int	fltlat = 0 )

FIR filter's frequency response calculation.

Function calculates frequency response of the specified FIR filter at the specified circular frequency. Phase can be calculated as atan2( im, re ).

Parameters

	flt	FIR filter's coefficients.
	fltlen	Number of coefficients (taps) in the filter.
	th	Circular frequency [0; pi].
[out]	re0	Resulting real part of the complex frequency response.
[out]	im0	Resulting imaginary part of the complex frequency response.
	fltlat	Filter's latency, in samples.

◆ calcMinPhaseTransform()

void r8b::calcMinPhaseTransform	(	double *const	Kernel,
		const int	KernelLen,
		const int	LenMult = 2,
		const bool	DoFinalMul = true,
		double *const	DCGroupDelay = R8B_NULL )

Calculates the minimum-phase transform of the filter kernel, using a discrete Hilbert transform in cepstrum domain.

For more details, see part III.B of http://www.hpl.hp.com/personal/Niranjan_Damera-Venkata/files/ComplexMinPhase.pdf

Parameters

[in,out]	Kernel	Filter kernel buffer.
	KernelLen	Filter kernel's length, in samples.
	LenMult	Kernel length multiplier. Used as a coefficient of oversampling in the frequency domain. Such oversampling is needed to improve the precision of the minimum-phase transform. If the filter's attenuation is high, this multiplier should be increased or otherwise the required attenuation will not be reached due to "smoothing" effect of this transform.
	DoFinalMul	`true` if the final multiplication after transform should be performed. Such multiplication returns the gain of the signal to its original value. This parameter can be set to `false` if normalization of the resulting filter kernel is planned to be used.
[out]	DCGroupDelay	If not `nullptr`, this variable receives group delay at DC offset, in samples (can be a non-integer value).

◆ calcSpline2p8Coeffs()

void r8b::calcSpline2p8Coeffs	(	double *const	c,
		const double	xm3,
		const double	xm2,
		const double	xm1,
		const double	x0,
		const double	x1,
		const double	x2,
		const double	x3,
		const double	x4 )

Calculates 2nd order spline coefficients, using 8 points.

Function calculates coefficients used to calculate 2nd order spline (polynomial) on the equidistant lattice, using 8 points. This function is based on the calcSpline3p8Coeffs() function, but without the 3rd order coefficient.

Parameters

[out]	c	Output coefficients buffer, length = 3.
	xm3	Point at x-3 position.
	xm2	Point at x-2 position.
	xm1	Point at x-1 position.
	x0	Point at x position.
	x1	Point at x+1 position.
	x2	Point at x+2 position.
	x3	Point at x+3 position.
	x4	Point at x+4 position.

◆ calcSpline3p4Coeffs()

void r8b::calcSpline3p4Coeffs	(	double *const	c,
		const double *const	y )

Calculates 3rd order spline coefficients, using 4 points.

Function calculates coefficients used to calculate 3rd order segment interpolation polynomial on the equidistant lattice, using 4 points.

Parameters

[out]	c	Output coefficients buffer, length = 4.
[in]	y	Equidistant point values. Value at offset 1 corresponds to x=0 point.

◆ calcSpline3p6Coeffs()

void r8b::calcSpline3p6Coeffs	(	double *const	c,
		const double *const	y )

Calculates 3rd order spline coefficients, using 6 points.

Function calculates coefficients used to calculate 3rd order segment interpolation polynomial on the equidistant lattice, using 6 points.

Parameters

[out]	c	Output coefficients buffer, length = 4.
[in]	y	Equidistant point values. Value at offset 2 corresponds to x=0 point.

◆ calcSpline3p8Coeffs()

void r8b::calcSpline3p8Coeffs	(	double *const	c,
		const double	xm3,
		const double	xm2,
		const double	xm1,
		const double	x0,
		const double	x1,
		const double	x2,
		const double	x3,
		const double	x4 )

Calculates 3rd order spline coefficients, using 8 points.

Function calculates coefficients used to calculate 3rd order spline (polynomial) on the equidistant lattice, using 8 points.

Parameters

[out]	c	Output coefficients buffer, length = 4.
	xm3	Point at x-3 position.
	xm2	Point at x-2 position.
	xm1	Point at x-1 position.
	x0	Point at x position.
	x1	Point at x+1 position.
	x2	Point at x+2 position.
	x3	Point at x+3 position.
	x4	Point at x+4 position.

◆ clampr()

double r8b::clampr	(	const double	Value,
		const double	minv,
		const double	maxv )

Clamps a value to be within the specified min-max range.

Function "clamps" (clips) the specified value so that it is not lesser than "minv", and not greater than "maxv".

Parameters

Value	Value to clamp.
minv	Minimal allowed value.
maxv	Maximal allowed value.

Returns: "Clamped" value.

◆ construct_ptr()

template<typename T>

T * r8b::construct_ptr	(	void *const	ptr,
		const size_t	c )

Performs placement new to turn a block of unoccupied memory into a "constructed" array of elements of type T.

Parameters

ptr	Pointer for placement of elements.
c	The number of elements to "construct".

Template Parameters

T	Required element type.

Returns: The pointer to a formally-constructed array of elements.

◆ findGCD()

bool r8b::findGCD	(	double	l,
		double	s,
		double &	GCD )

Interatively searches for a greatest common denominator (GCD) of 2 numbers.

Parameters

	l	Number 1.
	s	Number 2.
[out]	GCD	Resulting GCD.

Returns: true if the greatest common denominator of 2 numbers was found.

◆ getBitOccupancy()

int r8b::getBitOccupancy ( const int v )

Calculate the exact number of bits a value needs for representation.

Parameters

v	Input value.

Returns: Calculated bit occupancy of the specified input value. Bit occupancy means how many significant lower bits are necessary to store a specified value. Function treats the input value as unsigned.

◆ getWholeStepping()

bool r8b::getWholeStepping	(	const double	SSampleRate,
		const double	DSampleRate,
		int &	ResInStep,
		int &	ResOutStep )

Evaluates source and destination sample rate ratio and returns the required input and output stepping.

Function returns false if whole stepping cannot be used to perform interpolation using these sample rates.

Parameters

	SSampleRate	Source sample rate.
	DSampleRate	Destination sample rate.
[out]	ResInStep	Resulting input step.
[out]	ResOutStep	Resulting output step.

Returns: true if stepping was acquired.

◆ max()

template<typename T>

T r8b::max	(	const T &	v1,
		const T &	v2 )

Returns maximum of two values.

Parameters

v1	Value 1.
v2	Value 2.

Template Parameters

T	Values' type.

Returns: The maximum of 2 values.

◆ min()

template<typename T>

T r8b::min	(	const T &	v1,
		const T &	v2 )

Returns minimum of two values.

Parameters

v1	Value 1.
v2	Value 2.

Template Parameters

T	Values' type.

Returns: The minimum of 2 values.

◆ normalizeFIRFilter()

void r8b::normalizeFIRFilter	(	double *const	p,
		const int	l,
		const double	DCGain,
		const int	pstep = 1 )

FIR filter's gain normalization.

Function normalizes FIR filter so that its frequency response at DC is equal to DCGain.

Parameters

[in,out]	p	Filter coefficients.
	l	Filter length.
	DCGain	Filter's gain at DC (linear, non-decibel value).
	pstep	"p" array step.

◆ pow_a()

double r8b::pow_a	(	const double	v,
		const double	p )

Power of an absolute value.

Parameters

v	Input value.
p	Power factor.

Returns: Returns a precise, but generally approximate pow() function's value of absolute of input value.

◆ sqr()

double r8b::sqr ( const double x )

Returns square ot a value.

Parameters

x	Value to square.

Returns: Squared value of the argument.

Classes

Typedefs

Enumerations

Functions

Variables

Detailed Description

Enumeration Type Documentation

◆ EDSPFilterPhaseResponse

Function Documentation

◆ align_ptr()

◆ besselI0()

◆ calcFIRFilterGroupDelay()

◆ calcFIRFilterResponse()

◆ calcMinPhaseTransform()

◆ calcSpline2p8Coeffs()

◆ calcSpline3p4Coeffs()

◆ calcSpline3p6Coeffs()

◆ calcSpline3p8Coeffs()

◆ clampr()

◆ construct_ptr()

◆ findGCD()

◆ getBitOccupancy()

◆ getWholeStepping()

◆ max()

◆ min()

◆ normalizeFIRFilter()

◆ pow_a()

◆ sqr()