summaryrefslogtreecommitdiff
path: root/Drivers/CMSIS/DSP/DSP_Lib_TestSuite/RefLibs/src/BasicMathFunctions/scale.c
blob: 5ab655ce80f20c1d4efde1af2226ed2c811ebbdb (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
#include "ref.h"

void ref_scale_f32(
  float32_t * pSrc,
  float32_t scale,
  float32_t * pDst,
  uint32_t blockSize)
{
	uint32_t i;
	
	for(i=0;i<blockSize;i++)
	{
		pDst[i] = pSrc[i] * scale;
	}
}

void ref_scale_q31(
  q31_t * pSrc,
  q31_t scaleFract,
  int8_t shift,
  q31_t * pDst,
  uint32_t blockSize)
{
	uint32_t i;
  int8_t kShift = shift + 1;                     /* Shift to apply after scaling */
  int8_t sign = (kShift & 0x80);
	q63_t temp;
	
	for(i=0;i<blockSize;i++)
	{
		temp = ((q63_t) pSrc[i] * scaleFract) >> 32;
		if (sign)
			pDst[i] = temp >> -kShift;
		else
			pDst[i] = ref_sat_q31( (q63_t)temp << kShift );
	}
}

void ref_scale_q15(
  q15_t * pSrc,
  q15_t scaleFract,
  int8_t shift,
  q15_t * pDst,
  uint32_t blockSize)
{
	uint32_t i;
  int8_t kShift = 15 - shift;                      /* Shift to apply after scaling */
	
	for(i=0;i<blockSize;i++)
	{
		pDst[i] = ref_sat_q15(((q31_t) pSrc[i] * scaleFract) >> kShift);
	}
}

void ref_scale_q7(
  q7_t * pSrc,
  q7_t scaleFract,
  int8_t shift,
  q7_t * pDst,
  uint32_t blockSize)
{
	uint32_t i;
  int8_t kShift = 7 - shift;                      /* Shift to apply after scaling */
	
	for(i=0;i<blockSize;i++)
	{
		pDst[i] = ref_sat_q7(((q15_t) pSrc[i] * scaleFract) >> kShift);
	}
}