1 /*
2 * Copyright (c) 2022-2022 Huawei Device Co., Ltd. All rights reserved.
3 *
4 * Redistribution and use in source and binary forms, with or without modification,
5 * are permitted provided that the following conditions are met:
6 *
7 * 1. Redistributions of source code must retain the above copyright notice, this list of
8 * conditions and the following disclaimer.
9 *
10 * 2. Redistributions in binary form must reproduce the above copyright notice, this list
11 * of conditions and the following disclaimer in the documentation and/or other materials
12 * provided with the distribution.
13 *
14 * 3. Neither the name of the copyright holder nor the names of its contributors may be used
15 * to endorse or promote products derived from this software without specific prior written
16 * permission.
17 *
18 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
19 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
20 * THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
21 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR
22 * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
23 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
24 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
25 * OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
26 * WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
27 * OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
28 * ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
29 */
30
31 #include <stdio.h>
32 #include "arm_math.h"
33 #include "sys/time.h"
34
35 #define ROW 4
36 #define COW 4
37 #define LOOP_TIMES 1000
38 #define MS_PER_S 1000
39
40 float32_t F32MatrixTransA[ROW * COW];
41 float32_t F32MatrixATMA[ROW * COW];
42 float32_t F32MatrixATMAI[ROW * COW];
43
44 const float32_t F32MatrixA[ROW * COW] = {1.0, 32.0, 16.0, 64.0, 1.0, 32.0, 256.0, 256.0, 1.0, 64.0, 16.0, 512.0, 1.0, 64.0, 128.0, 1024.0};
45
test_dsp(void)46 int32_t test_dsp(void)
47 {
48 int index = LOOP_TIMES;
49 uint32_t srcRows = ROW;
50 uint32_t srcColumns = COW;
51 struct timeval start, end;
52 uint32_t timeCost;
53 arm_matrix_instance_f32 MatrixA;
54 arm_matrix_instance_f32 MatrixAT;
55 arm_matrix_instance_f32 MatrixATMA;
56 arm_matrix_instance_f32 MatrixATMAI;
57
58 gettimeofday(&start, NULL);
59 while (index--) {
60 arm_mat_init_f32(&MatrixA, srcRows, srcColumns, (float32_t *)F32MatrixA);
61 arm_mat_init_f32(&MatrixAT, srcRows, srcColumns, F32MatrixTransA);
62 arm_mat_trans_f32(&MatrixA, &MatrixAT);
63 arm_mat_init_f32(&MatrixATMA, srcRows, srcColumns, F32MatrixATMA);
64 arm_mat_mult_f32(&MatrixAT, &MatrixA, &MatrixATMA);
65 arm_mat_init_f32(&MatrixATMAI, srcRows, srcColumns, F32MatrixATMAI);
66 arm_mat_inverse_f32(&MatrixATMA, &MatrixATMAI);
67 arm_mat_mult_f32(&MatrixATMAI, &MatrixAT, &MatrixATMA);
68 }
69
70 gettimeofday(&end, NULL);
71 timeCost = (unsigned int)(end.tv_sec - start.tv_sec) * MS_PER_S + (unsigned int)(end.tv_usec - start.tv_usec) /
72 MS_PER_S;
73 printf("time cost: %d\n", timeCost);
74 return 0;
75 }
76
77