FFmpeg coverage

Directory: ../../../ffmpeg/
File: src/libavutil/float_dsp.c
Date: 2026-07-16 17:05:34
Exec Total Coverage
Lines: 73 74 98.6%
Functions: 12 12 100.0%
Branches: 23 24 95.8%
Line Branch Exec Source
1 /*
2 * Copyright 2005 Balatoni Denes
3 * Copyright 2006 Loren Merritt
4 *
5 * This file is part of FFmpeg.
6 *
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
11 *
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
16 *
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 */
21
22 #include "config.h"
23 #include "attributes.h"
24 #include "float_dsp.h"
25 #include "mem.h"
26
27 540435 static void vector_fmul_c(float *dst, const float *src0, const float *src1,
28 int len)
29 {
30 int i;
31
2/2
✓ Branch 0 taken 60451328 times.
✓ Branch 1 taken 540435 times.
 60991763 for (i = 0; i < len; i++)
32 60451328 dst[i] = src0[i] * src1[i];
33 540435 }
34
35 4 static void vector_dmul_c(double *dst, const double *src0, const double *src1,
36 int len)
37 {
38 int i;
39
2/2
✓ Branch 0 taken 1024 times.
✓ Branch 1 taken 4 times.
 1028 for (i = 0; i < len; i++)
40 1024 dst[i] = src0[i] * src1[i];
41 4 }
42
43 24583 static void vector_fmac_scalar_c(float *dst, const float *src, float mul,
44 int len)
45 {
46 int i;
47
2/2
✓ Branch 0 taken 11636384 times.
✓ Branch 1 taken 24583 times.
 11660967 for (i = 0; i < len; i++)
48 11636384 dst[i] += src[i] * mul;
49 24583 }
50
51 4 static void vector_dmac_scalar_c(double *dst, const double *src, double mul,
52 int len)
53 {
54 int i;
55
2/2
✓ Branch 0 taken 1024 times.
✓ Branch 1 taken 4 times.
 1028 for (i = 0; i < len; i++)
56 1024 dst[i] += src[i] * mul;
57 4 }
58
59 1284915 static void vector_fmul_scalar_c(float *dst, const float *src, float mul,
60 int len)
61 {
62 int i;
63
2/2
✓ Branch 0 taken 30426592 times.
✓ Branch 1 taken 1284915 times.
 31711507 for (i = 0; i < len; i++)
64 30426592 dst[i] = src[i] * mul;
65 1284915 }
66
67 4 static void vector_dmul_scalar_c(double *dst, const double *src, double mul,
68 int len)
69 {
70 int i;
71
2/2
✓ Branch 0 taken 1024 times.
✓ Branch 1 taken 4 times.
 1028 for (i = 0; i < len; i++)
72 1024 dst[i] = src[i] * mul;
73 4 }
74
75 235697 static void vector_fmul_window_c(float *dst, const float *src0,
76 const float *src1, const float *win, int len)
77 {
78 int i, j;
79
80 235697 dst += len;
81 235697 win += len;
82 235697 src0 += len;
83
84
2/2
✓ Branch 0 taken 48444996 times.
✓ Branch 1 taken 235697 times.
 48680693 for (i = -len, j = len - 1; i < 0; i++, j--) {
85 48444996 float s0 = src0[i];
86 48444996 float s1 = src1[j];
87 48444996 float wi = win[i];
88 48444996 float wj = win[j];
89 48444996 dst[i] = s0 * wj - s1 * wi;
90 48444996 dst[j] = s0 * wi + s1 * wj;
91 }
92 235697 }
93
94 3920290 static void vector_fmul_add_c(float *dst, const float *src0, const float *src1,
95 const float *src2, int len){
96 int i;
97
98
2/2
✓ Branch 0 taken 234706304 times.
✓ Branch 1 taken 3920290 times.
 238626594 for (i = 0; i < len; i++)
99 234706304 dst[i] = src0[i] * src1[i] + src2[i];
100 3920290 }
101
102 411438 static void vector_fmul_reverse_c(float *dst, const float *src0,
103 const float *src1, int len)
104 {
105 int i;
106
107 411438 src1 += len-1;
108
2/2
✓ Branch 0 taken 134833024 times.
✓ Branch 1 taken 411438 times.
 135244462 for (i = 0; i < len; i++)
109 134833024 dst[i] = src0[i] * src1[-i];
110 411438 }
111
112 595100 static void butterflies_float_c(float *restrict v1, float *restrict v2,
113 int len)
114 {
115 int i;
116
117
2/2
✓ Branch 0 taken 8787332 times.
✓ Branch 1 taken 595100 times.
 9382432 for (i = 0; i < len; i++) {
118 8787332 float t = v1[i] - v2[i];
119 8787332 v1[i] += v2[i];
120 8787332 v2[i] = t;
121 }
122 595100 }
123
124 836277 double ff_scalarproduct_double_c(const double *v1, const double *v2,
125 size_t len)
126 {
127 836277 double p = 0.0;
128
129
2/2
✓ Branch 0 taken 6691353 times.
✓ Branch 1 taken 836277 times.
 7527630 for (size_t i = 0; i < len; i++)
130 6691353 p += v1[i] * v2[i];
131
132 836277 return p;
133 }
134
135 1348 av_cold AVFloatDSPContext *avpriv_float_dsp_alloc(int bit_exact)
136 {
137 1348 AVFloatDSPContext *fdsp = av_mallocz(sizeof(AVFloatDSPContext));
138
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1348 times.
 1348 if (!fdsp)
139 ✗ return NULL;
140
141 1348 fdsp->vector_fmul = vector_fmul_c;
142 1348 fdsp->vector_dmul = vector_dmul_c;
143 1348 fdsp->vector_fmac_scalar = vector_fmac_scalar_c;
144 1348 fdsp->vector_fmul_scalar = vector_fmul_scalar_c;
145 1348 fdsp->vector_dmac_scalar = vector_dmac_scalar_c;
146 1348 fdsp->vector_dmul_scalar = vector_dmul_scalar_c;
147 1348 fdsp->vector_fmul_window = vector_fmul_window_c;
148 1348 fdsp->vector_fmul_add = vector_fmul_add_c;
149 1348 fdsp->vector_fmul_reverse = vector_fmul_reverse_c;
150 1348 fdsp->butterflies_float = butterflies_float_c;
151 1348 fdsp->scalarproduct_float = ff_scalarproduct_float_c;
152 1348 fdsp->scalarproduct_double = ff_scalarproduct_double_c;
153
154 #if ARCH_AARCH64
155 ff_float_dsp_init_aarch64(fdsp);
156 #elif ARCH_ARM
157 ff_float_dsp_init_arm(fdsp);
158 #elif ARCH_PPC
159 ff_float_dsp_init_ppc(fdsp, bit_exact);
160 #elif ARCH_RISCV
161 ff_float_dsp_init_riscv(fdsp);
162 #elif ARCH_X86 && HAVE_X86ASM
163 1348 ff_float_dsp_init_x86(fdsp);
164 #elif ARCH_MIPS
165 ff_float_dsp_init_mips(fdsp);
166 #endif
167 1348 return fdsp;
168 }
169