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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
/******************************************************************************
*
* Copyright 2022 Google LLC
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at:
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
******************************************************************************/
#include "neon.h"
#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>
/* -------------------------------------------------------------------------- */
#define TEST_NEON
#include <ltpf.c>
void lc3_put_bits_generic(lc3_bits_t *a, unsigned b, int c)
{ (void)a, (void)b, (void)c; }
unsigned lc3_get_bits_generic(struct lc3_bits *a, int b)
{ return (void)a, (void)b, 0; }
/* -------------------------------------------------------------------------- */
static int check_resampler()
{
int16_t __x[60+480], *x = __x + 60;
for (int i = -60; i < 480; i++)
x[i] = rand() & 0xffff;
struct lc3_ltpf_hp50_state hp50 = { 0 }, hp50_neon = { 0 };
int16_t y[128], y_neon[128];
resample_16k_12k8(&hp50, x, y, 128);
neon_resample_16k_12k8(&hp50_neon, x, y_neon, 128);
if (memcmp(y, y_neon, 128 * sizeof(*y)) != 0)
return printf("Error\n"), -1;
resample_32k_12k8(&hp50, x, y, 128);
neon_resample_32k_12k8(&hp50_neon, x, y_neon, 128);
if (memcmp(y, y_neon, 128 * sizeof(*y)) != 0)
return printf("Error\n"), -1;
resample_48k_12k8(&hp50, x, y, 128);
neon_resample_48k_12k8(&hp50_neon, x, y_neon, 128);
if (memcmp(y, y_neon, 128 * sizeof(*y)) != 0)
return -1;
return 0;
}
static int check_dot()
{
int16_t x[200];
for (int i = 0; i < 200; i++)
x[i] = rand() & 0xffff;
float y = dot(x, x+3, 128);
float y_neon = neon_dot(x, x+3, 128);
if (y != y_neon)
return -1;
return 0;
}
static int check_correlate()
{
int16_t alignas(4) a[500], b[500];
float y[100], y_neon[100];
for (int i = 0; i < 500; i++) {
a[i] = rand() & 0xffff;
b[i] = rand() & 0xffff;
}
correlate(a, b+200, 128, y, 100);
neon_correlate(a, b+200, 128, y_neon, 100);
if (memcmp(y, y_neon, 100 * sizeof(*y)) != 0)
return -1;
correlate(a, b+199, 128, y, 99);
neon_correlate(a, b+199, 128, y_neon, 99);
if (memcmp(y, y_neon, 99 * sizeof(*y)) != 0)
return -1;
return 0;
}
int check_ltpf(void)
{
int ret;
if ((ret = check_resampler()) < 0)
return ret;
if ((ret = check_dot()) < 0)
return ret;
if ((ret = check_correlate()) < 0)
return ret;
return 0;
}
|
